地下迷宫
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

思路:
首先对地图节点重新标号。假设E[i]表示DK从i点开始走出迷宫的期望值。
那么E[i]=(E[a1]+E[a2]+E[a3]+...+E[an])/n+1,其中a1...an是i的相邻节点。
那么对于每一个DK可达的节点来说,都可以为它建立这样的一个方程。现
在假设DK可达的点有N个,那么我们最终将会得到N元一次方程组。方程成
环所以利用高斯消元解出E[No[S]]。其中S是DK的起点,No[S]是重标号后的
起点这里要重点注意的是,我们联立方程的时候,一定要注意DK可达这个条
件,不然就会导致无解的情况。貌似zjutoj崩了。不能交题了。代码仅供参考。

详细见代码:
#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&高斯消元)的更多相关文章

  1. BZOJ_3143_[Hnoi2013]游走_期望DP+高斯消元

    BZOJ_3143_[Hnoi2013]游走_期望DP+高斯消元 题意: 一个无向连通图,顶点从1编号到N,边从1编号到M. 小Z在该图上进行随机游走,初始时小Z在1号顶点,每一步小Z以相等的概率随机 ...

  2. BZOJ.2707.[SDOI2012]走迷宫(期望 Tarjan 高斯消元)

    题目链接 一个点到达终点的期望步数 \(E_i=\sum_{(i,j)\in G}\frac{E_j+1}{out[i]}\),\(out[i]\)为点\(i\)的出度. 那么对于一个DAG可以直接在 ...

  3. 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 ...

  4. hdu4418 Time travel 【期望dp + 高斯消元】

    题目链接 BZOJ4418 题解 题意:从一个序列上某一点开始沿一个方向走,走到头返回,每次走的步长各有概率,问走到一点的期望步数,或者无解 我们先将序列倍长形成循环序列,\(n = (N - 1) ...

  5. 【noi2019集训题1】 脑部进食 期望dp+高斯消元

    题目大意:有n个点,m条有向边,每条边上有一个小写字母. 有一个人从1号点开始在这个图上随机游走,游走过程中他会按顺序记录下走过的边上的字符. 如果在某个时刻,他记录下的字符串中,存在一个子序列和S2 ...

  6. LightOJ 1151 Snakes and Ladders 期望dp+高斯消元

    题目传送门 题目大意:10*10的地图,不过可以直接看成1*100的,从1出发,要到达100,每次走的步数用一个大小为6的骰子决定.地图上有很多个通道 A可以直接到B,不过A和B大小不确定   而且 ...

  7. P4457-[BJOI2018]治疗之雨【期望dp,高斯消元】

    正题 题目链接:https://www.luogu.com.cn/problem/P4457 题目大意 开始一个人最大生命值为\(n\),剩余\(hp\)点生命,然后每个时刻如果生命值没有满那么有\( ...

  8. Codeforces.24D.Broken robot(期望DP 高斯消元)

    题目链接 可能这儿的会更易懂一些(表示不想再多写了). 令\(f[i][j]\)表示从\((i,j)\)到达最后一行的期望步数.那么有\(f[n][j]=0\). 若\(m=1\),答案是\(2(n- ...

  9. BZOJ2707: [SDOI2012]走迷宫(期望 tarjan 高斯消元)

    题意 题目链接 Sol 设\(f[i]\)表示从\(i\)走到\(T\)的期望步数 显然有\(f[x] = \sum_{y} \frac{f[y]}{deg[x]} + 1\) 证明可以用全期望公式. ...

随机推荐

  1. cygwin org/apache/zookeeper/KeeperException

    以前用cdh3-0.20的hbase,在windows下面直接启动就行了,但是最近安装0.94以上的,就不行了. 报标题的错误,搜遍网络,几乎都是要加HBASE_CLASSPATH的,后来看老外的文章 ...

  2. 跨平台生成GUID/UUID

    #ifndef XGUID_H#define XGUID_H #include <string>#include <stdio.h>#ifdef WIN32#include & ...

  3. Eddy's爱好(dfs+容斥)

    Eddy's爱好 Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total S ...

  4. easui tree载入时自动展开无子节点的节点

    利用loadFilter对后台返回的原始数据进行过滤处理,将数据中的state字段修改相应的值,若无子节点,则改成open,若有子节点,则改成closed. 由于一个node属性只有下面几项内容,因此 ...

  5. 2013 ACM/ICPC 成都网络赛解题报告

    第三题:HDU 4730 We Love MOE Girls 传送门:http://acm.hdu.edu.cn/showproblem.php?pid=4730 水题~~~ #include < ...

  6. LoaderManager使用具体解释(三)---实现Loaders

    这篇文字将介绍Loader<D>类,而且介绍自己定义Loader的实现.这是本系列的第三篇文章. 一:Loaders之前世界 二:了解LoaderManager 三:实现Loaders 四 ...

  7. UIBezierPath详解

    使用UIBezierPath类可以创建基于矢量的路径,这个类在UIKit中.此类是Core Graphics框架关于path的一个封装.使用此类可以定义简单的形状,如椭圆或者矩形,或者有多个直线和曲线 ...

  8. BZOJ 1925: [Sdoi2010]地精部落( dp )

    dp(i,j)表示1~i的排列中, 以1~j为开头且开头是下降的合法方案数 这种数列具有对称性, 即对于一个满足题意且开头是上升的n的排列{an}, 令bn = n-an+1, 那么{bn}就是一个满 ...

  9. Linux怎么设置PostgreSQL远程访问

    原文链接: Linux怎么设置PostgreSQL远程访问 安装PostgreSQL数据库之后,默认是只接受本地访问连接.如果想在其他主机上访问PostgreSQL数据库服务器,就需要进行相应的配置. ...

  10. Python 做过哪些有趣的项目

          1 icedx   241 天前 via Android   ♥ 1 考虑到Windows 下的类Alfred 软件都太傻逼 自己用PyQT 写了一个       2 crazyxin19 ...