拉格朗日插值模板 P4781 【模板】拉格朗日插值

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2e3+5;
const int mod=998244353;
ll x[N],y[N],n,k,ans;
ll quipow(ll a,ll n)
{
    ll ans=1%mod;
    while(n)
    {
        if(n&1)ans=ans*a%mod;
        n>>=1;
        a=a*a%mod;
    }
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    cin>>n>>k;
    for(int i=1;i<=n;i++)cin>>x[i]>>y[i];
    for(int i=1;i<=n;i++)
    {
        ll s1=(y[i]+mod)%mod,s2=1%mod;
        for(int j=1;j<=n;j++)
        {
            if(i==j)continue;
            s1=s1*(k-x[j]+mod)%mod,s2=s2*(x[i]-x[j]+mod)%mod; 
        }
        ans=(ans+s1*quipow(s2,mod-2)%mod)%mod;
    }
    cout<<ans<<"\n";
    return 0;
}