HDU 1565【状压DP】

今天脑袋不清楚了一天。。。TAT。。。这种题还调了那么久。。TAT。。。

之前有想过用贪心，因为每个数都是非负数，那么自然是每行拿的越多越好，就尝试0101或者1010这样的序列就好了，但是可能出现这种情况，上一行拿了0101的序列，这一行拿1010，有可能下一行1010的取值配上上一行0101的取值是最大的情况，或者是这一行01001更佳，所以还是要用状压DP。。

之前忘了初始化。。醉了。。刚开始的时候T，于是就把符合情况的二进制数先处理出来了，一共一万多个，这样复杂度也是可以忍受的。。

#include <stdio.h>
#include <iostream>
#include <string.h>
using namespace std;
#define maxn 1048578
int is[maxn],que[17714],vis[2][17714];
int a[21][21],top;
int n;
int dp[2][maxn];
int max(int x,int y){return x>y?x:y;}
void init()
{
    top=0;
    for(int i=0;i<maxn;i++) is[i]=1;
    for(int i=0;i<maxn;i++)
    {
        int j=i;
        while(j)
        {
            if(j&1)
            {
                if((j>>1)&1) {is[i]=0;break;}
                else j>>=1;
            }
            else
                j>>=1;
        }
        if(is[i])
        {
            que[top++]=i;
        }
    }
}
int getSum(int r,int x)
{
    int ret=0;
    for(int i=0;i<n&&x;i++)
    {
        if(x&1) ret+=a[r][i];
        x>>=1;
    }
    return ret;
}
void search()
{
    int maxnum=1<<n;
    int flag=0;
    memset(vis,1,sizeof(vis));
    memset(dp,0,sizeof(dp));
    for(int i=0;i<n;i++)
    {
        flag=i&1;
        for(int j=0;j<top&&que[j]<maxnum;j++)
        {
            int now=que[j];
            int sum=getSum(i,now);
            for(int k=0;k<top&&que[k]<maxnum;k++)
            {
                int pre=que[k];
                if(now&pre||!vis[!flag][k]) continue;
                dp[flag][now]=max(dp[flag][now],dp[!flag][pre]+sum);
                vis[flag][j]=1;
            }
        }
        for(int j=0;j<top&&que[j]<maxnum;j++) dp[!flag][que[j]]=0,vis[!flag][j]=0;
    }
    int ans=0;
    for(int j=0;j<top&&que[j]<maxnum;j++)
        ans=max(ans,dp[flag][que[j]]);
    cout<<ans<<endl;
}
int main()
{
    init();
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                scanf("%d",&a[i][j]);
        search();
    }
    return 0;
}