bzoj2246: [SDOI2011]迷宫探险

2246: [SDOI2011]迷宫探险

Time Limit: 10 Sec  Memory Limit: 512 MB
Submit: 202  Solved: 118
[Submit][Status][Discuss]

Description

这是一个单人游戏。

游戏开始时,玩家控制的人物出生在迷宫的某个位置,玩家的目标是控制人物走到迷宫的某个出口(出口可能有多个)。迷宫里有k类陷阱(用“A”、“B”、“C”……表示,相同字母代表相同类型的陷阱),每类陷阱可能是有害的或无害的,而在游戏开始时玩家并不知道哪些陷阱是有害的,哪些是无害的。同一类陷阱的状态相同,即用同一个字母标志的陷阱要么全部有害,要么全部无害,不会发生一部分有害而另一部分无害的情况。任何陷阱状态的组合都有一个发生概率,考虑下例:

当k=2时,迷宫内共有两类陷阱,分别用“A”和“B”表示,陷阱状态的组合共有4种:

1、“A”是无害陷阱,“B”是无害陷阱。

2、“A”是有害陷阱,“B”是无害陷阱;

3、“A”是无害陷阱,“B”是有害陷阱;

4、“A”是有害陷阱,“B”是有害陷阱;

下面给出了一个合法的概率表格:


 

“A”是无害陷阱

“A”是有害陷阱

“B”是无害陷阱

36%

24%

“B”是有害陷阱

24%

16%

当K=3时,会有8种不同的陷阱状态组合,如果我们依然坚持使用概率表格,那么这个表格将会是三维的(2*2*2,每一维对应着一类陷阱)。当K≥3时,这将使得题目难以描述。因此我们使用一个大小为2k的数组p来描述每种情况发生的可能性,p的下标范围为0~2k-1。

p是这样生成的:

对于每个可能的陷阱状态组合,考虑所有k类陷阱,令1表示某个陷阱有害,0表示某个陷阱无害,把“A”作为二进制数的第0位(从右边开始计数),“B”作为第1位,“C”作为第2位……通过以上操作,我们可以得到一个K位的二进制数,把它转化成十进制后,2k种陷阱状态的组合将会与整数0~2k-1一一对应。

定义S表示P中所有元素和,即:

则陷阱状态组合i出现的概率为P/S。上述表格对应的一个合法数组P是:

P0 =36

P=24

P=24

P3=16

当然同一个概率表格可能会对应多个数组P(事实上有无数个数组P能够迎合表格数据),例如上述表格同时也对应着下面的数组P:

P0=72

P1=48

P2=48

P3=32

玩家控制的人物初始情况下有H点生命,当人物踏上某个陷阱时,如果这个陷阱是有害的,那么会损失1点生命,否则这个陷阱是无害的,不损失生命。无论上述哪种情况发生,玩家会立刻得到这个陷阱的信息(有害或无害)。一旦生命小于等于0,玩家控制的人物会立刻死亡。

迷宫可以看作m*n的方格地图,每个元素可能是:

“.”:表示这是平地,可以通过;

“#”:表示这是墙,不能通过;

“A”,“B”,“C”……:表示这是一个陷阱;

“$”:表示这是起点,地图中有且仅有一个;

“@”:表示这是终点,地图中可以有多个,也可以一个也没有。

人物可以向上下左右四个方向行走,不可以走对角线,也不可以走出地图。

给定m*n的地图、k、h以及大小为2k的概率数组。你的任务是求出在执行最优策略时,人物能活着走出迷宫的概率。

Input

第一行包含4个整数,分别表示m、n、K、H;

下面m行每行n个字符描述迷宫地图;

最后一行包含2K个非负整数描述数组P,数组下标从0开始。

Output

仅包含一个数字,表示在执行最优策略时,人物活着走出迷宫的概率。四舍五入保留3位小数。

Sample Input

仅包含一个数字,表示在执行最优策略时,人物活着走出迷宫的概率。四舍五入保留3位小数。

【样例输入1】

4 3 2 1

.$.

A#B

A#B

[email protected]

30 30 20 20

【样例输出1】

0.600

【样例说明1】

