ZOJ_Magazine Delivery DP

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=1428

题意：用三辆车给n个地方送杂志，初始时三辆车都在位置1处，且位置1处能提供足够的杂志，没辆车能携带的杂志数是无限制的，杂志必须按照结点的顺序依次送到每个结点，任一时刻必须只有一辆车在送杂志。求最少耗时。

分析：状态DP，很容易就能想到用一个30位的二进制来表示车辆的停放情况，但是由于n<=30，这样显然不仅会超时，也会超内存。 这时我们换个角度考虑问题，因为每个状态需要记录的车辆信息只有三辆车，且每辆车最多只有n中选择，因此考虑用3为的30进制数记录车辆停放情况（考虑到运算方便，这里实际用了32进制），这样记录车辆情况的状态 j 就由原来的2^30 降低到32*32*32。时间和空间都可以接受。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>

#define min(a,b) (a>b?b:a)
using namespace std;

int n;
int map[31][31];
int dp[31][32*32*32] ;

int DP()
{
	int MA = n*32*32 + n*32+ n ;
	memset(dp,-1,sizeof(dp));	//-1表示该状态不合法 
	dp[1][1057] = 0 ;			//初始化 
	for(int i=2;i<=n;i++)
	{
		for(int j=1057;j<=MA;j++)
		{
			if(dp[i-1][j] == -1)	continue ;			//前一状态不合法，无法转移，跳过 
			int po1 = j%32,po2 = j/32%32 ,po3=j/(32*32);		//三辆车所在的位置。 
			if(po1>=i || po2>=i || po3 >=i)	continue ;
			//依次枚举三辆车，哪辆车负责来送杂志到i地 。 
			int new_s = po3*32*32 + po2*32 + i ;
			if(dp[i][new_s] == -1)
				dp[i][new_s] = dp[i-1][j] + map[po1][i] ;
			else
				dp[i][new_s] = min(dp[i][new_s] , dp[i-1][j] + map[po1][i]) ;
			new_s = po3*32*32 + i*32 + po1 ;
			if(dp[i][new_s] == -1)
				dp[i][new_s] = dp[i-1][j] + map[po2][i] ;
			else
				dp[i][new_s] = min(dp[i][new_s],dp[i-1][j] + map[po2][i]) ;
			new_s = i*32*32 + po2*32 + po1 ;
			if(dp[i][new_s] == -1)
				dp[i][new_s] = dp[i-1][j] + map[po3][i] ;
			else
				dp[i][new_s] = min(dp[i][new_s],dp[i-1][j] + map[po3][i]) ;
		}	
	}
	int min_ = 0x7fffffff ;
	for(int j=1057;j<=MA;j++)
	{
		if(dp[n][j] == -1)	continue ;
		min_ = min(min_,dp[n][j]) ;	
	}
	return min_ ;
}

int main()
{
    //freopen("1in.txt","r",stdin);
    //freopen("1out.txt","w",stdout);
	int T ;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&n);	
		for(int i=1;i<=n;i++)
		{
			for(int j=i+1;j<=n;j++)
			{
				int a;
				scanf("%d",&a);
				map[i][j] = a ;
				map[j][i] = a ;	
			}	
		}	
		int ans = DP();
		printf("%d\n",ans);
	} 



    return 0;
}