拉格朗日插值

在平面上给出$n$个点的坐标 { x i , y i } \{x_i,y_i\} {xi​,yi​},求唯一的$n-1$次多项式$f$,满足$\forall i\in [1,n],f(x_i)=y_i$,即经过这$n$个点

显然，这$n$个$x$互不相同

设多项式$l_j={\prod }_{i=1,i\ne j}^n\frac{x-x_i}{x_j-x_i}$

显然

$$l_j(i)=\{\begin{array}{rlrlrl}& & & 0& & x_i\ne x_j\\ & & & 1& & x_i=x_j\end{array}$$

$$f(x)=\sum _{i=1}{y}_i\times l_i(x)$$

$\forall i\in [1,n],x_i\in [1,n]$怎么快速求$f(qwq)$？

$$l_j=\prod _{i=1,i\ne j}^n\frac{x-i}{j-i}$$

上面用两个数组维护一下$\prod qwq-i$，下面阶乘+逆元

$$f(qwq)=\sum _{i=1}^n{y}_i\times l_i(qwq)$$
$$=\sum _{i=1}^n{y}_i\times \prod _{j=1}^{i-1}\frac{qwq-j}{i-j}\times \prod _{j=i+1}^n\frac{qwq-j}{i-j}$$
$$=\sum _{i=1}^n{y}_i\times \frac{{\prod }_{j=1}^{i-1}qwq-j}{(i-1)!}\times (-1{)}^{n-i}\frac{{\prod }_{j=i+1}^{n}qwq-j}{(n-i)!}$$
$$=\sum _{i=1}^n(-1{)}^{n-i}\times y_i\times \frac{{\prod }_{j=1}^{n}qwq-j}{(i-1)!\times (n-i)!\times (qwq-i)}$$