道路与航线acwing中级最短路

最新推荐文章于 2024-04-16 14:10:34 发布

mrgaohaihang

最新推荐文章于 2024-04-16 14:10:34 发布

阅读量161

点赞数

分类专栏：基础算法

本文链接：https://blog.csdn.net/mrgaohaihang/article/details/113061522

版权

基础算法专栏收录该内容

77 篇文章 0 订阅

订阅专栏

农夫约翰面临牛奶配送挑战，需要在多个城镇间寻找最经济的路线。文章介绍了一个图论问题，涉及道路和航线网络，其中航线费用可能为负，并且不允许回溯。解决方案利用了拓扑排序和优先队列，寻找从中心城镇出发到每个城镇的最小花费。程序实现了这一算法，输出每个城镇的最低配送成本，若无法到达则标记为'NOPATH'。

摘要由CSDN通过智能技术生成

农夫约翰正在一个新的销售区域对他的牛奶销售方案进行调查。

他想把牛奶送到T个城镇，编号为1~T。

这些城镇之间通过R条道路 (编号为1到R) 和P条航线 (编号为1到P) 连接。

每条道路 i
或者航线 i 连接城镇Ai到Bi，花费为Ci

。

对于道路，0≤Ci≤10,000
;然而航线的花费很神奇，花费Ci可能是负数(−10,000≤Ci≤10,000

)。

道路是双向的，可以从Ai
到Bi，也可以从Bi到Ai，花费都是Ci

。

然而航线与之不同，只可以从Ai
到Bi

。

事实上，由于最近恐怖主义太嚣张，为了社会和谐，出台了一些政策：保证如果有一条航线可以从Ai
到Bi，那么保证不可能通过一些道路和航线从Bi回到Ai

。

由于约翰的奶牛世界公认十分给力，他需要运送奶牛到每一个城镇。

他想找到从发送中心城镇S把奶牛送到每个城镇的最便宜的方案。
输入格式

第一行包含四个整数T,R,P,S。

接下来R行，每行包含三个整数（表示一个道路）Ai,Bi,Ci

。

接下来P行，每行包含三个整数（表示一条航线）Ai,Bi,Ci

。
输出格式

第1…T行：第i行输出从S到达城镇i的最小花费，如果不存在，则输出“NO PATH”。
数据范围

1≤T≤25000
,
1≤R,P≤50000,
1≤Ai,Bi,S≤T

,
输入样例：

6 3 3 4
1 2 5
3 4 5
5 6 10
3 5 -100
4 6 -100
1 3 -10

输出样例：

NO PATH
NO PATH
5
0
-95
-100


#include <iostream>
#include<queue>
#include<algorithm>
#pragma warning(disable:4996)
using namespace std;
int n, mt, nt, S;
const int N = 1e5 + 10;
int ne[N], h[N], e[N], d[N], idex;
bool sen[N];
int dist[N];
queue<int>q;
typedef pair<int, int>pii;
priority_queue<pii, vector<pii>, greater<pii> >qu;
int cntb = 1;
vector<int>block[N];
int cntblock[N];
int id[N];
void add(int a, int b, int c)
{
	d[idex] = b;
	e[idex] = c;
	ne[idex] = h[a];
	h[a] = idex++;
}
void dfs(int u, int cnt)
{
	id[u] = cnt;
	block[cnt].push_back(u);
	for (int i = h[u]; i != -1; i = ne[i])
	{
		int j = d[i];
		if (!id[j])
		{
			dfs(j, cnt);
		}
	}
}
void di(int u)
{
	for (int i = 0; i < block[u].size(); i++)
	{
		qu.push({ dist[block[u][i]],block[u][i] });
	}
	while (qu.size())
	{
		auto t = qu.top();
		qu.pop();
		int var = t.second, distance = t.first;
		if (sen[var])continue;
		sen[var] = true;
		for (int i = h[var]; i != -1; i = ne[i])
		{
			int j = d[i];
			if (id[j] != id[var] && --cntblock[id[j]] == 0)
			{
				q.push(id[j]);

			}
			if (dist[j] > dist[var] + e[i])
			{
				dist[j] = dist[var] + e[i];
				if (id[j] == id[var])
				{
					qu.push({ dist[j], j });
				}
			}
			
		}
	}
}

void topsort()
{
	memset(dist, 0x3f3f3f3f, sizeof(dist));
	dist[S] = 0;
	for (int i = 1; i <= cntb; i++)//按照块来分
	{
		if (!cntblock[i])
		{
			q.push(i);
		}
	}
	while (q.size())
	{
		int t = q.front();
		q.pop();
		di(t);
	}
}

int main()
{
	cin >> n >> mt >> nt >> S;
	memset(h, -1, sizeof(h));
	for (int i = 1; i <= mt; i++)
	{
		int a, b, c;
		scanf("%d%d%d", &a, &b, &c);
		add(a, b, c);
		add(b, a, c);
	}
	for (int i = 1; i <= n; i++)
	{
		if (!id[i])
		{
			dfs(i, cntb);
			cntb++;
		}
	}
	for (int i = 1; i <= nt; i++)
	{
		int a, b, c; scanf("%d%d%d", &a, &b, &c);
		add(a, b, c);
		cntblock[id[b]]++;

	}
	topsort();
	for (int i = 1; i <= n; i++)
	{
		if (dist[i] > 0x3f3f3f3f / 2)
		{
			puts("NO PATH");
		}
		else
		{
			printf("%d\n", dist[i]);
		}
	}
	return 0;

}