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

问题描述

• lagrange插值解决一种问题:给定 $n+1$个点和值（可以看作$(x,y)$）,用以描述一个$n$次多项式，需要求出其他$x$对应的多项式的值。

• 设前$n+1$个点为$(x_i,y_i)$，对于当前 的点$k$，有
  $$f(k)=\sum _{i=0}^n{y}_i\prod _{j\ne i}\frac{(k-x[j])}{(x[i]-x[j])}$$

• 证明：首先n+1个点可以确定一个n次多项式。
  对于给出的点$(x_i,y_i)$,带入 $x_i$可得结果为$y_i$，所以这$n+1$个点在此$n$次多项式上，即可确定有且仅有这一多项式。

• 应用：多点查值：
  已有$n+1$个点，要对于询问的$m$个点$k$查出多项式的值，要求做到$O(nm+n^2)$
  思路：我们把式子变形，分成与 给定点相关的项，$O(n^2)$预处理，再每次$O(n)$计算与$k$相关的项即可。

• 实现：
  $$f(k)=\sum _{i=0}^n{y}_i\prod _{j\ne i}\frac{(k-x[j])}{(x[i]-x[j])}$$

令$$t[i]=\frac{y_i}{{\prod }_{j=0}^n(x[i]-x[j])}$$

有

$$f(k)=\sum _{i=0}^{n}t[i]\prod _{i\ne j}(k-x[j])$$
令$$g=\prod _{i=0}^n(k-x[i])$$
有$$f(k)=g\sum _{i=0}^n\frac{t[i]}{(k-x[i])}$$

预处理出$t[i]$,$O(n)$查出$f(k)$。

---

• 模板：洛谷P4781
  注意$(x[i]-x[j])$不要改变其正负性