ZOJ 3777 Problem Arrangement (状态压缩 + 概率)

OJ题目：click here ~~

题目分析：给每个题目安排在每个位置的value。有一个程序随机选择安排的顺序，总的value值大于等于m时，就可以接受这个安排。问能够获得一次满足条件的安排的期望次数。

由于程序是每次随机选择一个排列，每次的选择之间不影响，而且每次选中的概率相等，都为dp[( 1 << n) - 1 ][ m ]/ fac[ n ] , 即满足条件的总数/ 总排列数。由于01分布的期望为p ， 这里p = dp[( 1 << n) - 1 ][ m ]/ fac[ n ]  ， p为试验一次事件发生的概率，那么，事件一定发生试验的次数为1/p 。 所以结果为  fac[ n ]/dp[( 1 << n) - 1 ][ m ]  ，结果用最简形式表示。

对于状态的表示，用状态压缩。 

AC_CODE

int n , m;
int dp[1<<12][502];
int x[13][13];
int fac[13];
int gcd(int x , int y)
{
    return y == 0 ? x : gcd(y , x%y);
}

void DP()
{
    int row , i , j , k;
    dp[0][0] = 1;
    for(i = 0;i < (1<<n);i++)//枚举所有的状态
    {
        for(row = 0 , j = 0;j < n;j++)//计算这个状态已经取的行数
            if(i & (1<<j)) row++;
        for(j = 0;j < n;j++)//枚举下一行所取的列数
        {
            if(i&(1<<j)) continue;//与上面的列数选择不能有冲突
            for(k = 0;k <= m;k++)//记录每种状态下的情况数
            {
                int t = k + x[row + 1][j];//计算新的和
                if(t >= m)//和大于等于m的算为一种
                {
                    dp[i|(1<<j)][m] += dp[i][k];
                }
                else
                    dp[i|(1<<j)][t] += dp[i][k];
            }
        }
    }
}
int main()
{
    //freopen("in.txt","r",stdin);
    int t, i , j;
    cin >> t;
    fac[0] = 1;
    for(i = 1;i <= 13;i++)
        fac[i] = fac[i - 1] * i;
    while(t--)
    {
        cin >> n >> m;
        memset(dp , 0 , sizeof(dp));
        for(i = 1;i <= n;i++)
            for(j = 0;j < n;j++)
            scanf("%d",&x[i][j]);
        DP();
        if(dp[(1<<n)-1][m] == 0) printf("No solution\n");
        else
        {
            int c = gcd(fac[n] , dp[(1<<n) - 1][m]);
            printf("%d/%d\n",fac[n]/c , dp[(1<<n)-1][m]/c);
        }

    }
    return 0;
}