【CF932E】Team Work

题目

luogu的Romtejudge挂了我就当我过了吧

求

\[\sum_{i=1}^n\binom{n}{i}i^k\]

其实是个思博套路题，但是我现在这个水平还是刷刷板子吧

处理\(x^k\)是一个套路了

\[x^k=\sum_{i=1}^k\begin{Bmatrix} k \\i \end{Bmatrix}\binom{x}{i}i!\]

于是

\[\sum_{i=1}^n\binom{n}{i}i^k=\sum_{i=1}^n\binom{n}{i}\sum_{j=1}^k\begin{Bmatrix} k \\j \end{Bmatrix}\binom{i}{j}j!\]

交换一下求和符号

\[\sum_{j=1}^n\begin{Bmatrix} k \\j \end{Bmatrix}j!\sum_{i=1}^n\binom{n}{i}\binom{i}{j}\]

后面这个\(\sum_{i=1}^n\binom{n}{i}\binom{i}{j}\)也是套路了，我们考虑一下组合意义，发现这个东西就是\(2^{n-j}\binom{n}{j}\)

我们再拆一下组合数，就得到了

\[\sum_{j=1}^k\begin{Bmatrix} k \\j \end{Bmatrix}2^{n-j}n^{\underline j}\]

这个题的\(k\leq5000\)，而且模数还不是ntt模数，于是直接递推就好了

递推的边界条件是\(\begin{Bmatrix} 0 \\0\end{Bmatrix}=1\)

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define re register
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
int n,k;
const int mod=1e9+7;
const int inv=500000004;
inline int ksm(int a,int b) {
    int S=1;
    while(b) {if(b&1) S=1ll*S*a%mod;b>>=1;a=1ll*a*a%mod;}
    return S;
}
int S[2][5005];
int main() {
    scanf("%d%d",&n,&k);
    S[0][1]=1;int o=0;
    for(re int i=2;i<=k;i++,o^=1) 
        for(re int j=1;j<=i;j++)
            S[o^1][j]=(S[o][j-1]+1ll*S[o][j]*j%mod)%mod;
    int ans=0,now=ksm(2,n-1),t=n,tot=n-1;
    for(re int i=1;i<=min(n,k);i++,--tot) {
        ans=(ans+1ll*S[o][i]*t%mod*now%mod)%mod;
        t=1ll*t*tot%mod;now=1ll*now*inv%mod;
    }
    printf("%d\n",ans);
    return 0;
}

转载于:https://www.cnblogs.com/asuldb/p/10864777.html