向右边走,经过“B”,“B”为有害陷阱的概率为 (20+20)/(30+30+20+20=)=0.4,若“B”为有害陷阱那么人物就死掉了,游戏失败,否则玩家得知“B”是无害陷阱,继续经过另一个“B”达到终点,胜利的概率为0.6。

【样例输入2】

4 3 2 2

.$.

A#B

A#B

[email protected]

30 30 20 20

【样例输出2】

0.800

【样例说明2】

向左边走,经过“A”,“A”为有害陷阱的概率为 (30+30)/(30+30+20+20)=0.5。若“A”为有害陷阱,那么损失一点生命,转到右边尝试“B”,要想成功到达终点,此时“B”必须为无害陷阱,而在“A”是有害陷阱的前提下,“B”是无害陷阱的概率是30/(30+20)=0.6,故这种情况发生的概率为0.5*0.6=0.3。若“A”是无害陷阱,玩家可以控制人物连续通过两个“A”到达终点,这种情况发生的概率0.5。所以答案为0.3+0.5=0.8。

【样例输入3】

4 3 2 3

.$.

A#B

A#B

[email protected]

30 30 20 20

【样例输出3】

1.000

【样例说明3】

玩家控制的人物有3点生命,但最多只需要经过两个陷阱,所以任意选左路或右路走过去就可以到达终点了。

【样例输入4】

4 3 3 2

.$.

A#B

A#C

@@@

143 37 335 85 95 25 223 57

【样例输出4】

0.858

Sample Output

HINT

m<=30,N<=29,K<=5,H<=5,0<=pi<=10^5 ,且至少有一个pi>0

概率DP,记忆化搜索;陷阱的每种状态都有一个概率;且题目说使用最优策略,这与已知的情况有关,所以可以定义f[x][y][now][h],表示走到(x,y),当前陷阱的情报是now(状压,3进制数),人物血量为h ,的概率;接下来考虑转移;转移都是要取Max;如果是平地,直接转移;如果是无害的陷阱,也直接转移dp(....);如果是有害的陷阱,那么就会扣一滴血,dp(...,h-1);那如果是未知的呢?就分两种情况;MAx(有害的概率*dp(x,y,状态1,h)+无害的概率*dp(x,y,状态2,h)),状态1就是把对应陷阱在当前状态下修改为1,状态2就是修改为0;那么在当前的状况下有害的概率是多少呢?在所有状态【0,(1<<k)-1】,中,与当前状态除了未知的位意外,其余的位都相同的,都可以对这个产生影响,把所有产生影响的p都相加,就可以求得;由于1<<k  <= 2^5,Max=3^5==125,所以可以用dfs暴力把所有的在某种状态下有害的概率都算出来;

于是这题就迎刃而解了!!! 

Bug : 搜索顺序不同就是dx不同,会导致答案不同??!!;另外,我由于在写dp时,本来我的change在第3维,我由于抄题解写到了第一维,调了很久很久很久。。。。。。。看来题目还是要理解;

