【[SCOI2008]奖励关】

又抄了一篇题解

要凉了要凉了，开学了我还什么都不会

文化课凉凉，NOIP还要面临爆零退役的历史进程

这道题挺神的，期望+状态压缩

我们设(dp[i][S])表示在第(i)天前，捡的宝物状态为(S)到第(K)天结束的期望收益是多少

于是我们的答案是(dp[1][0])，也就是第一天前（就是没开始）什么宝物也没有到结束的期望是多少

期望倒着推，我们的初始状态就是(dp[k+1]={0})

之后对于一个状态(dp[i][S])

在第(i)天(m)种宝物出现的概率都是(1/m)

所以枚举每一种宝物(j)，如果这个宝物的前提条件被(S)包含

就有

[dp[i][S]+=max(dp[i+1][S|(1<<(j-1))]+val[j],dp[i+1][S])*1/m ]

就是就算可以选择这个宝物的话，我们也两种选择拿或者不拿这个宝物

如过没有被包含，你只能不拿这个宝物

[dp[i][S]+=dp[i+1][S]*1/m ]

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#define re register
#define max(a,b) ((a)>(b)?(a):(b))
inline int read()
{
	char c=getchar();
	int x=0,r=1;
	while(c<'0'||c>'9') 
	{
		if(c=='-') r=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9')
	  x=(x<<3)+(x<<1)+c-48,c=getchar();
	return x*r;
}
double dp[102][32769];
int n,m,N;
int cost[16];
int f[16];
int main()
{
	n=read(),m=read();
	for(re int i=1;i<=m;i++)
	{
		cost[i]=read();
		int x=read(),y=0;
		while(x) y|=1<<(x-1),x=read();
		f[i]=y;
	}
	N=(1<<m)-1;
	for(re int i=0;i<=N;i++) dp[n+1][i]=0;
	for(re int i=n;i;i--)
	for(re int j=0;j<=N;j++)
	{
		for(re int k=1;k<=m;k++)
		if(((f[k]|j)==j))
		{
			dp[i][j]+=max(dp[i+1][j],dp[i+1][j|(1<<(k-1))]+cost[k])/double(m);
		}else dp[i][j]+=dp[i+1][j]/double(m);
	}
	printf("%.6lf",dp[1][0]);
	return 0;
}