社区讨论
求助 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 条回复,欢迎继续交流。
正在加载回复...