BZOJ 2179 [快速傅里叶变换高精度乘法]

2179: FFT快速傅立叶

Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 3108 Solved: 1599
[Submit][Status][Discuss]

Description

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

Input

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

数据范围：
n<=60000

扔个模板

注意读入字符转换成系数系数转换成整数

#include <iostream>

#include <cstdio>

#include <string>

#include <algorithm>

#include <cmath>

using namespace std;

const int N=3e5+;

inline int read(){

    char c=getchar();int x=,f=;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}

    return x*f;

}

const double PI=acos(-);

struct Vector{

    double x,y;

    Vector(double a=,double b=):x(a),y(b){}

};

typedef Vector CD;

Vector operator +(Vector a,Vector b){return Vector(a.x+b.x,a.y+b.y);}

Vector operator -(Vector a,Vector b){return Vector(a.x-b.x,a.y-b.y);}

Vector operator *(Vector a,Vector b){return Vector(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}

Vector conj(Vector a){return Vector(a.x,-a.y);}

struct FastFourierTransform{

    int n,rev[N];

    CD omega[N],omegaInv[N];

    void ini(int m){

        n=;

        while(n<m) n<<=;

        for(int k=;k<n;k++)

            omega[k]=CD(cos(*PI/n*k),sin(*PI/n*k)),

            omegaInv[k]=conj(omega[k]);

        int k=;

        while((<<k)<n) k++;

        for(int i=;i<n;i++){

            int t=;

            for(int j=;j<k;j++) if(i&(<<j)) t|=(<<(k-j-));

            rev[i]=t;

        }

    }

    void transform(CD *a,CD *omega){

        for(int i=;i<n;i++) if(i<rev[i]) swap(a[i],a[rev[i]]);

        for(int l=;l<=n;l<<=){

            int m=l>>;

            for(CD *p=a;p!=a+n;p+=l)

                for(int k=;k<m;k++){

                    CD t=omega[n/l*k]*p[k+m];

                    p[k+m]=p[k]-t;

                    p[k]=p[k]+t;

                }

        }

    }

    void DFT(CD *a,int flag){

        if(flag==) transform(a,omega);

        else{

            transform(a,omegaInv);

            for(int i=;i<n;i++) a[i].x/=(double)n;

        }

    }

    void FFT(CD *a,CD *b,int m){

        ini(m);

        DFT(a,);DFT(b,);

        for(int i=;i<n;i++) a[i]=a[i]*b[i];

        DFT(a,-);

    }

}fft;

CD A[N],B[N];

int n,m,c[N];

char s1[N],s2[N];

int main(){

    freopen("in","r",stdin);

    n=read();m=n+n-;

    scanf("%s%s",s1,s2);

    for(int i=;i<n;i++) A[i].x=s1[n-i-]-'',B[i].x=s2[n-i-]-'';

    fft.FFT(A,B,m);

    for(int i=;i<m;i++) c[i]=int(A[i].x+0.5);//printf("c %d\n",c[i]);

    for(int i=;i<m;i++) c[i+]+=c[i]/,c[i]%=;

    while(c[m]) m++;

    for(int i=m-;i>=;i--) printf("%d",c[i]);

    return ;

}

FFT 1880ms

#include <iostream>

#include <cstdio>

#include <string>

#include <algorithm>

#include <cmath>

using namespace std;

const int N=3e5+;

inline int read(){

    char c=getchar();int x=,f=;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}

    return x*f;

}

const double PI=acos(-);

struct Vector{

    double x,y;

    Vector(double a=,double b=):x(a),y(b){}

};

typedef Vector CD;

Vector operator +(Vector a,Vector b){return Vector(a.x+b.x,a.y+b.y);}

Vector operator -(Vector a,Vector b){return Vector(a.x-b.x,a.y-b.y);}

