SPFA入门

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]<<" ";
}