2019ICPC银川区域赛H BZOJ2200 AcWing 342. 道路与航线(连通块+dij+拓扑排序)(SLF优化的SPFA模板)

我以后再也不敢水题解了，就差这个题正确做法银川拿银！！

题干：

农夫约翰正在一个新的销售区域对他的牛奶销售方案进行调查。
他想把牛奶送到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

思路：

题意简化后就是从S点出发，求它到各个城市的最短路，道路表示双向边，航路表示单向边。因为有负边，所以不能用dijkstra算法，又因为不会有负环，所以还是可以跑最短路的。
（1）普通的SPFAtle，我们试下SLF优化的SPFA，居然氵过了。。。
(氵过你?呢（来自2019.10.21的我）)
（2）参考自LYD大佬的思路。
如果我们只看双向边，则会形成一个个的连通块，然后将每个连通块看作一个点，这样再加上那些单向边就会形成一个有向无环图，拓扑排序就可以o(n)求出单源最短路。然后对于每个连通块用dijkstra计算块内最短路。
流程：①输入双向边
②DFS划分连通块
③输入单向边并确定连通块的入度
④将起点和入度为零的连通块入队，准备开始拓扑
⑤得到队首，并将该连通块包含的所有点入堆（根据dis[i]的小根堆），开始dijkstra
⑥根据d[x]+z去更新d[y]
⑦若x和y在一个连通块中，且d[y]被更新，则将y入堆
⑧若x和y不在一个连通块中，则y所在的连通块入度-1，若减到零，则将该连通块插队。
⑨重复⑥~⑧直到堆空
⑩重复⑤~⑨直到队空

正确解法  
连通块+dij+拓扑排序
#include <cstdio>
#include <cmath>
#include <vector>
#include <queue>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
typedef pair<int,int> PII;
struct stu{
    int next;
    int cost;
};
vector<stu> tu[100100];
int dis[25010],vis[25010],cnt,v[25010],edg[25010];
void dfs(int x){
    v[x]=cnt;
    for(int i=0;i<tu[x].size();i++){
        int y=tu[x][i].next;
        if(v[y])
            continue;
        dfs(y);
    }
}
void add(int a,int b,int c){
    stu t;
    t.next=b;
    t.cost=c;
    tu[a].push_back(t);
}
int main()
{
    int r,n,p,a,b,c,s;
    scanf("%d%d%d%d",&n,&r,&p,&s);
    for(int i=0;i<r;i++){
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
        add(b,a,c);
    }
    for(int i=1;i<=n;i++){
        if(!v[i]){
            cnt++;
            dfs(i);
        }
    }
    for(int i=0;i<p;i++){
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
        edg[v[b]]++;
    }
    queue<int>q;
    q.push(v[s]);
    for(int i=1;i<=cnt;i++)
        if(!edg[i])
            q.push(i);
    memset(dis,0x3f,sizeof(dis));
    dis[s]=0;
    memset(vis,0,sizeof(vis));
    while(!q.empty()){
        int t=q.front();
        q.pop();
        priority_queue<PII,vector<PII>,greater<PII> >heap;
        for(int i=1;i<=n;i++){
            if(v[i]==t)
                heap.push({dis[i],i});
        }
        
        while(heap.size()){
            auto tt=heap.top();
            heap.pop();
            int x=tt.second;
            if(vis[x]) continue;
            vis[x]=1;
            for(int i=0;i<tu[x].size();i++){
                int y=tu[x][i].next,flag=0;
                if(dis[y]>dis[x]+tu[x][i].cost){
                    dis[y]=dis[x]+tu[x][i].cost;
                    flag=1;
                }
                if(v[x]==v[y]){
                    if(flag){
                        heap.push({dis[y],y});
                    }
                }
                else{
                    edg[v[y]]--;
                    if(!edg[v[y]]){
                        q.push(v[y]);
                    }
                }
            }
        }
    }
    for(int i=1;i<=n;i++){
        if(dis[i]>=1e9)
            printf("NO PATH\n");
        else
            printf("%d\n",dis[i]);
    }
    return 0;
}

//SLF优化的SPFA
#include <bits/stdc++.h>
using namespace std;
struct stu
{
    int next;
    int cost;
};
vector<stu> tu[100100];
int dis[25010],vis[25010];
void add(int a,int b,int c){
    stu t;
    t.next=b;
    t.cost=c;
    tu[a].push_back(t);
}
void spfa(int s)
{
    memset(vis,0,sizeof(vis));
    memset(dis,0x3f,sizeof(dis));
    deque<int> q;
    vis[s]=1;
    dis[s]=0;
    q.push_back(s);
    while(!q.empty()){
        int x=q.front();
        vis[x]=0;
        q.pop_front();
        for(int i=0;i<tu[x].size();i++){
            int v=tu[x][i].next;
            int c=tu[x][i].cost;
            if(dis[v]>dis[x]+c){
                dis[v]=dis[x]+c;
                if(!vis[v]){
                    vis[v]=1;
                    if(q.size()&&dis[v]<dis[q.front()])
                        q.push_front(v);
                    else
                        q.push_back(v);
                }
            }
        }
    }
}
int main()
{
    int t,r,p,s,a,b,c;
    scanf("%d%d%d%d",&t,&r,&p,&s);
    for(int i=0;i<r;i++){
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
        add(b,a,c);
    }
    for(int i=0;i<p;i++){
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
    }
    spfa(s);
    for(int i=1;i<=t;i++){
        if(dis[i]==0x3f3f3f3f)
            printf("NO PATH\n");
        else
            printf("%d\n",dis[i]);
    }

    return 0;
}