拉格朗日插值

【算法简介】

给定一些点值，更加快速的求出一个点对应的值

基础的拉格朗日插值可以在的时间算出来一个点对应的值

具体地，给定了点

假设要求的点为，那么

 

当x取值连续的时候，可以优化成

先预处理出来

    

那么

 

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

【代码】

#include<bits/stdc++.h>
using namespace std;
#define mp make_pair
#define fi first
#define se second
#define lson now<<1
#define rson now<<1|1
typedef long long ll;
const int maxn=2e3+5;
const int mod=998244353;
int n,k;
int x[maxn],y[maxn];
int qpow(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1) res=1LL*res*a%mod;
        a=1LL*a*a%mod;
        b>>=1;
    }
    return res;
}
int main()
{
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++)
        scanf("%d%d",&x[i],&y[i]);
    ll ans=0;
    for(int i=1;i<=n;i++)
    {
        ll res=y[i];
        for(int j=1;j<=n;j++)
        {
            if(i==j) continue;
            res=res*(k-x[j])%mod;
            res=res*qpow(x[i]-x[j],mod-2)%mod;
        }
        ans=(ans+res)%mod;
    }
    printf("%lld",(ans+mod)%mod);
    return 0;
}