博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
快速傅里叶变换(FFT)递归
阅读量:4500 次
发布时间:2019-06-08

本文共 934 字,大约阅读时间需要 3 分钟。

#include 
#include
#define Maxn 5000#define Pi 3.1415926535898using namespace std;int n,m;struct complex{complex (double xx=0,double yy=0){x=xx,y=yy;}double x,y;};complex operator + (complex a,complex b){ return complex(a.x+b.x,a.y+b.y);}complex operator - (complex a,complex b){ return complex(a.x-b.x,a.y-b.y);}complex operator * (complex a,complex b){ return complex(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}complex w[Maxn],b[Maxn],c[Maxn];void fft(complex *f,int len,short op){ if (!len)return ; complex fl[len+1],fr[len+1]; for (int k=0;k
>1,op); fft(fr,len>>1,op); complex tmp,buf; tmp=complex(cos(Pi/len),sin(Pi/len)*op); buf.x=1;buf.y=0; for (int k=0;k
>1,1); fft(c,n>>1,1); for(int i=0;i
>1,-1); for(int i=0;i<=m;++i)printf("%.0f ",fabs(b[i].x)/n); puts(""); return 0;}

转载于:https://www.cnblogs.com/ukcxrtjr/p/11194318.html

你可能感兴趣的文章
TCP/IP模型的一个简单解释
查看>>
解开最后期限的镣铐(转载)
查看>>
Kth Smallest Element in a BST
查看>>
ubuntu14.04利用aliyun安装docker
查看>>
iphone-命令行编译之--xcodebuild
查看>>
Plan : 破晓
查看>>
GNU make
查看>>
Visual Studio 2008 不能更改安装目录的原因
查看>>
threejs学习笔记04---相机动
查看>>
SAP Skill - How to search a field for which table it belongs
查看>>
parcel+vue入门
查看>>
基数排序
查看>>
Dell笔记本刷回低版本bios的方法
查看>>
HLP帮助文件源文件RTF文件的编写
查看>>
2.30模型字符串拷贝
查看>>
XPATH怎么获取TITLE中有中文的标签
查看>>
Tomcat中server.xml参数说明
查看>>
Wget下载终极用法和15个详细的例子
查看>>
JavaScript16进制颜色值和rgb的转换
查看>>
Laravel 输出Hellow World!
查看>>