原题链接也许没寄吧

题目描述

作为一位喜欢猫的人，小 B 家里养了 $n$ 只猫。这些猫从 $1$ 到 $n$ 标号。

他每天最喜欢的时候就是撸猫，可惜猫猫并不是每天都愿意被他撸。愿意被他撸的猫猫集合 $\subseteq [n]$ 根据某个概率分布 $P$ 随机而来。在知道了哪些猫猫愿意被撸之后，小B可以选择某一只愿意的来 rua。注意，如果 $S$ 集合为空，则小 B 不能撸任何一只猫。同时，小 B 的撸猫策略可以带有随机性。比如,在知道了 1,2 号猫猫都愿意被撸之后，他可以以 $50\%$ 的概率撸第一只猫，以另外 $50\%$ 的概率撸第二只猫。

作为一位公平的人，小 B 希望给每一只猫猫更多被 rua 的机会。他希望设计一个撸猫的策略来最大化 $\operatorname{Pr}[i_{\text{被撸}} | i\in S]$ 的最小值。换句话说，令 $p_i = \operatorname{Pr}[i\in S]$ 是第 $i$ 只猫猫愿意被撸的概率。我们需要设计一个撸猫的策略来最大化常数 $c$ ，使得每只猫被撸的概率至少是 $cp_i$ 。

输入格式

一行一个整数 $n$ 。
一行 $2^n$ 个浮点数，表示第 $i$ 种可能情况的概率。

输出格式

一行一个浮点数 $c$ ，表示答案。

样例

样例输入 1

2
0.25 0.25 0.25 0.25

样例输出 1

0.7500000000

样例输入 2

3
0.12 0.18 0.02 0.1 0.1 0.05 0.15 0.28

样例输出 2

0.5057471264

数据范围与提示

Subtask 1 [12 pts]： $1\leqslant n\leqslant 2$ 。
Subtask 2 [18 pts]： $1\leqslant n\leqslant 7$ 。
Subtask 3 [10 pts]： $1\leqslant n\leqslant 10$ 。
Subtask 4 [35 pts]： $1\leqslant n\leqslant 15$ 。
Subtask 5 [25 pts]： $1\leqslant n\leqslant 20$ 。

题解

听永神说是脑瘫题。

考虑二分答案，那么问题就转化为一个二分图网络流：左边是 $2^n$ 个点表示每种情况，入边容量为这种情况的概率，右边是 $n$ 个点表示每只猫，出边容量为 $c*p_i$ ，中间对应连上容量 $+\infty$ 的边。答案合法当且仅当图的右边满流。

我们可以直接套用 $\rm Hall$ 定理来判断是否可以满流，这样判断一次是 $O(2^n)$ ，二分答案可以直接通过。

当然你会发现完全没必要二分答案，直接 $\rm Hall$ 定理枚举时取个最小值即可。

代码

我的还是二分答案

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define lll __int128
#define uns unsigned
#define fi first
#define se second
#define IF (it->fi)
#define IS (it->se)
#define END putchar('\n')
#define lowbit(x) ((x)&-(x))
#define inline jzm
using namespace std;
const int MAXN=114514;
const ll INF=1e18;
ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[50],lpt;
void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt>0)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}
const double eps=1e-10;
int n;
double p[25],f[1<<20],g[1<<20],sum;
bool check(double c){
	for(int s=0;s<(1<<n);s++)g[s]=0;
	for(int i=1;i<=n;i++)
		for(int s=0,lim=1<<(i-1);s<lim;s++)g[s^lim]=g[s]+p[i]*c;
	for(int s=0,m=(1<<n)-1;s<=m;s++)if(g[s]>=sum-f[s^m]+eps)return 0;
	return 1;
}
int main()
{
	n=read();
	for(int s=0;s<(1<<n);s++)
		*new(int)=scanf("%lf",f+s),sum+=f[s];
	for(int s=0;s<(1<<n);s++)
		for(int i=1;i<=n;i++)if((s>>(i-1))&1)p[i]+=f[s];
	for(int k=1;k<(1<<n);k<<=1)
		for(int s=0;s<(1<<n);s++)if(s&k)f[s]+=f[s^k];
	double l=0,r=1,mid=0.5;
	for(int NND=50;NND--;){
		mid=(l+r)/2;
		if(check(mid))l=mid;
		else r=mid;
	}
	printf("%.10f\n",mid);
	return 0;
}