[bzoj 1076][SCOI2008]奖励关

传送门

Description

你正在玩你最喜欢的电子游戏,并且刚刚进入一个奖励关。在这个奖励关里,系统将依次随机抛出k次宝物, 每次你都可以选择吃或者不吃(必须在抛出下一个宝物之前做出选择,且现在决定不吃的宝物以后也不能再吃)。  宝物一共有n种,系统每次抛出这n种宝物的概率都相同且相互独立。也就是说,即使前k-1次系统都抛出宝物1( 这种情况是有可能出现的,尽管概率非常小),第k次抛出各个宝物的概率依然均为1/n。 获取第i种宝物将得到Pi 分,但并不是每种宝物都是可以随意获取的。第i种宝物有一个前提宝物集合Si。只有当Si中所有宝物都至少吃过 一次,才能吃第i种宝物(如果系统抛出了一个目前不能吃的宝物,相当于白白的损失了一次机会)。注意,Pi可 以是负数,但如果它是很多高分宝物的前提,损失短期利益而吃掉这个负分宝物将获得更大的长期利益。 假设你 采取最优策略,平均情况你一共能在奖励关得到多少分值?

Solution

\(n \leq 15\),可以对取到的种类集合状压

考虑逆推,\(f[i][S]\)表示到第\(i\)个位置,状态为\(S\),走到第\(k\)步的期望

\[f[i][S] +=max(f[i+1][S],f[i+1][S \cup j]+p[j]) \ \ \ j所需要的集合\subseteq S, \]

\[f[i][S]+ f[i+1][S] \ \ \ j所需要的集合 \nsubseteq S \]


Code 

#include<bits/stdc++.h>
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return x*f;
}
#define MN 105
int k,n,p[16],x,s[16];
double f[2][1<<15];
int main()
{
	k=read();n=read();
	register int i,S,j;
	for(i=0;i<n;++i)
	{
		p[i]=read();
		x=read();
		while(x)
		{
			s[i]|=1<<x-1;
			x=read();
		}
	}
	for(i=k;i>=1;--i)for(S=(1<<n)-1;~S;--S)
	{
		f[i&1][S]=0.;
		for(j=0;j<n;j++)
			f[i&1][S]+=(s[j]|S)==S?max(f[(i&1)^1][S],f[(i&1)^1][S|(1<<j)]+p[j]*1.):f[(i&1)^1][S];
		f[i&1][S]/=n*1.;
	}
	printf("%.6lf\n",f[1][0]);
}


Blog来自PaperCloud,未经允许,请勿转载,TKS!

posted @ 2019-01-03 13:21  PaperCloud  阅读(125)  评论(0编辑  收藏  举报