题目链接:http://poj.org/problem?id=1185
作为状态压缩的入门题的话,这一题的确是很基础,唯一需要注意的就是一些边界的范围。好吧,代码的注释写的比较详细,状态压缩的入门可以看看
http://blog.csdn.net/accry/article/details/6607703,博主绝对是业界良心
好吧,我英文不好,state是州,郡的意思,我给记成状态了,其实应该是statu,不管了,改起来略麻烦,就不管这些细节了╮(╯▽╰)╭
#include<stdio.h>
#include<string.h>
#define max(x,y) x<y?y:x
#define MEM(x,y) memset(x,y,sizeof(x))
int dixing[110]; //dixing[k]存储k行的状态
int dp[110][100][100]; //dp[i][p][q]表示第i行状态为state[p],第i-1行状态为state[q]时前i行的最大炮兵
int state[100],sum[100]; //state[]用来存储2的m次方里的合格状态,sum[i]存对应state[i]中的1的个数
char map[110][15];
int m,n,num; //num用来存2的m次方共有多少种合格的状态
bool ok(int x) //判断是否造成了行冲突
{
if(x & (x<<1))return false;
if(x & (x<<2))return false;
return true;
}
bool fit(int x,int k) //判断状态x与第k行的状态是否匹配
{
if(dixing[k]&x) return false; //不匹配
return true;
}
int count1(int x) //计算x的二进制中有多少个1
{
int cnt = 0;
while(x > 0 )
{
if(x&1) //如果最后一位是1,就加一
cnt++;
x = x>>1;; //然后x右移一位
}
return cnt;
}
void init()
{
num = 1;
for(int i = 0 ; i < (1<<m) ; i++)
{
if(ok(i))
{
state[num] = i;
sum[num] = count1(i);
num++;
}
}
}
int main()
{
while(scanf("%d%d",&n,&m) != EOF)
{
getchar();
MEM(map,0);
MEM(dp,-1);
MEM(dixing,0);
MEM(state,0);
for(int i = 0; i < n; i ++)
{
for(int j = 0; j < m ; j ++)
scanf("%c",&map[i][j]);
getchar();
}
init();
for(int i = 0; i < n;i++) //求每一行的地形状态
{
for(int j = 0; j < m;j++)
{
if(map[i][j] == 'H')
dixing[i] += 1<<(j); //把第i行倒过来作为二进制的数
}
}
for(int i = 1 ; i <= num;i++) //初始化dp
{
if(fit(state[i],0)) //state[i]不能与第0行的地形冲突
dp[0][i][1] = sum[i]; //第0行为i状态,第-1行必须为state[1](即全为0)状态,因为它不能与第0行的任意一种情况冲突
}
for(int i = 1; i < n; i++)
{
for(int p = 1 ;p <= num;p++) //枚举第i行的状态
{
if(!fit(state[p],i)) continue; //判断状态state[p]与第i行的地形状态是否匹配
for(int q = 1 ;q <= num; q++) //枚举第i-1行的状态
{
if(state[p] & state[q]) continue; //判断第i行的状态state[p]与第i-1行的状态state[q]是否匹配
for(int k = 1 ; k <= num; k ++) //枚举第i-2的状态
{
if(state[q]&state[k]) continue;
if(state[p]&state[k]) continue; 判断第i行的状态state[p]与第i-1行的状态state[q]是否匹配
if(dp[i-1][q][k] == -1) continue;
dp[i][p][q] = max(dp[i][p][q],dp[i-1][q][k]+sum[p]);
}
}
}
}
int ans = 0;
for(int j = 1 ; j <=num; j ++)
for(int k = 1; k <=num ; k++)
ans = max(ans,dp[n-1][j][k]); //共有n-1行
printf("%d\n",ans);
}
return 0;
}