码迷,mamicode.com
首页 > 其他好文 > 详细

2014牡丹江——Domination

时间:2014-10-13 16:42:19      阅读:151      评论:0      收藏:0      [点我收藏+]

标签:acm   dp   

题目链接

  • 题意:
    给一个n*m的矩阵,每天随机的在未放棋子的格子上放一个棋子。求每行至少有一个棋子,每列至少有一个棋子的天数的期望
     (1 <= N, M <= 50). 
  • 分析:
    比较明显的概率DP,难点在于如何设计状态。覆盖了多少行和列是必不可少的,之后比较关键的就是想到另一个属性:多少个交叉点(即放过的点)
double dp[55][55][2550];
bool vis[55][55][2550];
int tot;
int n, m;
inline double getp(int xx, int yy)
{
    return (xx * 1.0) / yy;
}
double dpf(int x, int y, int z)
{
    if (x > n || y > m || z > tot) return 0;
    if (x == n && y == m) return 0.0;

    if (vis[x][y][z]) return dp[x][y][z];
    vis[x][y][z] = true;

    int a = x * y - z;
    int b = tot - (x * m + y * n - x * y);
    int xc = x * m - x * y;
    int yc = y * n - x * y;

    int sum = a + b + xc + yc;
    dp[x][y][z] = 0;

    if (a)
        dp[x][y][z] += getp(a, sum) * dpf(x, y, z + 1);
    if (b)
        dp[x][y][z] += getp(b, sum) * dpf(x + 1, y + 1, z + 1);
    if (xc)
        dp[x][y][z] += getp(xc, sum) * dpf(x, y + 1, z + 1);
    if (yc)
        dp[x][y][z] += getp(yc, sum) * dpf(x + 1, y, z + 1);
    dp[x][y][z] += 1.0;
    return dp[x][y][z];
}


int main()
{
    int T;
    RI(T);
    while (T--)
    {
        RII(n, m); tot = n * m;
        memset(vis, 0, sizeof(vis));
        printf("%.12lf\n", dpf(0, 0, 0));
    }
    return 0;
}


2014牡丹江——Domination

标签:acm   dp   

原文地址:http://blog.csdn.net/wty__/article/details/40044019

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!