BZOJ-4827-FFT

先不考虑+c,那么只需要把原式展开,将a数组扩大一倍,求fft即可。

#include
#define N 262144
#define pi acos(-1)
#define ll long long
using namespace std;
typedef complex E;
int n,m,L;
int rev[N];
E f[N],_f[N],e[N];
void fft(E *a,int f)
{
	for(int i=0;irev[i])swap(a[i],a[rev[i]]);
	for(int i=1;i>1]>>1)|((i&1)<<(L-1));
	fft(f,1);fft(_f,1);
	for(int i=0;i

 

你可能感兴趣的:(FFT)