转载请注明出处: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;
}