Vector operator *(Vector a,Vector b){return Vector(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}

struct FastFourierTransform{

    int n,rev[N];

    void ini(int m){

        n=;

        while(n<m) n<<=;

        int k=;

        while((<<k)<n) k++;

        for(int i=;i<n;i++){

            int t=;

            for(int j=;j<k;j++) if(i&(<<j)) t|=(<<(k-j-));

            rev[i]=t;

        }

    }

    void DFT(CD *a,int flag){

        for(int i=;i<n;i++) if(i<rev[i]) swap(a[i],a[rev[i]]);

        for(int l=;l<=n;l<<=){

            int m=l>>;

            CD wn(cos(*PI/l),flag*sin(*PI/l));

            for(CD *p=a;p!=a+n;p+=l){

                CD w(,);

                for(int k=;k<m;k++){

                    CD t=w*p[k+m];

                    p[k+m]=p[k]-t;

                    p[k]=p[k]+t;

                    w=w*wn;

                }

            }

        }

        if(flag==-) for(int i=;i<n;i++) a[i].x/=n;

    }

    void FFT(CD *a,CD *b,int m){

        ini(m);

        DFT(a,);DFT(b,);

        for(int i=;i<n;i++) a[i]=a[i]*b[i];

        DFT(a,-);

    }

}fft;

CD A[N],B[N];

int n,m,c[N];

char s1[N],s2[N];

int main(){

    freopen("in","r",stdin);

    n=read();m=n+n-;

    scanf("%s%s",s1,s2);

    for(int i=;i<n;i++) A[i].x=s1[n-i-]-'',B[i].x=s2[n-i-]-'';

    fft.FFT(A,B,m);

    for(int i=;i<m;i++) c[i]=int(A[i].x+0.5);//printf("c %d\n",c[i]);

    for(int i=;i<m;i++) c[i+]+=c[i]/,c[i]%=;

    while(c[m]) m++;

    for(int i=m-;i>=;i--) printf("%d",c[i]);

    return ;

}

递推w的方法 FFT 1260ms

当然了，用NNT也可以，然而输给了常数

#include <iostream>

#include <cstdio>

#include <string>

#include <algorithm>

#include <cmath>

using namespace std;

typedef long long ll;

const int N=3e5+;

inline int read(){

    char c=getchar();int x=,f=;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}

    return x*f;

}

ll P=,MOD=P;

ll Pow(ll a,ll b,ll MOD){

    ll ans=;

    for(;b;b>>=,a=a*a%MOD)

        if(b&) ans=ans*a%MOD;

    return ans;

}

struct NumberTheoreticTransform{

    int n,rev[N];

    ll g;

    void ini(int m){

        n=;

        while(n<m) n<<=;

        int k=;

        while((<<k)<n) k++;

        for(int i=;i<n;i++){

            int t=;

            for(int j=;j<k;j++) if(i&(<<j)) t|=(<<(k-j-));

            rev[i]=t;

        }

        g=;

    }

    void DFT(ll *a,int flag){

        for(int i=;i<n;i++) if(i<rev[i]) swap(a[i],a[rev[i]]);

        for(int l=;l<=n;l<<=){

            int m=l>>;

            ll wn=Pow(g,flag==?(P-)/l:P--(P-)/l,P);

            for(ll *p=a;p!=a+n;p+=l){

                ll w=;

                for(int k=;k<m;k++){

                    ll t=w*p[k+m]%P;

                    p[k+m]=(p[k]-t+P)%P;

                    p[k]=(p[k]+t)%P;

                    w=w*wn%P;

                }

            }

        }

        if(flag==-){

            ll inv=Pow(n,P-,P);;

            for(int i=;i<n;i++) a[i]=a[i]*inv%P;

        }

    }

    void MUL(ll *A,ll *B){

        DFT(A,);DFT(B,);

        for(int i=;i<n;i++) A[i]=A[i]*B[i]%MOD;

        DFT(A,-);

    }

}fft;

int n,m,c[N];

char s1[N],s2[N];

ll A[N],B[N];

int main(){

    freopen("in","r",stdin);

    n=read();m=n+n-;

    scanf("%s%s",s1,s2);

    for(int i=;i<n;i++) A[i]=s1[n-i-]-'',B[i]=s2[n-i-]-'';

    fft.ini(m);

    fft.MUL(A,B);

    for(int i=;i<m;i++) c[i]=A[i];//printf("c %d\n",c[i]);

    for(int i=;i<m;i++) c[i+]+=c[i]/,c[i]%=;

    while(c[m]) m++;

    for(int i=m-;i>=;i--) printf("%d",c[i]);

}

NNT 3728ms

BZOJ 2179 [快速傅里叶变换高精度乘法]的更多相关文章

