ZJUT 1423 地下迷宫(期望DP&高斯消元)
地下迷宫
Time Limit:1000MS Memory Limit:32768K
Description:
由于山体滑坡,DK被困在了地下蜘蛛王国迷宫。为了抢在DH之前来到TFT,DK必须尽快走出此迷宫。此迷宫仅有一个出口,而由于大BOSS的力量减弱影响到了DK,使DK的记忆力严重下降,他甚至无法记得他上一步做了什么。所以他只能每次等概率随机的选取一个方向走。当然他不会选取周围有障碍的地方走。如DK周围只有两处空地,则每个都有1/2的概率。现在要求他平均要走多少步可以走出此迷宫。
Input:
先是一行两个整数N, M(1<=N, M<=10)表示迷宫为N*M大小,然后是N行,每行M个字符,'.'表示是空地,'E’表示出口,'D’表示DK,'X’表示障碍。
Output:
如果DK无法走出或要超过1000000步才能走出,输出tragedy!,否则输出一个实数表示平均情况下DK要走几步可以走出迷宫,四舍五入到小数点后两位。
Sample Input:
1 2
ED
3 3
D.X
.X.
X.E
Sample Output:
1.00
tragedy!
Source:
DK
那么对于每一个DK可达的节点来说,都可以为它建立这样的一个方程。现
在假设DK可达的点有N个,那么我们最终将会得到N元一次方程组。方程成
#include <iostream>
#include<stdio.h>
#include<math.h>
#include<string.h>
#include<queue>
using namespace std;
const int maxn=15;
const double eps=1e-9;
char maze[maxn][maxn];//记录地图
int pp[maxn][maxn];//重编号
int dx[4]={0,0,-1,1};
int dy[4]={1,-1,0,0};
double mat[maxn][maxn];//记录矩阵
int n,m,cnt,ptr;
struct node
{
int x,y;
node(int xx,int yy)
{
x=xx;
y=yy;
}
node(){}
} st,ed,t;
queue<node> q; bool isok(int x,int y)//判断是否越界
{
return x>=0&&x<n&&y>=0&&y>=0&&y<m&&maze[x][y]!='X';
}
void bfs()//宽搜。记录可到达点
{
int nx,ny,i;
while(!q.empty())
q.pop();
cnt=0;
nx=st.x;
ny=st.y;
pp[nx][ny]=cnt++;
q.push(st);
while(!q.empty())
{
t=q.front();
q.pop();
for(i=0;i<4;i++)
{
nx=t.x+dx[i];
ny=t.y+dy[i];
if(isok(nx,ny)&&pp[nx][ny]==-1)
{
q.push(node(nx,ny));
pp[nx][ny]=cnt++;//对可到达点编号
}
}
}
}
bool guass()//高斯消元
{
int row,i,j,id;
double maxx,var;
for(row=0;row<cnt;row++)//遍历行。重点在mat[row][row]先找此处最大系数。然后把以下方程的对应未知数消去
{
maxx=fabs(mat[row][row]);
id=row;//id记录位置
for(i=row+1;i<cnt;i++)
{
if(fabs(mat[i][row])>maxx)
{
maxx=fabs(mat[i][row]);//注意是绝对值大
id=i;
}
}
if(maxx<eps)
return false;
if(id!=row)//如果就是当前处理行就不用交换
{
for(i=row;i<=cnt;i++)//交换最大行和当前行
swap(mat[row][i],mat[id][i]);
}
for(i=row+1;i<cnt;i++)//遍历行。所以<cnt.把当前处理行以下的mat[row][row]变量消去。
{
if(fabs(mat[i][row])<eps)//本来就为0就不用处理了
continue;
var=mat[i][row]/mat[row][row];
for(j=row;j<=cnt;j++)//包括扩展矩阵所以c<=cnt。
mat[i][j]-=mat[row][j]*var;
}
}
for(i=cnt-1;i>=0;i--)//从最后一个系数开始
{
for(j=i+1;j<cnt;j++)
mat[i][cnt]-=mat[i][j]*mat[j][j];
mat[i][i]=mat[i][cnt]/mat[i][i];//现在系数矩阵的对角线用于记录答案。
}
return true;
}
int main()
{
int i,j,k,nx,ny,p; while(~scanf("%d%d",&n,&m))
{
for(i=0;i<n;i++)
{
scanf("%s",maze[i]);
for(j=0;j<m;j++)
{
if(maze[i][j]=='D')
st.x=i,st.y=j;
else if(maze[i][j]=='E')
ed.x=i,ed.y=j;
}
}
memset(pp,-1,sizeof pp);
bfs();
if(pp[ed.x][ed.y]==-1)
{
printf("tragedy!\n");
continue;
}
memset(mat,0,sizeof mat);
for(i=0;i<n;i++)
{
for(j=0;j<m;j++)
{
if(pp[i][j]!=-1)//以每个可到达点建立方程组
{
ptr=0;
p=pp[i][j];
for(k=0;k<4;k++)
{
nx=i+dx[k];
ny=j+dy[k];
if(isok(nx,ny))
{
mat[p][pp[nx][ny]]=-1;
ptr++;
}
}
mat[p][p]=ptr;
mat[p][cnt]=ptr;
}
}
}
p=pp[ed.x][ed.y];
memset(mat[p],0,sizeof mat[p]);
mat[p][p]=1;//在终点步数的期望为0.
if(guass())
{
p=pp[st.x][st.y];
if(mat[p][p]<=1000000)
printf("%.2lf\n",mat[p][p]);
else
printf("tragedy!\n");
}
else
printf("tragedy!\n");
}
return 0;
}
ZJUT 1423 地下迷宫(期望DP&高斯消元)的更多相关文章
- BZOJ_3143_[Hnoi2013]游走_期望DP+高斯消元
BZOJ_3143_[Hnoi2013]游走_期望DP+高斯消元 题意: 一个无向连通图,顶点从1编号到N,边从1编号到M. 小Z在该图上进行随机游走,初始时小Z在1号顶点,每一步小Z以相等的概率随机 ...
- BZOJ.2707.[SDOI2012]走迷宫(期望 Tarjan 高斯消元)
题目链接 一个点到达终点的期望步数 \(E_i=\sum_{(i,j)\in G}\frac{E_j+1}{out[i]}\),\(out[i]\)为点\(i\)的出度. 那么对于一个DAG可以直接在 ...
- HDU 2262 Where is the canteen 期望dp+高斯消元
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=2262 Where is the canteen Time Limit: 10000/5000 MS ...
- hdu4418 Time travel 【期望dp + 高斯消元】
题目链接 BZOJ4418 题解 题意:从一个序列上某一点开始沿一个方向走,走到头返回,每次走的步长各有概率,问走到一点的期望步数,或者无解 我们先将序列倍长形成循环序列,\(n = (N - 1) ...
- 【noi2019集训题1】 脑部进食 期望dp+高斯消元
题目大意:有n个点,m条有向边,每条边上有一个小写字母. 有一个人从1号点开始在这个图上随机游走,游走过程中他会按顺序记录下走过的边上的字符. 如果在某个时刻,他记录下的字符串中,存在一个子序列和S2 ...
- LightOJ 1151 Snakes and Ladders 期望dp+高斯消元
题目传送门 题目大意:10*10的地图,不过可以直接看成1*100的,从1出发,要到达100,每次走的步数用一个大小为6的骰子决定.地图上有很多个通道 A可以直接到B,不过A和B大小不确定 而且 ...
- P4457-[BJOI2018]治疗之雨【期望dp,高斯消元】
正题 题目链接:https://www.luogu.com.cn/problem/P4457 题目大意 开始一个人最大生命值为\(n\),剩余\(hp\)点生命,然后每个时刻如果生命值没有满那么有\( ...
- Codeforces.24D.Broken robot(期望DP 高斯消元)
题目链接 可能这儿的会更易懂一些(表示不想再多写了). 令\(f[i][j]\)表示从\((i,j)\)到达最后一行的期望步数.那么有\(f[n][j]=0\). 若\(m=1\),答案是\(2(n- ...
- BZOJ2707: [SDOI2012]走迷宫(期望 tarjan 高斯消元)
题意 题目链接 Sol 设\(f[i]\)表示从\(i\)走到\(T\)的期望步数 显然有\(f[x] = \sum_{y} \frac{f[y]}{deg[x]} + 1\) 证明可以用全期望公式. ...
随机推荐
- Windows Phone 8初学者开发—第6部分:设置应用程序的样式
原文 Windows Phone 8初学者开发—第6部分:设置应用程序的样式 Source Code: http://aka.ms/absbeginnerdevwp8 PDF Version: ht ...
- java学习之字符流与字节流的转换
package com.io; import java.io.File; import java.io.FileInputStream; import java.io.FileNotFoundExce ...
- ARPU_百度百科
ARPU_百度百科 ARPU
- 64位CentOS上编译 Hadoop 2.2.0
下载了Hadoop预编译好的二进制包,hadoop-2.2.0.tar.gz,启动起来后.总是出现这样的警告: WARN util.NativeCodeLoader: Unable to load n ...
- 你会用swift创建复杂的加载动画吗(1)
时至今日,iOS 应用商店已经拥有超过了140万 应用,让你自己的应用脱颖而出确实是个不小的挑战.不过,在你的应用掉入默默无闻的大黑洞之前,你拥有一个小小的机遇窗,它能帮你吸引用户的注意. AD: 时 ...
- 友盟iOS微信登陆没有回调的原因
1.在友盟文档中这样说: 链接 7.4 微信登录 添加配置文件参考文档:添加微信及朋友圈,添加相关库文件,配置URL schemes及添加系统回调 注意微信登录必须先在微信开放平台申请微信登录权限 在 ...
- Ural 1450 求最长路 SPFA
题意就是求S点到T点的有向无环图中的最长路. 用SPFA可以解决. 一开始一直RE的原因 QAQ 竟然是在开Edge 邻接表的时候开小了 改了一下4Y #include <stdio.h> ...
- Ural 1068 - Sum
Your task is to find the sum of all integer numbers lying between 1 and N inclusive. Input The input ...
- 设计模式 ( 十五 ) 中介者模式Mediator(对象行为型)
设计模式 ( 十五 ) 中介者模式Mediator(对象行为型) 1.概述 在面向对象的软件设计与开发过程中,根据“单一职责原则”,我们应该尽量将对象细化,使其只负责或呈现单一的职责,即将行为分布到各 ...
- oschina程序开发
程序开发 102脚本编程语言 36地图相关 7Epub电子图书工具 109UI组件库 16代码生成工具 25SVG开发包 17推荐引擎 3指纹识别 23拼音转换工具包 24蓝牙开发 295GUI开发框 ...