拉格朗日插值法：

①、CF622F The Sum of the k-th Powers：
求 sum of i^m ，i from 1 to n mod 1e9+7。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
#define ll long long
#define llu unsigned ll
using namespace std;
const int maxn=1001000;
const int mod=1e9+7;
ll pre[maxn],suf[maxn],fac[maxn],inv[maxn];
ll inv1;

ll mypow(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
ll lg(ll n,ll k)
{
    int m=k+2;
    n%=mod;
    inv1=mypow(mod-1,mod-2);
    fac[0]=pre[0]=1;
    for(int i=1;i<=m;i++)
        fac[i]=fac[i-1]*i%mod,pre[i]=pre[i-1]*(n-i+mod)%mod;

    suf[m+1]=1;
    suf[m]=suf[m+1]*(n-m+mod)%mod;
    inv[m]=mypow(fac[m],mod-2);
    for(int i=m-1;i>=0;i--)
        inv[i]=inv[i+1]*(i+1)%mod,suf[i]=suf[i+1]*(n-i+mod)%mod;

    ll ans=0,y=0;
    for(int i=1;i<=m;i++)
    {
        y=(y+mypow(i,k))%mod;
        ans=(ans+y*pre[i-1]%mod*suf[i+1]%mod*inv[i-1]%mod*inv[m-i]%mod*((m-i)&1?inv1:1)%mod)%mod;
    }
    return ans;
}

int main(void)
{
    ll n,m;
    scanf("%lld%lld",&n,&m);
    printf("%lld\n",lg(n,m));
    return 0;
}

②、P4593 [TJOI2018]教科书般的亵渎：
题目描述
小豆喜欢玩游戏，现在他在玩一个游戏遇到这样的场面,每个怪的血量为ai ，且每个怪物血量均不相同，小豆手里有无限张“亵渎”。亵渎的效果是对所有的怪造成1点伤害，如果有怪死亡，则再次施放该法术。我们认为血量为0怪物死亡。

小豆使用一张 “亵渎”会获得一定的分数，分数计算如下，在使用一张“亵渎”之后，每一个被亵渎造成伤害的怪会产生x^k ,其中x是造成伤害前怪的血量为x和需要杀死所有怪物所需的“亵渎”的张数k。

输入格式
第一行输入一个T(T≤10)，表示有多少组测试数据

每组组测试数据第一行为n，m，表示有当前怪物最高的血量n，和m种没有出现的血量

接下来m行，每行1个数ai ，表示场上没有血量为ai 的怪物

输出格式
一共T行，每行一个数， 第i行表示第i组测试数据中小豆的最后可以获得的分数, 因为这个分数会很大需要模10⁹+7

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
#define ll long long
#define llu unsigned ll
using namespace std;
const int maxn=1001000;
const int mod=1e9+7;
ll pre[maxn],suf[maxn],fac[maxn],inv[maxn],a[maxn];
ll inv1;

ll mypow(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
ll lg(ll n,ll k)
{
    int m=k+2;
    n%=mod;
    inv1=mypow(mod-1,mod-2);
    fac[0]=pre[0]=1;
    for(int i=1;i<=m;i++)
        fac[i]=fac[i-1]*i%mod,pre[i]=pre[i-1]*(n-i+mod)%mod;

    suf[m+1]=1;
    suf[m]=suf[m+1]*(n-m+mod)%mod;
    inv[m]=mypow(fac[m],mod-2);
    for(int i=m-1;i>=0;i--)
        inv[i]=inv[i+1]*(i+1)%mod,suf[i]=suf[i+1]*(n-i+mod)%mod;

    ll ans=0,y=0;
    for(int i=1;i<=m;i++)
    {
        y=(y+mypow(i,k))%mod;
        ans=(ans+y*pre[i-1]%mod*suf[i+1]%mod*inv[i-1]%mod*inv[m-i]%mod*((m-i)&1?inv1:1)%mod)%mod;
    }
    return ans;
}

int main(void)
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        ll n,m;
        scanf("%lld%lld",&n,&m);
        for(int i=1;i<=m;i++)
            scanf("%lld",&a[i]);
        sort(a+1,a+m+1);
        ll ans=0;
        for(int i=0;i<=m;i++)
        {
            ans=(ans+lg(n-a[i],m+1))%mod;
            for(int j=i+1;j<=m;j++)
            {
                ans=((ans-mypow(a[j]-a[i],m+1))%mod+mod)%mod;
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}