CF932E Team Work

luogu
{ } \{\} {}为第二类 s t r i n g string string
∑ i = 1 n C ( n , i ) ∗ i k \sum_{i=1}^{n}C(n,i)*i^k i=1nC(n,i)ik
∑ i = 1 n C ( n , i ) ∑ j = 0 k j ! ∗ { k , j } ∗ C ( i , j ) \sum_{i=1}^{n}C(n,i)\sum_{j=0}^{k}j!*\{k,j\}*C(i,j) i=1nC(n,i)j=0kj!{k,j}C(i,j)
∑ i = 1 n ∑ j = 0 k C ( n , i ) ∗ C ( i , j ) ∗ j ! ∗ { k , j } \sum_{i=1}^{n}\sum_{j=0}^{k}C(n,i)*C(i,j)*j!*\{k,j\} i=1nj=0kC(n,i)C(i,j)j!{k,j}
∑ j = 0 k j ! ∗ { k , j } ∑ i = 1 n C ( n , i ) ∗ C ( i , j ) \sum_{j=0}^{k}j!*\{k,j\}\sum_{i=1}^{n}C(n,i)*C(i,j) j=0kj!{k,j}i=1nC(n,i)C(i,j)
∑ j = 0 k j ! ∗ { k , j } ∑ i = 1 n C ( n , j ) ∗ C ( n − j , i − j ) \sum_{j=0}^{k}j!*\{k,j\}\sum_{i=1}^{n}C(n,j)*C(n-j,i-j) j=0kj!{k,j}i=1nC(n,j)C(nj,ij)
∑ j = 0 k j ! ∗ { k , j } ∑ i = 1 n C ( n , j ) ∗ C ( n − j , i − j ) \sum_{j=0}^{k}j!*\{k,j\}\sum_{i=1}^{n}C(n,j)*C(n-j,i-j) j=0kj!{k,j}i=1nC(n,j)C(nj,ij)
∑ j = 0 k j ! ∗ { k , j } ∗ C ( n , j ) ∑ i = 1 n C ( n − j , i − j ) \sum_{j=0}^{k}j!*\{k,j\}*C(n,j)\sum_{i=1}^{n}C(n-j,i-j) j=0kj!{k,j}C(n,j)i=1nC(nj,ij)
∑ j = 0 k j ! ∗ { k , j } ∗ C ( n , j ) ∗ 2 n − j \sum_{j=0}^{k}j!*\{k,j\}*C(n,j)*2^{n-j} j=0kj!{k,j}C(n,j)2nj

#include <bits/stdc++.h>
using namespace std;
const int mo = 1e9+7;
typedef long long LL;
#define int LL
const int maxk = 5e3+5;
int readint(){
	int x=0,f=1;char s=getchar();
	#define sc (s=getchar())
	while(s<'0'||s>'9'){
		if(s=='-')
			f=-1;
		sc;
	}
	while(s>='0'&&s<='9'){
		x=(x<<3)+(x<<1)+(s^48);
		sc;
	}
	#undef sc
	return x*f;
}
int C[maxk];
int fac[maxk];
int s[maxk][maxk];
int qkpow(int a,int b){
	int ans=1;
	while(b){
		if(b&1)
			ans=ans*a%mo;
		b>>=1;
		a=a*a%mo;
	}
	return ans;
}
void init(int n,int k){
	C[1]=n;C[0]=1;
	for(int i=2;i<=k;i++)
		C[i]=C[i-1]*(n-i+1)%mo*qkpow(i,mo-2)%mo;
	s[0][0]=1;
	for(int i=1;i<=k;i++){
		for(int j=1;j<=i;j++){
			s[i][j]=(s[i-1][j]*j%mo+s[i-1][j-1])%mo;
		}
	}
	fac[1]=1;
	for(int i=2;i<=k;i++){
		fac[i]=fac[i-1]*i%mo;
	}
}
signed main (){
	int n=readint(),k=readint();
	init(n,k);
	int pow_2=qkpow(2,n);
	int be=2;
	int ans=0;
	for(int j=1;j<=k;j++){
		if(j>n)
			break;
		ans=(ans+fac[j]*s[k][j]%mo*C[j]%mo*pow_2%mo*qkpow(be,mo-2))%mo;
		be=be*2%mo;
	}
	cout<<ans<<endl;
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值