FFT高精度乘法

众所周知，FFT(Fast Fourier Transformation)，即为快速傅氏变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的.
那么我们用这个nb的算法解决一些什么沙雕问题呢QwQ？？？
没错！就是大名鼎鼎的高精度乘法
FFT太难了我不会啊QAQ,NTT更难了哇怎么办啊

行吧，上代码

#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
#define N 200000
using namespace std;
const double DFT=2.0,IDFT=-2.0;
const double pi=acos(-1),eps=0.5;
struct Complex
{
    double x,y;
    Complex(double _=0,double __=0):x(_),y(__){}
    Complex operator + (const Complex &b) const{return Complex(x+b.x,y+b.y);}
    Complex operator - (const Complex &b) const{return Complex(x-b.x,y-b.y);}
    Complex operator * (const Complex &b) const{return Complex(x*b.x-y*b.y,x*b.y+y*b.x);}
}a[N],b[N],c[N];
char s1[N],s2[N];
int pos[N],n,top,len,ans[N];
inline void FFT(Complex *a,double mode)
{
    for(int i=0;i<len;i++)if(i<pos[i])swap(a[i],a[pos[i]]);
    for(int i=2,mid=1;i<=len;i<<=1,mid<<=1)
    {
        Complex wm(cos(2.0*pi/i),sin(mode*pi/i));
        for(int j=0;j<len;j+=i)
        {
            Complex w(1.0,0.0);
            for(int k=j;k<j+mid;k++,w=w*wm)
            {
                Complex l=a[k],r=w*a[k+mid];
                a[k]=l+r;a[k+mid]=l-r;
            }
        }
    }
    if(mode==IDFT) for(int i=0;i<len;i++) a[i].x/=len;
    return;
}
int main()
{
    scanf("%d%s%s",&n,s1,s2);
    for(len=1;len<=n<<1;len<<=1);
    for(int i=0;i<len;i++)
    {
        pos[i]=pos[i>>1]>>1;
        if(i&1) pos[i]|=len>>1;
    }
    int o=0;
    for(o=0;o<n;o++){a[o].x=s1[n-o-1]-'0';a[o].y=0.0;}
    while(o<len){a[o].x=0;a[o++].y=0;}
    for(o=0;o<n;o++){b[o].x=s2[n-o-1]-'0';b[o].y=0.0;}
    while(o<len){b[o].x=0;b[o++].y=0;}
    FFT(a,DFT);FFT(b,DFT);
    for(int i=0;i<len;i++) c[i]=a[i]*b[i];
    FFT(c,IDFT);
    for(int i=0;i<len;i++) ans[i]=int(c[i].x+0.5);
    for(int i=0;i<len;i++)
    {
        ans[i+1]+=ans[i]/10;
        ans[i]%=10;
    }
    for(int i=len-1;i>=0;i--) if(ans[i]){top=i;break;}
    for(int i=top;i>=0;i--) printf("%d",ans[i]);
    return 0;
}