题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4336
题目大意:每包中含有这些卡片的概率为p1,p2,……pN. 每包至多一张卡片,可能没有卡片。问需要买多少包才能集齐N张卡片,求次数的期望。
题目分析1:使用状态压缩表示卡片的收集情况,bin[i]=j表示在i状态的情况下,集齐卡片需要买东西次数的期望。
由于每包中有三种情况:
1. 没有卡片
2. 卡片是已经收集到的
3. 卡片是没有收集到的
得到公式:
bin[s] = 1 + (1-Σp[i])*bin[s] +Σ(p[j]*bin[s]) + Σ(p[k]*bin[s|(1<<k)])
->Σp[i]*bin[s] = 1 +Σ(p[i]*bin[s|(1<<i)
最后要求的是bin[0]。
代码参考1:
#include<queue>
#include<vector>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN = 1<<20;//最多20张卡片,所以会有1<<20种状态
double p[MAXN], bin[MAXN];
int main()
{
int i, j, k, n, m, t;
while(~scanf("%d", &n))
{
memset(bin, 0, sizeof(bin));
for(i=0; i<n; ++i)
{
scanf("%lf", &p[i]);
}
int all = (1<<n) - 1;//all表示全部集满时的状态
for(i=all-1; i>=0; --i)//枚举从哪个状态推过来的
{
double temp = 0.0;
bin[i] = 1;//如果是从此状态推过来的,那此状态必然存在,出现的概率为1
for(j=0; j<n; ++j)//枚举每一位上是否有1,即第j张卡片是否已经被收集
{
if(i & (1<<j)) continue;//如果第j张卡片已经收集到了则不再考虑
bin[i] += bin[i|(1<<j)] * p[j];//如果第j张卡片没有收集到过,则加入状态,第j位改为1,乘以第j张卡片出现的概率
temp += p[j];//所有有效概率之和
}
bin[i] /= temp;//除以概率总和则为从i状态推过来时所得期望
}
printf("%lf\n", bin[0]);
}
return 0;
}
====================================================================================================================================
题目分析2:
容斥原理。
代码参考2:
#include<cstdio>
const int N = 2000;
double p[N];
int main()
{
int i, j, k, n, m;
while(~scanf("%d", &n))
{
for(i=0; i<n; ++i)
{
scanf("%lf", &p[i]);
}
double sum = 0.0;
for(i=1; i<(1<<n); ++i)
{
int cnt = 0;
double cur = 0;
for(j=0; j<n; ++j)
{
if(i&(1<<j))
{
cnt++;
cur += p[j];
}
}
if(cnt&1) sum += 1.0/cur;//奇加偶减,根据公式计算即可
else sum -= 1.0/cur;
}
printf("%lf\n", sum);
}
return 0;
}