2020.02.29日常总结兼拉格朗日插值算法略讲

$$【前置知识】$$

• 快速幂($O(\log n)$)
• 利用快速幂求逆元

---

$$【算法略讲】$$

$【模板链接】：$ 洛谷P4781

$【算法略讲】：$ 所谓的拉格朗日插值算法，就是证明了，当直线$y=f(x)$经过了$n$个点$(x_1,y_1),(x_2,y_2)\cdots (x_n,y_n)$时，对于任意的$t$有：

$$f(t)=\sum _{i=1}^n({\Pi }_{j=1}^n\frac{t-x_j}{x_i-x_j}(i\ne j))$$

直接代入$t=k$求解，加以逆元辅助即可。时间复杂度$O(n^2)$。

$【代码】：$

const int mod=998244353;
int ksm(int a,int b){
	register int ret=1;
	while (b){
		if (b&1) ret=1ll*ret*a%mod;
		a=1ll*a*a%mod;b>>=1;
	}
	return ret;
}//快速幂算法 
const int N=2010;
int n,k,ans,x[N],y[N];
int main(){
	scanf("%d%d",&n,&k);
	for(int i=1;i<=n;i++)
		scanf("%d%d",&x[i],&y[i]);
//	输入n个点(xi,yi) 
	for(int i=1;i<=n;i++){
		int s1=y[i]%mod,s2=1;
		for(int j=1;j<=n;j++)
			if (i!=j){//注意边界条件，否则将RE 
				s1=1ll*s1*((k-x[j]+mod)%mod)%mod;
				s2=1ll*s2*((x[i]+mod-x[j])%mod)%mod;
			} 
			ans=(1ll*ans+(1ll*s1*ksm(s2,mod-2)%mod))%mod;
//			注意中间结果可能会爆int，需要中转long long求值 
	}
	printf("%d",ans);
	return 0;
}