【AcWing 342】 道路与航线 SPFA SLF优化

农夫约翰正在一个新的销售区域对他的牛奶销售方案进行调查。

他想把牛奶送到T个城镇，编号为1~T。

这些城镇之间通过R条道路 (编号为1到R) 和P条航线 (编号为1到P) 连接。

每条道路 i 或者航线 i 连接城镇Ai到Bi，花费为Ci。

对于道路，0≤Ci≤10,000;然而航线的花费很神奇，花费Ci可能是负数(−10,000≤Ci≤10,000)。

道路是双向的，可以从Ai到Bi，也可以从Bi到Ai，花费都是Ci。

然而航线与之不同，只可以从Ai到Bi。

事实上，由于最近恐怖主义太嚣张，为了社会和谐，出台了一些政策：保证如果有一条航线可以从Ai到Bi，那么保证不可能通过一些道路和航线从Bi回到Ai。

由于约翰的奶牛世界公认十分给力，他需要运送奶牛到每一个城镇。

他想找到从发送中心城镇S把奶牛送到每个城镇的最便宜的方案。

输入格式
第一行包含四个整数T,R,P,S。

接下来R行，每行包含三个整数（表示一个道路）Ai,Bi,Ci。

接下来P行，每行包含三个整数（表示一条航线）Ai,Bi,Ci。

输出格式
第1…T行：第i行输出从S到达城镇i的最小花费，如果不存在，则输出“NO PATH”。

数据范围
1≤T≤25000,
1≤R,P≤50000,
1≤Ai,Bi,S≤T,

输入样例：
6 3 3 4
1 2 5
3 4 5
5 6 10
3 5 -100
4 6 -100
1 3 -10
输出样例：
NO PATH
NO PATH
5
0
-95
-100

题意：n个点，给r条双向边和p条单向边，判断并输出单源s到其他点的最短路

思路（SPFA + SLF优化）：

很明显的单源最短路问题了，本来想直接SPFA完事，但是这个题专门卡常，所以要用SLF优化，借助双端队列，在原来SPFA算法基础上，d[v]比队头小的进队头，否则进队尾。每次取当前点的时候仍然取队头

AC代码：

#include<iostream>
#include<string>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<map>
#include <queue>
#include<sstream>
#include <stack>
#include <set>
#include<vector>
#define FAST ios::sync_with_stdio(false)
#define abs(a) ((a)>=0?(a):-(a))
#define sz(x) ((int)(x).size())
#define all(x) (x).begin(),(x).end()
#define mem(a,b) memset(a,b,sizeof(a))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define rep(i,a,n) for(int i=a;i<=n;++i)
#define per(i,n,a) for(int i=n;i>=a;--i)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
using namespace std;
typedef pair<int,int> PII;
typedef long long ll;
const int maxn = 1e6+5;
const int inf=0x3f3f3f3f;
const double eps = 1e-7;
const double pi=acos(-1.0);
const int mod = 1e9+7;
inline int lowbit(int x){return x&(-x);}
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
void ex_gcd(ll a,ll b,ll &d,ll &x,ll &y){if(!b){d=a,x=1,y=0;}else{ex_gcd(b,a%b,d,y,x);y-=x*(a/b);}}//x=(x%(b/d)+(b/d))%(b/d);
inline ll qpow(ll a,ll b,ll MOD=mod){ll res=1;a%=MOD;while(b>0){if(b&1)res=res*a%MOD;a=a*a%MOD;b>>=1;}return res;}
inline ll inv(ll x,ll p){return qpow(x,p-2,p);}
inline ll Jos(ll n,ll k,ll s=1){ll res=0;rep(i,1,n+1) res=(res+k)%i;return (res+s)%n;}
inline ll read(){ ll f = 1; ll x = 0;char ch = getchar();while(ch>'9'||ch<'0') {if(ch=='-') f=-1; ch = getchar();}while(ch>='0'&&ch<='9') x = (x<<3) + (x<<1) + ch - '0',  ch = getchar();return x*f; }
int dir[4][2] = { {1,0}, {-1,0},{0,1},{0,-1} };
const int V = 1e5+5, E=5e5+5;
ll d[V],cost[E];
ll head[V],pnt[E],nxt[E],e=0;
ll vis[V];
ll n,r,p,s;
void addedge(ll u,ll v,ll c)
{
    pnt[e]=v;       //当前以u为顶点，c为边长的，到v的一条边
    cost[e]=c;      //存入当前边权值
    nxt[e]=head[u];     //下一个其实是前一个
    head[u]=e++;        //当前边编号
}

void SPFA()
{
    deque<ll > q;
    d[s] = 0;
    q.push_back(s);
    vis[s] = 1;
    while(!q.empty())
    {
        ll x = q.front();
        q.pop_front();
        vis[x] = 0;
        for(int i=head[x]; i!=-1; i = nxt[i])
        {
            ll v = pnt[i];
            if(d[v] > d[x] + cost[i] )
            {
                d[v] = d[x] + cost[i] ;
                if(vis[v]) continue;
                vis[v] = 1;
                if(q.size()&&d[v]<d[q.front()])     //SLF，小的进队头，大的到队尾
                q.push_front(v);
                else
                q.push_back(v);
            }
        }
    }
}


int main()
{
   // freopen("DATA.txt","r",stdin);
        FAST;
       // cin>>n>>r>>p>>s;
        n = read();
        r = read();p = read(); s = read();
        mem(vis,0);
        mem(head,-1);
        mem(cost,0);
        rep(i,1,n)
        d[i] = inf;
        rep(i,1,r)
        {
            ll x, y,z;
            x = read(); y = read(); z = read();
            addedge(x,y,z);
            addedge(y,x,z);
        }
        rep(i,1,p)
        {
            ll x,y,z;
            x = read();
            y = read();
             z = read();
            addedge(x,y,z);
        }
        SPFA();
        rep(i,1,n)
        if(d[i]==inf) cout<<"NO PATH"<<endl;
        else cout<<d[i]<<endl;
        return 0;
}