ZOJ 3822 概率DP

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3822

记得比赛之前，看到WALK那道概率DP的题，然后就觉得应该可能出概率DP，十一还专门看了以前做的，

现场的时候，真的出了，然后好激动的样子，开始推得二维，然后感觉好难推，发现n只有50，所以就去推三维，然后发现k<max(i,j)的时候，有无用状态，无用状态初始化不会处理，然后觉得好像也用不到，但是感觉更稳的还是去推二维，然后就陷入二维--------悲剧了  最后打铁，真他妈不爽

也发现了自己的问题：1、平时学的东西掌握的不够扎实

2、其实如果队友没明白自己的意思，完全自己推吧，自信点，别想着必须队友确认才去做

想要保研牛校或者出国，单单小的竞赛还是不够的，ACM如果不出结果会有很大压力，估计保牛校无望，我校牛逼的学长学姐，也就是清华或者交大的硕博连读或者博士，这对于92年出生的我实在不能忍......

还有西安一站，这几天我会全力最后一搏，下次比赛绝不犯此次比赛的错误，如果自己的思路队友这边难以沟通，必要的话，就自己坚持自己想自己做，尤其自信点，有很多人其实不像他们想象的那样强，我自己也不想我想的那么弱

E[i][j][k]  已经占据i行j列，走了k步的时候，还需要走的步数的期望，其他看代码就好

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
#define CL(a,b) memset(a,b,sizeof(a))
#define IN(s) freopen(s,"r",stdin)
const int MAXN = 55;
double dp[MAXN][MAXN][MAXN*MAXN];

int main()
{
    //IN("D.txt");
    int ncase;
    scanf("%d",&ncase);
    int n,m;
    while(ncase--)
    {
        scanf("%d%d",&n,&m);
        CL(dp,0);
        for(int i=n;i>=0;i--)
            for(int j=m;j>=0;j--)
            {
                if(i==n&&j==m)continue;
                for(int k=i*j;k>=max(i,j);k--)
                {
                    dp[i][j][k]+=1.0*j*(n-i)/(1.0*n*m-k)*dp[i+1][j][k+1];
                    dp[i][j][k]+=1.0*i*(m-j)/(1.0*n*m-k)*dp[i][j+1][k+1];
                    dp[i][j][k]+=1.0*(i*j-k)/(1.0*n*m-k)*dp[i][j][k+1];
                    dp[i][j][k]+=1.0*(n-i)*(m-j)/(1.0*n*m-k)*dp[i+1][j+1][k+1];
                    dp[i][j][k]+=1.0;
                }
            }
        printf("%.12lf\n",dp[0][0][0]);
    }
    return 0;
}