第五周学习

问题

农夫约翰正在针对一个新区域的牛奶配送合同进行研究。他打算分发牛奶到T个城镇（标号为1…T），这些城镇通过R条标号为（1…R）的道路和P条标号为（1…P）的航路相连。

每一条公路i或者航路i表示成连接城镇Ai（1<=A_i<=T）和Bi（1<=Bi<=T）代价为Ci。每一条公路，Ci的范围为0<=Ci<=10,000；由于奇怪的运营策略，每一条航路的Ci可能为负的，也就是-10,000<=Ci<=10,000。

每一条公路都是双向的，正向和反向的花费是一样的，都是非负的。

每一条航路都根据输入的Ai和Bi进行从Ai->Bi的单向通行。实际上，如果现在有一条航路是从Ai到Bi的话，那么意味着肯定没有通行方案从Bi回到Ai。

农夫约翰想把他那优良的牛奶从配送中心送到各个城镇，当然希望代价越小越好，你可以帮助他嘛？配送中心位于城镇S中（1<=S<=T）。

Input

输入的第一行包含四个用空格隔开的整数T，R，P，S。

接下来R行，描述公路信息，每行包含三个整数，分别表示Ai，Bi和Ci。

Output

输出T行，分别表示从城镇S到每个城市的最小花费，如果到不了的话输出NO PATH。

示例

输入

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

代码

#include<iostream>
#include<cstring>
#include<queue>
#include<vector>
using namespace std;

typedef pair<int,int> PII;

const int N = 25010, M = 150010, INF = 0x3f3f3f3f;
int h[N], e[M], ne[M], w[M], idx;
int dist[N];
bool st[N];

int n, mr, mp, s;
//id:点所在的连通块号, bcnt:连通块数量
int id[N], bcnt;
//连通块的入度数组
int deg[N];
vector<int> block[N];

int q[N], hh, tt = -1;

void add(int a, int b, int c)
{
    e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx++;
}

void dfs(int u)
{
    block[bcnt].push_back(u);
    id[u] = bcnt;

    for(int i = h[u]; i != -1; i = ne[i])
    {
        int j = e[i];
        if(!id[j]) dfs(j);
    }
}

void dijkstra(int block_id)
{
    priority_queue<PII, vector<PII>, greater<PII>> heap;

    //取出该连通块所有的点, 入堆
    for(int u : block[block_id])
    {
        heap.push({dist[u], u});
    }

    while(heap.size())
    {
        PII t = heap.top();
        heap.pop();
        int u = t.second;

        if(st[u]) continue;
        st[u] = true;

        for(int i = h[u]; i != -1; i = ne[i])
        {
            int j = e[i];
            if(dist[j] > dist[u] + w[i])
            {
                dist[j] = dist[u] + w[i];
                if(id[j] == block_id) heap.push({dist[j], j});
            }

            if(id[j] != block_id && --deg[id[j]] == 0) q[++tt] = id[j];
        }

    }
}

void topoSort()
{
    //dist赋初值
    memset(dist, 0x3f, sizeof dist);
    dist[s] = 0;

    //入度为0的连通块入队列
    for(int i = 1; i <= bcnt; i++)
    {
        if(!deg[i]) q[++tt] = i;
    }

    while(hh <= tt)
    {
        int t = q[hh++];
        dijkstra(t);
    }
}

int main()
{
    scanf("%d%d%d%d", &n, &mr, &mp, &s);
    memset(h, -1, sizeof h);
    for(int i = 0; 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);
        }
    }
    //在连通块之间建边
    for(int i = 0; i < mp; i++)
    {
        int a, b, c;
        scanf("%d%d%d", &a, &b, &c);
        add(a, b, c), deg[id[b]]++;
    }

    //算法开始
    topoSort();

    for(int i = 1; i <= n; i++)
    {
        if(dist[i] > INF / 2) puts("NO PATH");
        else printf("%d\n", dist[i]);
    }

    return 0;
}