HDU 3853 LOOPS (简单概率DP)

题目： LINK

题意：迷宫是一个R*C的布局，每个格子中给出停留在原地，往右走一个，往下走一格的概率，起点在(1,1)，终点在(R,C)，每走一格消耗两点能量，求出最后所需要的能量期望.
基础的概率期望DP吧，E[i][j] 表示从(i, j)这个位置到达(n, m) 这个位置花费的期望。
可以得到 E[i][j] = p1 * E[i][j] + p2 * E[i][j+1] + p3 * E[i+1][j] + 2 
注意一下p1 == 0 的情况，之后可以到这推，E[1][1]就是要求的结果了

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <string>
#include <vector>
#include <cmath>
#include <queue>
#include <map>
#include <set>
using namespace std; 
#define INF 1000000000
//typedef __int64 LL; 
#define N 1111
double dp[N][N], p1[N][N], p2[N][N], p3[N][N]; 
int n, m; 
int main() {
#ifndef ONLINE_JUDGE
	freopen("in.txt", "r", stdin); 
#endif // ONLINE_JUDGE
	while(scanf("%d%d", &n, &m) != EOF) {
		for(int i = 1; i <= n; i++) {
			for(int j = 1;j <= m; j++) {
				scanf("%lf%lf%lf",&p1[i][j], &p2[i][j], &p3[i][j]); 
			}
		}
		dp[n][m] = 0; 
		for(int i = n; i >= 1; i--) {
			for(int j = m; j >= 1; j --) {
				if(i == n && j == m) continue; 
				if(p1[i][j] >=1)  dp[i][j] = 0;  
				else dp[i][j] = (p3[i][j] * dp[i+1][j] + p2[i][j] * dp[i][j+1] + 2) / (1.0 - p1[i][j]); 
			}
		}
		printf("%.3lf\n", dp[1][1]);  
	}	
	return 0; 
}