jzxx1853道路和航線
題目描述
Farmer John正在一個新的銷售區域對他的牛奶銷售方案進行調查。他想把牛奶送到T個城鎮(1 <= T <= 25,000)，編號為1…T。這些城鎮之間通過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) 把奶牛送到每個城鎮的最便宜的方案，或者知道這是不可能的。
輸入
第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
輸出
第1到T行：從S到達城鎮i的最小花費，如果不存在輸出"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
來源
USACO 2011 January Gold
滿分程式碼： #include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <vector> #include <queue> #include <stack> #define MAXN 55555 #define INF 1000010080 using namespace std; typedef pair<int, int> P; priority_queue<P, vector<P>, greater<P> > q; vector<P> g[MAXN]; int dis[MAXN], vis[MAXN], cango[MAXN]; int src, n, R, m; int scc, ind; int dfn[MAXN], instack[MAXN], low[MAXN], fa[MAXN], in[MAXN], last[MAXN], num[MAXN], c[MAXN]; stack<int>st; vector<int>bl[MAXN], got[MAXN]; void tarjan(int u) { dfn[u] = low[u] = ++ind; st.push(u); instack[u] = 1; int sz = g[u].size(); for(int i = 0; i < sz; i++) { int v = g[u][i].first; if(!dfn[v]) { tarjan(v); low[u] = min(low[u], low[v]); } else if(instack[v]) low[u] = min(low[u], dfn[v]); } if(dfn[u] == low[u]) { scc++; while(true) { int tmp = st.top(); st.pop(); instack[tmp] = 0; fa[tmp] = scc; bl[scc].push_back(tmp); if(tmp == u) break; } } } int que[MAXN]; int l, r; void dij(int now) { while(!q.empty()) q.pop(); for(int i = 0; i < got[now].size(); i++) { int u = got[now][i]; q.push(P(dis[u], u)); } while(!q.empty()) { P p = q.top(); q.pop(); int u = p.second, d = p.first; if(vis[u]) continue; vis[u] = 1; int sz = g[u].size(); for(int i = 0; i < sz; i++) { int v = g[u][i].first; int w = g[u][i].second; if(fa[v] == fa[u]) { if(d + w < dis[v]) { dis[v] = d + w; q.push(P(dis[v], v)); } } } } } void dfs(int u) { cango[u] = 1; int sz = g[u].size(); for(int i = 0; i < sz; i++) { int v = g[u][i].first; if(!cango[v]) dfs(v); } } void gao() { que[r++] = fa[src]; while(l < r) { int now = que[l++]; dij(now); int sz = bl[now].size(); for(int i = 0; i < sz; i++) { int u = bl[now][i]; for(int j = 0; j < g[u].size(); j++) { int v = g[u][j].first; int w = g[u][j].second; if(fa[v] != now) { --in[fa[v]]; if(dis[u] + w < dis[v]) dis[v] = dis[u] + w, got[fa[v]].push_back(v); if(in[fa[v]] == 0) que[r++] = fa[v]; } } } } } int main() { int u, v, w; scanf("%d%d%d%d", &n, &R, &m, &src); for(int i = 1; i <= n; i++) dis[i] = INF; dis[src] = 0; for(int i = 0; i < R; i++) { scanf("%d%d%d", &u, &v, &w); g[u].push_back(P(v, w)); g[v].push_back(P(u, w)); } for(int i = 0; i < m; i++) { scanf("%d%d%d", &u, &v, &w); g[u].push_back(P(v, w)); } for(int i = 1; i <= n; i++) if(!dfn[i]) tarjan(i); got[fa[src]].push_back(src); dfs(src); for(int i = 1; i <= n; i++) { if(!cango[i]) continue; int sz = g[i].size(); for(int j = 0; j < sz; j++) { v = g[i][j].first; if(fa[v] != fa[i]) in[fa[v]]++; } } gao(); for(int i = 1; i <= n; i++) { if(dis[i] == INF) printf("NO PATH\n"); else printf("%d\n", dis[i]); } return 0; }