SPFA
今天看看SPFA。
至于SPFA的来历,在此我不再赘述,详情查询百度百科。
这一次我们借一个单源最短路的题来看看SPFA,这次只讲朴素SPFA,不讲SLF与LLL:
(摘自Luogu模板)
题目描述
如题,给出一个有向图,请输出从某一点出发到所有点的最短路径长度。
输入输出格式
输入格式:
第一行包含三个整数N、M、S,分别表示点的个数、有向边的个数、出发点的编号。
接下来M行每行包含三个整数Fi、Gi、Wi,分别表示第i条有向边的出发点、目标点和长度。
输出格式:
一行,包含N个用空格分隔的整数,其中第i个整数表示从点S出发到点i的最短路径长度(若S=i则最短路径长度为0,若从点S无法到达点i,则最短路径长度为2147483647)
输入输出样例
输入样例#1:
4 6 1
1 2 2
2 3 2
2 4 1
1 3 5
3 4 3
1 4 4
输出样例#1:
0 2 4 3
说明
时空限制:1000ms,128M
数据规模:
对于20%的数据:N<=5,M<=15
对于40%的数据:N<=100,M<=10000
对于70%的数据:N<=1000,M<=100000
对于100%的数据:N<=10000,M<=500000
样例说明:
面对这样的单源最短路问题,我们可以选择SPFA,下面讲讲怎么操作。
初始化
做过图论题的应该都会处理,这里我们用链表。
n=r(),m=r(),s=r();
for(i=1;i<=m;i++)
{
xa=r(),xb=r(),xc=r();
a[i].nxt=&a[hd[xa]];
a[i].v=xc;
a[i].y=xb;
hd[xa]=i;
}
for(i=1;i<=n;i++)
mul[i]=2147483647;
mul[s]=0;
解释一下,r()是读入优化,mul[]是一个记录开始点到该点距离的一个数组,我们要求最短路,所以当然就要把它设置为MAX(注意题中已给MAX取值),并把开始点的mul设为0。
存储节点
SPFA这里需要存储已经修改数据后的节点,用Queue或Deque就可以,二者区别是后者可以在前后插入或弹出(可供SPFA中的SLF优化),而前者只能是前面出后面进。
和其他的最短路求法一样,我们只不过是把更新后的节点存入队列中。
现在我来模拟一下SPFA的过程,来更好的理解这件事情。
每一步里mul[]的变化如下:
这里一共有4步。
它们对应的队列如下图:
这里,如果被更新的节点已经在队列中存在,那么我们就不需要再push了,结束条件就是当Q非空。
代码
C++
下面的代码有几点缺陷:
1.忘记检查需要更新的节点是否已经在队列中了。这里需要一个数组来确定节点是否已经入队,并在pop与push时及时更新,虽然不会影响最终结果的正确性,但还是会拖慢速度。
2.这个代码里用到了SLF,其实就是一个判断来尽量减少操作的次数,以判断的方式防止更新次数更多。
if(mul[p->y]<mul[x])//SLF
q.push_front(p->y);
else
q.push_back(p->y);
接下来就是这个模板题的代码。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<deque>
#include<cstring>
using namespace std;
int i,j,m,n,t,s,e;
int xa,xb;
int mul[10005];
int hd[10005];
double xc;
deque<int> q;
struct data
{
int y;
int v;
struct data *nxt;
}a[500005];
int r()
{
int ans=0,f=1;
char ch=getchar();
while(ch>'9'||ch<'0')
{
if(ch=='-')
f=-1;
ch=getchar();
}
while(ch<='9'&&ch>='0')
{
ans*=10;
ans+=ch-'0';
ch=getchar();
}
return ans;
}
void spfa()
{
struct data *p;
while(!q.empty())
{
int x=q.front();
p=&a[hd[x]];
q.pop_front();
while(p!=NULL)
{
if(mul[p->y]>mul[x]+(p->v))
{
mul[p->y]=mul[x]+(p->v);
if(mul[p->y]<mul[x])
q.push_front(p->y);
else
q.push_back(p->y);
}
p=p->nxt;
}}
}
int main()
{
n=r(),m=r(),s=r();
for(i=1;i<=m;i++)
{
xa=r(),xb=r(),xc=r();
a[i].nxt=&a[hd[xa]];
a[i].v=xc;
a[i].y=xb;
hd[xa]=i;
}
for(i=1;i<=n;i++)
mul[i]=2147483647;
q.push_front(s);
mul[s]=0;
spfa();
for(i=1;i<=n;i++)
cout<<mul[i]<<" ";
}