SPFA的SLF与LLL优化

    SPFA算法的优化有两个：SLF与LLL，这两个优化都是利用deque进行的。

    1、SLF(Small Label First)：为了方便，把新加入的顶点记为v，队首记为u。每次进队时如果dist[v]<dist[u]则插入队首，否则插入队尾。

    2、LLL(Large Label Last)：把队伍中顶点的dist平均值记为x。如果dist[v]>x则插入队尾，然后查找下一顶点直到dist[v]<=x再进行松弛操作。

    虽然传说 SLF + LLL 可以使速度提高50%，但是根据我自己测试的情况，以及各位巨巨的说法 LLL 不是很给力啊，很多时候加上 LLL 反而没有只用 SLF 快。也许这就是很多人说 LLL 不常用，一般只用 SLF 的原因吧。

    用HDU 4725来实现 SLF 与 LLL：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <vector>
#include <queue>
#include <map>
#include <cmath>
#include <stack>
#include <deque>
using namespace std;
#define INF 0x3f3f3f3f
#define fi first
#define se second
#define mem(a,b) memset((a),(b),sizeof(a))

const int MAXV=300000+3;
const int MAXE=600000+3;

struct Edge
{
    int to,cost,next;
    Edge(int t=0,int c=0,int n=0):to(t),cost(c),next(n){}
}edge[MAXE];

int N,M,C,V,head[MAXV],cnt;
int sum;//队列中顶点dist之和（用于lll优化）
bool vis[MAXV];
int dist[MAXV];
deque<int> que;

void spfa()
{
    sum=0;
    mem(vis,0);
    mem(dist,0x3f);
    vis[1]=true;
    dist[1]=0;
    while(!que.empty())
        que.pop_back();
    que.push_back(1);
    while(!que.empty())
    {
        int u=que.front(); que.pop_front();
        sum-=dist[u];
        if(dist[u]*que.size()>sum)//lll优化
        {
            que.push_back(u);
            sum+=dist[u];
            continue;
        }
        vis[u]=false;
        for(int i=head[u];i!=-1;i=edge[i].next)
        {
            int v=edge[i].to;
            if(dist[v]>dist[u]+edge[i].cost)
            {
                dist[v]=dist[u]+edge[i].cost;
                if(!vis[v])
                {
                    vis[v]=true;
                    if(!que.empty()&&dist[v]<dist[que.front()])//slf优化
                        que.push_front(v);
                    else que.push_back(v);
                    sum+=dist[v];
                }
            }
        }
    }
}

void init()
{
    for(int i=1;i<=V;++i)
        head[i]=-1;
    cnt=0;
}

void add_edge(int from,int to,int cost)
{
    edge[cnt].to=to;
    edge[cnt].cost=cost;
    edge[cnt].next=head[from];
    head[from]=cnt++;
}

int main()
{
    int T_T;
    scanf("%d",&T_T);
    for(int cas=1;cas<=T_T;++cas)
    {
        // 1 ~ N        点
        // N+1 ~ N*2    层（入）
        // N*2+1 ~ N*3  层（出）
        scanf("%d%d%d",&N,&M,&C);
        V=N*3;
        init();
        for(int u=1;u<=N;++u)
        {
            int tmp;
            scanf("%d",&tmp);
            add_edge(u,N+tmp,0);
            add_edge(N*2+tmp,u,0);
        }
        for(int i=1;i<N;++i)
        {
            add_edge(N+i,N*2+i+1,C);
            add_edge(N+i+1,N*2+i,C);
        }
        for(int i=0;i<M;++i)
        {
            int u,v,c;
            scanf("%d%d%d",&u,&v,&c);
            add_edge(u,v,c);
            add_edge(v,u,c);
        }
        spfa();
        printf("Case #%d: %d\n",cas,dist[N]==INF?-1:dist[N]);
    }
    
    return 0;
}