[SCOI2008]奖励关 状压dp+期望

原题:https://www.luogu.org/problemnew/show/P2473

题解:很容易想到一个方程，f[i][S]，第i轮 选S的期望值，转移应该是第i-1轮能转移的和的平均值，不过不好判断转移的条件，所以我们考虑倒序枚举。即f[i][S]表示 1-i-1轮S状态，i-k的期望值。这样就能用&运算来判断。

即当能转移时 即a[j]&S ==a[j] S中包含a[j](a[j] 表示第i种宝物所需要的集合压成二进制的值)

f[i][S]=sum{ max{f[i+1][s]  , f[i+1][S|1<<(j-1)]}  }/n   S∈[0,(1<<n)-1] j∈[1,n]

不能转移时

f[i][S]=sum{ f[i+1][S]}/n

答案就是 f[1][0]

#include<bits/stdc++.h>
using namespace std;
const int N=110;
const int M=55000;
int n,k,maxn;
int a[N],p[N];
double f[N][M];//表示在1-i-1轮里是否S i-k的期望 
inline int rd(){
	int x=0;int f=1;char s=getchar();
	while(!isdigit(s)) f=(s=='-'?-1:f),s=getchar();
	while(isdigit(s)) x=(x<<1)+(x<<3)+s-'0',s=getchar();
	return x*f;
}
int main(){
	k=rd();n=rd();
	for(int i=1,x;i<=n;i++){
		p[i]=rd();
		while(x=rd())	a[i]|=(1<<(x-1));
	}
	maxn=(1<<n)-1;
	for(int i=k;i;i--){//倒序枚举 
		for(int s=0;s<=maxn;s++){
			for(int j=1;j<=n;j++){//枚举宝藏 
				if((a[j]&s)==a[j]) f[i][s]+=max(f[i+1][s],f[i+1][s|(1<<(j-1))]+p[j]);
				else f[i][s]+=f[i+1][s]; 
			}
			f[i][s]/=n;
		}
	}
	printf("%.6lf\n",f[1][0]);
	
	return 0;
}