拉格朗日插值

拉格朗日插值

给$n$个点$(x_i,y_i)$，要唯一确定一个$n-1$次多项式$f(x)$，使得$\forall i,f(x_i)=y_i$。

有简单的结论：
$$f(x)=\sum _{i=1}^n{y}_i\prod _{j\ne i}\frac{x-x_j}{x_i-x_j}$$

任意条件$O(n^2)$插值

从这个式子可以简单的得到一个$O(n^2)$插值的做法。

令$h(x)={\prod }_{i=1}^n(x-x_i)$，$f(x)$可以简写为${\sum }_{i=1}^n\frac{y_i}{{\prod }_{j\ne i}(x_i-x_j)}\frac{h(x)}{x-x_i}$，其中前面的分式可以暴力$O(n^2)$的算，$h(x)$也可以暴力$O(n^2)$卷，后面的$\frac{h(x)}{x-x_i}$每次可以$O(n)$递推算出来，故总时间复杂度$O(n^2)$。

一定条件下$O(n{\log }^{2}n)$插值

观察一下上面$f(x)$的简写形式，发现可以分治NTT做到$O(n{\log }^{2}n)$，但是难点在怎么快速的计算${\prod }_{j\ne i}(x_i-x_j)$。

令$g_i(x)=\frac{h(x)}{x-x_i}$，发现就是要求$g_i(x_i)$，把$x_i$带入并用洛必达法则可以得到：
$$\underset{x\to x_i}{\lim }\frac{h(x)}{x-x_i}=\underset{x\to x_i}{\lim }{h}^{\prime }(x)=h^{\prime }(x_i)$$
所以就是要求$h^{\prime }(x)$的$n$个点值，直接$O(n{\log }^{2}n)$多点求值即可。

但是这个做法只有在模数比较好的情况下才能实现。