拉格朗日插值法学习笔记

对于一个n次函数$P(x)=a_0+a_1x+a_2x^2+....a_nx^n$
如果要求这个函数的系数，有两种通用的方法，一种是$n^3\ $的高斯消元，一种是$n^2logn \ $的拉格朗日插值
但如果要求这个函数在$x=k\ $处的点值，拉格朗日插值可以做到$O(n)$（或者$O(nlogn)$）

这个东西主要是一个构造
一个n次函数要至少n+1个点才能唯一确定，我们先取$(x_0,y_0),(x_1,y_1).....(x_n,y_n)$这n+1个点
然后我们构造出一些多项式
$l(i)=\prod_{j=0,j\neq i}^n \dfrac{(x-x_j)}{(x_i-x_j)}$，可以发现，$l(i)$只有在$x_i$处取到值1，在$x_j(j \neq i)\ $处都取到0
那么$y_il(i)\ $这个多项式就可以取到点$(x_i,y_i)\ $且不影响其他n个点

$P(x)=\sum_{i=0}^nl(i)y_i$

然后这个东西用fft加速卷积的话，可以$n^2logn\ \ $找到一个n次多项式，如果是取一个点的值的话，发现分子那里可以$O(n)\ $算出n+1个数的乘积再乘上对应的逆元，分子那里是两个阶乘乘积，注意一下符号，如果逆元那里可以做到$O(n)$就可以$O(n)$求点值

一道例题：bzoj3453
先对g(i)做插值，卷积可以直接暴力卷，n^3求出一个k+2次的多项式，然后要求g(A)+g(A+d)+g(A+2d)….+g(A+nd)，可以枚举我们插出来的g(i)的多项式的第i项，相当于要计算$A^i+(A+d)^i+....+(A+nd)^i$，这个东西是个i+1次的多项式，我们可以用上面提到的那个方法，$O(nlogn)$得到他在n处的值，带进g(i)的多项式