FFT
今天把之前一直没有系统整理的FFT弄了一下,考虑到分治版本的效率太低,于是还是写了蝴蝶变换的板子,感觉理解起来有些困难但是还是能够接受。
要记得 i n v = 1 inv=1 inv=1是DFT, i n v = − 1 inv = -1 inv=−1是IDFT,之前比赛就这边有问题,要记住。
还有在计算答案的时候,需要对于结果除以整个的长度进行修正!!
我终于把之前没有看懂的复数根看懂了QAQ
找了一道十分裸的板子写了一下,luogu上的A*B问题作为自己FFT的板子。
#include <bits/stdc++.h>
#define rep( i , l , r ) for( int i = (l) ; i <= (r) ; ++i )
#define per( i , r , l ) for( int i = (r) ; i >= (l) ; --i )
#define erep( i , u ) for( int i = head[(u)] ; ~i ; i = e[i].nxt )
using namespace std;
struct Complex{
double x , y;
Complex( double _x = 0, double _y = 0. ){
x = _x;
y = _y;
}
Complex operator +(const Complex &b) {
return Complex( x + b.x , y + b.y );
}
Complex operator -(const Complex &b) {
return Complex( x - b.x , y - b.y