【学习笔记】FFT

1、内容

由于noble_太懒 不想写了

非常好的博客：

https://www.cnblogs.com/rvalue/p/7351400.html

http://www.cnblogs.com/candy99/p/6641972.html

http://www.gatevin.moe/acm/fft%E7%AE%97%E6%B3%95%E5%AD%A6%E4%B9%A0%E7%AC%94%E8%AE%B0/

http://hzwer.com/6896.html 黄学长模板

https://oi.men.ci/fft-notes/

https://blog.csdn.net/ggn_2015/article/details/68922404

http://www.cnblogs.com/zwfymqz/p/8244902.html

http://www.cnblogs.com/19992147orz/p/8041323.html

2、模板

洛谷A了，maxn要开大一点

 #include <bits/stdc++.h>
 using namespace std;
 typedef complex<double> com;
 const int maxn=3e7;
 const double PI=acos(-);
 com a[maxn], b[maxn];
 int rev[maxn];

 void FFT(com* a,int n,int type){
     for(int i=;i<n;i++){
         if(rev[i]>i) swap(a[i],a[rev[i]]);
     }

     for(int step=;step<n;step<<=){ //待合并区域中点
         com wn(cos(PI/step),type*sin(PI/step));
         for(int j=;j<n;j+=(step<<)){ //step<<1是区间右端点
             com w(,); //幂
             for(int k=j;k<j+step;k++,w*=wn){//枚举左半部分
                 com x=a[k], y=w*a[k+step];
                 a[k]=x+y; a[k+step]=x-y;
             }
         }
     }
 //    if(type==-1) for(int i=0;i<n;i++) a[i]/=n;
 }
 int main()
 {
     int n1,n2,n,x,L=;
     scanf("%d%d",&n1,&n2);
     for(int i=;i<=n1;i++){
         scanf("%d",&x); a[i]=x;
     }
     for(int i=;i<=n2;i++){
         scanf("%d",&x); b[i]=x;
     }
     for(n=;n<=n1+n2;n*=) L++;
     for(int i=;i<n;i++){
         rev[i]=(rev[i>>]>>)|((i&)<<(L-));
     }
     FFT(a,n,); FFT(b,n,);
     for(int i=;i<=n;i++) a[i]*=b[i];
     FFT(a,n,-);
     for(int i=;i<=n1+n2;i++){
         printf("%d ",(int)(a[i].real()/n+0.5));
     }
     return ;
 }