BZOJ 2839: 集合计数

最新推荐文章于 2019-03-02 10:29:42 发布

weixin_33743703

最新推荐文章于 2019-03-02 10:29:42 发布

阅读量68

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/silenty/p/9898570.html

本文探讨了容斥原理在组合数学中的应用，通过具体的算法实现，展示了如何使用容斥原理解决复杂的问题。代码示例中，利用了C++语言进行编程，详细解释了mi、ni、tmp数组的作用，以及如何通过递推公式计算特定条件下可能的组合数量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

容斥

#include<cstdio>
using namespace std;
const int mod=1e9+7;
int mi[1000005],ni[1000005],tmp[1000005];
int C(int n,int m){
	return 1ll*mi[n]*ni[m]%mod*ni[n-m]%mod;
}
int main(){
	int n,k;
	scanf("%d%d",&n,&k);
	tmp[n]=1;
	for (int i=n-1; i>=k; i--) tmp[i]=1ll*tmp[i+1]*(tmp[i+1]+2)%mod;
	mi[0]=1;
	for (int i=1; i<=n; i++) mi[i]=1ll*mi[i-1]*i%mod;
	ni[0]=1;
	ni[1]=1;
	for (int i=2; i<=n; i++) ni[i]=1ll*ni[mod%i]*(mod-mod/i)%mod;
	for (int i=1; i<=n; i++) ni[i]=1ll*ni[i-1]*ni[i]%mod;
	int ans=0,f=1;
	for (int i=k; i<=n; i++){
		(ans+=1ll*f*C(i,k)%mod*C(n,i)%mod*tmp[i]%mod)%=mod;
		(ans+=mod)%=mod;
		f=-f;
	}
	printf("%d\n",ans);
	return 0;
}

转载于:https://www.cnblogs.com/silenty/p/9898570.html