仙人NTT的入门

仙人NTT的入门

更详细的NTT在这里

在快速傅里叶变换中,我们利用$w_n$单位复数根实现了消去引理和折半引理
但是由于复数运算的关系,导致精度问题,使人十分捉鸡
那么,有没有什么整数也满足消去引理和折半引理来代替$w_n$单位复数根呢?
这就是所谓的原根

那么什么是原根呢?
定义$p$的原根为满足$g^{\phi (p)}\equiv 1 \pmod p$的整数$g$,其中,$\phi(p)为满足条件的最小指数$
并且要求$g^i\not \equiv g^j \pmod p$
那么我们用$g^i$生成的横坐标就不会相同,以此来生成点集表示法

现在我们用$g^{\frac{\phi(p)}{n}}$来代替$e^{\frac{2\pi i}{n}}$
因为我们想利用整数进行计算,那么$^{\frac{\phi(p)}{n}}$应该为整数,因为当$p$为质数时,$\phi(p)=p-1$,那么$p-1\equiv 0\pmod {2^q},n\leq 2^q$

消去引理: $w_{dn}^{dk}=w_{n}^{k}$
* 证明:$w_{dn}^{dk}=g^{\frac{dk\phi(p)}{dn}}=g^{\frac{k\phi(p)}{n}}=w_{n}^{k}$

折半引理: 若$n$为偶数,则$n$次单位复数根的平方的集合就是$\frac{n}{2}$次单位复数根的集合,特别的,每个$\frac{n}{2}$次单位复数根出现两次
* 证明:$(w_{n}^{k})^2=w_{\frac{n}{2}}^k,k\in[0,n-1]$ [消去引理]

非常nice

所以就和FFT一样一样的啦!

板子!
多项式乘法

#include <bits/stdc++.h>
using namespace std;
const int N = 3000010;
const int G=3,mod=(119<<23)+1;
int n,m,l;
int a[N],b[N],r[N];
int read() {
    int ans=0,flag=1;
    char ch=getchar();
    while((ch>'9' || ch<'0') && ch!='-') ch=getchar();
    if(ch=='-') flag=-1,ch=getchar();
    while(ch>='0' && ch<='9') ans=ans*10+ch-'0',ch=getchar();
    return ans*flag;
}
int ksm(int a,int b) {
    int ans=1;
    while(b) {
        if(b&1) ans=1ll*ans*a%mod;
        a=1ll*a*a%mod;
        b>>=1;
    }
    return ans;
}
void ntt(int *now,int f) {
    for(int i=0;i<n;++i)
        if(i<r[i]) swap(now[i],now[r[i]]);
    for(int i=1;i<n;i<<=1) {
        int gn=ksm(G,(mod-1)/(i<<1));
        for(int j=0;j<n;j+=(i<<1)) {
            int x,y,g=1;
            for(int k=0;k<i;++k,g=1ll*g*gn%mod) {
                x=now[j+k],y=1ll*g*now[j+k+i]%mod;
                now[j+k]=(x+y)%mod;
                now[j+k+i]=(x-y+mod)%mod;
            }
        }
    }
    if(f!=1) {
        int ny=ksm(n,mod-2);
        reverse(now+1,now+n);
        for(int i=0;i<n;++i) now[i]=1ll*now[i]*ny%mod;
    }
}
int main() {
    n=read();m=read();
    for(int i=0;i<=n;++i) a[i]=read();
    for(int i=0;i<=m;++i) b[i]=read();
    m+=n;
    for(n=1;n<=m;n<<=1) ++l;
    for(int i=0;i<n;++i)
        r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));
    ntt(a,1);ntt(b,1);
    for(int i=0;i<n;++i) a[i]=1ll*a[i]*b[i]%mod;
    ntt(a,-1);
    for(int i=0;i<=m;++i)
        printf("%d ",a[i]);
    return 0;
}