萌新初学FFT，想来简单的试一下，代码也没有完全搞懂，照着其他博客的思路打了一遍，结果40pts，有RE,WA

#include<bits/stdc++.h>
using namespace std;
typedef struct node{
    double a,b;
    node(double x=0,double y=0){
        a=x,b=y;
    }
    node operator+(const node&t)const{
        return node(a+t.a,b+t.b);
    }
    node operator-(const node&t)const{
        return node(a-t.a,b-t.b);
    }
    node operator*(const node&t)const{
        return node(a*t.a-b*t.b,a*t.b+b*t.a);
    }
}aplusb;
const double pi=2*acos(-1);
char f[20],g[20],c[20];
int n,m,r[20];
aplusb a[20],w[20];
int log_2(int x){
    int ans=0;
    if(x&0xffff0000){
    	ans+=16,x>>=16;
	}
    if(x&0xff00){
    	ans+=8,x>>=8;
	}
    if(x&0xf0){
    	ans+=4,x>>=4;
	}
    if(x&0xc){
    	ans+=2,x>>=2;
	}
    if(x&2){
    	++ans;
	}
    return ans;
}
void init(){
    m+=n;
    int Log_2=log_2(m);
    n=1<<Log_2+1;
    for(int i=0;i<n;++i){
    	r[i]=r[i>>1]>>1|(i&1)<<Log_2;
	} 
    aplusb t(cos(pi/n),sin(pi/n));
    w[0].a=1;
    for(int i=1;i<n;++i){
    	w[i]=w[i-1]*t;
	}
}
void swap(aplusb& a,aplusb& b){
    static aplusb t;
    t=a,a=b,b=t;
}
void fft(bool type){
    for(int i=0;i<n;++i){
    	if(i<r[i]){
			swap(a[i],a[r[i]]);
		}
	}
    static aplusb x,y;
    for(int len=2;len<=n;len<<=1){
    	for(int l=0;l<n;l+=len){
    		for (int i=0,j=len>>1;j<len;++i,++j){
                x=a[l+i],y=w[n/len*i];
                if(type){
                	y=y*a[l+j];
				}else{
                	y=aplusb(y.a,-y.b)*a[l+j];
				}
                a[l+i]=x+y,a[l+j]=x-y;
            }
		}
	}
}
int main(){
    cin>>f>>g;
    n=strlen(f)-1;
	m=strlen(g)-1;
    for(int i=n;~i;--i){
    	a[i].a=(double)(f[n-i]^48);
	}
    for(int i=m;~i;--i){
    	a[i].b=(double)(g[m-i]^48);
	} 
    init();
    fft(1);
    for(int i=0;i<n;++i){
    	a[i]=a[i]*a[i];
	}
    fft(0);
    int t=0;
    for(int i=0;i<=m;++i){
        t+=(int)(a[i].b/n/2.0+0.5);
        c[i]=t%10^48;
		t/=10;
    }
    int ans=0;
    if(t!=0){
    	ans=ans*10+t;
	}
    for(int i=m;i>=0;--i){
    	ans=ans*10+(c[i]-'0');
	}
    n=strlen(f)-1;
	m=strlen(g)-1;
    int a=0,b=0;
    for(int i=0;i<=n;++i){
    	a=a*10+(f[i]-'0');
	}
    for(int i=0;i<=m;++i){
    	b=b*10+(g[i]-'0');
	}
	cout<<ans-a*b+a+b;
    return 0;
}