HDU4336：Card Collector(min-max容斥)

题面

传送门

Sol

方法一

直接状压就好了

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;

int n;
double p[21], f[1 << 20];

int main(RG int argc, RG char *argv[]){
	while(scanf("%d", &n) != EOF){
		for(RG int i = 1; i <= n; ++i) scanf("%lf", &p[i]);
		RG int S = 1 << n; Fill(f, 0);
		for(RG int i = S - 2; ~i; --i){
			RG double s = 1.0;
			for(RG int j = 1; j <= n; ++j) if(~i & (1 << (j - 1))) s -= p[j];
			s = 1.0 - s, f[i] = 1.0 / s;
			for(RG int j = 1; j <= n; ++j)
				if(~i & (1 << (j - 1))) f[i] += p[j] * f[i | (1 << (j - 1))] / s;
		}
		printf("%.6lf\n", f[0]);
	}
    return 0;
}

方法二

方法一实在太水了，显然不是重点
下面介绍一种容斥方法

min-max容斥

\(E[max(S)]=\sum(-1)^{k+1}E[min(S')]\)
其中集合\(S'\subseteq S\)，\(k=|S'|\)

\(max(S)\)指的是这个集合内最后出现的元素
\(min(S)\)指的是这个集合内最先出现的元素
\(E\)表示期望第几步出现

具体到这个题
就是要求\(E[max(\)全集\()]\)
\(E[min(S')]\)就是指\(S'\)任意出现一个的期望步数
举例：
每步出现\(x_i\)的概率\(p_i\)
\(min\{x_1, x_2, x_3 \}\)的概率就是\(p_1+p_2+p_3\)
期望步数就是\(\frac{1}{p_1+p_2+p_3}\)

然后就容斥一下这个题

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;

int n;
double p[21], ans;

IL void Dfs(RG int x, RG double E, RG int op){
	if(x > n){
		if(E > 1e-7) ans += 1.0 * op / E;
		return;
	}
	Dfs(x + 1, E, op);
	Dfs(x + 1, E + p[x], -op);
}

int main(RG int argc, RG char *argv[]){
	while(scanf("%d", &n) != EOF){
		for(RG int i = 1; i <= n; ++i) scanf("%lf", &p[i]);
		ans = 0, Dfs(1, 0, -1);
		printf("%.6lf\n", ans);
	}
    return 0;
}