题意:n道题目放在n个位置,任意题目放在任意位置会得到一个分数,求总分>=m的方案数;也可以看成在n*n矩阵中选n个不同行不同列的数。
分析:dp[i][s][k]表示前i个数放在s集合的这些位置得分为k的方案数,s是现在各位置的状态,可以省略i这一维,s已经包含前多少个位置的信息了,而且可以由s递推之后的状态,这是因为由s递推其他状态时,s已经被递推完成了。
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
int n,m;
int d[15][15];
int dp[(1<<12)][505];
int num[(1<<12)];
long long gcd(long long a,long long b)
{
return b==0?a:gcd(b,a%b);
}
int main()
{
memset(num,0,sizeof(num));
for(int i=0;i<(1<<12);i++)
{
int l=0;
for(int j=0;j<12;j++)
{
if(i&(1<<j))
l++;
}
num[i]=l;
}
int t;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&n,&m);
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
scanf("%d",&d[i][j]);
memset(dp,0,sizeof(dp));
dp[0][0]=1;
for(int i=0;i<(1<<n);i++)
{
for(int j=0;j<n;j++)
{
if((1<<j)&i)
continue;
for(int k=0;k<=m;k++)
{
if(k+d[num[i]][j]>m)
dp[i|(1<<j)][m]+=dp[i][k];
else
dp[i|(1<<j)][k+d[num[i]][j]]+=dp[i][k];
}
}
}
long long ans=1;
for(int i=n;i>=2;i--)
ans*=i;
long long u=dp[(1<<n)-1][m];
if(u==0)
{
printf("No solution\n");
continue;
}
long long b=gcd(ans,u);
printf("%lld/%lld\n",ans/b,u/b);
}
return 0;
}