退役前寫的東西
令\(F(x)\)為\(n\)次項多項式
拉格朗日插值：\(f(x)=\sum\limits_{k=0}^n f(x_k)l_k(x)=\sum\limits_{k=0}^n f(x_k)\prod\limits_{i\neq k}^n \frac{x-x_i}{x_k-x_i}\)
因為很簡單記憶，在OI中應用廣泛
缺點：在增加或減少次項時需要重新全部計算
為實現在增加或減少次項時快速計算，我們構造：
這種形式的插值多項式稱為牛頓插值多項式。記為\(N_n(x)\)，需要滿足\(N_n(x_i)=y_i\)
定義1：設函數\(f(x)\)在點\(x_0,x_1,\cdots\)的值為\(f(x_0),f(x_1),\cdots\)，稱\(\frac{f(x_j)-f(x_i)}{x_j-x_i}(i\neq j)\)為\(f(x)\)在點\(x_i,x_j\)處的一階差商，記為\(f[x_i,x_j]\)。稱一階差商的差商\(\frac{f[x_j,x_k]-f[x_i,x_j]}{x_k-x_i}(i,j,k互異)\)為\(f(x)\)在\(x_i,x_j,x_k\)處的二階差商，記為\(f[x_i,x_j,x_k]\)。
一般地，稱m-1階差商的差商：\(f[x_0,x_1,\cdots,x_m]=\frac{f[x_1,x_2,\cdots ,x_m]-f[x_0,x_1,\cdots ,x_{m-1}]}{x_m-x_0}\)，為\(f(x)\)在點\(x_0,x_1,\cdots, x_m\)處的\(m\)階差商。特別地，零階差商\(f[x_i]=f(x_i)\)
性質1：\(k\)階差商\(f[x_0,x_1,\cdots ,x_k]\)是由\(f(x_0),f(x_1),\cdots ,f(x_k)\)線性組合成的。有\(f[x_0,x_1,\cdots,x_k]=\sum\limits_{i=0}^k \frac{f(x_i)}{(x_i-x_0)\cdots (x_j-x_{j-1})(x_j-x_{j+1})\cdots (x_j-x_k)}\)
證明：
靠歸納顯然成立
性質2：差商具有對稱性，即在\(k\)階差商\(f[x_0,x_1,\cdots ,x_k]\)中任意調換\(2\)個節點\(x_i,x_j\)的順序，其值不變。
證明：
根據結論1顯然成立
性質3：\(k\)階差商\(f[x_0,x_1,\cdots ,x_k]\)和\(k\)階導數\(f^{(k)}(x)\)之間有如下關系：
\(f[x_0,x_1,\cdots ,x_k]=\frac{f^{(k)}(\delta)}{k!}(\delta \in (min\{x_0,x_1,\cdots,x_k\},max\{x_0,x_1,\cdots ,x_k\})\)
證明：
挖坑。下面用不到這條性質。
考慮一個一個添加項：
由\(N_0(x_0)=f(x_0)\)，可得
\(a_0=f(x_0)=f[x_0]\)
由\(N_1(x_1)=f(x_1)\)，可得
\(a_1=\frac{f(x_1)-f(x_0)}{x_1-x_0}=f[x_0,x_1]\)
由\(N_n(x_2)=f(x_2)\)，可得
\(\begin{aligned} a_2&=\frac{f(x_2)-f(x_0)-f[x_0,x_1](x_2-x_0)}{(x_2-_0)(x_2-x_1)}\\ &=\frac{\frac{f(x_2)-f(x_0)}{x_2-x_0}-f[x_0,x_1]}{x_2-x_1}\\ &=\frac{f[x_0,x_2]-f[x_0,x_1]}{x_2-x_1}\\ &=f[x_1,x_0,x_2]\\ &=f[x_0,x_1,x_2]\\ \end{aligned}\)
一般地，可以證明有\(a_k=f[x_0,x_1,\cdots ,x_k]\)
故\(N_n(x_i)=f(x_i)(i=0,1,2,\cdots ,n)\)的\(n\)次牛頓插值多項式為：
\(\begin{aligned}\\ N_n(x)&=f[x_0]+f[x_0,x_1](x-x_0)+f[x_0,x_1,x_2](x-x_0)(x-x_1)\\ &+\cdots +f[x_0,x_1,\cdots ,x_n](x-x_0)(x-x_1)\cdots (x-x_{n-1})\\ \end{aligned}\)
另一種簡單的證明方式：
我們來證明\([y_1,\cdots,y_n](x_n-x_1)\cdots (x_n-x_{n-1})=y_n-P(x_n)\)（\(P(x)\)是最高度為\(n-2\)的通過點\((x_1,y_1),\cdots,(x_{n-1},y_{n-1})\)）。
若成立，容易得到牛頓多項式就為原多項式。
歸納法：
邊界：\([y_1,y_2](x_2-x_1)=\frac{y_2-y_1}{x_2-x_1}(x_2-x_1)=y_2-y_1=y_2-P(x_2)\)（\(P(x)=y_1\)）
\(\begin{aligned} &[y_1,\cdots,y_{n+1}](x_{n+1}-x_1)\cdots (x_{n+1}-x_{n})\\ &=\frac{[y_2,\cdots,y_{n+1}]-[y_1,\cdots,y_n]}{x_{n+1}-x_1}(x_{n+1}-x_1)\cdots (x_{n+1}-x_{n})\\ &=([y_2,\cdots,y_{n+1}]-[y_1,\cdots,y_n])(x_{n+1}-x_2)\cdots (x_{n+1}-x_{n})\\ &=(y_{n+1}-Q(x_{n+1}))-[y_1,\cdots,y_n](x_{n+1}-x_2)\cdots (x_{n+1}-x_{n})\\ &=y_{n+1}-(Q(x_{n+1})+[y_1,\cdots,y_n](x_{n+1}-x_2)\cdots (x_{n+1}-x_{n}))\\ \end{aligned}\)
（其中\(Q(x)\)是最高度為\(n-2\)通過點\((x_2,y_2),\cdots,(x_n,y_n)\)）
令\(P(x)=Q(x)+[y_1,\cdots,y_n](x-x_2)\cdots(x-x_n)\)，下面證明\(P(x)\)通過點\((x_1,y_1)\cdots(x_{n},y_n)\)且最高度為\(n-1\)。下面證明\(P(x)\)通過\((x_1,y_1)\)，根據性質2剩下點也容易證明。
\(\begin{aligned} &Q(x_1)+[y_1,\cdots,y_n](x_1-x_2)\cdots(x_1-x_n)\\ &=Q(x_1)+[y_n,\cdots,y_1](x_1-x_n)\cdots(x_1-x_2)\\ &=Q(x_1)+y_1-Q(x_1)\\ &=y_1\\ \end{aligned}\)