多项式小结

多项式

一 简介

**多项式的度：**对于一个多项式$f(x)$，称其最高次项的次数为该多项式的度，记作$degf$。

多项式的乘法：

可以通过快速傅里叶变换在$O(nlogn)$下计算。

给定多项式$f(x)$和$g(x)$：
$$\begin{gathered} f(x)=a_0+a_{1}x+\cdot \cdot \cdot +a_n{x}^n \\ g(x)=b_0+b_{1}x+\cdot \cdot \cdot +b_m{x}^m \end{gathered}$$
计算多项式$Q(x)=f(x)\cdot g(x)$：
$$Q(x)=\sum _{i=0}^n\sum _{j=0}^m{a}_i{b}_j{x}^{i+j}=c_0+c_{1}x+\cdot \cdot \cdot +c_{n+m}{x}^{n+m}$$
多项式的逆元：

对于多项式$f(x)$，若存在$g(x)$满足：
$$\begin{gathered} f(x)g(x)\equiv 1(mod{x}^n) \\ degg\le degf \end{gathered}$$
则称$g(x)$为$f(x)$在模$x^n$意义下的逆元，记作$f^{-1}(x)$。

多项式的余数和商：

对于多项式$f(x),g(x)$，存在唯一的$Q(x),R(x)$满足：
$$\begin{gathered} f(x)=Q(x)g(x)+R(x) \\ degQ=degf-degg \\ degR<degg \end{gathered}$$
称$Q(x)$为$g(x)$除$f(x)$的商，$R(x)$为$g(x)$除$f(x)$的余数，可记作：
$$f(x)=R(x)(modg(x))$$
多项式的对数函数与指数函数：

对于一个多项式$f(x)$，可以将其对数函数看作其与麦克劳林级数的复合：
$$\begin{gathered} ln(1-f(x))=-\sum _{i=1}^{+\infty }\frac{f^i(x)}{i} \\ ln(1+f(x))=\sum _{i=1}^{+\infty }\frac{(-1{)}^{i-1}{f}^i(x)}{i} \end{gathered}$$
其指数函数可以定义：
$$expf(x)=e^{f(x)}=\sum _{i=0}^{+\infty }\frac{f^i(x)}{i!}$$
多项式的多点求值和插值：

多项式的多点求值：给出一个多项式$f(x)$和$n$个点$x_1,x_2,...,x_n$，求$f(x_1),f(x_2),...,f(x_n)$。

多项式的插值：给出$n+1$个点$(x_0,y_0),(x_1,y_1),...,(x_n,y_n)$，求一个$n$次多项式$f(x)$使得这$n+1$个点都在$f(x)$上。

两种操作实质是将多项式在系数表示和点值表示间转化。

二 拉格朗日插值

2.1 知识点

时间复杂度：$O(n^2)$

**问题：**给出$n$个点$P_i(x_i,y_i)$，将过这$n$个点的最多$n-1$次的多项式记为$f(x)$，求$f(k)$的值。

**差分法：**适用于$x_i=i$的情况。

利用差分法求多项式$f(x)={\sum }_{i=0}^n{a}_i{x}^i$的多项式形式，时间复杂度$O(n^2)$。

已知$f(1)$至$f(6)$的值分别为$1,5,14,30,55,91$
KaTeX parse error: Undefined control sequence: \ at position 140: …\ \ 2\ \ \ \ \ \̲ ̲
第一行为$f(x)$的连续的前$n$项，之后每一行为上一行相邻两项的差值，最终会变成定值行，即可求解。

待定系数法：

设$f(x)={\sum }_{i=0}^{n-1}{a}_i{x}^i$，将每个$x_i$代入$f(x)$，得到一个由$n$条$n$元一次方程所组成的方程组，然后使用高斯消元求解

拉格朗日插值法：

有：
$$f(x)\equiv f(a)(mod(x-a))$$
因为$f(x)-f(a)=(a_0-a_0)+a_1(x^1-a^1)+a_2(x^2-a^2)+\cdot \cdot \cdot +a_n(x^n-a^n)$有因子$(x-a)$。

列出多项式线性同余方程组：
$$

\begin{cases} f(x)≡y_1\ (mod\ (x-x_1)) \ f(x)≡y_2\ (mod\ (x-x_2)) \ ··· \ f(x)≡y_n\ (mod\ (x-x_n)) \end{cases}
$$根据中国剩余定理：$$
M=\prod_{i=1}^{n}(x-x_i) \
m_i=\frac{M}{x-x_i}=\prod_{j≠i}(x-x_j)
KaTeX parse error: Can't use function '$' in math mode at position 2: $̲m_i$模$(x-x_i)$意…
m_i^{-1}=\prod_{j≠i}\frac{1}{x_i-x_j}
$$所以：$$
f(x)=\sum_{i=1}^{n}y_im_im_i{-1}≡\sum_{i=1}^ny_i\prod_{j≠i}\frac{x-x_i}{x_i-x_j}\ (mod\ M)
$$

