FFT NNT

算算劳资已经多久没学新算法了，又要重新开始学辣。直接扔板子，跑。。。话说FFT算法导论里讲的真不错，去看下就懂了。

//FFT
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
typedef __int64 ll;
const double pi = acos(-1.0);
const int N = +;
const double eps = 1e-;

struct Complex {
        double x,y;
        Complex(double _x=,double _y=)
                        :x(_x),y(_y) {}
        Complex operator + (const Complex& tt)const { return Complex(x+tt.x,y+tt.y); }
        Complex operator - (const Complex& tt)const { return Complex(x-tt.x,y-tt.y); }
        Complex operator * (const Complex& tt)const { return Complex(x*tt.x-y*tt.y,x*tt.y+y*tt.x); }
};
void FFT(Complex *a, int n, int rev) {
        // n是大于等于相乘的两个数组长度的2的幂次
        // 从0开始表示长度，对a进行操作
        // rev==1进行DFT，==-1进行IDFT
        for (int i = ,j = ; i < n; ++ i) {
                for (int k = n>>; k > (j^=k); k >>= );
                if (i<j) swap(a[i],a[j]);
        }
        for (int m = ; m <= n; m <<= ) {
                Complex wm(cos(*pi*rev/m),sin(*pi*rev/m));
                for (int i = ; i < n; i += m) {
                        Complex w(1.0,0.0);
                        for (int j = i; j < i+m/; ++ j) {
                                Complex t = w*a[j+m/];
                                a[j+m/] = a[j] - t;
                                a[j] = a[j] + t;
                                w = w * wm;
                        }
                }
        }
        if (rev==-) {
                for (int i = ; i < n; ++ i) a[i].x = (a[i].x/n+eps);
        }
}

//NNT
#include <cstdio>
#include <cmath>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
typedef __int64 ll;
const double pi = acos(-1.0);
const int N = +;
const double eps = 1e-;
const int MOD = ; 

int f[N<<];

int pow_mod(int x, int n) {
    int ret = ;
    while(n) {
        if(n&) ret = (ll)ret*x%MOD;
        x = (ll)x*x%MOD;
        n >>= ;
    }
    return ret;
}

int Top, m1[N<<], m2[N<<];

void FFT(int *a, int n, int rev) {
    for(int i = , j = ;i < n; i++) {
        for(int k = n>>;k > (j^=k); k >>= ) ;
        if(i < j)   swap(a[i], a[j]);
    }
    for(int i = ;i < n; i++)   m2[i] = m1[i*(Top/n)];
    for(int m = ;m <= n; m <<= ) {
        int wm = rev ==  ? n-n/m : n/m;
        for(int i = ;i < n; i += m) {
            int w = ;
            for(int j = i;j < i+m/; j++) {
                int t = (ll)m2[w]*a[j+m/]%MOD;
                a[j+m/] = a[j] - t;
                if(a[j+m/] < )    a[j+m/] += MOD;
                a[j] = a[j] + t;
                if(a[j] >= MOD) a[j] -= MOD;
                w += wm;
                if(w >= n)  w -= n;
            }
        }
    }
    if(rev == -) {
        int inv = pow_mod(n, MOD-);
        for(int i = ;i < n; i++) a[i] = (ll)a[i]*inv%MOD;
    }
}