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

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

#include
#include
#define Maxn 1000500const double Pi=acos(-1);using namespace std;int n,m,r[Maxn<<2];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 b[Maxn<<2],c[Maxn<<2];void fft(complex *f,short op){ for (int i=0;i
>1]>>1)|((i&1)?n>>1:0); fft(b,1); fft(c,1);//DFT for(int i=0;i

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

你可能感兴趣的文章
SQL Server--导入和导出向导
查看>>
python 数据类型
查看>>
05-linux文件属性-硬链接-时间戳
查看>>
2015-2016 ACM-ICPC, Central Europe Regional Contest (CERC 15)
查看>>
malloc 实现二维数组
查看>>
P2661 信息传递
查看>>
[HDU] 1025 Constructing Roads In JGShining's Kingdom - 二分的求最大递增非连续子序列
查看>>
decode函数
查看>>
SCP注意事项
查看>>
英国NHS
查看>>
ScrollView嵌套GridView和ListView行高问题
查看>>
测试秒杀新版本3.5 stieserver cms
查看>>
Lua获取当前时间
查看>>
redis5.0主从配置
查看>>
JavaScript严谨模式(Strict Mode)提升开发效率和质量
查看>>
[洛谷P4092][HEOI2016/TJOI2016]树
查看>>
nginx配置比较杂乱的总结
查看>>
docker 真实---安装基本映像 (一)
查看>>
php中使用array_slice将数组中的元素分类
查看>>
关于C#的partial修饰符
查看>>