SPFA算法模版+邻接表实现

SPFA即shotest path faster algorithm,由意思就可以看出该算法效率比较高。

其实SPFA就是bellman-ford算法的一个优化。

具体做法是用一个队列保存待松弛的点,然后对于每个出队的点依次遍历每个与他有边相邻的点(用邻接表效率较高),如果该点可以松弛并且队列中没有该点则将它加入队列中,如此迭代直到队列为空。

据说平均效率是O(E),可见对边稀疏的图用此算法效果是相当可观的。

 

若要判负环路,则记录一个点的入队次数,若超过边数,则有负权环。

 

#include  < iostream >
#include 
< queue >
using   namespace  std;

const   long  MAXN = 10000 ;
const   long  lmax = 0x7FFFFFFF ;

typedef 
struct   
ExpandedBlockStart.gifContractedBlock.gif
{
    
long v;
    
long next;
    
long cost;
}
Edge;


Edge e[MAXN];
long  p[MAXN];
long  Dis[MAXN];
bool  vist[MAXN];

queue
< long >  q;

long  m,n; // 点,边
void  init()
ExpandedBlockStart.gifContractedBlock.gif
{
    
long i;
    
long eid=0;

    memset(vist,
0,sizeof(vist));
    memset(p,
-1,sizeof(p));
    fill(Dis,Dis
+MAXN,lmax);

    
while (!q.empty())
ExpandedSubBlockStart.gifContractedSubBlock.gif    
{
        q.pop();
    }


    
for (i=0;i<n;++i)
ExpandedSubBlockStart.gifContractedSubBlock.gif    
{
        
long from,to,cost;
        scanf(
"%ld %ld %ld",&from,&to,&cost);

        e[eid].next
=p[from];
        e[eid].v
=to;
        e[eid].cost
=cost;
        p[from]
=eid++;

        
//以下适用于无向图
        swap(from,to);
        
        e[eid].next
=p[from];
        e[eid].v
=to;
        e[eid].cost
=cost;
        p[from]
=eid++;

    }

}


void  print( long  End)
ExpandedBlockStart.gifContractedBlock.gif
{
    
//若为lmax 则不可达
    printf("%ld\n",Dis[End]);    
}


void  SPF()
ExpandedBlockStart.gifContractedBlock.gif
{

    init();

    
long Start,End;
    scanf(
"%ld %ld",&Start,&End);
    Dis[Start]
=0;
    vist[Start]
=true;
    q.push(Start);

    
while (!q.empty())
ExpandedSubBlockStart.gifContractedSubBlock.gif    
{
        
long t=q.front();
        q.pop();
        vist[t]
=false;
        
long j;
        
for (j=p[t];j!=-1;j=e[j].next)
ExpandedSubBlockStart.gifContractedSubBlock.gif        
{
            
long w=e[j].cost;
            
if (w+Dis[t]<Dis[e[j].v])
ExpandedSubBlockStart.gifContractedSubBlock.gif            
{
                Dis[e[j].v]
=w+Dis[t];
                
if (!vist[e[j].v])
ExpandedSubBlockStart.gifContractedSubBlock.gif                
{
                    vist[e[j].v]
=true;
                    q.push(e[j].v);
                }

            }

        }

    }


    print(End);

}


int  main()
ExpandedBlockStart.gifContractedBlock.gif
{
    
while (scanf("%ld %ld",&m,&n)!=EOF)
ExpandedSubBlockStart.gifContractedSubBlock.gif    
{
        SPF();
    }

    
return 0;
}

转载于:https://www.cnblogs.com/zhuangli/archive/2008/07/26/1251869.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值