社区讨论

求助 90pts

P4781【模板】拉格朗日插值参与者 1已保存回复 1

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
1 条
当前快照
1 份
快照标识符
@lo18m7t4
此快照首次捕获于
2023/10/22 16:59
2 年前
此快照最后确认于
2023/11/02 16:49
2 年前
查看原帖
CPP
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int mod=998244353;
const int N=2e3;
int n,k,x[N],y[N];
int quickly_pow(int a,int p,int mod){
	int tmp=1;
	while(p){
		if(p&1) tmp=tmp*a%mod;
		a=a*a%mod;
		p>>=1;
	} 
	return tmp%mod;
}
int inv(int x){
	return quickly_pow(x,mod-2,mod);
}
signed main(){
	cin>>n>>k;
	for(int i=1;i<=n;i++){
		cin>>x[i]>>y[i];
	}
	int ans=0;
	for(int i=1;i<=n;i++){
		int s1=y[i]%mod,s2=1;
		for(int j=1;j<=n;j++){
			if(i!=j){
				s1=s1*(k-x[j])%mod;
				s2=s2*(x[i]-x[j])%mod;
			}
		}
		ans=(ans+inv(s2)*s1%mod)%mod;
	}
	cout<<(ans+mod)%mod;
	return 0;
} 
最后一个点WA了

回复

1 条回复,欢迎继续交流。

正在加载回复...