网上多是一行一行处理的,其实因为距离为2,可以把图分为奇数点和偶数点,分成两张图,这两张图的点互不影响
用插头dp dp[i][j] 表示前i个点能到达j的状态,判断i+1个点时根据这个点是否能放士兵,
以及上方位置,左上方位置,右上方位置,左方位置有无士兵来决定这个位置是否可以放置士兵。
奇偶图分别dp,每次dp复杂度 = (n*m*2^m) 因为有效的状态并不多(100-),
还可以删除一些无用状态,优化的复杂度可以到(有效状态)*n*m
优化方法:map[i]表示i这个状态是否可达,用-1表示不可达,可达的话记录在state[]中的位置,
用state[]记录有效的状态,转移的时候如果转移的状态u可达则更新
插头dp示意图
0 1
0 1 0 0
0 0 X
00010010 表示现在对x所在的位置进行处理,在x之前的位置记录的i行,或者i-1行信息,
在x以及x之后的位置记录的是i-1,i-2行的信息,(体现奇偶图的优势了),然后判断x是否可以为1进行转移到
0 1
0 1 0 0
0 0 0
和
0 1
0 1 0 0
0 0 1
两个状态,但是第二个状态不合法,不用记录。。
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
int map[101][11];
int dp[2][1077];
int n,m;
int work(int f){
int p = 0 ,q = 1,ans = 0,u,v,i,j,k;
memset(dp[p],-1,sizeof(dp[p]));
dp[p][0] = 0;
for( i = 0;i < n; i++){
for( j = 0;j < m; j++){
if((i+j)%2 != f) continue;
for( k = (1<<m)-1;k >= 0;k--)
dp[q][k] = -1;
for( k = (1<<m)-1;k >= 0;k--){
if(dp[p][k] < 0) continue;
if(map[i][j]){
v = (1<<j);
if((!(k&v)) && (!(k&(v<<1))) && (!(k&(v>>1))) && (!(k&(v>>2)))){
u = k|v;
if(dp[q][u] < dp[p][k]+1)
dp[q][u] = dp[p][k]+1;
}
}
u = k;
if(u&(1<<j)) u -= (1<<j);
if(dp[q][u] < dp[p][k]) dp[q][u] = dp[p][k];
}
p^=1,q^=1;
}
}
for( k = (1<<m)-1; k >= 0;k--)
if(dp[p][k] > ans) ans = dp[p][k];
return ans;
}
int main(){
while(scanf("%d%d",&n,&m)!=EOF){
for(int i = 0;i < n; i++)
for(int j = 0;j < m; j++)
scanf("%d",&map[i][j]);
int ans = work(0) + work(1);
printf("%d\n",ans);
}
return 0;
}
/*
4 4
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1
*/