[BZOJ1076][SCOI2008]奖励关(概率DP)

最新推荐文章于 2019-10-01 10:26:38 发布

转载最新推荐文章于 2019-10-01 10:26:38 发布 · 59 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/void-f/p/9129442.html

本文介绍了一种使用动态规划算法解决特定概率问题的方法，通过计算每个状态转移的期望值来达到最终目标。代码中详细展示了如何初始化数据结构、读取输入参数，并递归地更新动态规划表。

Code

#include <cstdio>
#include <algorithm>
#include <cstring>
#define N 110
#define db double 
using namespace std;

int n,m,nd[N],v[N];
db f[N][1<<16];

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*10+ch-'0';ch=getchar();}
    return x*f;
}

db DP(int i,int S){
	if(i>m) return 0;
	db &tmp=f[i][S];
	if(tmp+1e9>1e-9) return tmp;
	tmp=0;
	for(int j=1;j<=n;++j)
		if((S|nd[j])==S)
			tmp+=max(DP(i+1,S),DP(i+1,S|(1<<(j-1)))+v[j]);
		else tmp+=DP(i+1,S);	
	tmp/=n;
	return tmp;
}

int main(){
	m=read(),n=read();
	for(int i=1;i<=n;++i){
		v[i]=read();
		for(int x=read();x!=0;x=read()) nd[i]|=(1<<(x-1));
	}
	for(int i=0;i<=m;++i)for(int j=0;j<(1<<n);++j)f[i][j]=-2e9;
	printf("%.6lf\n",DP(1,0));
}

转载于:https://www.cnblogs.com/void-f/p/9129442.html