①、CF622F The Sum of the k-th Powers:
求 sum of im ,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怪物死亡。
小豆使用一张 “亵渎”会获得一定的分数,分数计算如下,在使用一张“亵渎”之后,每一个被亵渎造成伤害的怪会产生xk ,其中x是造成伤害前怪的血量为x和需要杀死所有怪物所需的“亵渎”的张数k。
输入格式
第一行输入一个T(T≤10),表示有多少组测试数据
每组组测试数据第一行为n,m,表示有当前怪物最高的血量n,和m种没有出现的血量
接下来m行,每行1个数ai ,表示场上没有血量为ai 的怪物
输出格式
一共T行,每行一个数, 第i行表示第i组测试数据中小豆的最后可以获得的分数, 因为这个分数会很大需要模109+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;
}