蒟蒻的内心独白

什么蝴蝶操作什么复数单位根。。老子统统不会

那么怎么办呢，，，当然是什么板子好记背什么板子咯。。。

所以就来个递归版本的吧，，，

#include<bits/stdc++.h>
#define pi acos(-1)
using namespace std;
typedef complex<double> E;
int n,m;
E a[262145],b[262145];
void fft(E *x,int n,int type)
{
if(n==1)return;
E l[n>>1],r[n>>1];
for(int i=0;i<n;i+=2)
l[i>>1]=x[i],r[i>>1]=x[i+1];
fft(l,n>>1,type);fft(r,n>>1,type);
E wn(cos(2*pi/n),sin(type*2*pi/n)),w(1,0),t;
for(int i=0;i<n>>1;i++,w*=wn)
t=w*r[i],x[i]=l[i]+t,x[i+(n>>1)]=l[i]-t;
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=0,x;i<=n;i++)scanf("%d",&x),a[i]=x;
for(int i=0,x;i<=m;i++)scanf("%d",&x),b[i]=x;
m=n+m;for(n=1;n<=m;n<<=1);
fft(a,n,1);fft(b,n,1);
for(int i=0;i<=n;i++)a[i]=a[i]*b[i];
fft(a,n,-1);
for(int i=0;i<=m;i++)
printf("%d ",(int)(a[i].real()/n+0.5));
return 0;
}卡死算我认了，，，，