【NTT】hdu1402 A * B Problem Plus

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
#define N ((1<<18)+5)
#define MOD 1004535809ll
ll Quick_Pow(ll a,ll p){
	if(p==0){
		return 1ll;
	}
    ll res=Quick_Pow(a,p>>1);
    res=res*res%MOD;
    if((p&1ll)==1ll){
    	res=(a%MOD*res)%MOD;
    }
    return res;
}
struct NTT{
	int n,rev[N];
	ll g;
	void ini(int lim) {
		g=3;//1004535809,998244353的原根都是3
		n=1;
		int k=0;
		while(n<lim){
			n<<=1;
			++k;
		}
		for(int i=0;i<n;++i){
			rev[i]=((rev[i>>1]>>1)|((i&1)<<(k-1)));
		}
	}
	void dft(ll a[],int DFT) {
		for(int i=0;i<n;++i){
			if(i<rev[i]){
				swap(a[i],a[rev[i]]);
			}
		}
		for(int l=2;l<=n;l<<=1){
			int m=l>>1;
			ll wn=Quick_Pow(g,DFT==1 ? (MOD-1ll)/(ll)l : MOD-1ll-(MOD-1ll)/(ll)l);
			for(int i=0;i<n;i+=l){
				ll w=1;
				for(int k=0;k<m;++k){
					ll t=w*a[i+k+m]%MOD;
					a[i+k+m]=(a[i+k]-t+MOD)%MOD;
					a[i+k]=(a[i+k]+t)%MOD;
					w=w*wn%MOD;
				}
			}
		}
		if(DFT==-1){
			ll inv=Quick_Pow(n,MOD-2ll);
			for(int i=0;i<n;++i){
				a[i]=a[i]*inv%MOD;
			}
		}
	}
	void mul(ll a[],ll b[],int len) {
		ini(len);
		dft(a,1);
		dft(b,1);
		for(int i=0;i<n;++i){
			a[i]=a[i]*b[i];
		}
		dft(a,-1);
	}
}ntt;
int len1,len2,len,c[N];
ll a[N],b[N];
char s1[N],s2[N];
int main() {
//	freopen("ntt.in","r",stdin);
	while(scanf("%s%s",s1,s2)!=EOF){
		memset(c,0,sizeof(c));
		memset(a,0,sizeof(a));
		memset(b,0,sizeof(b));
		len1=strlen(s1);
		len2=strlen(s2);
		for(int i=0;i<len1;++i){
			a[i]=s1[len1-i-1]-‘0‘;
		}
		for(int i=0;i<len2;++i){
			b[i]=s2[len2-i-1]-‘0‘;
		}
		len=len1+len2-1;
		ntt.mul(a,b,len);
		for(int i=0;i<len;++i){
			c[i]=a[i];
		}
		for(int i=0;i<len;++i){
			c[i+1]+=c[i]/10;
			c[i]%=10;
		}
//		if(c[len]){
//			++len;
//		}//两个数乘积的长度要么是A+B-1，要么是A+B。
//		for(int i=len-1;i>=0;--i){
//			printf("%d",c[i]);
//		}
//		puts("");
		for(int i=len;i>=0;--i){
			if(c[i]!=0 || i==0){
				for(int j=i;j>=0;--j){
					printf("%d",c[j]);
				}
				puts("");
				break;
			}
		}
	}
	return 0;
}