首页 > 其他 > 详细

zoj 3777

时间:2014-04-23 07:06:12      阅读:578      评论:0      收藏:0      [点我收藏+]

状态压缩dp

bubuko.com,布布扣
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 4100
using namespace std;

int map[13][13];
int dp[maxn][510];

int gcd(int x,int y)
{
    return y==0?x:gcd(y,x%y);
}

int main()
{
    int t;
    int n,m;
    scanf("%d",&t);
    while(t--)
    {
        memset(dp,0,sizeof dp);
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                scanf("%d",&map[i][j]);
        dp[0][0]=1;
        int cnt,tmp;
        for(int i=0;i<(1<<n);i++)
        {
            cnt=0;
            for(int j=0;j<n;j++)
                if(i&(1<<j))cnt++;//更新到第cnt+1列了;
            for(int j=0;j<n;j++)
            {
                if(i&(1<<j))continue;
                for(int k=0;k<=m;k++)
                {
                    tmp=min(m,k+map[j][cnt]);
                    dp[i|(1<<j)][tmp]+=dp[i][k];
                }
            }
        }
        tmp=1;
        for(int i=1;i<=n;i++)
            tmp*=i;
        int ans=dp[(1<<n)-1][m];
        if(ans==0)
        {
            puts("No solution");
            continue;
        }
        cnt=gcd(tmp,ans);
        printf("%d/%d\n",tmp/cnt,ans/cnt);
    }
    return 0;
}
View Code

 

zoj 3777,布布扣,bubuko.com

zoj 3777

原文:http://www.cnblogs.com/yours1103/p/3680883.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!