1076: [SCOI2008]奖励关

Description

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

Input

  第一行为两个正整数k和n,即宝物的数量和种类。以下n行分别描述一种宝物,其中第一个整数代表分值,随
后的整数依次代表该宝物的各个前提宝物(各宝物编号为1到n),以0结尾。

Output

  输出一个实数,保留六位小数,即在最优策略下平均情况的得分。

Sample Input

1 2
1 0
2 0

Sample Output

1.500000

HINT

 

【数据规模】

1<=k<=100,1<=n<=15,分值为[-10^6,10^6]内的整数。

 
 
状压DP。。。
http://lvat2000.is-programmer.com/posts/196742.html
 1 #include<iostream>
 2 #include<cstdlib>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<cstdio>
 6 #include<algorithm>
 7 #include<string>
 8 #include<map>
 9 #include<queue>
10 #include<vector>
11 #include<set>
12 #define inf 1000000000
13 #define maxn 10000+5
14 #define maxm 10000+5
15 #define eps 1e-10
16 #define ll long long
17 #define for0(i,n) for(int i=0;i<=(n);i++)
18 #define for1(i,n) for(int i=1;i<=(n);i++)
19 #define for2(i,x,y) for(int i=(x);i<=(y);i++)
20 #define for3(i,x,y) for(int i=(x);i>=(y);i--)
21 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)
22 using namespace std;
23 int read(){
24     int x=0,f=1;char ch=getchar();
25     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
26     while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
27     return x*f;
28 }
29 double f[101][65536];
30 int n,k,t;
31 int v[20],d[20],p[20];
32 int main(){
33     //freopen("input.txt","r",stdin);
34     //freopen("output.txt","w",stdout);
35     for1(i,16)p[i]=1<<(i-1);
36     n=read();k=read();
37     for1(i,k){
38         v[i]=read();t=read();
39         while(t){
40             d[i]+=p[t];
41             t=read();
42         }
43     }
44     for(int i=n;i;i--)
45         for0(j,p[k+1]-1){
46             for1(l,k)
47                 if((d[l]&j)==d[l])
48                     f[i][j]+=max(f[i+1][j],f[i+1][j|p[l]]+v[l]);
49                 else f[i][j]+=f[i+1][j];
50             f[i][j]/=k;
51         }
52     printf("%.6lf",f[1][0]);
53     return 0;
54 }
View Code

 

posted @ 2016-05-25 21:08  HTWX  阅读(152)  评论(0编辑  收藏  举报