【学习笔记】拉格朗日插值法

已知这样一个方程。
$$\{\begin{array}{c}f(1)=y_1\\ f(2)=y_2\\ f(3)=y_3\\ \cdots \\ f(n)=y_n\end{array}$$
求$f(x)$的式子。

---

首先由小奥知识知道，上面的这个方程可唯一确定一个n次多项式。然后就可以用高斯消元来搞了。
然而时间复杂度是$O(n^3)$的，在这里你应该过不去（小常数选手、暴力碾标算选手等神仙可以在下方教♂导本萌新交流）。
于是我们开始介绍拉格朗日插值法。

---

拉格朗日插值法可以将上面的那个方程组（称作点值表达）转化为系数表达（即多项式）。
那么这个插值的过程是怎样的呢？
我们不妨重新定义$f(x)$：
$$\begin{gathered} f(x)=\sum _{i=1}^n{l}_i(x_i)y_i \\ 其中l_i(x)=\prod _{j=1\&j\ne i}^n\frac{x-x_j}{x_i-x_j} \end{gathered}$$
发现将$y_i,i\in [1,n]$代入这个式子都成立。
这就完了？

---

好像真完了
那。。。接下来就贴代码吧qwq

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int mn = 2005, mod = 998244353;
ll a[mn], x[mn], y[mn];
inline ll ksm(ll a, int b)
{
    if(a < 0) a += mod;
    ll ret = 1;
    while(b)
    {
        if(b & 1)
            (ret *= a) %= mod;
        (a *= a) %= mod, b >>= 1;
    }
    return ret;
}
int main()
{
    int n, k;
    cin >> n >> k;
    for(int i = 1; i <= n; i++)
        cin >> x[i] >> y[i];
    ll ans = 0;
    for(int i = 1; i <= n; i++)
    {
        ll tmp = y[i];
        for(int j = 1; j <= n; j++)
            if(i != j)
                (tmp *= (k - x[j] + mod) % mod * ksm(x[i] - x[j], mod - 2) % mod) %= mod;
        (ans += tmp) %= mod;
    }
    printf("%lld\n", ans);
}