在做概率dp突然发现一道容斥的水体:
简答考虑两个集合的时候,E(1) = 1/p1 , E(2) = 1/p2 ,直接相加会出现一个问题,在没抽到1的时候抽到了2,只要减一个E(1,2) = 1/(p1+p2); 这成立了就可以地推n个集合也是如此的.
#include<cstdio>
#include<iostream>
const double eps = 1e-10;
const int maxn = 3e2;
using namespace std;
double p[maxn];
int main(){
int n;
while(~scanf("%d",&n)){
bool flag = 0;
for(int i = 0; i < n; i++){
scanf("%lf",&p[i]);
}
double ans = 0;
double temp;int cc;
for(int i = 1; i < (1<<n); i++){
temp = 0;
cc = 0;
for(int j = 0; j < n; j++){
if(i&(1<<j)){
temp+=p[j];
cc++;
}
}
if(cc&1)ans+=1.0/temp;
else ans -= 1.0/temp;
}
printf("%.6lf\n",ans);
}
return 0;
}