hdu 2167 Pebbles （状态压缩DP）

题目点我点我点我

题目大意：给你一个矩阵，让你取其中的数使得这些数的和尽可能地大，并且取的数周围的8个都不能取。

解题思路：状态压缩DP，预处理出来哪些状态可以用的，然后枚举每种可行状态，最后max一下求得最大值即可。

dp[i][j]表示前i行状态为j的和。

转移方程为：dp[i][j] = max(dp[i][j],dp[i-1][k]+sum[i][j])，其中k为前一行的可行状态。

sum[i][j]表示第i行状态为j的数和。

/* ***********************************************
┆  ┏┓　　　┏┓ ┆
┆┏┛┻━━━┛┻┓ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　━　　　┃ ┆
┆┃　┳┛　┗┳　┃ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　┻　　　┃ ┆
┆┗━┓　马　┏━┛ ┆
┆　　┃　勒　┃　　┆　　　　　　
┆　　┃　戈　┗━━━┓ ┆
┆　　┃　壁　　　　　┣┓┆
┆　　┃　的草泥马　　┏┛┆
┆　　┗┓┓┏━┳┓┏┛ ┆
┆　　　┃┫┫　┃┫┫ ┆
┆　　　┗┻┛　┗┻┛ ┆
************************************************ */

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
using namespace std;

#define rep(i,a,b) for (int i=(a),_ed=(b);i<=_ed;i++)
#define per(i,a,b) for (int i=(b),_ed=(a);i>=_ed;i--)
#define pb push_back
#define mp make_pair
const int inf_int = 2e9;
const long long inf_ll = 2e18;
#define inf_add 0x3f3f3f3f
#define mod 1000000007
#define LL long long
#define ULL unsigned long long
#define MS0(X) memset((X), 0, sizeof((X)))
#define SelfType int
SelfType Gcd(SelfType p,SelfType q){return q==0?p:Gcd(q,p%q);}
SelfType Pow(SelfType p,SelfType q){SelfType ans=1;while(q){if(q&1)ans=ans*p;p=p*p;q>>=1;}return ans;}
#define Sd(X) int (X); scanf("%d", &X)
#define Sdd(X, Y) int X, Y; scanf("%d%d", &X, &Y)
#define Sddd(X, Y, Z) int X, Y, Z; scanf("%d%d%d", &X, &Y, &Z)
inline int read(){int ra,fh;char rx;rx=getchar(),ra=0,fh=1;while((rx<'0'||rx>'9')&&rx!='-')rx=getchar();if(rx=='-')fh=-1,rx=getchar();while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra*fh;}
//#pragma comment(linker, "/STACK:102400000,102400000")

int n,m;
int a[20][20];
int dp[20][1<<15];
int sum[20][1<<15];
int s[1<<15];
char str[100];


int solve()
{
    MS0(dp);
    MS0(sum);
    int cnt = 0;
    for(int i=0;i<(1<<n);i++)  //选出可行状态
    {
        if(i&(i<<1))continue;
        s[cnt++] = i;
    }
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<cnt;j++)
        {
            for(int k=0;k<n;k++)
            {
                if(s[j]&(1<<k))sum[i][j] += a[i][k];
            }
        }
    }
    for(int i=0;i<cnt;i++)dp[0][i] = sum[0][i];
    for(int i=1;i<n;i++)
    {
        for(int j=0;j<cnt;j++)      //枚举当前行可行状态
        {
            for(int k=0;k<cnt;k++)  //枚举前一行可行状态
            {
                if(s[j]&s[k])continue;
                if(s[j]&(s[k]>>1))continue;
                if(s[j]&(s[k]<<1))continue;
                dp[i][j] = max(dp[i][j],dp[i-1][k]+sum[i][j]);
            }
        }
    }
    int res = 0;
    for(int i=0;i<cnt;i++)
        res = max(res,dp[n-1][i]);
    return res;
}


int main()
{
	//freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	ios::sync_with_stdio(0);
	cin.tie(0);
	while(gets(str))
    {
        int len = strlen(str);
        n=0;
        for(int i=0; i<len; i+=3)
        {
            a[0][n++]=(str[i]-'0')*10+(str[i+1]-'0');
        }
        for(int i=1; i<n; i++)
        {
            int m=0;
            gets(str);
            for(int j=0; j<len; j+=3)
                a[i][m++]=(str[j]-'0')*10+(str[j+1]-'0');
        }
        printf("%d\n",solve());
        getchar();
    }


	return 0;
}