BZOJ 2194 [快速傅里叶变换卷积]
题意:请计算C[k]=sigma(a[i]*b[i-k]) 其中 k < = i < n ,并且有 n < = 10 ^ 5. a,b中的元素均为小于等于100的非负整数. 卷积 ( ...
BZOJ 2179 FFT快速傅立叶题解
bzoj 2179 Description 给出两个n位10进制整数x和y,你需要计算x*y. [题目分析] 高精裸题.练手. [代码] 1.手动高精 #include<cstdio> # ...
BZOJ 2179 FFT快速傅立叶 ——FFT
[题目分析] 快速傅里叶变换用于高精度乘法. 其实本质就是循环卷积的计算,也就是多项式的乘法. 两次蝴蝶变换. 二进制取反化递归为迭代. 单位根的巧妙取值,是的复杂度成为了nlogn 范德蒙矩阵计算逆 ...
【POJ 1001】Exponentiation (高精度乘法+快速幂)
BUPT2017 wintertraining(15) #6A 题意求$R^n$ ( 0.0 < R < 99.999 )(0 < n <= 25) 题解将R用字符串读 ...
Algorithm: 多项式乘法 Polynomial Multiplication: 快速傅里叶变换 FFT / 快速数论变换 NTT
Intro: 本篇博客将会从朴素乘法讲起,经过分治乘法,到达FFT和NTT 旨在能够让读者(也让自己)充分理解其思想模板题入口:洛谷 P3803 [模板]多项式乘法(FFT) 朴素乘法约定:两个多 ...
[学习笔记] 多项式与快速傅里叶变换(FFT)基础
引入可能有不少OIer都知道FFT这个神奇的算法, 通过一系列玄学的变化就可以在 $O(nlog(n))$ 的总时间复杂度内计算出两个向量的卷积, 而代码量却非常小. 博主一年半前曾经因COGS的一 ...
快速傅里叶变换FFT& 数论变换NTT
相关知识时间域上的函数f(t)经过傅里叶变换(Fourier Transform)变成频率域上的F(w),也就是用一些不同频率正弦曲线的加权叠加得到时间域上的信号. \[ F(\omega)=\m ...
多项式之快速傅里叶变换(FFT)/数论变换(NTT)/常用套路【入门】
原文链接https://www.cnblogs.com/zhouzhendong/p/Fast-Fourier-Transform.html 多项式之快速傅里叶变换(FFT)/数论变换(NTT)/ ...
HDU 1402 A * B Problem Plus 快速傅里叶变换 FFT 多项式
http://acm.hdu.edu.cn/showproblem.php?pid=1402 快速傅里叶变换优化的高精度乘法. https://blog.csdn.net/ggn_2015/artic ...

随机推荐

SQL强化(二) 在Oracle 中写代码
一 : 关于查询中的转换 -- 字符串转换一 : decode 函数转换 SELECT DECODE ( PROTYPE.PRO_TYPE_DATE, 'L', '长', 'm', '短', ' ...
kafka producer生产数据到kafka异常：Got error produce response with correlation id 16 on topic-partition...Error: NETWORK_EXCEPTION
kafka producer生产数据到kafka异常:Got error produce response with correlation id 16 on topic-partition... ...
腾讯云服务器php+mysq+nginx配置出现的问题及解决方法（亲测）
http://blog.csdn.net/hfdmv/article/details/50900043 删除文件命令 sudo rm -f /usr/share/nginx/html/home.php ...
git gui提交无法获知你的身份 20
刚刚学习,请说的详细一些,谢谢 callct | 浏览 3382 次我有更好的答案 1条回答你没有定义你的名字和邮箱.你打开git console/shell, #输入下面两句,并且替换成你的名字 ...
Python3 的列表
1:列表: Python的列表比C语言的数组强大的多,数组只能存放相同类型的数据,而列表则像一个大集装箱可以存放整形.浮点型.字符串.对象等 2:创建列表的方法 #创建一个普通列表 list1=[1, ...
Apache Commons介绍（转载）
一.Commons BeanUtils说明:针对Bean的一个工具集.由于Bean往往是有一堆get和set组成,所以BeanUtils也是在此基础上进行一些包装. 二.Commons CLI说明:这 ...
io利用率100%问题
iostat -mx 1 dm-60 dm-61 dm-62 dm-63 dm-64 dm-65 dm-66 dm-67 Device: rrqm/s wrqm/s r/s ...
docker 安装NexusRepository Manager
今天学习了一下docker 感觉这东西要学习好多的命令,但是自己又是喜欢这种命令,感觉linux总是高一个等级的东西,这几天学习使用docker安装各种东西,下面记录一些我安装nexus的步骤,还是不 ...
一个自己稍作修改了的美赛论文LaTeX模板
原模板(5.0)来自LaTeX工作室(latexstudio.net),我按照比赛规范做了一点小小的修改(5.0y),并加上了比原来更详细一些的注释,方便使用. 仅仅分享一下方便大家使用,模板的原创者 ...
jdk源码->集合->HashSet
类的属性 public class HashSet<E> extends AbstractSet<E> implements Set<E>, Cloneable, ...

BZOJ 2179 [快速傅里叶变换 高精度乘法]