高精度乘法 普通(n^2)+fft(nlogn)

高精度乘法核心为

ci=j=1iajbij+1

普通算法时间复杂度为 O(n2) .
又由于是卷积形式,可用fft优化为 O(nlogn)

普通版

#include
#include
#include
#include
#include
using namespace std;

const int N=405;
int n,m,a[N],b[N],c[N];
char s[N];

int main()
{
    //freopen("lx.in","r",stdin);

    scanf("%s",s+1);
    n=strlen(s+1);
    reverse(s+1,s+n+1);
    for(int i=1;i<=n;i++)a[i]=s[i]-'0';

    scanf("%s",s+1);
    m=strlen(s+1);
    reverse(s+1,s+m+1);
    for(int i=1;i<=m;i++)b[i]=s[i]-'0';

    for(int i=1;i<=m+n;i++)
        for(int j=1;j<=n&&i-j+1>0;j++)
            c[i]+=a[j]*b[i-j+1];
    for(int i=1;i<=m+n;i++)
        c[i+1]+=c[i]/10,c[i]%=10;
    int len=m+n;
    while(!c[len])len--;
    for(int i=len;i>=1;i--)
        cout<return 0;
}

FFT版

#include
#include
#include
#include
#include
#include
#include
using namespace std;

const int N=513;
const double PI=acos(-1);
int n,m,pos[N],ans[N];
char s[N];
complex<double>f1[N],f2[N];

void rev(int k)
{
    for(int i=1;i1)?(pos[i>>1]>>1)^(k>>1):pos[i>>1]>>1;
}

void fft(complex<double>f[],int len,int on)
{
    for(int i=0;iif(ifor(int i=1;i1)
    {
        complex<double>wn(cos(on*PI/i),sin(on*PI/i));
        for(int j=0;j1))
        {
            complex<double>wi(1,0);
            for(int k=j;kcomplex<double>u=f[k],v=f[k+i]*wi;
                f[k]=u+v;
                f[k+i]=u-v;
                wi*=wn;
            }
        }
    }

    if(on==-1)
    {
        for(int i=0;ivoid multi(complex<double>f1[],complex<double>f2[])
{
    int len=1;
    while(len1;

    rev(len);

    fft(f1,len,1),fft(f2,len,1);
    for(int i=0;i1);

    for(int i=0;iint(f1[i].real()+0.5);
    for(int i=0;iif(ans[i]>=10)ans[i+1]+=ans[i]/10,ans[i]%=10;
    while(!ans[len])len--;
    for(int i=len;i>=0;i--)cout<int main()
{
    //freopen("lx.in","r",stdin);

    scanf("%s",s);
    n=strlen(s);
    reverse(s,s+n);
    for(int i=0;i'0';

    scanf("%s",s);
    m=strlen(s);
    reverse(s,s+m);
    for(int i=0;i'0';

    multi(f1,f2);

    return 0;
}

你可能感兴趣的:(高精度算法,多项式运算)