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

时间：2023-03-09 16:10:11
[状压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

相关文章