TOT 傅立叶变换 FFT 入门
HDU 1402,计算很大的两个数相乘。
FFT 只要78ms,这里;
一些FFT 入门资料:http://wenku.baidu.com/view/8bfb0bd476a20029bd642d85.html (讲解的很详细
http://blog.csdn.net/iamzky/article/details/22712347 (这个也不错
另外算导的其实也蛮好,只是怕公式的看前面的也可。
IDFT只是FFT的逆变换,这里想了很久原来只要在FFT 变换后的结果后/N 即可,算实数部分即可。
前面的一份模板 :
1 /*
2 algorithm : High-Precision FFT
3
4 */
5 #include <cstdio>
6 #include <cstring>
7 #include <cmath>
8 #include <algorithm>
9 #define N 200005
10 #define pi acos(-1.0) // PI值
11 using namespace std;
12 struct complex
13 {
14 double r,i;
15 complex(double real=0.0,double image=0.0){
16 r=real; i=image;
17 }
18 // 以下为三种虚数运算的定义
19 complex operator + (const complex o){
20 return complex(r+o.r,i+o.i);
21 }
22 complex operator - (const complex o){
23 return complex(r-o.r,i-o.i);
24 }
25 complex operator * (const complex o){
26 return complex(r*o.r-i*o.i,r*o.i+i*o.r);
27 }
28 }x1[N],x2[N];
29 char a[N/],b[N/];
30 int sum[N]; // 结果存在sum里
31 void brc(complex *y,int l) // 二进制平摊反转置换 O(logn)
32 {
33 register int i,j,k;
34 for(i=,j=l/;i<l-;i++)
35 {
36 if(i<j) swap(y[i],y[j]); // 交换互为下标反转的元素
37 // i<j保证只交换一次
38 k=l/;
39 while(j>=k) // 由最高位检索,遇1变0,遇0变1,跳出
40 {
41 j-=k;
42 k/=;
43 }
44 if(j<k) j+=k;
45 }
46 }
47 void fft(complex *y,int l,double on) // FFT O(nlogn)
48 // 其中on==1时为DFT,on==-1为IDFT
49 {
50 register int h,i,j,k;
51 complex u,t;
52 brc(y,l); // 调用反转置换
53 for(h=;h<=l;h<<=) // 控制层数
54 {
55 // 初始化单位复根
56 complex wn(cos(on**pi/h),sin(on**pi/h));
57 for(j=;j<l;j+=h) // 控制起始下标
58 {
59 complex w(,); // 初始化螺旋因子
60 for(k=j;k<j+h/;k++) // 配对
61 {
62 u=y[k];
63 t=w*y[k+h/];
64 y[k]=u+t;
65 y[k+h/]=u-t;
66 w=w*wn; // 更新螺旋因子
67 } // 据说上面的操作叫蝴蝶操作…
68 }
69 }
70 if(on==-) for(i=;i<l;i++) y[i].r/=l; // IDFT
71 }
72 int main(void)
73 {
74 int l1,l2,l;
75 register int i;
76 while(scanf("%s%s",a,b)!=EOF)
77 {
78 l1=strlen(a);
79 l2=strlen(b);
80 l=;
81 while(l<l1* || l<l2*) l<<=; // 将次数界变成2^n
82 // 配合二分与反转置换
83 for(i=;i<l1;i++) // 倒置存入
84 {
85 x1[i].r=a[l1-i-]-'';
86 x1[i].i=0.0;
87 }
88 for(;i<l;i++) x1[i].r=x1[i].i=0.0;
89 // 将多余次数界初始化为0
90 for(i=;i<l2;i++)
91 {
92 x2[i].r=b[l2-i-]-'';
93 x2[i].i=0.0;
94 }
95 for(;i<l;i++) x2[i].r=x2[i].i=0.0;
96 fft(x1,l,); // DFT(a)
97 fft(x2,l,); // DFT(b)
98 for(i=;i<l;i++) x1[i]=x1[i]*x2[i]; // 点乘结果存入a
99 fft(x1,l,-); // IDFT(a*b)
for(i=;i<l;i++) sum[i]=x1[i].r+0.5; // 四舍五入
for(i=;i<l;i++) // 进位
{
sum[i+]+=sum[i]/;
sum[i]%=;
}
l=l1+l2-;
while(sum[l]<= && l>) l--; // 检索最高位
for(i=l;i>=;i--) putchar(sum[i]+''); // 倒序输出
putchar('\n');
}
return ;
112 }
UOJ#34
多项式也是常用FFT的,似乎代码可以更短
1 #include <cstdio>
2 #include <cstring>
3 #include <cmath>
4 #include <algorithm>
5 #define N 600005
6 #define pi acos(-1.0) // PI值
7 using namespace std;
8 struct complex
9 {
10 double r,i;
11 complex(double real=0.0,double image=0.0){
12 r=real; i=image;
13 }
14 // 以下为三种虚数运算的定义
15 complex operator + (const complex o){
16 return complex(r+o.r,i+o.i);
17 }
18 complex operator - (const complex o){
19 return complex(r-o.r,i-o.i);
20 }
21 complex operator * (const complex o){
22 return complex(r*o.r-i*o.i,r*o.i+i*o.r);
23 }
24 }x1[N],x2[N];
25 char a[N/],b[N/];
26 int sum[N]; // 结果存在sum里
27 void brc(complex *y,int l) // 二进制平摊反转置换 O(logn)
28 {
29 register int i,j,k;
30 for(i=,j=l/;i<l-;i++)
31 {
32 if(i<j) swap(y[i],y[j]); // 交换互为下标反转的元素
33 // i<j保证只交换一次
34 k=l/;
35 while(j>=k) // 由最高位检索,遇1变0,遇0变1,跳出
36 {
37 j-=k;
38 k/=;
39 }
40 if(j<k) j+=k;
41 }
42 }
43
44
45 void fft(complex *y,int l,double on) // FFT O(nlogn)
46 // 其中on==1时为DFT,on==-1为IDFT
47 {
48 register int h,i,j,k;
49 complex u,t;
50 brc(y,l); // 调用反转置换
51 for(h=;h<=l;h<<=) // 控制层数
52 {
53 // 初始化单位复根
54 complex wn(cos(on**pi/h),sin(on**pi/h));
55 for(j=;j<l;j+=h) // 控制起始下标
56 {
57 complex w(,); // 初始化螺旋因子
58 for(k=j;k<j+h/;k++) // 配对
59 {
60 u=y[k];
61 t=w*y[k+h/];
62 y[k]=u+t;
63 y[k+h/]=u-t;
64 w=w*wn; // 更新螺旋因子
65 } // 据说上面的操作叫蝴蝶操作…
66 }
67 }
68 if(on==-) for(i=;i<l;i++) y[i].r/=l; // IDFT
69 }
70
71 int main(void)
72 {
73 int n,m,i;
74 scanf("%d",&n);
75 scanf("%d",&m);
76 n++;
77 m++;
78 int l=;
79 while(l<n* || l<m*) l<<=; // 将次数界变成2^n
80 // 配合二分与反转置换
81 for(i=;i<n;i++) // 倒置存入
82 {
83 scanf("%lf",&x1[i].r);
84 x1[i].i=0.0;
85 }
86 for(;i<l;i++) x1[i].r=x1[i].i=0.0;
87 // 将多余次数界初始化为0
88 for(i=;i<m;i++)
89 {
90
91 scanf("%lf",&x2[i].r);
92 x2[i].i=0.0;
93 }
94
95 //for (int i=0;i<m/2;i++) swap(x1[i],x1[m-i-1]);
96 for(;i<l;i++) x2[i].r=x2[i].i=0.0;
97
98 fft(x1,l,); // DFT(a)
99 fft(x2,l,); // DFT(b)
for(i=;i<l;i++) x1[i]=x1[i]*x2[i]; // 点乘结果存入a
fft(x1,l,-); // IDFT(a*b)
for(i=;i<l;i++) sum[i]=x1[i].r+0.5; // 四舍五入
/*
for(i=0;i<l;i++) // 进位
{
sum[i+1]+=sum[i]/10;
sum[i]%=10;
}
*/
l=n+m-;
//cout<<l<endl;
// while(sum[l]<=0 && l>0) l--; // 检索最高位
for(i=;i<l-;i++) printf("%d ",sum[i]); // 倒序输出
printf("%d\n",sum[l-]);
return ;
116 }
BZOJ 2194
http://hzwer.com/6902.html 参考这篇blog
卷积啥的现在还是晕晕的。
BZOJ 3527
有了前面的卷积,那么这道题就是化简了
http://www.cnblogs.com/iwtwiioi/p/4126284.html
http://blog.csdn.net/regina8023/article/details/44910225
TOT 傅立叶变换 FFT 入门的更多相关文章
- 快速傅立叶变换(FFT)算法
已知多项式f(x)=a0+a1x+a2x2+...+am-1xm-1, g(x)=b0+b1x+b2x2+...+bn-1xn-1.利用卷积的蛮力算法,得到h(x)=f(x)g(x),这一过程的时间复 ...
- 傅立叶变换—FFT
FFT(快速傅立叶变换)使用“分而治之”的策略来计算一个n阶多项式的n阶DFT系数的值.定义n为2的整数幂数,为了计算一个n阶多项式f(x),算法定义了连个新的n/2阶多项式,函数f[0](x)包含了 ...
- 快速傅立叶变换FFT模板
递归版 UOJ34多项式乘法 //容易暴栈,但是很好理解 #include <cmath> #include <iostream> #include <cstdio> ...
- 傅立叶变换—FFT(cuda实现)
背景: 无意间看到cuda解决FFT有一个cufft函数库,大体查看了有关cufft有关知识,写了一个解决一维情况的cuda代码,据调查知道cufft在解决1D,2D,3D的情况时间复杂度都为O(nl ...
- 为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换
写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,绝大部分内容非我所原创.在此向多位原创作者致敬!!!一.傅立叶变换的由来关于傅立叶变换,无论是书本还是在网上可以很容易找到关于傅立叶 ...
- Matlab图像处理系列4———傅立叶变换和反变换的图像
注意:这一系列实验的图像处理程序,使用Matlab实现最重要的图像处理算法 1.Fourier兑换 (1)频域增强 除了在空间域内能够加工处理图像以外,我们还能够将图像变换到其它空间后进行处理.这些方 ...
- Matlab图像处理系列4———图像傅立叶变换与反变换
注:本系列来自于图像处理课程实验.用Matlab实现最主要的图像处理算法 1.Fourier变换 (1)频域增强 除了在空间域内能够加工处理图像以外.我们还能够将图像变换到其它空间后进行处理.这些方法 ...
- 离散傅立叶变换与快速傅立叶变换(DFT与FFT)
自从去年下半年接触三维重构以来,听得最多的词就是傅立叶变换,后来了解到这个变换在图像处理里面也是重点中的重点. 本身自己基于高数知识的理解是傅立叶变换是将一个函数变为一堆正余弦函数的和的变换.而图像处 ...
- 傅立叶变换系列(五)快速傅立叶变换(FFT)
说明: 傅里叶级数.傅里叶变换.离散傅里叶变换.短时傅里叶变换...这些理解和应用都非常难,网上的文章有两个极端:“Esay” Or “Boring”!如果单独看一两篇文章就弄懂傅里叶,那说明你真 ...
随机推荐
- C#值类型和引用类型与Equals方法
1. C#的值类型和引用类型 C#的对象里面有两种类型,一个是引用类型,一个是值类型,值类型和引用类型的具体分类可以看下面的分类. 在C#中,不管是引用类型还是值类型,他们都隐式继承Object类 ...
- ios开发--常用的高效开发的宏
本次在做项目的时候使用了下面的一些宏定义 以及 建立宏定义的一些规则.虽然只用了其中的一点点,但是还是极大的提高了开发效率.. 将这些宏放到一个头文件里然后再放到工程中,在需要使用这些宏定义的地方体检 ...
- js截屏
<html><head> <meta name="layout" content="main"> <meta http ...
- C#图形学习笔记
绘图常用控件.类和结构 颜色 使用System.Drawing.Color结构表示 设置颜色的方法 调用静态函数:Color.FromArgb() public static Color FromAr ...
- tkinter学习-滚动条
阅读目录 Listbox 以列表的形式显示 Scrollbar 滚动条 Scale 更滚动条很相似,但更精准 Listbox: 说明:列表框控件,在Listbox窗口小部件是用来显示一个字符串列表给 ...
- python项目开发视频
精品Python项目开发学习视频 所属网站分类: 资源下载 > python视频教程 作者:乐天派 链接:http://www.pythonheidong.com/blog/article/44 ...
- Python三元表达式和列表生成式
三元表达式 取代 if …… else……的简单表达方式 # 常规写法 x = 1 y = 2 if x>y: print(x) else: print(y) #三元表达式写法 res ...
- mysql 面安装配置
解压MySQL压缩包 将以下载的MySQL压缩包解压到自定义目录下,我的解压目录是: "D:\Program Files\MySQL\mysql-5.6.13-win32&quo ...
- luogu2805 [NOI2009]植物大战僵尸
想象一下,要搞掉一个植物,必须先搞掉另一些植物--我们可以发现这是一个最大权闭合子图的问题. 最大权闭合子图的话,太空飞行计划问题是一个入门题,可以一看. 然而我们手玩一下样例就会惊恐地发现,保护关系 ...
- 自定义iOS上双击Home键图切换
如果双击Home,会来到iOS App的switcher页面,在这儿列出了当前系统挂起的App, 上面有每个App的切屏,相信大家都熟悉这个东东了.它其实是每个App在挂起前,对App后个载屏. 那么 ...