ACM多项式算法入门

在摸索了很久很久，还是没找到多项式的正确入门方式，所以尝试着自己写一下

内容包含

• 多项式乘法
• 多项式求逆
• 多项式除法
• 多项式取模
• 多项式开根（待更新）
• 多项式求ln（待更新）
• 多项式求exp（待更新）

多项式乘法

• 对于多项式$A(x),B(x)$ $$A(x)=a_n{x}^n+a_{n-1}{x}^{n-1}+a_{n-2}{x}^{n-2}+a_{n-3}{x}^{n-3}+...+a_{1}x+a_0$$$$B(x)=b_n{x}^n+b_{n-1}{x}^{n-1}+b_{n-2}{x}^{n-2}+b_{n-3}{x}^{n-3}+...+b_{1}x+b_0$$
  $G(x)=A(x)B(x)$ 中第 $x$ 项的系数 $$g_x=\sum _{i+j=x}{a}_i{b}_j$$
  这是一个卷积的形式，因此我们可以考虑使用$FFT$ 或者 $NTT$ 来优化这个卷积，这样就可以在 $O(nlog(n))$ 时间内实现多项式乘法

void mul(LL* a,int n, LL* b,int m){
   
    for(int i=n+1; i<=Lim; i++) a[i] = 0;
    for(int i=m+1; i<=Lim; i++) b[i] = 0;
    for(int i=0; i<=n; i++) a[i] = (a[i]%P + P) % P; // 第一个多项式系数
    for(int i=0; i<=m; i++) b[i] = (b[i]%P + P) % P; // 第二个多项式系数
    NTT(a,1); NTT(b,1);
    for(int i=0; i<= Lim; i++) a[i] = a[i]*b[i] % P; // O（n）求点值表示法
    NTT(a,-1); // 逆变换变成系数表示法
}

多项式求逆

• 定义

  • 对于多项式 $A(x)$ 若存在多项式 $B(x)$ 满足 $A(x)B(x)=1\mathrm{m}\mathrm{o}\mathrm{d}(x^n)$，那么称