p1919 A*B Problem升级版

传送门

题目

给出两个n位10进制整数x和y，你需要计算x*y。

输入格式：

第一行一个正整数n。 第二行描述一个位数为n的正整数x。 第三行描述一个位数为n的正整数y。

输出格式：

输出一行，即x*y的结果。（注意判断前导0）

数据范围：

n<=60000

分析

将数字倒序读入，即个位数是1次方系数，十位数是2次方系数，以此类推。然后进行朴素的FFT，最后答案同样也倒序输出即可，注意输出地数应截止至2次方系数，因为1次方乘1次方得到的是2次方。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<ctime>
#include<vector>
#include<set>
#include<map>
#include<stack>
using namespace std;
#define pi acos(-1.0)
#define ri register int
int n,m,len,r[];
inline int read(){
      int x=,f=;char s=getchar();
      while(s<''||s>''){if(s=='-')f=-;s=getchar();}
      while(s>=''&&s<=''){x=(x<<)+(x<<)+(s-'');s=getchar();}
      return x*f;
}
inline int read2(){
      int f=;char s=getchar();
      while(s<''||s>''){if(s=='-')f=-;s=getchar();}
      return (s-'')*f;
}
struct node {
      double x,y;
      node(){};
      node(double a,double b){x=a;y=b;}
}a[],b[];
inline node operator + (const node a,const node b){
      return node(a.x+b.x,a.y+b.y);
}
inline node operator - (const node a,const node b){
      return node(a.x-b.x,a.y-b.y);
}
inline node operator * (const node a,const node b){
      return node(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);
}
inline void fft(node *a,int f){
      for(ri i=;i<n;i++)
        if(i<r[i])
          swap(a[i],a[r[i]]);
      for(ri k=;k<n;k<<=){
          node wn(cos(pi/k),sin(pi/k)*f);
          for(ri i=;i<n;i+=(k<<)){
            node w(,),p,q;
            for(ri j=;j<k;j++,w=w*wn){
                p=a[i+j],q=a[i+j+k]*w;
                a[i+j]=p+q,a[i+j+k]=p-q;
            }
        }
      }
      if(f==-)
        for(ri i=;i<n;i++)
          a[i].x=a[i].x/n;
}
int ans[];
int main()
{     n=read();m=n;
      a[].x=b[].x=;
      for(ri i=n;i>;i--)a[i].x=read2();
      for(ri i=m;i>;i--)b[i].x=read2();
      m+=n;
      for(n=;n<=m;n<<=)len++;
      for(ri i=;i<n;i++)r[i]=((r[i>>]>>)|((i&)<<(len-)));
      fft(a,),fft(b,);
      for(ri i=;i<n;i++)a[i]=a[i]*b[i];
      fft(a,-);
      for(ri i=;i<=m;i++)ans[i]=(int)(a[i].x+0.5);
      for(ri i=;i<=m;i++)
         if(ans[i]>=){
           ans[i+]+=(ans[i]/);
           ans[i]%=;
         }
      int ok=;
      for(ri i=m+;i>;i--){
         if(ans[i])ok=;
         if(ok)printf("%d",ans[i]);
      }
      return ;
}