请柬（双向SPFA及SLF LLL优化法模板题）

请柬

（双向SPFA模板题）

（题目来源：Luogu.org）
P1342 请柬
难度 普及+/提高
时空限制 1s / 128MB

题目描述

在电视时代,没有多少人观看戏剧表演。Malidinesia古董喜剧演员意识到这一事实，他们想宣传剧院,尤其是古色古香的喜剧片。他们已经打印请帖和所有必要的信息和计划。许多学生被雇来分发这些请柬。每个学生志愿者被指定一个确切的公共汽车站，他或她将留在那里一整天,邀请人们参与。

这里的公交系统是非常特殊的:所有的线路都是单向的，连接两个站点。公共汽车离开起始点，到达目的地之后又空车返回起始点。学生每天早上从总部出发，乘公交车到一个预定的站点邀请乘客。每个站点都被安排了一名学生。在一天结束的时候,所有的学生都回到总部。现在需要知道的是，学生所需的公交费用的总和最小是多少。
输入输出格式
输入格式：

第1行有两个整数n、m(1<=n,m<=1000000)，n是站点的个数，m是线路的个数。

然后有m行，每行描述一个线路，包括3个整数，起始点，目的地和价格。

总部在第1个站点，价钱都是整数，且小于1000000000。

输出格式：

输出一行，表示最小费用。

输入输出样例
输入样例#1：

4 6
1 2 10
2 1 60
1 3 20
3 4 10
2 4 5
4 1 50

输出样例#1：

210

说明

【注意】

此题数据规模较大，需要使用较为高效的算法，此题不设小规模数据分数。

思路

题目分析

这道题需要求两种最短路，一个是从1点到其余各点的最短路，单源多目标的最短路我们选用SPFA，仍然是与链表结合在一起来讲。第二种最短路就是各点到1点的最短路，也就是从1点到其余各点的逆向最短路。也就是下面我们要看到的双向SPFA。

双向SPFA讲解

顾名思义，双向SPFA就是向两个方向进行SPFA，所以我们很自然的可以想到开另外一个结构体数组，把输入数据反着存，这里要记录上一个存的是哪个点，还有r()是读入优化。

数据存储：

for(i=1;i<=m;i++)
    {
        st=r(),en=r(),w=r();
        a[i].y=en;
        a[i].v=w;
        a[i].next=&a[head[st]];
        head[st]=i;
        e[i].v=w;
        e[i].y=st;
        e[i].next=&e[hed2[en]];
        hed2[en]=i;
    }

与正常的SPFA相同，我们只不过就是把SPFA的操作对象由一个变为两个。

SPFA优化法–SLF与LLL

接下来看看SLF和LLL，这里我们重点看SLF，因为LLL不太好理解，（蒟蒻说：）背下来就行了。

SLF

SLF优化法就是把SPFA原来的函数中加了一个判断，就是如果修改后从1到y的距离len[y]>len[q.front()]，就把它推到用来存点的队列q的末尾，反之存在开始处。
要达到这个操作目的，需要用到双向队列deque（不会的查一下C++ Reference Book），在这里主要讲讲它的几种操作：
与queue相似，如果你想在队列前面插入元素，就是q.push_ front()，末尾就是q.push_ back()，弹出元素pop同理。
这种优化法的原因就是尽量减少更新判断成功后的更新步骤，也就是说要尽量让短路靠前，那么我们就把长路放后面，短路放前面（蒟蒻的解释）。
这样的速度会增加50%。

LLL

LLL优化法也是加了一个判断len[yy]*tot<=sum，不大理解，打算背下来（有dalao可以讲解的话欢迎回复赐教）。

结果

因为这个都是大数据嘛，所以dijikstra可能只能过一个点，然后这个题是专门用来卡朴素SPFA的，所以上述优化法很有必要。
带读入优化代码长136 Lines。
还有两位dalao张逸凡和李子晗只用了SLF，
大佬zyf的成绩：
1 AC 844ms/51703kB
2 AC 793ms/51703kB
3 AC 516ms/29253kB
4 AC 27ms/9597kB
只用了SLF都比我快，ORZ。

大佬lzh的成绩
1 AC 970ms/116492kB
2 AC 937ms/88281kB
3 AC 982ms/116492kB
4 AC 12ms/56683kB
这就显得有点卡时间了，还是ORZ。

附上蒟蒻的成绩：

代码

（C++）

#include<iostream>
#include<deque>
#include<cstring>
using namespace std;
int i,j,m,n,st,en,w,tot,sum;

struct node
{
    int y;
    int v;
    struct node *next;
};
struct node a[1000005];
struct node e[1000005];
int head[1000005];
int hed2[1000005];
long long len[1000005],ans;
bool b[1000005];
deque<int> q;
int r()
{
    int ans=0,f=0;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch='-')
        f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        ans*=10;
        ans+=ch-'0';
        ch=getchar();
    }
    return ans;
}

int spfa1(int x)
{
    tot=1;sum=0;
    memset(len,0x7f7f7f,sizeof(len));
    len[x]=0;
    q.push_front(x);
    b[x]=1;
    struct node *p;
    while(!q.empty())
    {
    p=&a[head[q.front()]];
    x=q.front();
    q.pop_front();
    b[x]=0;
    tot--;
    sum-=len[p->y];
    while(p!=NULL)
    {
        int yy=p->y;
        if(len[yy]>p->v+len[x])
        {
            len[yy]=p->v+len[x];
            if(!b[yy])
            {
                b[yy]=1;
                if(q.empty()||len[yy]*tot<=sum||len[yy]>len[q.front()])
                q.push_back(yy);
                else
                q.push_front(yy);
                tot++;sum+=len[yy];
            }
        }
        p=p->next;
    }} 
}

int spfa2(int x)
{
    tot=1;sum=0;
    memset(len,0x7f7f7f,sizeof(len));
    len[x]=0;
    q.push_front(x);
    b[x]=1;
    struct node *p;
    while(!q.empty())
    {
    p=&e[hed2[q.front()]];
    x=q.front();
    q.pop_front();
    b[x]=0;
    tot--;
    sum-=len[p->y];
    while(p!=NULL)
    {
        int yy=p->y;
        if(len[yy]>p->v+len[x])
        {
            len[yy]=p->v+len[x];
            if(!b[yy])
            {
                b[yy]=1;
                if(q.empty()||len[yy]*tot<=sum||len[yy]>len[q.front()])
                q.push_back(yy);
                else
                q.push_front(yy);
                tot++;sum+=len[yy];
            }
        }
        p=p->next;
    }} 
}


int main()
{
    n=r(),m=r();
    for(i=1;i<=m;i++)
    {
        st=r(),en=r(),w=r();
        a[i].y=en;
        a[i].v=w;
        a[i].next=&a[head[st]];
        head[st]=i;
        e[i].v=w;
        e[i].y=st;
        e[i].next=&e[hed2[en]];
        hed2[en]=i;
    }
    spfa1(1);
    for(i=2;i<=n;i++)
    ans+=len[i];
    spfa2(1);
    for(i=2;i<=n;i++)
    ans+=len[i];

    cout<<ans;

}