[ 斯特林数 ] Codeforces961G Partitions

显然最终答案中 wi w i 的系数是相等的。我们尝试计算 wi w i 的系数。
因为 i i j 在同一集合中会使 wiwj w i , w j 的系数加 1 1 ,所以 wi 的系数可以表示成

j=1ngi,j ∑ j = 1 n g i , j

其中 gi,j g i , j 表示 i i j 在同一集合内的概率。
利用第二类斯特林数, gi,i=S(n,k),gi,j,ij=S(n1,k) g i , i = S ( n , k ) , g i , j , i ≠ j = S ( n − 1 , k )
求出 S(n,k) S ( n , k ) S(n1,k) S ( n − 1 , k ) 就好了。

#include<bits/stdc++.h>
using namespace std;
const int N=200010;
const int P=1000000007;
int k,n,m,x;
int s1[N],s2[N];
int Ans;
int Pow(int x,int y) {
    int Ans=1;
    for(;y;y>>=1,x=1ll*x*x%P)
        if(y&1) Ans=1ll*Ans*x%P;
    return Ans;
}
void Add(int& x,int y) {
    x=(x+y)%P;
}
void init() {
    s1[0]=s1[1]=s2[0]=s2[1]=1;
    for(int i=2;i<=k;i++) s2[i]=1ll*s2[P%i]*(P-P/i)%P;
    for(int i=2;i<=k;i++) s1[i]=1ll*s1[i-1]*i%P,s2[i]=1ll*s2[i-1]*s2[i]%P;
}
int C(int n,int m) {
    return 1ll*s1[n]*s2[m]%P*s2[n-m]%P;
}
int Calc() {
    int s=0,t=1;
    for(int i=0;i<k;i++,t*=-1) Add(s,1ll*t*C(k,i)*Pow(k-i,n)%P);
    s=1ll*s*s2[k]%P;
    int S=0;t=1;
    for(int i=0;i<k;i++,t*=-1) Add(S,1ll*t*C(k,i)*Pow(k-i,n-1)%P);
    S=1ll*S*s2[k]%P;
    Add(s,1ll*(n-1)*S%P);
    return (s+P)%P;
}
int main() {
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++) scanf("%d",&x),Add(Ans,x);
    init();
    cout<<1ll*Ans*Calc()%P<<endl;
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值