从多项式乘法到快速傅里叶变换,Miskcoo
FFT用到的各种素数,Miskcoo
直接上三模数NTT
#pragma GCC optimize(2)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const int MAXN=2100000;//两倍空间
int n,m,GP,r[MAXN],nn,mm;
LL a[MAXN],b[MAXN],c[MAXN],d[MAXN],ans[MAXN][4];
inline int read() {
char ch=getchar();int sum=0;
while (ch>'9' || ch<'0') ch=getchar();
while (ch>='0' && ch<='9') sum=sum*10+ch-'0',ch=getchar();
return sum;
}
void print(LL x) {
if (x>9) print(x/10);
putchar(x%10+'0');
}
//2009国家集训队论文:骆可强:《论程序底层优化的一些方法与技巧》
inline LL fast_mul(LL a,LL b,