ZOJ 3822 可能性DP

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

本场比赛之前,我记得。见WALK概率路DP称号。那么它应该是可以考虑的概率DP,十一还特意看了碍着,

当场景。真的OUT了。然后好激动的样子,開始推得二维。然后感觉好难推。发现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;
}


版权声明:本文博客原创文章,博客,未经同意,不得转载。

posted @ 2015-07-25 19:00  phlsheji  阅读(210)  评论(0编辑  收藏  举报