LOOPS HDU - 3853(概率dp,期望)

题意:

有一个R*C的方格。一个人想从(1,1)走到(r,c)。在每个格子都有三种选择,向下,向右,或者原地不动。每个格子里的每个选择都有一定的概率。而每次移动都需要消耗2点的能量,问期望消耗的能量是多少。

题目:

Akemi Homura is a Mahou Shoujo (Puella Magi/Magical Girl).

Homura wants to help her friend Madoka save the world. But because of the plot of the Boss Incubator, she is trapped in a labyrinth called LOOPS.
在这里插入图片描述

The planform of the LOOPS is a rectangle of R*C grids. There is a portal in each grid except the exit grid. It costs Homura 2 magic power to use a portal once. The portal in a grid G(r, c) will send Homura to the grid below G (grid(r+1, c)), the grid on the right of G (grid(r, c+1)), or even G itself at respective probability (How evil the Boss Incubator is)!
At the beginning Homura is in the top left corner of the LOOPS ((1, 1)), and the exit of the labyrinth is in the bottom right corner ((R, C)). Given the probability of transmissions of each portal, your task is help poor Homura calculate the EXPECT magic power she need to escape from the LOOPS.

Input

The first line contains two integers R and C (2 <= R, C <= 1000).

The following R lines, each contains C*3 real numbers, at 2 decimal places. Every three numbers make a group. The first, second and third number of the cth group of line r represent the probability of transportation to grid (r, c), grid (r, c+1), grid (r+1, c) of the portal in grid (r, c) respectively. Two groups of numbers are separated by 4 spaces.

It is ensured that the sum of three numbers in each group is 1, and the second numbers of the rightmost groups are 0 (as there are no grids on the right of them) while the third numbers of the downmost groups are 0 (as there are no grids below them).

You may ignore the last three numbers of the input data. They are printed just for looking neat.

The answer is ensured no greater than 1000000.

Terminal at EOF

Output

A real number at 3 decimal places (round to), representing the expect magic power Homura need to escape from the LOOPS.

Sample Input

2 2
0.00 0.50 0.50 0.50 0.00 0.50
0.50 0.50 0.00 1.00 0.00 0.00

Sample Output

6.000

分析:

f[i][j]为从(i,j)到(r,c)的期望消耗。G[i][j][k]记录每次走动三种选择的概率从(i,j)有三种转移方法向下。
在原地不动的概率是G[i][j][1];
向右的概率为G[i][j][2];
向下的概率为G[i][j][3]。
但是在原地不动是不消耗能量的.所以转移我们只考虑转移到右边和下边的情况。
f[i][j]=f[i][j]*G[i][j][1]+f[i+1][j]*G[i][j][2]+f[i][j+1]*G[i][j][3]+2
将f[i][j]*G[i][j][1]移到左边,得到状态转移方程
f[i][j]=(f[i+1][j]*G[i][j][2]+f[i][j+1]*G[i][j][3]+2)/(1-G[i][j][1]).

