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;
}