HDU 4539 状压DP

这道题是腾讯马拉松的题，中文题，就不解释题意了。

状压DP，同炮兵阵地，写着练练手。

攻击范围是曼哈顿距离2，画个图能知道是这样的

       x

    x    x

x     o     x

下半部分就不画了。其实就是裸题。

 

#include <set>

#include <map>

#include <stack>

#include <cmath>

#include <queue>

#include <cstdio>

#include <string>

#include <vector>

#include <iomanip>

#include <cstring>

#include <iostream>

#include <algorithm>

#define Max 2505

#define FI first

#define SE second

#define ll long long

#define PI acos(-1.0)

#define inf 0x3fffffff

#define LL(x) ( x << 1 )

#define bug puts("here")

#define PII pair<int,int>

#define RR(x) ( x << 1 | 1 )

#define mp(a,b) make_pair(a,b)

#define mem(a,b) memset(a,b,sizeof(a))

#define REP(i,s,t) for( int i = ( s ) ; i <= ( t ) ; ++ i )



using namespace std;



inline void RD(int &ret) {

    char c;

    int flag = 1 ;

    do {

        c = getchar();

        if(c == '-')flag = -1 ;

    } while(c < '0' || c > '9') ;

    ret = c - '0';

    while((c=getchar()) >= '0' && c <= '9')

        ret = ret * 10 + ( c - '0' );

    ret *= flag ;

}



inline void OT(int a) {

    if(a >= 10)OT(a / 10) ;

    putchar(a % 10 + '0') ;

}



/*********************************************/

#define N 1005

#define M 11

int n , m ;

int Map[111][11] ;

int MM[111] ;

int st[170] ;

int cnt[170] ;

int num = 0 ;

int dp[105][170][170] ;//j上一状态,k当前状态.

void ok(){

    for (int i = 0 ; i < 1 << m ; i ++ ){

        if((i >> 2) & i) continue ;

        if((i << 2) & i) continue ;

        int x = i ;

        int nn = 0 ;

        while(x){

            nn += x & 1 ;

            x >>= 1 ;

        }

        st[num] = i ;

        cnt[num ++] = nn ;

    }

}

int main() {

//    for (int i = 1 ; i <= 10 ;i ++ ){

//        num = 0 ;

//        m = i ;

//        ok() ;

//        cout << num << endl ;

//    }

    while(cin >> n >> m){

        mem(MM ,0) ;

        mem(dp ,0) ;

        num = 0 ;

        for (int i = 1 ; i <= n ; i ++ ){

            for (int j = 1 ; j <= m ; j ++ ){

                RD(Map[i][j]) ;

                if(!Map[i][j])MM[i] += 1 << ( j - 1 ) ;

            }

        }

        ok() ;

        int ans = 0 ;

        //预处理第一行

        for (int i = 0 ; i < num ; i ++ ){

            if(st[i] & MM[1])continue ;

            dp[1][0][i] = max(dp[1][0][i] , cnt[i]) ;

            ans = max(ans , dp[1][0][i]) ;

        }

        for (int i = 2 ; i <= n ; i ++ ){

            for (int j = 0 ; j < num ; j ++ ){

                if(MM[i - 1] & st[j])continue ;

                for (int k = 0 ; k < num ; k ++ ){

                    if(MM[i] & st[k])continue ;

                    if(st[k] >> 1 & st[j])continue ;

                    if(st[k] << 1 & st[j])continue ;

                    for (int l = 0 ; l < num ; l ++ ){

                        if(MM[i - 2] & st[l])continue ;

                        if(st[k] & st[l])continue ;

                        dp[i][j][k] = max(dp[i][j][k] , dp[i - 1][l][j] + cnt[k]) ;

                        ans = max(dp[i][j][k] , ans) ;

                    }

                }

            }

        }

        OT(ans) ;puts("") ;

    }

    return 0 ;

}