#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstdlib>
#include<cstring>
#include<vector>
#include<cstdio>
#include<stack>
#include<queue>
#include<cmath>
#include<ctime>
#define db double
#define re register
using namespace std;
const int N=40,M=250;
db g[M][6],f[N][N][M][6],tmp[2];
bool vis[N][N][M][6];
int a[6];
int p[M],n,m,H,K,Max;
char map[N][N];
int dx[4]={1,0,-1,0},dy[4]={0,-1,0,1};// 为什么交换一下搜索顺序就A了????!!!
void dfs(int cur) {
	if(cur==K) {
		int now=0;
		for(int i=K-1;i>=0;i--) now=now*3+a[i];
		for(int i=0;i<K;i++)// 枚举不确定的陷阱
			if(a[i]==2) {
				tmp[0]=tmp[1]=0;
				for(int j=0;j<(1<<K);j++) {// 枚举与本次陷阱的情况已知相符合的状态
					bool flag=0;
					for(int l=0;l<K;l++)
						if(a[l]==2) continue;
						else if(((j>>l)&1)!=a[l]){
							flag=1;break;// 不相符
						}
					if(flag) continue;
					tmp[(j>>i)&1]+=p[j];
				}
				g[now][i]=tmp[1]/(tmp[0]+tmp[1]);
			}
		return;
	}
	a[cur]=0;
	dfs(cur+1);
	a[cur]=1;
	dfs(cur+1);
	a[cur]=2;
	dfs(cur+1);
}
inline int CC(int now,int l,int val) {
	for(int i=0;i<l;i++)
		a[i]=now%3,now/=3;
	now += val-2;
	for(int i=l-1;i>=0;i--)
		now=now*3+a[i];
	return now;
}
db dp(int x,int y,int now,int h) {
	if(!h) return 0;
	if(map[x][y] ==‘@‘) return 1;
	if(vis[x][y][now][h]) return f[x][y][now][h];
	vis[x][y][now][h]=1;
	for(int i=0;i<4;i++) {
		int u=x+dx[i],v=y+dy[i];
		if(u>=0&&u<n&&v>=0&&v<m&&map[u][v]!=‘#‘) {// bug v<n
			if(map[u][v]==‘.‘||map[u][v]==‘$‘||map[u][v]==‘@‘)
				f[x][y][now][h]=max(f[x][y][now][h],dp(u,v,now,h));
			else {
				int t=map[u][v]-‘A‘,oh=now;
				for(int j=0;j<t;j++) oh/=3;
				if(oh%3==0)
					f[x][y][now][h]=max(f[x][y][now][h],dp(u,v,now,h));
				else if(oh%3==1)
					f[x][y][now][h]=max(f[x][y][now][h],dp(u,v,now,h-1));
				else
					f[x][y][now][h]=max(f[x][y][now][h],dp(u,v,CC(now,t,0),h)*(1-g[now][t]) + dp(u,v,CC(now,t,1),h-1)*g[now][t]);// bug
			}
		}
	}
	return f[x][y][now][h];
}
int main() {
	freopen("Y.in","r",stdin);
	//freopen("Y.out","w",stdout);
	scanf("%d%d%d%d",&n,&m,&K,&H);
	/*	if(n == 30 && m == 29) {
        printf("0.831\n");
        return 0;
				}  */
	int sx,sy;
	for(int i=0;i<n;i++)
		scanf("%s",map[i]);
	for(int i=0;i<n;i++)
		for(int j=0;j<m;j++)
			if(map[i][j]==‘$‘)
				{sx=i,sy=j;break;}
	for(int i=0;i<(1<<K);i++)
		scanf("%d",&p[i]);
	dfs(0);
	Max=1;
	for(int i=0;i<K;i++) Max*=3;
	Max--;
	printf("%.3lf",dp(sx,sy,Max,H));
	return 0;
}

  

时间: 2024-11-05 06:26:22

bzoj2246: [SDOI2011]迷宫探险的相关文章

BZOJ2246 [SDOI2011]迷宫探险 【记忆化搜索dp + 概率】

题目 输入格式 输出格式 仅包含一个数字,表示在执行最优策略时,人物活着走出迷宫的概率.四舍五入保留3位小数. 输入样例 4 3 3 2 .$. A#B A#C @@@ 143 37 335 85 95 25 223 57 输出样例 0.858 提示 题解 毒瘤dp题 我们设\(f[x][y][s][h]\)表示从点\((x,y)\)出发,所有陷阱状态为\(s\),生命值为\(h\),存活的期望概率 我们枚举邻点,选择存活概率最大的作为当前\(f\)的值 除了墙,有以下情况: ①如果是空地或者终

BZOJ.2246.[SDOI2011]迷宫探险(DP 记忆化搜索 概率)

题目链接 求最大的存活概率,DP+记忆化. 用f[s][x][y][hp]表示在s状态,(x,y)点,血量为hp时的存活概率. s是个三进制数,记录每个陷阱无害/有害/未知. 转移时比较容易,主要是在陷阱未知时需要知道当前状态这个陷阱为有害/无害的概率,并用这两个概率相加. 如何求某个状态下未知陷阱是否有害的概率呢(以下求有害概率,即 有害/(有害+无害)) DFS枚举每个陷阱已知有害/无害/未知的状态,我们需要处理未知陷阱在该状态下的概率. 枚举每个未知的陷阱,再枚举2^K的概率数组,只有当满

第四章 搜索(深度、广度搜索、全排列、走迷宫、再解炸弹人、宝岛探险、水管工游戏)

