hdu 3853 概率dp

/*
题目大意：一个n*c的网格，求从(1,1)到(n,c)花费魔法值的期望。
每次开启一个通道需要花费2魔法值，有三个方向可以走
在(x,y)位置时,分别转向
1.(x,y),概率是p(x,y*3-2)
2.(x,y+1),概率是p(x,y*3-1)
3.(x+1,y),概率是p(x,y*3)

dp[x][y]=p[x][y*3-2]*d[x][y]+p[x][y*3-1]*dp[x][y+1]+p[x][y*3]*dp[x+1][y]+2;
*/
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

const double eps=1e-8;
const int maxn=1005;
double p[maxn][3*maxn],dp[maxn][maxn];

int main()
{
    int n,m,i,j;
    while(~scanf("%d%d",&n,&m))
    {
        memset(p,0,sizeof(p));
        memset(dp,0,sizeof(dp));
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=3*m;j++)
                scanf("%lf",&p[i][j]);
        }
        for(i=n;i>0;i--)
        {
            for(j=m;j>0;j--)
            {
                if(i==n && j==m) continue;
                if(fabs(1-p[i][j*3-2])<eps) //呆在原地的概率为1
                    continue;
                dp[i][j]=(p[i][j*3-1]*dp[i][j+1]+p[i][j*3]*dp[i+1][j]+2)/(1-p[i][j*3-2]);
            }
        }
        printf("%.3lf\n",dp[1][1]);
    }
    return 0;
}