[poj1185][NOI2001]炮兵阵地 状压dp

                                                       炮兵阵地

Time Limit: 2000MS Memory Limit: 65536KTotal Submissions:0 Accepted:0

Description

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成，地图的每一格可能是山地（用"H" 表示），也可能是平原（用"P"表示），如下图。在每一格平原地形上最多可以布置一支炮兵部队（山地上不能够部署炮兵部队）；一支炮兵部队在地图上的攻击范围如图中黑色区域所示：

如果在地图中的灰色所标识的平原上部署一支炮兵部队，则图中的黑色的网格表示它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。
现在，将军们规划如何部署炮兵部队，在防止误伤的前提下（保证任何两支炮兵部队之间不能互相攻击，即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内），在整个地图区域内最多能够摆放多少我军的炮兵部队。

Input

第一行包含两个由空格分割开的正整数，分别表示N和M；
接下来的N行，每一行含有连续的M个字符('P'或者'H')，中间没有空格。按顺序表示地图中每一行的数据。N <= 100；M <= 10。

Output

仅一行，包含一个整数K，表示最多能摆放的炮兵部队的数量。

Sample Input

5 4PHPPPPHHPPPPPHPPPHHP

Sample Output

6

Source

预处理每一行可以有的状态

dp[i][j][k]:前i行，i行状态为j，i-1行状态为k的最大解（注意状态是编号不是二进制押位，因为这道题就算全是平原也没有多少种所以要开编号不然目测会mle（然而我开的很大））

然后这不是水题吗？（开始看成求方案数想问怎么没有模数写完样例输出400多的zz）

#include <algorithm>#include <iostream>#include <cstring>#include <cstdlib>#include <cstdio>#include <bitset>#include <queue>#include <ctime>#include <map>#include <set>using namespace std;typedef long long ll;int q[105][1025], a[105], n, m, ans, dp[105][505][505];char ch[105][15];void init(){for( int j = 1; j <= n; j++ )for( int i = 0; i < (1<<m); i++ )if( !(i&(i<<1)) && !(i&(i<<2)) )if( !(i&a[j]) )q[j][++q[j][0]] = i;}int count( int x ){ int cnt = 0; while(x){ cnt++; x -= x&-x; } return cnt; }int main(){scanf( "%d%d", &n, &m );for( int i = 1; i <= n; i++ ) scanf( "%s", ch[i] );for( int i = 1; i <= n; i++ ) for( int j = 0; j < m; j++ ) if( ch[i][j] == 'H' ) a[i] += (1<<j);init();if( n == 1 ){for( int i = 1; i <= q[1][0]; i++ ) ans = max( ans, count(q[1][i]) );cout<<ans<<endl;return 0;}for( int i = 1; i <= q[2][0]; i++ )for( int j = 1; j <= q[1][0]; j++ )if( (q[1][j]&q[2][i]) == 0 )dp[2][i][j] = count(q[1][j]) + count(q[2][i]);for( int i = 3; i <= n; i++ )for( int j = 1; j <= q[i][0]; j++ )for( int k = 1; k <= q[i - 1][0]; k++ )if( (q[i][j]&q[i - 1][k]) == 0 )for( int l = 1; l <= q[i - 2][0]; l++ )if( (q[i][j]&q[i - 2][l]) == 0 && (q[i - 1][k]&q[i - 2][l]) == 0 )dp[i][j][k] = max( dp[i][j][k], dp[i - 1][k][l] + count(q[i][j]) );for( int i = 1; i <= q[n][0]; i++ )for( int j = 1; j <= q[n - 1][0]; j++ )ans = max( ans, dp[n][i][j] );cout<<ans<<endl;return 0;}/*5 4PHPPPPHHPPPPPHPPPHHP*/