道路与航线 AcWing.342

题目描述:

评注: 这道题有一定的思路难度，并且具有一定的实现难度，是一道很好的题，值得细细品味

思路: 首先看懂题面，很容易找到“入手点”是航道，航道是单向边，并且权值可能为负数，并且航道不被包含于图中的环。我们考虑“去边法”，先将图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;
}