其实一个多项式f(x)被表示成广泛一点的形式应该是这样的:
f(x)=i∑gi(x)×ci
平常我们写的多项式是自然数幂的形式,也就是gi(x)=xi
插值的定义是在离散数据的基础上补插连续函数,使得这条连续曲线通过全部给定的离散数据点。
前置技能: 组合数交错和求和
这里组合数意义用小学教的方法表示是因为这样第一项为负数的时候依然有意义(无组合意义)
i=0∑n(−1)i(im)
=i=0∑n(−1)ii!m×(m−1)×⋯×(m−i+1)
=i=0∑ni!−m×(−m+1)×⋯×(−m+i−1)
=i=0∑n(i−m+i−1)
这一步证明合集有,就是前面加一项(−1−m−1)
=(nn−m)
线性插值
多项式有其他的表示方式,这里用到的是这样的:
f(k)=i∑(ik)ci
现在我们知道f(1⋯n)的值是a1⋯n(n足够),求f(k)。
当n比较小的时候,我们可以通过偏序函数的容斥把每一项ci都求出来:
ci=j=0∑i(−1)i−j(ji)aj
复杂都O(N2),没有更好的解决方法,于是只能将式子带入。
f(k)=i∑(ik)ci
=i=0∑n(ik)j=0∑i(−1)i−j(ji)aj
=i=0∑nj=0∑i(−1)i−j(ik)(ji)aj
(ca)(bc)的组合意义可以理解为在一个大范围内确定大圈,然后再里面找一个小圈。
这样等价于找一个小圈,再枚举可以贡献到多少个大圈上。
顺便交换一下i,j的枚举顺序。
=j=0∑ni=j∑n(−1)i−j(jk)(i−jk−i)aj
=j=0∑n(jk)aji=j∑n(−1)i−j(i−jk−i)
有点乱,重定义一下i=j,j=i-j(等号后面都是原来的i,j)
=i=0∑n(ik)aij=0∑n−i(−1)j(jk−i)
发现后面的sum式子其实是一个前面说的交错和的形式,所以:
=i=0∑n(ik)ai(n−in−k)
注意到n特别小,所以每次求的时候就可以快速预处理一下,就O(n)了。