概率DP/数学期望

　　kuangbin总结中的第7题

　　其实跟UVA 11762 Race To 1 那道题差不多……直接推下公式，然后倒推即可

 

Trick：有的点可能是p1[i][j]==1……这样的点是永远不会走出去的……所以也不能走到……遇到这样的点直接跳过就好了TAT

　　　　但是！！浮点数不能直接判定相等……应该写成 fabs(1-p[i][j][0])<eps 来判定！！！sad……

【HDOJ】【3853】LOOPS_概率DP【HDOJ】【3853】LOOPS_概率_02
 1 //HDOJ 3853
 2 #include<cmath>
 3 #include<cstdio>
 4 #define F(i,j,n) for(int i=j;i<=n;++i)
 5 #define D(i,j,n) for(int i=j;i>=n;--i)
 6 
 7 const int N=1100,INF=~0u>>2;
 8 const double eps=1e-5;
 9 double f[N][N],p[N][N][3];
10 
11 int main(){
12     int n,m;
13     while(scanf("%d%d",&n,&m)!=EOF){
14         F(i,1,n) F(j,1,m) F(k,0,2)
15             scanf("%lf",&p[i][j][k]);
16         f[n][m]=0;
17         D(i,n,1)
18             D(j,m,1){
19                 if (i==n && j==m) continue;
20                 if (fabs(1-p[i][j][0])<eps) continue;
21                 double *P=p[i][j];
22                 f[i][j]=2/(1-P[0]) + f[i][j+1]*P[1]/(1-P[0]) + f[i+1][j]*P[2]/(1-P[0]);
23             }
24         printf("%.3lf\n",f[1][1]);
25     }
26     return 0;
27 }
View Code