bzoj2179 FFT快速傅立叶

来源：互联网发布：程序员眼里只有两种人编辑：程序博客网时间：2024/04/30 19:35

题目

先写模板，再理解吧。

FFT高精度乘法模板。

//FFT 高精度乘法模板 #include<bits/stdc++.h>#define N 131075 using namespace std;int n,c[N],i,j;complex<double> a[N],b[N],p[N];const double pi=acos(-1);void FFT(complex<double> *x,int n,int type){    int i,j,k,t;    for(i=0,t=0;i<n;i++)    {        if(i>t)swap(x[i],x[t]);        for(int j=n>>1;(t^=j)<j;j>>=1);    }    for(int k=2;k<=n;k<<=1)    {        complex<double> wn(cos(type*2*pi/k),sin(type*2*pi/k));        for(int i=0;i<n;i+=k)        {            complex<double> w(1,0),t;            for(int j=0;j<k>>1;w*=wn,j++)            {                t=w*x[i+j+(k>>1)];                x[i+j+(k>>1)]=x[i+j]-t,x[i+j]+=t;            }        }    }}int main(){    //freopen("in.txt","r",stdin);    cin>>n;    for(getchar(),i=n-1;i>=0;i--)a[i]=getchar()-'0';    for(getchar(),i=n-1;i>=0;i--)b[i]=getchar()-'0';    for(j=n,i=1;i>>2<j;i<<=1)n=i;    FFT(a,n,1),FFT(b,n,1);    for(int i=0;i<n;i++)p[i]=a[i]*b[i];    FFT(p,n,-1);    int len=0;    for(int i=0;i<n;i++)c[i]=p[i].real()/n+0.1;    for(int i=0;i<n;i++)if(c[i])        len=i,c[i+1]+=c[i]/10,c[i]%=10;    for(int i=len;i>=0;i--)printf("%d",c[i]);     return 0;}

实在不行背板吧，反正noip不考（flag

阅读全文

0 0