zoj3822 Domination 概率dp --- 2014 ACM-ICPC Asia Mudanjiang Regional Contest

一个n行m列的棋盘，每次可以放一个棋子，问要使得棋盘的每行每列都至少有一个棋子 需要的放棋子次数的期望。

dp[i][j][k]表示用了k个棋子共能占领棋盘的i行j列的概率。

那么对于每一颗棋子，在现有的棋盘上，它可能有四种影响：新占了一行，新占了一列，既占了新的一行又占了新的一列，无影响。

对于每一种情况，dp[i][j][k]=原始状态的概率×选到这样的位置的概率

最后算答案的时候注意到，题目问的是到第k个刚好放完n行m列，也就是最后一个棋子一定是有作用的，

所以用dp[i][j][k]-dp[i][j][k-1]得到是第k个棋子恰好使得每行每列都占领的概率。

#include<cstdio>
#include<cstring>
double dp[55][55][2550];
int n,m;
int main()
{
    int T,i,j,k;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        int sum=n*m;
        for(i=0;i<=n;i++)
            for(j=0;j<=m;j++)
                for(k=0;k<=sum;k++) dp[i][j][k]=0;
        dp[0][0][0]=1.0;
        for(k=1;k<=sum;k++)
            for(i=1;i<=n;i++)
                for(j=1;j<=m;j++)
                {
                    dp[i][j][k]+=(dp[i][j][k-1]*(i*j-k+1)*1.0/(sum-k+1));//添加的位置没有新增新行或新列
                    dp[i][j][k]+=(dp[i-1][j][k-1]*((n-i+1)*j)*1.0/(sum-k+1));//增加行不增加列
                    dp[i][j][k]+=(dp[i][j-1][k-1]*(m-j+1)*i*1.0/(sum-k+1));//增加列不增加行
                    dp[i][j][k]+=(dp[i-1][j-1][k-1]*(n-i+1)*(m-j+1)*1.0/(sum-k+1));//既增加列也增加行
                 //   printf("i:%d j;%d k;%d dp:%lf\n",i,j,k,dp[i][j][k]);
                }
        double ans=0;
        for(k=1;k<=sum;k++) ans+=(dp[n][m][k]-dp[n][m][k-1])*k;
        printf("%.10lf\n",ans);
    }
    return 0;
}