ZOJ 3288 Domination

这个题只要设出数组来， 这个题就比较好做了。

设d【i】【j】【k】 为  i 行 j列 k个石子的时候概率。

d【i】【j】【k】 = p1*d【i-1】【j】【k-1】 + p2 * d【i】【j-1】【k-1】 + p3 * d【i-1】【j-1】【k-1】 + p4*d【i】【j】【k-1】；

一个状态由四个可能的状态转化而来。

p1   p2  p3 p4 都是很好算的。

当 i == m   j == n时   也就是说 当石子放的正好 每一列 每一行都有石子的时候 就不放了。 所以此时的转化方程应该没有 p4*d【i】【j】【k-1】  这一项

#include <cstdio>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <string>
#include <map>
#include <vector>
#include <set>
#include <queue>
#include <stack>
#include <cctype>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define MAXN 3
#define INF 1<<30
#define mod 1000007
double  d[51][51][2501];
int main (){
    int t;
    scanf("%d",&t);
    while(t--){
        int m,n;
        scanf("%d%d",&m,&n);
        memset(d,0,sizeof(d));
        d[0][0][0] = 1;
        for(int i = 1; i <= m; i++){
            for(int j = 1; j <= n; j++){
                for(int k = 1; k <= m*n; k++){
                    if(i == m && j == n)
                    d[i][j][k] = d[i-1][j][k-1]*(j)*(m-i+1)/(m*n-k+1)+d[i][j-1][k-1]*(i)*(n-j+1)/(m*n-k+1)+d[i-1][j-1][k-1]*(m-i+1)*(n-j+1)/(m*n-k+1);
                    else
                    d[i][j][k] = d[i-1][j][k-1]*(j)*(m-i+1)/(m*n-k+1)+d[i][j-1][k-1]*(i)*(n-j+1)/(m*n-k+1)+d[i-1][j-1][k-1]*(m-i+1)*(n-j+1)/(m*n-k+1) + d[i][j][k-1]*(i*j-(k-1))/(m*n-k+1);
                }
            }
        }
        double sum = 0;
        for(int i = 1; i <= m*n; i++){
            sum += d[m][n][i]*i;
        }
        printf("%.12lf\n",sum);
    }
    return 0;
}