Marriage Ceremonies LightOJ

Marriage Ceremonies LightOJ - 1011

常规状压dp。popcount(S)表示S集合中元素数量。ans[S]表示S中的女性与前popcount(S)个男性结婚的最大收益。

那么，$ans[S]=max{ans[S-p]+a[popcount(S)][p]}$

老代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int ans[17][70000],a[17][17],T,TT,n;
int main()
{
    int i,j,k;
    scanf("%d",&T);
    for(TT=1;TT<=T;TT++)
    {
        scanf("%d",&n);
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
                scanf("%d",&a[i][j]);
        memset(ans,0,sizeof(ans));
        for(i=1;i<=n;i++)
            for(j=0;j<(1<<n);j++)
            {
                if(__builtin_popcount(j)!=i)    continue;
                for(k=1;k<=n;k++)
                    if(j&(1<<(k-1)))
                        ans[i][j]=max(ans[i][j],ans[i-1][j^(1<<(k-1))]+a[i][k]);
            }
        printf("Case %d: %d
",TT,ans[n][(1<<n)-1]);
    }
    return 0;
}

新代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int ans[70000],a[17][17],T,TT,n;
int main()
{
    int i,j,k;
    scanf("%d",&T);
    for(TT=1;TT<=T;TT++)
    {
        scanf("%d",&n);
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
                scanf("%d",&a[i][j]);
        memset(ans,0,sizeof(ans));
        for(i=1;i<(1<<n);i++)
            for(j=1;j<=n;j++)
                if(i&(1<<(j-1)))
                    ans[i]=max(ans[i],ans[i^(1<<(j-1))]+a[__builtin_popcount(i)][j]);
        printf("Case %d: %d
",TT,ans[(1<<n)-1]);
    }
}

另：貌似可以用二分图最大完美匹配做