ZOJ 3822 Domination

题目：

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

题意：

随机往一个n*m的棋盘里放棋子，求每行每列都有棋子所需的棋子数的期望。

思路：

概率dp，然而是强行暴力所有情况出来的跑了2s多，dp[i][x][y]表示放了i枚棋子有x行有棋子y列有棋子的概率，然后从dp[i-1][x][y],dp[i-1][x-1][y],dp[i-1][x][y-1],dp[i-1][x-1][y-1],4种情况中强行暴力求解就行了，，，这么暴力难怪变成全场题，

代码：

#define N 51
int n,m;
double dp[N*N][N][N],ans;
int main()
{
    int i,j,k,kk,t,x,y,z;
    scanf("%d",&k);
    while(k--)
    {
        scanf("%d%d",&n,&m);
        memset(dp,0,sizeof(dp));
        dp[1][1][1]=1.0;
        for(i=2;i<=n*m;i++)
        for(x=1;x<=n&&x<=i;x++)
        for(y=1;y<=m&&y<=i;y++)
        {
            if((x<n||y<m)&&x*y>=i)dp[i][x][y]+=dp[i-1][x][y]*(x*y-i+1)/(n*m-i+1);
            dp[i][x][y]+=dp[i-1][x-1][y]*((n-x+1)*y)/(n*m-i+1);
            dp[i][x][y]+=dp[i-1][x][y-1]*((m-y+1)*x)/(n*m-i+1);
            dp[i][x][y]+=dp[i-1][x-1][y-1]*(n-x+1)*(m-y+1)/(n*m-i+1);
        }
        for(ans=0,i=0;i<=n*m;i++) ans+=i*dp[i][n][m];
        printf("%.12lf\n",ans);
    }
    return 0;
}