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)<