FFT模板

我终于下定决心学习FFT了。

orzCHX，得出模板：

#include<cstdio>
#include<cctype>
#include<queue>
#include<cmath>
#include<cstring>
#include<algorithm>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define ren for(int i=first[x];i!=-1;i=next[i])
using namespace std;
inline int read() {
    int x=,f=;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-;
    for(;isdigit(c);c=getchar()) x=x*+c-'';
    return x*f;
}
const int maxn=;
const double PI=acos(-1.0);
struct FFT {
    struct cox {
        double r,i;
        cox(double _r=0.0,double _i=0.0) {r=_r;i=_i;}
        cox operator + (const cox& b) const {return cox(r+b.r,i+b.i);}
        cox operator - (const cox& b) const {return cox(r-b.r,i-b.i);}
        cox operator * (const cox& b) const {return cox(r*b.r-i*b.i,r*b.i+i*b.r);}
    }f[maxn];
    int len;
    void init(int* A,int Len,int L) {
        len=L;rep(i,,Len-) f[i]=cox(A[Len-i-],);
    }
    void cal(int tp) {
        int j=len>>;
        rep(i,,len-) {
            if(i<j) swap(f[i],f[j]);int k=len>>;
            while(j>=k) j-=k,k>>=;j+=k;
        }
        double lm=-*tp*PI;
        for(int i=;i<=len;i<<=) {
            cox wn(cos(lm/i),sin(lm/i));
            for(int j=;j<len;j+=i) {
                cox w(,);
                for(int k=j;k<j+(i>>);k++) {
                    cox u=f[k],v=w*f[k+(i>>)];
                    f[k]=u+v;f[k+(i>>)]=u-v;w=w*wn;
                }
            }
        }
        if(tp<) rep(i,,len-) f[i].r/=len;
    }
}a,b;
void mul(int* A,int* B,int L1,int L2,int& L,int* ans) {
    L=;while(L<L1<<||L<L2<<) L<<=;
    a.init(A,L1,L);b.init(B,L2,L);
    a.cal();b.cal();
    rep(i,,L-) a.f[i]=a.f[i]*b.f[i];
    a.cal(-);rep(i,,L-) ans[i]=int(a.f[i].r+0.5);
}
int A[maxn],B[maxn],ans[maxn];
int main() {
    int L1=read()+,L2=read()+,L;
    rep(i,,L1-) A[i]=read();
    rep(i,,L2-) B[i]=read();
    mul(A,B,L1,L2,L,ans);
    while(L>L1+L2-&&!ans[L-]) L--;
    dwn(i,L-,) printf("%d ",ans[i]);
    return ;
}