#include<iostream>
#include<bits/stdc++.h>
const int N=505;
#define inf 0x3f3f3f3f
using namespace std;
int mp[N][N];
int maxx[N]; // 记录救援数量
int num[N]; // 记录同一个路径的个数
int dis[N]; //最短路径
int visit[N]; // 记录是否被访问过
int pre[N]; // 记录前驱
int a[N]; //记录每个点的救援队数量
int n,m,s,d;
void dj(int start)
{
// 进行初始化
memset(visit,0,sizeof(visit));
for(int i=0; i<n; i++)
{
if(i==start)
{
dis[i]=0;
maxx[i]=a[i];
}
else
{
dis[i]=mp[start][i]; //初始为起点到他的长度 没有直接通道初始化为inf
maxx[i]=a[i]+a[start]; // 都要加上起点的救援队数量
}
if(i!=start&&mp[start][i]!=inf) //不是起点 且 起点到i点有直接通道
{
pre[i]=start;
}
else
{
pre[i]= -1;
}
}
visit[start]=1;
for(int i=0; i<n; i++) // 开始遍历
{
int t=inf,k;
for(int j=0; j<n; j++) // 找到最小的那个点
{
if(visit[j]==0&&dis[j]<t)
{
t=dis[j];
k=j; // 这个最小点的坐标为k,到他的最小距离为t
}
}
if(t==inf ) break; //再也找不到没有被访问过的点了
visit[k]=1;
for(int j=0;j<n;j++)
{
if(!visit[j])
{
if(dis[j]>dis[k]+mp[j][k])
{
dis[j]=dis[k]+mp[j][k];
maxx[j]=maxx[k]+a[j];
pre[j]=k;
num[j]=num[k];
}
else if(dis[j]==dis[k]+mp[k][j])
{
num[j]+=num[k];
if(maxx[j]<maxx[k]+a[j]) //第二个判断条件,当距离相同时,找到最多的救援队
{
maxx[j]=maxx[k]+a[j];
pre[j]=k;
}
}
}
}
}
}
void print_path(){
stack<int >x;
// cout <<"打印路径"<<endl;
for(int i=d;i!=-1;i=pre[i]){
// cout <<i<<endl;
x.push(i);
}
// cout <<"完全进栈"<<endl;
// cout <<x.size()<<endl;
int flag=0;
while(!x.empty())
{
if(flag) cout <<" ";
else flag=1;
cout <<x.top();
x.pop();
}
}
int main()
{
cin >>n>>m>>s>>d;
for(int i=0; i<n; i++)
{
cin >>a[i];
num[i]=1; //每一个点的最短路径初始化为1
}
memset(mp,inf,sizeof(mp));
for(int i=0; i<m; i++)
{
int x,y,len;
cin >>x>>y>>len;
mp[x][y]=min(mp[x][y],len);
mp[y][x]=min(mp[y][x],len);
}
dj(s);
cout <<num[d]<<" "<<maxx[d]<<endl;
// cout <<pre[s];
print_path();
return 0;
}
紧急救援(迪杰斯特拉模板)
最新推荐文章于 2024-08-24 09:12:23 发布