bzoj2194: 快速傅立叶之二

了解到了FFT求卷积，但是还是感性的认识。。

取反就可以了。输出一定要加int！！！！

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
const double pi=acos(-1.0);struct Complex
{double r,i;Complex(){}Complex(double _r,double _i){r=_r, i=_i;}friend Complex operator +(Complex x,Complex y){return Complex(x.r+y.r,x.i+y.i);}friend Complex operator -(Complex x,Complex y){return Complex(x.r-y.r,x.i-y.i);}friend Complex operator *(Complex x,Complex y){return Complex(x.r*y.r-x.i*y.i,x.r*y.i+x.i*y.r);}
}A[410000],B[410000],C[410000];int R[410000];
void fft(Complex *a,int n,int op)
{for(int i=0;i<n;i++)if(i<R[i])swap(a[i],a[R[i]]);for(int i=1;i<n;i*=2){Complex wn(cos(pi/i),sin(pi*op/i));for(int j=0;j<n;j+=(i<<1)){Complex w(1,0);for(int k=0;k<i;k++,w=w*wn){Complex a1=a[j+k],a2=a[j+k+i];a[j+k]  =a1+w*a2;a[j+k+i]=a1-w*a2;}}}if(op==-1)for(int i=0;i<n;i++)a[i].r/=n;
}
int main()
{int n,m,L=0;scanf("%d",&m);m--;for(int i=0;i<=m;i++)scanf("%lf%lf",&A[i].r,&B[m-i].r);m*=2;for(n=1;n<=m;n*=2)L++;for(int i=0;i<n;i++)R[i]=( R[i>>1]>>1 )|( (i&1) << (L-1) );fft(A,n,1);fft(B,n,1);for(int i=0;i<n;i++)C[i]=A[i]*B[i];fft(C,n,-1);for(int i=m/2;i<=m;i++)printf("%d\n",int(C[i].r+0.5));return 0;
}

转载于:https://www.cnblogs.com/AKCqhzdy/p/7976703.html

bzoj2194: 快速傅立叶之二相关推荐

bzoj2194 快速傅立叶之二
2194: 快速傅立叶之二 Time Limit: 10 Sec Memory Limit: 259 MB Submit: 1730 Solved: 1026 [Submit][Status][D ...
BZOJ2194 快速傅立叶之二【fft】
题目请计算C[k]=sigma(a[i]*b[i-k]) 其中 k < = i < n ,并且有 n < = 10 ^ 5. a,b中的元素均为小于等于100的非负整数. 输入格式 ...
BZOJ-2194 快速傅立叶之二
FFT模版题. 观察题目,我们可以发现,只要把序列b倒过来,再联想一下乘法运算... 我们会发现,将序列a和序列b当作100进制数,做一次乘法,然后从低到高每一位便是答案了(乘完无需进位) #incl ...
解题报告（二）C、（darkBZOJ 2194）快速傅立叶之二（FFT、卷积的概念、常用变换）
繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量题解和代码,题目难度不一定按照题号排序,我会在每道题后面加上题目难度指数(1∼51 \sim 51∼5),以模板题难度 11 ...
[BZOJ]2194: 快速傅立叶之二
题目大意:给定序列a,b,求序列c满足c[k]=sigma(a[i]*b[i-k]) (k<=i<n).(n<=10^5) 思路:观察发现就是普通的卷积反一反(翻转ab其中一个后做卷 ...
BZOJ_2194_快速傅立叶之二_(FFT+卷积)
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=2194 给出序列\(a[0],a[1],...,a[n-1]\)和\(b[0],b[1],... ...
【刷题】BZOJ 2194 快速傅立叶之二
Description 请计算C[k]=sigma(a[i]*b[i-k]) 其中 k < = i < n ,并且有 n < = 10 ^ 5. a,b中的元素均为小于等于100的非 ...
bzoj 2194 快速傅立叶之二
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2194 因为卷积的第 k 项是 sigma(i=0~k)a[ i ]*b[ k-i ] ,也就 ...
c语言二维图像傅立叶变换,C语言数据结构算法之实现快速傅立叶变换.pdf
C语语言言数数据据结结构构算算法法之之实实现现快快速速傅傅立立叶叶变变换换 C语语言言数数据据结结构构算算法法之之实实现现快快速速傅傅立立叶叶变变换换本实例将实现二维快速傅立叶变换,同时也将借此实例 ...

bzoj2194: 快速傅立叶之二

bzoj2194: 快速傅立叶之二相关推荐

最新文章

热门文章