题目描述:
评注: 这道题有一定的思路难度,并且具有一定的实现难度,是一道很好的题,值得细细品味
思路: 首先看懂题面,很容易找到“入手点”是航道,航道是单向边,并且权值可能为负数,并且航道不被包含于图中的环。我们考虑“去边法”,先将图G中所有航道全部去掉,我们来分析此时的图G'会有什么性质:
此时图G'只含有双向边,那么图G'一定由若干个连通块组成(我们将一个连通块称为“团”)。此时我们再将航道加入图G',分析不难得到: 航道必然只能出现在团与团之间,而不能出现在一个团的内部.
Lemma1: 航道不能出现在图G'的团的内部
Proof: 假设一条航道e出现在了团1的内部,考虑e所连接的两个点u,v,不妨假设u->v. 由团的定义,在加入航道e之前,u和v本身就连通,i.e. 存在一条路径l, 从v出发,到达u.
此时,u通过航道到达v,并且v可以通过一条路径到达u,与题目条件矛盾。End of proof.
因此,航道只能出现在团与团之间,进一步分析不难得到以下结论:
Lemma2: 若把团视为一个节点,团与团之间的航道作为有向边,那么图G就是一个有向无环图
Proof: 取一个有向环,简单分析即可,这里不再赘述。
有了以上的引理,不难想到以下思路:
在一个团内,跑堆优化dijkstra算法;在团与团之间,跑拓扑排序算法。时间复杂度O(mlogn)
一个团,可以被拓扑排序遍历,当且仅当这个团的所有节点的入度均为0(否则统计到的答案就是不正确的,容易举反例)
具体实现的时候需要注意很多细节,实现的方法也很多,可以自己实现一下,看看会踩到哪些坑
AC代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
const int N=25010;
const int inf=1e9+10;
struct node{
int v;
int w;
bool friend operator<(node x,node y){
return x.w > y.w;
}
};
vector<node> rout[N];
vector<node> nav[N];
priority_queue<node> q;
queue<int> topo;
int ins[N];
int ins_t[N]; // ins[i] indicate the number of input of the tuan i
int n,r,p,st;
int dis[N];
int cnt=0;
int vis[N];
vector<int> tuan[N];
int tuan_d[N];
void dfs(int x){
vis[x]=1;
ins_t[cnt]+=ins[x];
tuan[cnt].push_back(x);
tuan_d[x]=cnt;
for(int i=0;i<rout[x].size();i++){
int v=rout[x][i].v;
if(vis[v]) continue;
dfs(v);
}
return ;
}
void dijkstra(int s){ // s is the number of tuan
while(!q.empty()) q.pop();
for(int i=0;i<tuan[s].size();i++){
int u=tuan[s][i];
q.push(node{u,dis[u]});
}
while(!q.empty()){
node cur=q.top();
q.pop();
int u=cur.v;
if(vis[u]) continue;
vis[u]=1;
for(int i=0;i<rout[u].size();i++){
int v=rout[u][i].v;
int w=rout[u][i].w;
if(vis[v]) continue;
if(dis[v]>dis[u]+w){
dis[v]=dis[u]+w;
q.push(node{v,dis[v]});
}
}
for(int i=0;i<nav[u].size();i++){
int v=nav[u][i].v;
int w=nav[u][i].w;
if(vis[v]) continue;
if(dis[u]!=inf&&dis[v]>dis[u]+w){
dis[v]=dis[u]+w;
}
ins_t[tuan_d[v]]--;
if(ins_t[tuan_d[v]]==0){
topo.push(tuan_d[v]);
}
}
}
return ;
}
int main(void){
scanf("%d%d%d%d",&n,&r,&p,&st);
for(int i=1;i<=r;i++){
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
rout[x].push_back(node{y,z});
rout[y].push_back(node{x,z});
}
for(int i=1;i<=p;i++){
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
nav[x].push_back(node{y,z});
ins[y]++;
}
memset(vis,0,sizeof(vis));
for(int i=1;i<=n;i++){
if(vis[i]) continue;
++cnt;
dfs(i);
}
for(int i=1;i<=n;i++) dis[i]=inf;
memset(vis,0,sizeof(vis));
dis[st]=0;
topo.push(tuan_d[st]);
for(int i=1;i<=cnt;i++){
if(i==tuan_d[st]) continue;
if(ins_t[i]==0) topo.push(i);
}
while(!topo.empty()){
int u=topo.front();
topo.pop();
dijkstra(u);
}
for(int i=1;i<=n;i++){
if(dis[i]==inf) printf("NO PATH\n");
else printf("%d\n",dis[i]);
}
return 0;
}