码迷,mamicode.com
首页 > 其他好文 > 详细

LOOPS

时间:2016-08-23 18:55:20      阅读:128      评论:0      收藏:0      [点我收藏+]

标签:

LOOPS

题目链接:http://acm.split.hdu.edu.cn/showproblem.php?pid=3853

递推

设dp[i][j]为(i,j)到终点期望的使用魔力值,mp[i][j][k]为(i,j)到三个方向的概率。

那么,dp[i][j]=(2+dp[i][j+1])*mp[i][j][1]+(2+dp[i+1][j])*mp[i][j][2]+(2+dp[i][j])*mp[i][j][0],

移向得,dp[i][j]*(1-mp[i][j][0])=(2+dp[i][j+1])*mp[i][j][1]+(2+dp[i+1][j])*mp[i][j][2]+2*mp[i][j][0],

化简得,dp[i][j]=(2+dp[i][j+1]*mp[i][j][1]+dp[i+1][j]*mp[i][j][2])/(1-mp[i][j][0]).

注意当mp[i][j][0]==1时,此时将陷入LOOPS.

代码如下:

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #define EPS 1e-8
 5 #define N 1005
 6 #define INF 1000005
 7 using namespace std;
 8 int n,m;
 9 double mp[N][N][3];
10 double dp[N][N];
11 int main(void){
12     while(~scanf("%d%d",&n,&m)){
13         memset(dp,0,sizeof(dp));
14         for(int i=0;i<n;++i)
15         for(int j=0;j<m;++j)
16         for(int k=0;k<3;++k)
17             scanf("%lf",&mp[i][j][k]);
18         for(int i=n-1;i>=0;--i)
19         for(int j=m-1;j>=0;--j){
20             if(i==n-1&&j==m-1)continue;
21             if(fabs(1-mp[i][j][0])<EPS)dp[i][j]=INF;
22             else dp[i][j]=(2+dp[i][j+1]*mp[i][j][1]+dp[i+1][j]*mp[i][j][2])/(1-mp[i][j][0]);
23         }
24         printf("%.3lf\n",dp[0][0]);
25     }
26 }

 

LOOPS

标签:

原文地址:http://www.cnblogs.com/barrier/p/5800401.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!