如果给定的点为连续的点，那么可以使用前缀后缀优化。

例如给定$(0,f_0),(1,f_1),(2,f_2),...,(n,f_n)$

则记：$pre[i]={\prod }_{j=0}^i(k-j),suff[i]={\prod }_{j=i}^n(k-j),fact[i]=i!$

则$ f(k)=\sum_{i=0}^{n} y_i\frac{pre[i-1]*suff[i+1]}{fact[i]*fact[n-i]}$,但是分母$fact[n-i]$可能出现负数，需要去判断，当 $n-i$ 为奇数时，分母应该取负号。这样就可以预处理得到$fact$数组，然后根据输入的k来$O(n)$计算得到$pre、suff$数组，从而$O(n)$完成$f(k)$计算。

2.2 模板

2.1 随机给定n个点插值：$O(n^{2}lo{g}_n)$

#include <bits/stdc++.h>

#define int long long
using namespace std;

const int MAXN = 2010;
int mod = 998244353;
int n, k, x[MAXN], y[MAXN], res, s1, s2;

int qmi(int a, int k, int p) {
   
    int res = 1 % p;  // res记录答案, 模上p是为了防止k为0，p为1的特殊情况
    while(k) {
     // 只要还有剩下位数
        if (k & 1) res = res * a % p;  // 判断最后一位是否为1，如果为1就乘上a，模上p, 乘法时可能爆int，所以变成long long
        k >>= 1;  // 右移一位
        a = a * a % p;  // 当前a等于上一次的a平方，取模，平方时可能爆int，所以变成long long
    }
    return res;
}

int get_inv(int a, int p) {
   
    return a % p == 0? -1: qmi(a, p - 2, p);
}

// n个点(x[i], y[i])，求f(k)
int Lagrange(int x[], int y[], int n, int k) {
   
    int res = 0;
    for (int i = 1; i <= n; i++) {
   
        s1 = y[i] % mod;  // s1计算分子
        s2 = 1;  // s2计算分母
        for (int j = 1; j <= n; j++)
            if (i != j)
                s1 = s1 * (k - x[j]) % mod, s2 = s2 * (x[i] - x[j]) % mod;  // 分别计算分子和分母
        res += s1 * get_inv(s2, mod) % mod;
    }
    return (res % mod + mod) % mod;
}

signed main() {
   
    scanf("%lld%lld", &n, &k);
    for (int i = 1; i <= n; i++) scanf("%lld%lld", x + i, y + i);  // 读入n个点
    cout << Lagrange(x, y, n, k);
    return 0;
}

2.2 给定连续n个点插值:$O(n)$

void init(int n, int mod) {
     // 预处理逆元
    LL res = 1;
    for (int i = 1; i <= n + 2; ++i) res = res * i % mod;  // 计算阶乘
    fact[n + 2] = qmi(res, mod - 2, mod);  
    for (int i = n + 1; i >= 0; --i) fact[i] = (i + 1) * fact[i + 1] % mod;  // 计算逆元
}

//给定最高次为n的多项式的n+1个点分别为：(0,f0),(1,f1),(2,f2)...(n,fn),求f(k)的值
LL Lagrange(LL f[], int n, int k) {
   
    if (k <= n) return f[k];
    pre[0] = suff[n + 1] = 1;
    for (int i = 0; i <= n; ++i) pre[i + 1] = pre[i] * (k - i) % mod;
    for (int i = n; i >= 0; --i) suff[i] = suff[i + 1] * (k - i) % mod;
    LL fk = 0;
    for (int i = 0; i <= n; ++i) {
   
        LL tmp = f[i] * pre[i] % mod * suff[i + 1] % mod * fact[i] % mod * fact[n - i] % mod;  // 计算 (pre[i-1]*suff[i+1]) / (fact[i]*fact[n-i])
        if ((n - i) & 1)  // 当 n-i 为奇数时，分母应该取负号。
            fk = (fk - tmp + mod) % mod;
        else
            fk = (fk + tmp) % mod;
    }
    return fk;
}

int main() {
   
    ...
    init(MAXN - 10, mod);  // 预处理逆元
    for (int i = 0; i <= n; ++i) cin >> f[i];  // 读入n+1个点
    f[n + 1] = Lagrange(f, n, n + 1);  // 得到f[n+1]
    ...
}

2.3 例题

P4781 【模板】拉格朗日插值

**题意：**给定$n$个点$(x_i,y_i)$确定唯一的多项式，并求出$f(k)mod998244353$。

数据：$1\le n\le 2\times 10^3$，$1\le x_i,y_i,$