SPFA

SPFA

题目链接

Time Limit: 10000ms

Case Time Limit: 1000ms

Memory Limit: 256MB

描述

万圣节的晚上，小Hi和小Ho在吃过晚饭之后，来到了一个巨大的鬼屋！

鬼屋中一共有N个地点，分别编号为1..N，这N个地点之间互相有一些道路连通，两个地点之间可能有多条道路连通，但是并不存在一条两端都是同一个地点的道路。

不过这个鬼屋虽然很大，但是其中的道路并不算多，所以小Hi还是希望能够知道从入口到出口的最短距离是多少？

  

提示：Super Programming Festival Algorithm。

“唔……地点很多，道路很少，这个鬼屋是一个稀疏图，既然这一点被特地标注出来，那么想来有其作用的咯？”小Ho道。

“是的，正好有一种最短路径算法，它的时间复杂度只和边的条数有关，所以特别适合用来解决这种边的数量很少的最短路问题！”小Hi点了点头道：“它就是SPFA算法，即Shortest Path Faster Algorithm。”

“听上去很厉害的样子，但是实际上怎么做的呢？”小Ho问道。

“你会用宽度优先搜索写这道题么？”小Hi反问道。

“这个当然会啊，构造一个队列，最开始队列里只有(S, 0)——表示当前处于点S，从点S到达该点的距离为0，然后每次从队首取出一个节点(i, L)——表示当前处于点i，从点S到达该点的距离为L，接下来遍历所有从这个节点出发的边(i, j, l)——表示i和j之间有一条长度为l的边，将(j, L+l)加入到队尾，最后看所有遍历的(T, X)节点中X的最小值就是答案咯~”小Ho对于搜索已经是熟稔于心，张口便道。

“SPFA算法呢，其实某种意义上就是宽度优先搜索的优化——如果你在尝试将(p, q)加入到队尾的时候，发现队列中已经存在一个(p, q')了，那么你就可以比较q和q'：如果q>=q'，那么(p, q)这个节点实际上是没有继续搜索下去的必要的——算是一种最优化剪枝吧。而如果q<q'，那么(p, q')也是没有必要继续搜索下去的——但是它已经存在于队列里了怎么办呢？很简单，将队列中的(p, q')改成(p, q)就可以了！”

“那我该怎么知道队列中是不是存在一个(p, q')呢？”<额，维护一个position[1..N]的数组就可以了，如果不在队列里就是-1，否则就是所在的位置！”

“所以说这本质上就是宽度优先搜索的剪枝咯？”小Ho问道。

小Hi笑道：“怎么可能！SPFA算法其实是BELLMAN-FORD算法的一种优化版本，只不过在成型之后可以被理解成为宽度优先搜索的！这个问题，我们会在之后好好讲一讲的！”

    

输入

每个测试点（输入文件）有且仅有一组测试数据。

在一组测试数据中：

第1行为4个整数N、M、S、T，分别表示鬼屋中地点的个数和道路的条数，入口（也是一个地点）的编号，出口（同样也是一个地点）的编号。

接下来的M行，每行描述一条道路：其中的第i行为三个整数u_i, v_i, length_i，表明在编号为u_i的地点和编号为v_i的地点之间有一条长度为length_i的道路。

对于100%的数据，满足N<=10^5，M<=10^6, 1 <= length_i <= 10^3, 1 <= S, T <= N, 且S不等于T。

对于100%的数据，满足小Hi和小Ho总是有办法从入口通过地图上标注出来的道路到达出口。

输出

对于每组测试数据，输出一个整数Ans，表示那么小Hi和小Ho为了走出鬼屋至少要走的路程。

AC代码：

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
int n = 100010;  //点
int m = 1000010;  //路径
long long len[1000010];
int vis[100010];
int head[100010];
queue<int>Q;
struct node {
	int v;
	int next;
	int weight;
}point[2000010];


void run(int num,int s,int t)
{
	Q.push(s);
	while (!Q.empty())
	{
		int h=Q.front();
		Q.pop();
		vis[h] = 0;
		for(int i=head[h];i!=-1;i=point[i].next)
		{
			int v = point[i].v;
			if (len[h] + point[i].weight < len[v])
			{
				len[v] = len[h] + point[i].weight;
				if (!vis[v])
				{
					vis[v] = 1;
					Q.push(v);
				}
			}
		}
	}
	cout <<len[t] << endl;
}

int main()
{
	memset(len, 0x37, sizeof(len));
	memset(head, -1, sizeof(head));
	memset(vis, 0, sizeof(vis));
	int s, t; //入口，出口
	while (scanf("%d%d%d%d", &n, &m, &s, &t)!=EOF)
	{
		len[s] = 0;
		int num = 1;
		for (int i = 1; i <= m; i++)
		{
			int a, b, c;
			scanf("%d%d%d", &a, &b, &c);
			point[num].v = b;
			point[num].weight = c;
			point[num].next = head[a];
			head[a] = num++;
 			point[num].v = a;
			point[num].weight = c;
			point[num].next = head[b];
			head[b] = num++;
		}
		run(num,s,t);
	}
	return 0;
}