标签：动态规划   状态压缩   dp   算法   

方格取数(1)

3
75 15 21 
75 15 28 
34 70 5 

188

Poj - 3254 Corn Fields

用dp[i][j]表示前i行，第i行选第j种状态时的最优解，

首先找出所有本行不冲突的状态（即这一行中没有相邻的情况），存入state数组

计算出第i行取第j种状态时可得到的数值stn[i][j]

那么dp[i][j]=max{dp[i][j],dp[i-1][k]+stn[i][j]} (k表示第i-1行取第k种状态

最终答案即为dp[n][j]中的最大值。

代码：

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int hpn=18000;
int state[hpn],stn[25][hpn],dp[25][hpn];
//dp[i][j]:前i行，第i行选第j种状态时的最优解
int mst,map[25][25];//第i行选第j种状态时的值

inline int bet(int x,int y)
{
	if(x>y)	return x;
	return y;
}

void find_all_state(int n)
{
	memset(state,0,sizeof(state));
	mst=0;//最多有多少种状态
	int lin=(1<<n),index=1;
	for(int i=0;i<lin;++i)
	{
		if((i&(i<<1))==0)
		{
			state[index]=i;
			++mst;
			++index;
		}
	}
}

int main()
{
	int n;
	while(cin>>n)
	{
		if(n==0)
		{
			cout<<0<<endl;
			continue;
		}
		find_all_state(n);

		for(int i=1;i<=n;++i)
		{
			for(int j=1;j<=n;++j)
			{
				scanf("%d",&map[i][j]);
			}
		}

		memset(stn,0,sizeof(stn));
		for(int i=1;i<=n;++i)
		{
			for(int j=1;j<=mst;++j)
			{
				int b=n;
				for(int t=1;t<(1<<n);t=(t<<1))
				{
					if((t&state[j])!=0)
					{
						stn[i][j]+=map[i][b];
					}
					--b;
				}
			}
		}//统计第i行选第j种状态时可得到的值

		memset(dp,0,sizeof(dp));
		for(int i=1;i<=mst;++i)//第一行取第i种状态时
		{
			dp[1][i]=stn[1][i];
		}
		for(int i=2;i<=n;++i)//对于第i行
		{
			for(int j=1;j<=mst;++j)//取第j种状态
			{
				for(int k=1;k<=mst;++k)//第i-1行取第k种状态
				{
					if((state[j]&state[k])!=0)	continue;
					dp[i][j]=bet(dp[i][j],dp[i-1][k]+stn[i][j]);
				}
			}
		}

		int ans=0;
		for(int i=1;i<=mst;++i)
		{
			ans=bet(ans,dp[n][i]);
		}

		cout<<ans<<endl;
	}
	return 0;
}

Hdu-1565 方格取数(1) （状态压缩dp入门题,布布扣,bubuko.com

Hdu-1565 方格取数(1) （状态压缩dp入门题

标签：动态规划   状态压缩   dp   算法   

原文地址：http://blog.csdn.net/harrypoirot/article/details/26023519