https://www.lydsy.com/JudgeOnline/problem.php?id=5339
https://www.luogu.org/problemnew/show/P4593
小豆喜欢玩游戏, 现在他在玩一个游戏遇到这样的场面,每个怪的血量为ai,且每个怪物血量均不相同, 小豆手里有无限张“亵渎”。亵渎的效果是对所有的怪造成1点伤害,如果有怪死亡,则再次施放该法术。我们认为血量为0时怪物死亡。小豆使用一张“亵渎”会获得一定的分数,分数计算如下,在使用一张“亵渎”之后,每一个被亵渎造成伤害的怪会产生x^k,其中x是造成伤害前怪的血量为x和需要杀死所有怪物所需的“亵渎”的张数k。
参考:https://www.luogu.org/blog/user44829/solution-p4593 ,你可以在这个博客里面找到各种各样的本题做法。
题意很乱,但是整理整理后发现实际上你只需要知道如何求出\\(S(n,k)=\\sum_{i=1}^na_i^k\\)即可。
直接给公式\\(S(n,k)=\\frac{1}{k+1}\\sum_{i=1}^{k+1}C^i_{k+1}B_{k+1-i}(n+1)^i\\)
组合数递推:
\\(C^m_n=C^m_{n-1}+C^{m-1}_{n-1}\\)
伯努利数递推:
\\(B_n=[m=0]-\\sum_{k=0}^{m-1}C_m^k\\frac{B_k}{m-k+1}\\)
\\(m=n-1\\)
#include<cmath>
#include<queue>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll p=1e9+7;
const int N=65;
inline ll read(){
ll X=0,w=0;char ch=0;
while(!isdigit(ch)){w|=ch==\'-\';ch=getchar();}
while(isdigit(ch))X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
return w?-X:X;
}
ll n,a[N],inv[N],c[N][N],b[N];
int m;
ll qpow(ll x,ll y){
ll res=1;
while(y){
if(y&1)res=res*x%p;
x=x*x%p;y>>=1;
}
return res;
}
ll f(ll x,ll y){
ll res=0;
for(int i=1;i<=y+1;i++){
res=(res+c[y+1][i]*b[y+1-i]%p*qpow(x+1,i)%p)%p;
}
res=res*inv[y+1]%p;
return res;
}
inline void init(){
for(int i=1;i<N;i++)inv[i]=qpow(i,p-2);
for(int i=0;i<N;i++){
c[i][0]=1;
for(int j=1;j<=i;j++){
c[i][j]=(c[i-1][j]+c[i-1][j-1])%p;
}
}
b[0]=1;
for(int i=1;i<N;i++){
b[i]=0;
for(int j=0;j<=i-1;j++)b[i]=(b[i]+c[i+1][j]*b[j]%p)%p;
b[i]=(p-b[i])*inv[i+1]%p;
}
}
int main(){
init();
int t=read();
while(t--){
n=read(),m=read();
for(int i=1;i<=m;i++)a[i]=read();
a[++m]=++n;
sort(a+1,a+m+1);
ll ans=0;
for(int i=1;i<=m;i++){
for(int j=i;j<=m;j++){
ans=(ans+f(a[j]-1,m)-f(a[j-1],m)+p)%p;
}
for(int j=i+1;j<=m;j++)a[j]-=a[i];
a[i]=0;
}
printf("%lld\\n",ans);
}
return 0;
}
+++++++++++++++++++++++++++++++++++++++++++
+本文作者:luyouqi233。 +
+欢迎访问我的博客:http://www.cnblogs.com/luyouqi233/ +
+++++++++++++++++++++++++++++++++++++++++++