FFT一周目开坑！

先来一段非递归！

#include<bits/stdc++.h>
using namespace std;
#define N ((1<<18)+3)
const double eps=0.5,pi=acos(-1);
struct vec{
    double r,i;
    vec operator +(const vec& w){return (vec){r+w.r,i+w.i};}
    vec operator -(const vec& w){return (vec){r-w.r,i-w.i};}
    vec operator *(const vec& w){return (vec){r*w.r-i*w.i,i*w.r+r*w.i};}
}a[N],b[N];
int n,m,len,rev[N],t;
bool che[N];
int read(){
    int f=1,a=getchar(),tmp=0;
    while(a<'0' || a>'9') {if(a=='-') f=-1; a=getchar();}
    while(a>='0'&& a<='9') tmp=tmp*10+a-'0',a=getchar();
    return tmp*=f;
}
void fft(vec *x,int t){
    memset(che,0,sizeof(che));
    for(int i=0;i<len;i++) if(!che[i]) swap(x[i],x[rev[i]]),che[rev[i]]=1;
    for(int lnow=2;lnow<=len;lnow<<=1){   //从最底层开始往上推，一开始每一小坨只有2个，慢慢往上到最后达到len 
        vec w0=(vec){cos(t*2*pi/lnow),sin(t*2*pi/lnow)},w,t1,t2;
        for(int i=0;i<len;i+=lnow){ //每一坨
            vec w=(vec){1,0};
            for(int j=0;j<lnow/2;j++,w=w*w0){ 
                t1=x[i+j]; t2=w*x[i+j+lnow/2];
                x[i+j]=t1+t2; x[i+j+lnow/2]=t1-t2;
            }
        }
    }
}
int main(){
    n=read(); m=read();
    for(len=1;len<(n+m+1);len<<=1,t++); t--;
    for(int i=1;i<=len;i++) rev[i]=(rev[i>>1]>>1)|(i&1?(1<<t):0); //二进制反转，注意之前的t-- 
    for(int i=0;i<=n;i++) a[i].r=read();
    for(int i=0;i<=m;i++) b[i].r=read();
    fft(a,1); fft(b,1);
    for(int i=0;i<=len;i++) a[i]=a[i]*b[i];
    fft(a,-1);
    for(int i=0;i<=n+m;i++) printf("%d ",(int)(a[i].r/len+eps));
    return 0;
}