Broken robot CodeForces - 24D 期望dp

题意:n*m的棋盘,一个机器人在(i,j)处,每次等概率地停在原地,向左移动一格,向右移动一格,向下移动一格(不能移出棋盘).求走到最后一行所需期望步数.n<=1000,m<=1000

先讲坑点，1原地停留算一步，2要考虑机器人在同一层楼无限走路，3，m=1情况特殊考虑；

思路：因为是期望dp所以是逆推；从底层推到机器人所在层数；

先最外面套一层关于层数的循环；然后在套一层代表机器人可以在同一层游荡的循环，（循环次数随意，只要不超时，且使得结果在误差范围内都行，66好看）；

然后就是同一层的循环；

若机器人不在两端，就是有个地方可以呆

dp[i][j]=dp[i+1][j](1.0/4)+dp[i][j-1](1.0/4)+dp[i][j+1](1.0/4)+dp[i][j](1.0/4)+1;

如果在两端就是三个地方可以呆；

公式就是上面的稍改；

#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
double dp[1010][1010];
int main()
{
int n,m,a,b;
scanf("%d%d%d%d",&n,&m,&a,&b);
if(n==a)
printf("0.0000000000\n");
else
{
for(int i=n-1;i>=a;i--)
{
for(int k=1;k<=66;k++)
{
for(int j=1;j<=m;j++)
{
if(m==1)
dp[i][1]=(dp[1+i][1]*1.0/2+1)/(1.0/2);
else if(j==1)
dp[i][j]=(dp[i+1][j]*(1.0/3)+dp[i][j+1]*      (1.0/3)+1)/(2.0/3);
else if(j==m)
dp[i][j]=(dp[i+1][j]*(1.0/3)+dp[i][j-1]*(1.0/3)+1)/(2.0/3);
else
dp[i][j]=(dp[i+1][j]*(1.0/4)+dp[i][j-1]*(1.0/4)+dp[i][j+1]*(1.0/4)+1)/(3.0/4);
}
}
}
printf("%.10f\n",dp[a][b]);
}
}