题目传送门:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3777
题意:给出n道题目以及每一道题目在不同顺序做的时候的兴趣值,让你求出所有的做题顺序与所有做题顺序中兴趣值不小于m的比例。按一个最简分数表示。
分析:
全排列肯定会超时,所以只能用状压。但是想状压的时候,我想到的只是记录已经安排好了几个做题,但不知道怎么确定所有的n的题目中是哪几个被安排好了。所有就被卡主了。其实把安排了几个题当做状态就可以,不需要去管到底安排好的到底是前面的题还是后面的题,然后用所有待选题目中的前几个去满足上面的那几个安排好了的题就行,状态的转移就相当于了一个全排列的过程。
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
typedef long long ll;
#define pii pair<int, int>
int dp[1<<13][510];
int f[15], a[15][15];
bool vis[1<<13];
void init()
{
f[1] = 1;
for(int i = 2; i <= 12; i++)
f[i] = f[i-1]*i;
}
int main()
{
int T, n, m;
init();
scanf("%d", &T);
while(T--)
{
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i++)
for(int j = 1; j <= n; j++)
scanf("%d", &a[i][j]);
memset(dp, 0, sizeof(dp));
memset(vis, false, sizeof(vis));
dp[0][0] = 1;
queue<pii> q;
q.push(pii(0, 0));
while(!q.empty())
{
pii now = q.front();
q.pop();
int s = now.first;
int cnt = now.second;
for(int i = 1; i <= n; i++)
{
if(s>>(i-1) & 1)
continue;
for(int v = 0; v <= m; v++)
{
if(dp[s][v] == 0)
continue;
if((v + a[cnt+1][i]) >= m)
continue;
dp[s|(1<<(i-1))][v+a[cnt+1][i]] += dp[s][v];
if(!vis[s|(1<<(i-1))] && cnt + 1 != n)
{
q.push(pii(s|(1<<(i-1)), cnt+1));
vis[s|(1<<(i-1))] = true;
}
}
}
}
int ans = 0, all = f[n];
for(int i = 0; i < m; i++)
ans += dp[(1<<n)-1][i];
ans = all - ans;
int gcd = __gcd(all, ans);
all /= gcd;ans /= gcd;
if(ans == 0)
printf("No solution\n");
else
printf("%d/%d\n", all, ans);
}
return 0;
}