POJ 1185 炮兵阵地

题目链接：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;
}