[状压dp]HDOJ1565 方格取数(1)

中文题~~ 题意略

$n\le 20$ ! 很明显是状压!

 #include <cstdio>
 #include <cstdlib>
 #include <cstring>
 #include <climits>
 #include <cctype>
 #include <cmath>
 #include <string>
 #include <sstream>
 #include <iostream>
 #include <algorithm>
 #include <iomanip>
 using namespace std;
 #include <queue>
 #include <stack>
 #include <vector>
 #include <deque>
 #include <set>
 #include <map>
 typedef long long LL;
 typedef long double LD;
 #define pi acos(-1.0)
 #define lson l, m, rt<<1
 #define rson m+1, r, rt<<1|1
 typedef pair<int, int> PI;
 typedef pair<int, PI> PP;
 #ifdef _WIN32
 #define LLD "%I64d"
 #else
 #define LLD "%lld"
 #endif
 //#pragma comment(linker, "/STACK:1024000000,1024000000")
 //LL quick(LL a, LL b){LL ans=1;while(b){if(b & 1)ans*=a;a=a*a;b>>=1;}return ans;}
 //inline int read(){char ch=' ';int ans=0;while(ch<'0' || ch>'9')ch=getchar();while(ch<='9' && ch>='0'){ans=ans*10+ch-'0';ch=getchar();}return ans;}
 //inline void print(LL x){printf(LLD, x);puts("");}
 //inline void read(double &x){char c = getchar();while(c < '0') c = getchar();x = c - '0'; c = getchar();while(c >= '0'){x = x * 10 + (c - '0'); c = getchar();}}

 int mp[][];
 int dp[][<<];   // n格为一个状态
 int p[<<], d;
 void pre()       //先记录所有的合法状态 即不相邻
 {
     d=;
     for(int i=;i<(<<);i++)
         if((i & (i<<))==)
             p[d++]=i;
 }
 int main()
 {
 #ifndef ONLINE_JUDGE
     freopen("in.txt", "r", stdin);
     freopen("out.txt", "w", stdout);
 #endif
     pre();
     int n;
     while(~scanf("%d", &n))
     {
         for(int i=;i<n;i++)
             for(int j=;j<n;j++)
                 scanf("%d", &mp[i][j]);
         memset(dp, , sizeof(dp));
         int cur=;
         for(int i=;i<n;i++)
         {
             for(int k=;k<d;k++)
             {
                 int sum=;
                 if(p[k]>=(<<n))    // 算是优化吧 没有就tle了
                         break;
                 for(int j=;j<n;j++)
                     if(p[k] & (<<j))
                         sum+=mp[i][j];
                 for(int j=;j<d;j++)
                 {
                     if(p[j]>=(<<n))   // 算是优化吧 没有就tle了
                         break;
                     if((p[k] & p[j])==)
                         dp[cur][p[k]]=max(dp[cur][p[k]], dp[cur^][p[j]]+sum);
                 }
             }
             cur^=;
         }
         cur^=;
         int ans=;
         for(int i=;i<(<<n);i++)
             ans=max(ans, dp[cur][i]);
         printf("%d\n", ans);
     }
     return ;
 }

HDOJ 1565