【UOJ34】—多项式乘法（FFT）

传送门

学了$2$个小时才把$fft$搞懂

我真是太菜了（膜拜$ldx$神仙$20$分钟学完$fft$）

结果发现自己根本看不懂代码（233…）

对着板子敲了一发，等有时间再回来补学习笔记吧

#include<bits/stdc++.h>
using namespace std;
inline int read(){
    char ch=getchar();
    int res=0;
    while(!isdigit(ch))ch=getchar();
    while(isdigit(ch))res=(res<<3)+(res<<1)+(ch^48),ch=getchar();
    return res;
}
const int N=400005;
const double pi=acos(-1.0);
int n,m,lim=1,tim=0,pos[N];
struct plx{
    double x,y;
    friend inline plx operator+(const plx&a,const plx&b){
        return (plx){a.x+b.x,a.y+b.y};
    }
    friend inline plx operator-(const plx&a,const plx&b){
        return (plx){a.x-b.x,a.y-b.y};
    }
    friend inline plx operator*(const plx&a,const plx&b){
        return (plx){a.x*b.x-a.y*b.y,a.y*b.x+a.x*b.y};
    }
}a[N],b[N];
inline void fft(plx a[],int kd){
    for(int i=0;i<lim;i++)if(i<pos[i])swap(a[i],a[pos[i]]);
    for(int mid=1;mid<lim;mid<<=1){
        plx nw=(plx){cos(pi/mid),kd*sin(pi/mid)};
        for(int j=0,len=mid<<1;j<lim;j+=len){
            plx w=(plx){1,0};
            for(int k=0;k<mid;k++,w=w*nw){
                plx a0=a[j+k],a1=w*a[j+k+mid];
                a[j+k]=a0+a1,a[j+k+mid]=a0-a1;
            }
        }
    }
}
int main(){
    n=read(),m=read();
    for(int i=0;i<=n;i++)a[i].x=read();
    for(int i=0;i<=m;i++)b[i].x=read();
    while(lim<=n+m)lim<<=1,tim++;
    for(int i=0;i<lim;i++)pos[i]=(pos[i>>1]>>1)|((i&1)<<(tim-1));
    fft(a,1),fft(b,1);
    for(int i=0;i<lim;i++)a[i]=a[i]*b[i];
    fft(a,-1);
    for(int i=0;i<=n+m;i++){
        cout<<(int)(a[i].x/lim+0.5)<<" ";
    }
}

转载于:https://www.cnblogs.com/stargazer-cyk/p/10366367.html