【BZOJ2339】卡农(递推,容斥)

【BZOJ2339】卡农(递推,容斥)

题面

BZOJ

题解

先简化一下题意:
\([1,2^n-1]\)中选择不重复的\(m\)个数,使得他们异或和为\(0\)的方案数。

我们设\(f[i]\)表示选择\(i\)个数异或和为\(0\)的方案数。
直接算是很麻烦的,所以我们反过来,总数减去不合法的。
因为确定了前\(i-1\)个数最后一个数就已经知道了。
所以总方案数是\(A_{2^n-1}^{i-1}\),不合法的有两种,一种是选择了\(0\),一种是有重复。
选择了\(0\),意味着前\(i-1\)个数的异或和为\(0\),所以方案数是\(f[i-1]\)种。
有重复,我们枚举哪个数重复了,那么剩下的\(i-2\)个数的异或和仍然为\(0\)
所以方案数是\(f[i-2]\times (2^n-1-(i-2))\),题目没有考虑顺序,但是我们计算的时候先考虑了顺序,所以这里的方案数还需要考虑在哪个位置,也就是再乘上一个\((i-1)\)
所以$$f[i]=A_{2n-1}{i-1}-f[i-1]-(i-1)\times f[i-2]\times(2^n-1-(i-2))$$
最终的答案再把顺序的问题处理一下就好了。

#include<iostream>
#include<cstdio>
using namespace std;
#define MOD 100000007
#define MAX 1000100
inline int read()
{
	int x=0,t=1;char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}
int jv=1,A[MAX],n,m,p,f[MAX];
int fpow(int a,int b)
{
	int s=1;
	while(b){if(b&1)s=1ll*s*a%MOD;a=1ll*a*a%MOD;b>>=1;}
	return s;
}
int main()
{
	n=read();m=read();p=fpow(2,n);f[0]=A[0]=1;f[1]=0;
	for(int i=1;i<=m;++i)jv=1ll*jv*i%MOD;jv=fpow(jv,MOD-2);
	for(int i=1;i<=m;++i)A[i]=1ll*A[i-1]*(p-i+MOD)%MOD;
	for(int i=2;i<=m;++i)f[i]=((A[i-1]-f[i-1]+MOD)%MOD-1ll*f[i-2]*(i-1)%MOD*(p-1-(i-2)+MOD)%MOD+MOD)%MOD;
	printf("%lld\n",1ll*jv*f[m]%MOD);
	return 0;
}

posted @ 2018-07-02 10:49  小蒟蒻yyb  阅读(406)  评论(0编辑  收藏  举报