博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
洛谷 3803 【模板】多项式乘法(FFT)
阅读量:6657 次
发布时间:2019-06-25

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

题目:

第一道FFT!

就是把系数转化为2*n个点值,点值相乘一下,再转化回2*n个系数的过程。

转化为点值的过程就是倍增一样,第一步是w_{1,0},也就是说x都是1,所以一开始2*n个位置上的点值都是原来的系数;然后变成两个一组取w_{2,0},w_{2,1}的点值,最后变成2*n个分别取w_{2*n,0},w_{2*n,1},......,w_{2*n,2*n-1}的点值。过程就是DFT,证明可见上面博客。

从点值转化回系数的方法和DFT差不多,似乎只要把 x 都变成倒数、做一边刚才的就行。变成倒数的方法就是那个Wn方向变成负的,这样 x^k 就是倒着转的,上面那个角标就一直是和原来相反的了。

最后算答案的时候似乎直接把虚数的部分舍弃了。

对于iDFT的证明还有点不明白。

#include
#include
#include
#include
#include
#define db doubleusing namespace std;const int N=1e6+5;const db pi=acos(-1.0);int n,m,len,r[N<<2];//<<2! for (n+m)<<1struct cpl{ db x,y;}I,a[N<<2],b[N<<2];cpl operator+ (cpl a,cpl b){
return (cpl){a.x+b.x,a.y+b.y};}cpl operator- (cpl a,cpl b){
return (cpl){a.x-b.x,a.y-b.y};}cpl operator* (cpl a,cpl b){
return (cpl){a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x};}int rdn(){ int ret=0;bool fx=1;char ch=getchar(); while(ch>'9'||ch<'0'){
if(ch=='-')fx=0;ch=getchar();} while(ch>='0'&&ch<='9') ret=(ret<<3)+(ret<<1)+ch-'0',ch=getchar(); return fx?ret:-ret;}void fft(cpl *a,bool fx){ for(int i=0;i
>1; cpl Wn=(cpl){ cos(pi/m),(fx?-1:1)*sin(pi/m) }; for(int i=0;i
>1]>>1)+((i&1)?len>>1:0); fft(a,0); fft(b,0); for(int i=0;i

 

转载于:https://www.cnblogs.com/Narh/p/10019574.html

你可能感兴趣的文章
matlab中如何求某一个矩阵的标准差和均值
查看>>
Leetcode: Find Peak Element
查看>>
Android 抓包,监控流量工具之 mitmproxy
查看>>
hOAuth2.0认证和授权原理
查看>>
Leetcode: Sentence Screen Fitting
查看>>
How To PLAY_SOUND in Oracle Forms
查看>>
atitit.提升开发效率---使用server控件生命周期 asp.net 11个阶段 java jsf 的6个阶段比較...
查看>>
资源下载南方cass视频教程,包括文档,数据,很全的
查看>>
CoreAudio实现录音播音和扬声器听筒模式的切换
查看>>
CSV
查看>>
Atiti.ui原理与gui理论
查看>>
六:二叉树中第k层节点个数与二叉树叶子节点个数
查看>>
编写可测试的代码
查看>>
jQuery中的编程范式
查看>>
Mac OS X 下安装MySQL 5.7
查看>>
爆料!传阿里人靠脸吃饭……
查看>>
马云亲口确认明年退休:转型当老师,张勇将接任
查看>>
布局产业互联网下半场,腾讯云开发者平台扮演何种角色?
查看>>
华为Mate 10来袭 适配HUAWEI VR 2带来全新视听体验
查看>>
乐视网今日股票开盘大涨,5分钟密集成交2亿!
查看>>