2339: [HNOI2011]卡农
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 842 Solved: 510
[Submit][Status][Discuss]
Description
可以把集合视作有序的,当做排列做,最后再 /m!
设f[i]表示选出i个集合的合法方案
选出了(i-1)个集合后,最后一个集合是唯一确定的
总数就是A(2^n - 1,i-1)
但是最后确定的集合可能使方案不合法,有两种情况
1.最后确定的集合为空,这种情况的方案数=f[i-1]
2.最后确定的集合和之前确定的集合重复,因为有重复,所以删去这两个重复的集合,
依旧满足所有元素出现偶数次的性质, 这种情况的方案数 =f[i-2]*(2^n-1-(i-2))
ans就可以计算了
还有一种理解方式,理解成无序的,用组合搞
推荐blog http://blog.csdn.net/dflasher/article/details/51615325
#include<iostream>
#include<cstdio>
#define N 1000010
#define P 100000007
using namespace std;
long long n,m,p[N],f[N],temp;
long long power(long long a,long long b)
{
long long ans(1);
for(long long i=b;i;i>>=1,(a*=a)%=P) if(i&1)(ans*=a)%=P;
return ans;
}
void pre()
{
p[0]=1;
for (long long i=1;i<=m;i++) p[i]=(p[i-1]*((temp-i+1+P)%P))%P;
}
int main()
{
scanf("%lld%lld",&n,&m);
temp=power(2,n);temp--;
if (temp<0) temp+=P;
pre();
for (long long i=3;i<=m;i++)
f[i]=((p[i-1]-f[i-1]-f[i-2]*(i-1)%P*(temp-(i-2))%P)+P)%P;
temp=1;
for (long long i=1;i<=m;i++) (temp*=i)%=P;
(f[m]*=power(temp,P-2))%=P;
cout<<f[m];
}