SPFA ---可以解决dijsktar出现负权的

算法流程

算法大致流程是用一个队列来进行维护。初始时将源加入队列。每次从队列中取出一个元素,并对所有与他相邻的点进行松弛,若某个相邻的点松弛成功,则将其入队。直到队列为空时算法结束。

这个算法,简单的说就是队列优化的bellman-ford,利用了每个点不会更新次数太多的特点发明的此算法

SPFA——Shortest Path Faster Algorithm,它可以在O(kE)的时间复杂度内求出源点到其他所有点的最短路径,可以处理负边。SPFA的实现甚至比Dijkstra或者Bellman_Ford还要简单:

设Dist代表S到I点的当前最短距离,Fa代表S到I的当前最短路径中I点之前的一个点的编号。开始时Dist全部为+∞,只有Dist[S]=0,Fa全部为0。

维护一个队列,里面存放所有需要进行迭代的点。初始时队列中只有一个点S。用一个布尔数组记录每个点是否处在队列中。

每次迭代,取出队头的点v,依次枚举从v出发的边v->u,设边的长度为len,判断Dist[v]+len是否小于 Dist[u],若小于则改进Dist[u],将Fa[u]记为v,并且由于S到u的最短距离变小了,有可能u可以改进其它的点,所以若u不在队列中,就将它放入队尾。这样一直迭代下去直到队列变空,也就是S到所有的最短距离都确定下来,结束算法。若一个点入队次数超过n,则有负权环

SPFA 在形式上和宽度优先搜索非常类似,不同的是宽度优先搜索中一个点出了队列就不可能重新进入队列,但是SPFA中一个点可能在出队列之后再次被放入队列,也就是一个点改进过其它的点之后,过了一段时间可能本身被改进,于是再次用来改进其它的点,这样反复迭代下去。设一个点用来作为迭代点对其它点进行改进的平均次数为k,有办法证明对于通常的情况,k在2左右。


由于一般不存在负权环,所以我就不写出现负权的问题了。

本代码采用邻接矩阵的方法储存图,方便理解。

其实我也不知道fa数组是干嘛的。。。。。。。。。

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
#define maxn 1000
#define INF 1000000

using namespace std;

int dist[maxn+2],fa[maxn+2],bool_n[maxn+2],use[maxn+2];
int grap[maxn+2][maxn+2];
int n;
queue<int>q;

void clear()
{
    memset(fa,0,sizeof(fa));
    memset(bool_n,0,sizeof(bool_n));
    memset(dist,INF,sizeof(dist));
    memset(grap,0,sizeof(grap));
    dist[1]=0;
}

void spfa(int x)
{
    int m,i;
    q.push(x);
    use[x]=1;
    while(!q.empty())
    {
        m=q.front();
        bool_n[m]=0;
        q.pop();
        for(i=1;i<=n;i++)
        {
            if(i!=m && grap[m][i])
            {
                if(dist[i]>dist[m]+grap[m][i])
                {
                    dist[i]=dist[m]+grap[m][i];
                    fa[i]=m;
                    if(!bool_n[i])
                        q.push(i),bool_n[i]=1;
                }
            }
        }
    }
}

int main()
{
    int m,x,y,z,i;
    while(~scanf("%d%d",&n,&m))
    {
        clear();
        for(i=0;i<m;i++)
        {
            scanf("%d%d%d",&x,&y,&z);
            grap[x][y]=z;
            grap[y][x]=z;
        }
        spfa(1);
        while(scanf("%d",&x),x)
            printf("%d\n",dist[x]);
    }
    return 0;
}



                
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值