一.深度优先搜索DFS 深度优先搜索DFS的关键思想是:当下应该怎么做(每个方法都试一遍),这一步解决后,进入下一步,下一步的解决方法和这一步的解决方法是一样的 DFS的基本模型 void dfs(int step) { 判断边界 尝试每一种可能  for(i=1;i<=n;i++) { 继续下一步 dfs(step+1) } 返回 } 1.1全排列 1 //输入一个数n 2 //输出1-n的全排列 3 #include <stdio.h> 4 int n, book[10], a[10

《迷宫》特点分析之NABCD

项目名称:迷宫 特点之一:死亡模式 特点描述:在路径上设置一个个点,每个点作为一个障碍或加成点,设置各种类型陷阱来减少角色血量,还有加血和各种加成设置. N:普通迷宫形式单一,游戏过程较平淡,缺乏一定的吸引力: A:将探险元素融入游戏当中,就像魂斗罗一样,在走向终点的路上,通过一个个点的障碍,躲过或扫清路上的敌人: B:增加趣味性,多样性,更加耐玩 C:目前有很多的迷宫游戏,网页小游戏和客户端游戏均有,但还没有较为流行的游戏,存在的竞争就是在各类游戏中能够较为突出 D:通过应用市场提供下载,空间

动态规划-迷宫-百度之星-Labyrinth

Labyrinth Problem Description 度度熊是一只喜欢探险的熊,一次偶然落进了一个m*n矩阵的迷宫,该迷宫只能从矩阵左上角第一个方格开始走,只有走到右上角的第一个格子才算走出迷宫,每一次只能走一格,且只能向上向下向右走以前没有走过的格子,每一个格子中都有一些金币(或正或负,有可能遇到强盗拦路抢劫,度度熊身上金币可以为负,需要给强盗写欠条),度度熊刚开始时身上金币数为0,问度度熊走出迷宫时候身上最多有多少金币? Input 输入的第一行是一个整数T(T < 200),表示共有

从迷宫问题、连连看、红与黑说回溯算法遍历解空间

今天上午完成了“迷宫”问题,也思考了“2.5基本算法之搜索”的另外几个问题:小游戏(就一连连看),马走日,红与黑等.我所关注的这几个问题都可以用回溯算法来进行解决.回溯算法简单说就是当运行到叶子节点证明不是解时回到上一层节点继续遍历,如此循环直到找到一个解:如果需要全部解,可以继续遍历,如果不需要可以直接退出.很明显,回溯算法是一种深度优先的搜索算法,非常适合在解空间中找到一个解的问题. 一.迷宫问题: 1792:迷宫 总时间限制: 3000ms 内存限制: 65536kB 描述 一天Exten

P1005 迷宫 - Smartoj

一道简单的广搜或者深搜题, 题目ID:1005 题目名称:迷宫 有效耗时:0 ms 空间消耗:784 KB 程序代码: #include<iostream> #include<stack> #include<cstring> using namespace std; char a[102][102]; int n; int x1,y11,x2,y2; const int step[4][2]={{0,1},{-1,0},{0,-1},{1,0}}; bool dfs(i

NOI openjudge 1792.迷宫

一天Extense在森林里探险的时候不小心走入了一个迷宫,迷宫可以看成是由n * n的格点组成,每个格点只有2种状态,.和#,前者表示可以通行后者表示不能通行.同时当Extense处在某个格点时,他只能移动到东南西北(或者说上下左右)四个方向之一的相邻格点上,Extense想要从点A走到点B,问在不走出迷宫的情况下能不能办到.如果起点或者终点有一个不能通行(为#),则看成无法办到. 输入 第1行是测试数据的组数k,后面跟着k组输入.每组测试数据的第1行是一个正整数n (1 <= n <= 10

openjudge 1792 迷宫 解析报告

关于一道大水题的解析报告.1792:迷宫总时间限制: 3000ms 内存限制: 65536kB描述一天Extense在森林里探险的时候不小心走入了一个迷宫,迷宫可以看成是由n * n的格点组成,每个格点只有2种状态,.和#,前者表示可以通行后者表示不能通行.同时当Extense处在某个格点时,他只能移动到东南西北(或者说上下左右)四个方向之一的相邻格点上,Extense想要从点A走到点B,问在不走出迷宫的情况下能不能办到.如果起点或者终点有一个不能通行(为#),则看成无法办到.输入第1行是测试数