AC代码

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;
const int M=1e3+10;
const int eps=1e-5;
double G[M][M][5];
double f[M][M];
int r,c;
int main()
{while(~scanf("%d%d",&r,&c)&&r&&c){for(int i=1; i<=r; i++)for(int j=1; j<=c; j++)for(int l=1; l<=3; l++)scanf("%lf",&G[i][j][l]);memset(f,0,sizeof(f));for(int i=r; i>=1; i--)for(int j=c; j>=1; j--){if(i==r&&j==c)continue;if(fabs(1-G[i][j][1])<=eps)continue;f[i][j]=(f[i][j+1]*G[i][j][2]+f[i+1][j]*G[i][j][3]+2)/(1-G[i][j][1]);}printf("%.3f\n",f[1][1]);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/309806.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Sql Server之旅——第四站 你必须知道的非聚集索引扫描

非聚集索引&#xff0c;这个是大家都非常熟悉的一个东西&#xff0c;有时候我们由于业务原因&#xff0c;sql写的非常复杂&#xff0c;需要join很多张表&#xff0c;然后就泪流满面了。。。这时候就有DBA或者资深的开发给你看这个猥琐的sql&#xff0c;通过执行计划一分析,或许…

[C++11]可调用对象绑定器

std::bind用来将可调用对象与其参数一起进行绑定。绑定后的结果可以使用std::function进行保存&#xff0c;并延迟调用到任何我们需要的时候。通俗来说&#xff0c;它主要有两个作用: 1.将可调用对象与其参数一起绑定成一个仿函数。 2.将多元(参数个数为n&#xff0c;n > …

电视android已停止运行是什么意思,智能电视提示应用停止运行怎么办?当贝市场三招解决...

智能电视提示应用停止运行怎么办&#xff1f;当贝市场三招解决2019年11月28日 17:53作者&#xff1a;网络编辑&#xff1a;王动分享智能电视使用久了之后,电视页面会提示我们应用停止运行,这是怎么回事?当贝小编针对这个问题,整理了一份解决教程,大家可以看看有没有什么帮助。…

重磅!2020年微软开发者大会落幕,.NET迎来新机遇!

两天前微软举行了首个线上Build大会&#xff0c;而开发者成为大会里唯一的主角。和所有技术公司一样&#xff0c;开发者对于微软来说&#xff0c;同样也越来越重要了。如同血肉相依的关系&#xff0c;谁也离不开谁。在这次大会上&#xff0c;开发者是最大的宠儿&#xff0c;成了…

番茄时间有感之关于在疫情期间我与ACM不得不说的故事

哼哼~首先声明&#xff0c;我不是来讲故事的&#xff0c;我来总结一下在疫情这段时间&#xff0c;在ACM训练过程中的自我情况的总结和反思&#xff0c;嘻嘻&#xff0c;我是一个标题党&#xff0c;如果是被标题骗进来哒&#xff0c;抱歉啦&#xff0c;有句话说的好&#xff0c;…

[PAT乙级]1042 字符统计

请编写程序&#xff0c;找出一段给定文字中出现最频繁的那个英文字母。 输入格式&#xff1a; 输入在一行中给出一个长度不超过 1000 的字符串。字符串由 ASCII 码表中任意可见字符及空格组成&#xff0c;至少包含 1 个英文字母&#xff0c;以回车结束&#xff08;回车不算在内…

获取壁纸设置背景android,【Android学习】获取Bing 15天前到明天的壁纸,并设置为背景...

//设置壁纸public void setWallpaper() {new Thread() {public void run() {try {WebImage web new WebImage(imgRealPath);Bitmap bmp web.getBitmap(MainActivity.this);Message msg new Message();msg.obj bmp;msg.what 2;mHandler.sendMessage(msg);} catch (Exceptio…

基于 abp vNext 和 .NET Core 开发博客项目 - 自定义仓储之增删改查

上一篇文章我们用Code-First的方式创建了博客所需的实体类&#xff0c;生成了数据库表&#xff0c;完成了对EF Core的封装。本篇说一下自定义仓储的实现方式&#xff0c;其实在abp框架中已经默认给我们实现了默认的通用(泛型)仓储&#xff0c;IRepository<TEntity, TKey>…

计算机操作系统第四章作业

计算机操作系统第四章作业 1.何为静态链接&#xff1f;静态链接时需要解决两个什么问题? 答&#xff1a;静态链接是指在程序运行之前&#xff0c;先将各自目标模块及它们所需的库函数&#xff0c;链接成一个完整的装入模块&#xff0c;以后不再拆开的链接方式。   将几个目…

[PAT乙级]1046 划拳

划拳是古老中国酒文化的一个有趣的组成部分。酒桌上两人划拳的方法为&#xff1a;每人口中喊出一个数字&#xff0c;同时用手比划出一个数字。如果谁比划出的数字正好等于两人喊出的数字之和&#xff0c;谁就赢了&#xff0c;输家罚一杯酒。两人同赢或两人同输则继续下一轮&…

走进WebApiClientCore的设计

WebApiClientWebApiClient是NCC开源社区的一个项目&#xff0c;是目前微服务里http接口调用的一把锋利尖刀&#xff0c;项目早期设计与开发的时候&#xff0c;是基于.netframework的&#xff0c;然后慢慢加入netstandard和netcoreapp多个框架的支持&#xff0c;设计能力出众&am…

android 片段,android – 将片段添加到片段中(嵌套片段)

我想动态地将youtube片段添加到我已经存在的片段中.我使用的代码如下&#xff1a;// setting the Youtube Player Dynamicallyprivate int setYoutubePlayer(String desc,View view,int prevID,Bundle input) {if (desc.indexOf("") ! -1) {desc desc.substring(des…

HDU 3062 Party(2-sat题模板+tarjan )

题目: 有n对夫妻被邀请参加一个聚会&#xff0c;因为场地的问题&#xff0c;每对夫妻中只有1人可以列席。在2n 个人中&#xff0c;某些人之间有着很大的矛盾&#xff08;当然夫妻之间是没有矛盾的&#xff09;&#xff0c;有矛盾的2个人是不会同时出现在聚会上的。有没有可能会…

[PAT乙级]1043 输出PATest

给定一个长度不超过 10​4​​ 的、仅由英文字母构成的字符串。请将字符重新调整顺序&#xff0c;按 PATestPATest… 这样的顺序输出&#xff0c;并忽略其它字符。当然&#xff0c;六种字符的个数不一定是一样多的&#xff0c;若某种字符已经输出完&#xff0c;则余下的字符仍按…

Blazor WebAssembly 3.2 正式发布

5月 20日&#xff0c;微软 发布了 Blazor WebAssembly 3.2(https://devblogs.microsoft.com/aspnet/blazor-webassembly-3-2-0-now-available/) 。Blazor 是 ASP.NET Core 中的一个新框架&#xff0c;支持使用 C#和 HTML 创建交互式 Web 应用程序。Blazor WebAssembly 使用基于…

android 语音助手官网,breeno语音助手最新版

breeno语音助手最新版是一款手机中的导航软件&#xff0c;在这款软件中你能享受到非常方便的导航体验&#xff0c;这款软件中的指令不需要用户手动去输入&#xff0c;现在只需要你使用的语音就直接能对其进行操控了。感兴趣的用户就来去我下载网进行下载使用吧&#xff01;bree…

[PAT乙级]1047 编程团体赛

编程团体赛的规则为&#xff1a;每个参赛队由若干队员组成&#xff1b;所有队员独立比赛&#xff1b;参赛队的成绩为所有队员的成绩和&#xff1b;成绩最高的队获胜。 现给定所有队员的比赛成绩&#xff0c;请你编写程序找出冠军队。 输入格式&#xff1a; 输入第一行给出一个…

2-SAT适定性(Satisfiability)问题知识点详解

SAT是适定性(Satisfiability)问题的简称。一般形式为k-适定性问题&#xff0c;简称 k-SAT。而当k>2时该问题为NP完全的&#xff0c;所以我们只研究k2时情况。 2-SAT问题 现有一个由N个布尔值组成的序列A&#xff0c;给出一些限制关系&#xff0c;比如A[x] AND A[y]0、A[x]…

温故知新:Docker基础知识知多少?

【云原生】| 作者/Edison Zhou这是恰童鞋骚年的第233篇原创文章记得之前曾经粗略的写过一篇Docker的基础及ASP.NET Core部署Docker示例的入门文章&#xff0c;但那个时候刚刚学习对Docker的认知还比较浅&#xff0c;现在重新来温故知新一下。本文预计阅读时间为10min。1容器的用…

自动备份html文件,windows下定期自动备份本地文件(文件夹)

虽然网上有一些免费的文件自动备份软件&#xff0c;但是没有自己编写一段批处理来完成备份任务来的放心&#xff0c;而且不用占用系统资源。就给大家讲一下如何利用批处理完成本地文件或者文件夹的备份。1、批处理脚本该方法可把某文件夹下的文件同步到另外的文件夹&#xff0c…