HDU 4571——Travel in time（dp）

转载请注明出处：http://blog.csdn.net/u010734277

Bob厌倦了游戏，离开了Alice，独自一人去长沙，Yuelu山，世界之窗等是Bob想要游览的地方。但是他的时间很有限，他不能浏览它们全部。假设长沙有N个景点，Bob给每一个景点Si定义了一个满意值。如果他游览了这个经典，他的总满意值就会增加Si。Bob希望在有限的时间T里，他可以从景点S开始，选择性的游览一些景点，最终到达景点E，使得总满意值可以达到尽可能的最大。很明显游览一个景点同样需要一些时间，假设游览景点i(0<=i<N)需要花掉Ci个单位的时间。

记住，Bob可以选择经过一个景点而不游览它们（包括S和E），可能是他只想走更近的路线来省时间。

Bob同样有一个特殊的需求，他只会游览那些满意值比他上一个游览的景点大的景点。例如，如果他游览了一个满意值为50的经典，他就只会游览满意值是51或者更大的。当然，景点之间的路是双向的。

最开始的时候不会做，参考了别人的博客。首先用弗洛伊德算法求出任意两点之间的最短路径，然后动态规划出游览的最优解。动态规划方程是dp(i,k)=max{dp(i,k),dp(j,k-a[i].cost-map[i][j])+a[i].value}。其中i代表了当前所在的景点，k代表当前的时间。dp(i,k)及在k时刻到达地点i时总的最大的满意值，是当前已求得的dp(i,k)与从景点j到达i得到的满意值之中较大的一个。j是在i之前访问的景点。

其中要注意的是开dp数组的时候第二维最大值不是100而是300，开始开错了，然后WA了好多次。

代码如下：

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=100000000;

typedef struct point{
	int cost;
	int value;
	int num;
	bool operator < (const point &ii) const{
		if(value!=ii.value)return value<ii.value;
		return num<ii.num;
	}
}Point;

int map[105][105];
Point a[105];
int dp[105][305];
int hash[105];

int main(){
	ios::sync_with_stdio(false);
//	freopen("data.txt","r",stdin);
	int w;
	cin>>w;
	for(int time=1;time<=w;++time){
		int n,m,t,s,e;
		cin>>n>>m>>t>>s>>e;
		for(int i=0;i<n;++i){
			cin>>a[i].cost;
		}
		for(int i=0;i<n;++i){
			cin>>a[i].value;
			a[i].num=i;
		}
		sort(a,a+n);
		for(int i=0;i<n;++i){
			hash[a[i].num]=i;
		}
		s=hash[s];
		e=hash[e];
		for(int i=0;i<n;++i){
			for(int j=0;j<n;++j){
				if(i==j)map[i][j]=0;
				else map[i][j]=maxn;
			}
		}
		for(int i=0;i<m;++i){
			int u,v,l;
			cin>>u>>v>>l;
			u=hash[u];
			v=hash[v];
			map[u][v]=min(l,map[u][v]);
			map[v][u]=min(l,map[v][u]);
		}
		for(int u=0;u<n;++u){
			for(int i=0;i<n;++i){
				for(int j=0;j<n;++j){
					map[i][j]=min(map[i][j],map[i][u]+map[u][j]);
				}
			}
		}
		memset(dp,-1,sizeof(dp));
		for(int i=0;i<n;++i){
			for(int j=a[i].cost+map[i][s];j<=t;++j){
				dp[i][j]=a[i].value;
			}
		}
		for(int i=0;i<n;++i){
			for(int j=0;j<i;++j){
				if(a[j].value==a[i].value)break;
				for(int k=0;k<=t;++k){
					if(k-a[i].cost-map[i][j]<0)continue;
					if(dp[j][k-a[i].cost-map[i][j]]==-1)continue;
					dp[i][k]=max(dp[i][k],dp[j][k-a[i].cost-map[i][j]]+a[i].value);
				}
			}
		}
		int ans=0;
		for(int i=0;i<n;++i){
			for(int k=0;k<=t;++k){
				if(k+map[i][e]>t)break;
				ans=max(ans,dp[i][k]);
			}
		}
		cout<<"Case #"<<time<<':'<<endl;
		cout<<ans<<endl;
	}
	return 0;
}