BZOJ 2200 [Usaco2011 Jan]道路和航线

Description

Farmer John正在一个新的销售区域对他的牛奶销售方案进行调查。他想把牛奶送到T个城镇 (1 <= T <= 25,000)，编号为1T。这些城镇之间通过R条道路 (1 <= R <= 50,000，编号为1到R) 和P条航线 (1 <= P <= 50,000，编号为1到P) 连接。每条道路i或者航线i连接城镇A_i (1 <= A_i <= T)到B_i (1 <= B_i <= T)，花费为C_i。对于道路，0 <= C_i <= 10,000;然而航线的花费很神奇，花费C_i可能是负数(-10,000 <= C_i <= 10,000)。道路是双向的，可以从A_i到B_i，也可以从B_i到A_i，花费都是C_i。然而航线与之不同，只可以从A_i到B_i。事实上，由于最近恐怖主义太嚣张，为了社会和谐，出台 了一些政策保证：如果有一条航线可以从A_i到B_i，那么保证不可能通过一些道路和航线从B_i回到A_i。由于FJ的奶牛世界公认十分给力，他需要运送奶牛到每一个城镇。他想找到从发送中心城镇S(1 <= S <= T) 把奶牛送到每个城镇的最便宜的方案，或者知道这是不可能的。

Input

第1行：四个空格隔开的整数: T, R, P, and S 第2到R+1行：三个空格隔开的整数（表示一条道路）：A_i, B_i 和 C_i 第R+2到R+P+1行：三个空格隔开的整数（表示一条航线）：A_i, B_i 和 C_i

Output

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

Sample Input

 6 3 3 4
 1 2 5
 3 4 5
 5 6 10
 3 5 -100
 4 6 -100
 1 3 -10

样例输入解释：

一共六个城镇。在1-2，3-4，5-6之间有道路，花费分别是5，5，10。同时有三条航线：3->5，

4->6和1->3，花费分别是-100，-100，-10。FJ的中心城镇在城镇4。

Sample Output

 NO PATH
 NO PATH
 5
 0
 -95
 -100

样例输出解释：

FJ的奶牛从4号城镇开始，可以通过道路到达3号城镇。然后他们会通过航线达到5和6号城镇。
但是不可能到达1和2号城镇。

---

联通块+拓扑排序+Dijkstra

先把道路加进图中，划分联通块，把每个联通块看做一个整体，这样整个图是DAG，可以统计入度进行拓扑排序，个人认为按照拓扑序进行扫描也就是在没有环时扫描所有的边，所以不受边权正负影响

计算时联通块内部用dijkstra求最短路，如果道路通到另一个联通块就将另一个块的编号加进拓扑序中

#include <queue>
#include <cstdio>
#include <cstring>
using namespace std;

using paii = pair<int, int>;
const int maxn = 25005;
const int maxm = 200005;
int n, m, p, st, tot, tun, cnt;    	   //点,道路,航路,起点,边邻接表,联通块邻接表,联通块
int head1[maxn], Next1[maxm], ver1[maxm], edge[maxm]; //边邻接表
int head2[maxn], Next2[maxn], ver2[maxn];             //联通块邻接表
int c[maxn], vis[maxn], deg[maxn], d[maxn];           //联通块编号,访问数组,入度,距离

void add(int x, int y, int z) //边
{
    Next1[++tot] = head1[x], head1[x] = tot;
    ver1[tot] = y, edge[tot] = z;
}

void addun(int x, int y) //编号为cnt的点的邻接表
{
    Next2[++tun] = head2[x], head2[x] = tun;
    ver2[tun] = y;
}

void dfs(int x)
{
    c[x] = cnt;
    addun(cnt, x); //编号为cnt的点的邻接表
    for (int i = head1[x]; i; i = Next1[i])
        if (!c[ver1[i]])
            dfs(ver1[i]);
}

queue<int> top; //拓扑队列
priority_queue<paii> q;//dikstra堆

void dijkstra(int s)
{
    for (int i = head2[s]; i; i = Next2[i])
        q.push(paii(-d[ver2[i]], ver2[i]));
    while (!q.empty())
    {
        int x = q.top().second;
        q.pop();
        if (vis[x])
            continue;
        vis[x] = 1;
        for (int i = head1[x]; i; i = Next1[i])
        {
            int y = ver1[i], z = edge[i];
            if (d[y] > d[x] + z)//更新距离
                d[y] = d[x] + z;
            if (c[y] == s)	//还再当前联通块内，不严谨，应该d[y]被更新才加进去
                q.push(paii(-d[y], y));
            else	//可以到另外一个联通块，如果入度为0，就加入拓扑序
            {
                if (--deg[c[y]] == 0)
                    top.push(c[y]);
            }
        }
    }
}

void topsort()
{
    for (int i = 1; i <= cnt; ++i)
        if (deg[i] == 0)
            top.push(i);
    while (!top.empty())
    {
        int x = top.front();
        top.pop();
        dijkstra(x);
    }
}

int main()
{
    scanf("%d %d %d %d", &n, &m, &p, &st);
    int x, y, z;
    for (int i = 0; i < m; ++i)
    {
        scanf("%d %d %d", &x, &y, &z);
        add(x, y, z);
        add(y, x, z);
    }
    for (int i = 0; i <= n; ++i)
        if (!c[i])
        {
            ++cnt;
            dfs(i);
        }
    for (int i = 0; i < p; ++i)
    {
        scanf("%d %d %d", &x, &y, &z);
        add(x, y, z);
        ++deg[c[y]]; //联通块整体入度加
    }
    memset(d, 0x3f, sizeof(d));//初始化
    d[st] = 0;
    topsort();
    for (int i = 1; i <= n; ++i)
        if (d[i] > 1000000000)
            printf("NO PATH\n");
        else 
            printf("%d\n", d[i]);
    return 0;
}