单源最短路——道路与航线

道路与航线

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

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

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

每条道路 i 或者航线 i 连接城镇$A_i$到$B_i$，花费为$C_i$。

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

道路是双向的，可以从$A_i$到$B_i$，也可以从Bi到$A_i$，花费都是$C_i$。

然而航线与之不同，只可以从$A_i$到$B_i$。

事实上，由于最近恐怖主义太嚣张，为了社会和谐，出台了一些政策：保证如果有一条航线可以从$A_i$到$B_i$，那么保证不可能通过一些道路和航线从$B_i$回到$A_i$。
由于约翰的奶牛世界公认十分给力，他需要运送奶牛到每一个城镇。
他想找到从发送中心城镇S把奶牛送到每个城镇的最便宜的方案。
输入格式
第一行包含四个整数T,R,P,S。
接下来R行，每行包含三个整数（表示一个道路）$A_i,B_i,C_i$。
接下来P行，每行包含三个整数（表示一条航线）$A_i,B_i,C_i$。
输出格式
第1…T行：第i行输出从S到达城镇i的最小花费，如果不存在，则输出“NO PATH”。

数据范围
$1\le T\le 25000$,
$1\le R,P\le 50000$,
$1\le A_i,B_i,S\le 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

题意：

说这张图上面有单向边也有双向边，道路中是双向图并且保证边权为正，但是航线是单相边并且边权可正可负。

题解：

这道题蛮好的，对自己的图论知识提升挺大的，开始看到有正有负直接就去跑spfa了，但是直接就被卡掉了。然后y老师说过如果我们在一张无环的拓扑图中我们直接跑一边拓扑序是可以找出最短路的。所以我们就可以针对航线跑拓扑序，道路跑地杰斯特拉求最短路。跑拓扑序是线性的，地杰斯特拉是mlog所以最后跑过去了。

#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#define x first
#define y second
using namespace std;
typedef pair<int, int> PII;
const int N = 25010, M = 150010, inf = 0x3f3f3f3f;
int n, mr, mp, S;
int id[N];
int head[N], e[M], w[M], ne[M], idx;
int dist[N], din[N];
vector<int> block[N];
int bcnt;
bool vis[N];
queue<int> q;
void add(int a, int b, int c)
{
    e[idx] = b, w[idx] = c, ne[idx] = head[a], head[a] = idx ++ ;
}
void dij(int x)
{
    priority_queue<PII,vector<PII>,greater<PII> >p;
    for(auto i:block[x]){
        p.push({dist[i],i});
    }
    while(p.size()){
        auto u=p.top(); p.pop();
        if(vis[u.y]) continue;
        vis[u.y]=1;
        for(int i=head[u.y];i!=-1;i=ne[i]){
            int wi=w[i];
            int j=e[i];
            if(id[j]!=id[u.y]&&--din[id[j]]==0) q.push(id[j]);
            if(dist[j]>dist[u.y]+wi){
                dist[j]=dist[u.y]+wi;
                if(id[j]==id[u.y]) p.push({dist[j],j});
            }
        }
    }
}
void top_sort()
{
    memset(dist,inf,sizeof dist);
    memset(vis,0,sizeof vis);
    dist[S]=0;
    for(int i=1;i<=bcnt;i++){
        if(!din[i]) q.push(i);
    }
    while(q.size()){
        int u=q.front(); q.pop();
        dij(u);
    }
}
void dfs(int x,int cnt)
{
    id[x]=cnt; block[cnt].push_back(x);
    for(int i=head[x];i!=-1;i=ne[i]){
        int j=e[i];
        if(id[j]) continue;
        dfs(j,cnt);
    }
}
int main()
{
    cin>>n>>mr>>mp>>S;
    memset(head,-1,sizeof head);
    for(int i=1;i<=mr;i++){
        int a,b,c; scanf("%d%d%d",&a,&b,&c);
        add(a,b,c),add(b,a,c);
    }
    for(int i=1;i<=n;i++){
        if(!id[i]) bcnt++,dfs(i,bcnt);
    }
    for(int i=1;i<=mp;i++){
        int a,b,c; scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
        din[id[b]]++;
    }
    top_sort();
    for (int i = 1; i <= n; i ++ )
        if (dist[i] > inf / 2) cout << "NO PATH" << endl;
        else cout << dist[i] << endl;

    return 0;
}