作为一个城市的应急救援队伍的负责人,你有一张特殊的全国地图。在地图上显示有多个分散的城市和一些连接城市的快速道路。每个城市的救援队数量和每一条连接两个城市的快速道路长度都标在地图上。当其他城市有紧急求助电话给你的时候,你的任务是带领你的救援队尽快赶往事发地,同时,一路上召集尽可能多的救援队。
输入格式:
输入第一行给出4个正整数N、M、S、D,其中N(2≤N≤500)是城市的个数,顺便假设城市的编号为0 ~ (N−1);M是快速道路的条数;S是出发地的城市编号;D是目的地的城市编号。
第二行给出N个正整数,其中第i个数是第i个城市的救援队的数目,数字间以空格分隔。随后的M行中,每行给出一条快速道路的信息,分别是:城市1、城市2、快速道路的长度,中间用空格分开,数字均为整数且不超过500。输入保证救援可行且最优解唯一。
输出格式:
第一行输出最短路径的条数和能够召集的最多的救援队数量。第二行输出从S到D的路径中经过的城市编号。数字间以空格分隔,输出结尾不能有多余空格。
输入样例:
4 5 0 3
20 30 40 10
0 1 1
1 3 2
0 3 3
0 2 2
2 3 2
输出样例:
2 60
0 1 3
知识点
迪杰斯特拉算法(Dijkstra):是从一个顶点到其余各顶点的最短路径算法,解决的是有权图中最短路径问题。迪杰斯特拉算法主要特点是从起始点开始,采用贪心算法的策略,每次遍历到始点距离最近且未访问过的顶点的邻接节点,直到扩展到终点为止。
DP策略:DP算法的思想与分治算法有点类似, 也是将待求解的问题分解为若干个子问题(称为阶段), 按顺序求解子阶段, 前一子问题的解为后一子问题的求解提供了有用的信息。在求解任一子问题时, 列出各种可能的局部解, 通过决策保留那些有可能达到最优的局部解, 丢弃其他局部解。依次解决各子问题, 最后一个子问题就是初始问题的解。
思路
早知道先做这道题,很经典。此题借鉴,写在这是多打了注释方便自己这个渣渣看懂。(唉,学渣就是学渣,大佬的看不懂啊)
代码
#include<bits/stdc++.h>
using namespace std;
const int N=510,INF=1e8;
int n,m,s,d,city[N],team[N],xx[N][N];
int dist[N],visit[N],pre[N],r[N];
void DIJ(int start,int end){
memset(visit,0,sizeof(visit));
memset(pre,-1,sizeof(pre));
for(int i=0;i<=n;i++){
dist[i]=xx[start][i];
team[i]=city[i];
if(dist[i]<INF){
team[i]=team[start]+city[i];
pre[i]=start;
r[i]=1;
}
}
visit[start]=1;
pre[start]=-1;
team[start]=city[start];
for(int i=1;i<n;i++){
int minpos=n;
for(int j=0;j<n;j++){
if(!visit[j]&&dist[j]<dist[minpos]) minpos=j;
}
if(minpos==n) break;
visit[minpos]=1;
for(int j=0;j<n;j++){
if(!visit[j]){
if(dist[j]>dist[minpos]+xx[minpos][j]){
dist[j]=dist[minpos]+xx[minpos][j];
pre[j]=minpos;
team[j]=team[minpos]+city[j];
r[j]=r[minpos];
}
else if(dist[j]==dist[minpos]+xx[minpos][j]){
r[j]+=r[minpos];
if(team[j]<team[minpos]+city[j]){
team[j]=team[minpos]+city[j];
pre[j]=minpos;
}
}
}
}
}
}
int main(){
for(int i=0;i<=N;i++){
for(int j=0;j<=N;j++) xx[i][j]=INF;
}
cin>>n>>m>>s>>d;
for(int i=0;i<n;i++) cin>>city[i];
for(int i=0;i<m;i++){
int a,b,c;
cin>>a>>b>>c;
xx[a][b]=xx[b][a]=c;
}
DIJ(s,d);
cout<<r[d]<<" "<<team[d]<<endl;
vector<int> v;
for(int i=d;i!=-1;i=pre[i]){
v.push_back(i);
}
for(int i=v.size()-1;i>=0;i--){
cout<<v[i];
if(i) cout<<" ";
}
return 0;
}
学了数据结构以后再研究,好像就没那么难了
上面写的太乱了,当时不会做,看得别人的博客
下面代码是按照课本上写的。。。。。
第37行总是很容易敲错
迪杰斯特拉算法
#include<bits/stdc++.h>
#define x first
#define y second
#define max1 1e8+5
using namespace std;
const int N=505,INF=1e8;
int n,m,s,d;
int num[N],team[N],dist[N],pre[N];
// 路径条数,救援队,
int xx[N][N],city[N],visit[N];
void DIJ(){
for(int i=0;i<=n;i++){
dist[i]=xx[s][i];
team[i]=city[s]+city[i];
num[i]=1;
pre[i]=s;
}
visit[s]=1;
pre[s]=-1;
for(int i=0;i<n-1;i++){
int min1=n;
for(int j=0;j<n;j++){
if(!visit[j]){
if(dist[min1]>dist[j]) min1=j;
}
}
visit[min1]=1;
for(int j=0;j<n;j++){
if(!visit[j]){
if(dist[j]>dist[min1]+xx[min1][j]){
pre[j]=min1;
dist[j]=dist[min1]+xx[min1][j];
team[j]=team[min1]+city[j];
num[j]=num[min1];
}
else if(dist[j]==dist[min1]+xx[min1][j]){
num[j]+=num[min1];
if(team[min1]+city[j]>team[j]){
pre[j]=min1;
team[j]=team[min1]+city[j];
}
}
}
}
}
}
int main(){
for(int i=0;i<N;i++){
for(int j=0;j<N;j++) xx[i][j]=INF;
}
memset(pre,-1,sizeof(pre));
cin>>n>>m>>s>>d;
for(int i=0;i<n;i++) cin>>city[i];
for(int i=0;i<m;i++){
int a,b,c;
cin>>a>>b>>c;
xx[a][b]=xx[b][a]=c;
}
DIJ();
cout<<num[d]<<" "<<team[d]<<endl;
vector<int> v1;
for(int i=d;i!=-1;i=pre[i]){
v1.push_back(i);
}
for(int i=v1.size()-1;i>=0;i--){
cout<<v1[i];
if(i) cout<<" ";
}
return 0;
}