正题
题目链接:https://www.luogu.com.cn/problem/P4781
题目大意
给出nnn个点确定一个多项式y=f(x)y=f(x)y=f(x),求f(k)f(k)f(k)的值
解题思路
拉格朗日插值:
f(k)=∑i=1nyi∏i≠jk−xjxi−xjf(k)=\sum_{i=1}^ny_i\prod_{i\neq j}\frac{k-x_j}{x_i-x_j}f(k)=i=1∑nyii=j∏xi−xjk−xj
时间复杂度O(n2)O(n^2)O(n2)
codecodecode
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=2e3+10,XJQ=998244353;
ll n,k,x[N],y[N],ans;
ll power(ll x,ll b){x%=XJQ;ll ans=1;while(b){if(b&1)ans=ans*x%XJQ;x=x*x%XJQ;b>>=1;}return ans;
}
int main()
{scanf("%lld%lld",&n,&k);for(ll i=1;i<=n;i++)scanf("%lld%lld",&x[i],&y[i]);for(ll i=1;i<=n;i++){for(ll j=1;j<=n;j++)if(i!=j)(y[i]*=(k-x[j]+XJQ)%XJQ*power(x[i]-x[j]+XJQ,XJQ-2)%XJQ)%=XJQ;(ans+=y[i])%=XJQ; }printf("%lld",ans);
}