BZOJ 2179 [快速傅里叶变换 高精度乘法]
2179: FFT快速傅立叶
Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 3108 Solved: 1599
[Submit][Status][Discuss]
Description
Input
n<=60000
#include <iostream> #include <cstdio> #include <string> #include <algorithm> #include <cmath> using namespace std; const int N=3e5+5; inline int read(){char c=getchar();int x=0,f=1;while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x*f; } const double PI=acos(-1); struct Vector{double x,y;Vector(double a=0,double b=0):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=1;while(n<m) n<<=1;for(int k=0;k<n;k++) omega[k]=CD(cos(2*PI/n*k),sin(2*PI/n*k)),omegaInv[k]=conj(omega[k]);int k=0;while((1<<k)<n) k++;for(int i=0;i<n;i++){int t=0;for(int j=0;j<k;j++) if(i&(1<<j)) t|=(1<<(k-j-1));rev[i]=t;}}void transform(CD *a,CD *omega){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;for(CD *p=a;p!=a+n;p+=l)for(int k=0;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==1) transform(a,omega);else{transform(a,omegaInv);for(int i=0;i<n;i++) a[i].x/=(double)n;}}void FFT(CD *a,CD *b,int m){ini(m);DFT(a,1);DFT(b,1);for(int i=0;i<n;i++) a[i]=a[i]*b[i];DFT(a,-1);} }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-1;scanf("%s%s",s1,s2);for(int i=0;i<n;i++) A[i].x=s1[n-i-1]-'0',B[i].x=s2[n-i-1]-'0';fft.FFT(A,B,m);for(int i=0;i<m;i++) c[i]=int(A[i].x+0.5);//printf("c %d\n",c[i]);for(int i=0;i<m;i++) c[i+1]+=c[i]/10,c[i]%=10;while(c[m]) m++;for(int i=m-1;i>=0;i--) printf("%d",c[i]);return 0; }
FFT 1880ms
#include <iostream> #include <cstdio> #include <string> #include <algorithm> #include <cmath> using namespace std; const int N=3e5+5; inline int read(){char c=getchar();int x=0,f=1;while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x*f; } const double PI=acos(-1); struct Vector{double x,y;Vector(double a=0,double b=0):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=1;while(n<m) n<<=1;int k=0;while((1<<k)<n) k++;for(int i=0;i<n;i++){int t=0;for(int j=0;j<k;j++) if(i&(1<<j)) t|=(1<<(k-j-1));rev[i]=t;}}void DFT(CD *a,int flag){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;CD wn(cos(2*PI/l),flag*sin(2*PI/l));for(CD *p=a;p!=a+n;p+=l){CD w(1,0);for(int k=0;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==-1) for(int i=0;i<n;i++) a[i].x/=n;}void FFT(CD *a,CD *b,int m){ini(m);DFT(a,1);DFT(b,1);for(int i=0;i<n;i++) a[i]=a[i]*b[i];DFT(a,-1);} }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-1;scanf("%s%s",s1,s2);for(int i=0;i<n;i++) A[i].x=s1[n-i-1]-'0',B[i].x=s2[n-i-1]-'0';fft.FFT(A,B,m);for(int i=0;i<m;i++) c[i]=int(A[i].x+0.5);//printf("c %d\n",c[i]);for(int i=0;i<m;i++) c[i+1]+=c[i]/10,c[i]%=10;while(c[m]) m++;for(int i=m-1;i>=0;i--) printf("%d",c[i]);return 0; }
递推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+5; inline int read(){char c=getchar();int x=0,f=1;while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x*f; } ll P=1004535809,MOD=P; ll Pow(ll a,ll b,ll MOD){ll ans=1;for(;b;b>>=1,a=a*a%MOD)if(b&1) ans=ans*a%MOD;return ans; } struct NumberTheoreticTransform{int n,rev[N];ll g;void ini(int m){n=1;while(n<m) n<<=1;int k=0;while((1<<k)<n) k++;for(int i=0;i<n;i++){int t=0;for(int j=0;j<k;j++) if(i&(1<<j)) t|=(1<<(k-j-1));rev[i]=t;}g=3;}void DFT(ll *a,int flag){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=Pow(g,flag==1?(P-1)/l:P-1-(P-1)/l,P);for(ll *p=a;p!=a+n;p+=l){ll w=1;for(int k=0;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==-1){ll inv=Pow(n,P-2,P);;for(int i=0;i<n;i++) a[i]=a[i]*inv%P;}}void MUL(ll *A,ll *B){DFT(A,1);DFT(B,1);for(int i=0;i<n;i++) A[i]=A[i]*B[i]%MOD;DFT(A,-1);} }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-1;scanf("%s%s",s1,s2);for(int i=0;i<n;i++) A[i]=s1[n-i-1]-'0',B[i]=s2[n-i-1]-'0';fft.ini(m);fft.MUL(A,B);for(int i=0;i<m;i++) c[i]=A[i];//printf("c %d\n",c[i]);for(int i=0;i<m;i++) c[i+1]+=c[i]/10,c[i]%=10;while(c[m]) m++;for(int i=m-1;i>=0;i--) printf("%d",c[i]); }
NNT 3728ms
转载于:https://www.cnblogs.com/candy99/p/6388403.html
BZOJ 2179 [快速傅里叶变换 高精度乘法]相关推荐
- 使用快速傅里叶变换计算大整数乘法-代码
在上一篇随笔"使用快速傅里叶变换计算大整数乘法"中,已经讲述了使用快速傅里叶变换计算大整数乘法的原理.在这一篇随笔中,我们就使用快速傅里叶变换来实现一个提供任意精度的算术运算的静态 ...
- 十五 多项式乘法与快速傅里叶变换
分享一下我老师大神的人工智能教程!零基础,通俗易懂!http://blog.csdn.net/jiangjunshow 也欢迎大家转载本篇文章.分享知识,造福人民,实现我们中华民族伟大复兴! 十五.多 ...
- 多项式乘法与快速傅里叶变换
前言 经典算法研究系列,已经写到第十五章了,本章,咱们来介绍多项式的乘法以及快速傅里叶变换算法.本博客之前也已详细介绍过离散傅里叶变换(请参考:十.从头到尾彻底理解傅里叶变换算法.上,及十.从头到尾彻 ...
- c语言fft乘法步骤,C语言实现FFT(快速傅里叶变换).doc
C语言实现FFT(快速傅里叶变换) 择蚁牙幸帆揣邓淌港烬粹甩滋整维含兔忿茂慨渔下餐随扼哇房坏鹅穆礼围引介害芝共茨恿把喜恤寇杖除冕嗓停揍猫调锚遭傀个碱晓频斌硕宾撕坪莱哩腊养掘蹄轴国繁蔬虞靡砖焙倍勾呸怀怒 ...
- FFT快速傅里叶变换 超详细的入门学习总结
FFT快速傅里叶变换 说明 本文创作的目的是为自己巩固该算法,加深印象并深入理解,同时也为FFT入门学者提供一份可鉴的学习总结. 原文链接:https://blog.csdn.net/qq_39565 ...
- [转]十分简明易懂的FFT(快速傅里叶变换)
FFT前言 快速傅里叶变换 (fast Fourier transform),即利用计算机计算离散傅里叶变换(DFT)的高效.快速计算方法的统称,简称FFT.快速傅里叶变换是1965年由J.W.库利和 ...
- FFT(快速傅里叶变换)学习笔记
简介 FFTFFTFFT (法法塔)是个什么玩意?他的全名叫快速傅里叶变换(然而貌似和傅里叶并没有太大关系),用来快速求出多项式的点值表示,这个东西一般用来解决多项式相乘的问题. 一般的高精度乘法,我 ...
- 【算法竞赛学习笔记】快速傅里叶变换FFT-数学提高计划
tilte : 快速傅里叶变换FFT学习笔记 tags : ACM,数论 date : 2021-7-18 简介 FFT(Fast Fourier Transformation),中文名快速傅里叶变换 ...
- 超详细易懂FFT(快速傅里叶变换)及代码实现
前言 昨天学了一晚上,终于搞懂了FFT.希望能写一篇清楚易懂的题解分享给大家,也进一步加深自己的理解. FFT算是数论中比较重要的东西,听起来就很高深的亚子.但其实学会了(哪怕并不能完全理解),会实现 ...
最新文章
- Swift 异常处理
- HTTP中Accept与Content-Type区别
- android java和c混合编程_C/C++在Java、Android和Objective-C三大平台下实现混合编程
- 文件夹里面照片自动分成子文件夹_Windows居然自带这个功能,自动整理你硬盘里的照片...
- java多线程之:SynchronousQueue队列
- 微抖动,繁忙的等待和绑定CPU
- 华为入局 VR 眼镜能让 VR 早普及几年?|CSDN博文精选
- NetApp F3020 盘柜报警升级修复全过程
- SQL server 数据导入导出BCP工具使用详解
- 电脑硬件故障排除经验
- 2017 4 自考java_自考00831英语语法2017年4月真题及答案【解析版】
- LeetCode 707. Design Linked List
- 错误:Redis----(error) MISCONF Redis is configured to save RDB snapshots
- 手机usb共享计算机网络连接,手机如何通过USB共享电脑网络
- Python爬虫框架Scrapy入门(三)爬虫实战:爬取长沙链家二手房
- Thinkphp5+JWT开发 api接口
- 华中科技大学2021计算机学院,2021年华中科技大学计算机考研科目
- 大数据与云安全专题-1
- java判断php的emoji,php 过滤emoji
- 数据分析中缺失值的处理方法
热门文章
- [coco2d]pageView:addPage时,page无法对齐
- 转---谈谈HTTP协议中的短轮询、长轮询、长连接和短连接
- Design Patterns in Java
- timestamp 数据类型在 sql_mode 主从不一致引起的不同步问题解决
- RAC测试及命令示例
- (android实战)Service 生命周期和使用注意项
- android removeRule的使用
- Android滑屏 mScrollX mScrollY scrollTo() scrollBy()
- 从open系统调用的源码看文件的打开过程
- Android 自定义View,自定义属性--自定义圆形进度条(整理)