P3214 [HNOI2011]卡农

题目

P3214 [HNOI2011]卡农

在被一题容斥(dp)完虐之后，打算做一做集合容斥这类的题了

第一次深感HNOI的毒瘤(题做得太少了！！)

做法

求([1,n])组成的集合中选(m)个不同集合且每个元素出现偶数的组合方案

无序(打乱顺序仍记为一种)通常我们对于有序的做法更简单，怎么转换呢

C组合数的公式是怎么得来的？别说你是背来的(emmm)(那也没有做这题的必要了)

有序(m!)就得到了无序的

我们考虑(dp)，数组(dp_i)表示选i个不同集合的排列方案

异或和为(0)，则，确定前(i-1)个集合则第(i)个集合自然也出来了，方案数为(A_{2^n-1}^{i-1})

如果前面(i-1)个集合异或和已为(0)，那第(i)个集合为空集，不符题意，这部分的方案数就是(dp_{i-1})

保证所选集合不重复，若(i)与前(i-1)任意重复，去掉这个重复的集合，为(dp_{i-2})，可能的位置有((i-1))个，重复集合个数有((2^n-1-(i-2)))

(dp_i=A_{2^n-1}^{i-1}-dp_{i-1}-dp_{i-2}*(i-1)*(2^n-i+1))

最后再乘下逆元就好了

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
const LL p=100000007;
const int maxn=1e6+9;
inline LL Pow(LL base,LL b){
	LL ret(1);
	while(b){
		if(b&1)
		    ret=ret*base%p;
		base=base*base%p,
		b>>=1;
	}
	return ret;
}
LL n,m,a,Up,A,ans;
LL dp[maxn];
int main(){
	scanf("%lld%lld",&n,&m);
	dp[1]=dp[2]=0,
	Up=(Pow(2ll,n)-1ll+p)%p,
	A=Up;
	for(LL i=3;i<=m;++i)
		A=A*(Up-i+2)%p,
		dp[i]=((A-dp[i-1]+p)%p-dp[i-2]*(i-1)%p*((Up-(i-2)+p)%p)%p +p)%p;
	a=1;
	for(LL i=2;i<=m;++i)
	    a=a*i%p;
	ans=dp[m]*Pow(a,p-2)%p;
	printf("%lld
",ans);
	return 0;
}/*
100 1000
*/