bzoj2339[HNOI2011]卡农 dp+容斥

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];
}