BZOJ5339：[TJOI2018]教科书般的亵渎——题解

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/ +

+++++++++++++++++++++++++++++++++++++++++++