zoj 3777 Problem Arrangement

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=5264

题意：给出n道题目以及每一道题目不同时间做的兴趣值，让你求出所有做题顺序中兴趣值大于等于m的比例。用一个分数表示。

状压dp。 枚举每一个状态，用二进制表示。dp[i][j]表示第i个题目，兴趣值为j的个数。

转移方程 dp[i|(1<<j)][k+a[num][j]]+=dp[i][k];兴趣值大于m的为 dp[1|(1<<j)][m]+=dp[i][k];

#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL long long
using namespace std;

int t,n,m;
int p[15][15];
int dp[1<<13][1000];
int f[10000];

LL gcd(LL a,LL b)
{
    return b==0?a:gcd(b,a%b);
}

int main()
{
    scanf("%d",&t);
    f[1]=1;
    for(int i=2; i<=12; i++)
    {
        f[i]=f[i-1]*i;
    }
    while(t--)
    {
        scanf("%d%d",&n,&m);
        memset(dp,0,sizeof(dp));
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<n; j++)
            {
                scanf("%d",&p[i][j]);
            }
        }
        dp[0][0]=1;
        for(int i=0; i<=(1<<n); i++)
        {
            int num=0;
            for(int j=0; j<n; j++)
            {
                if(i&(1<<j)) num++;
            }
            for(int j=0; j<n; j++)
            {
                if(i&(1<<j)) continue;
                for(int k=0; k<=m; k++)
                {
                    if(k+p[num][j]>=m)
                    {
                        dp[i|(1<<j)][m]+=dp[i][k];
                    }
                    else
                    {
                        dp[i|(1<<j)][k+p[num][j]]+=dp[i][k];
                    }
                }
            }
        }
        if(dp[(1<<n)-1][m]==0)
        {
            printf("No solution
");
        }
        else
        {
            int g=gcd(f[n],dp[(1<<n)-1][m]);
            printf("%d/%d
",f[n]/g,dp[(1<<n)-1][m]/g);
        }
    }
    return 0;
}