[USACO 2011Jan Gold]道路与航线

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) 把奶牛送到每个城镇的最便宜的方案，或者知道这是不可能的。 

**********************************************************************************************************************

Input Format
* 第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 Format

* 第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

**********************************************************************************************************************

Sample Output

NO PATH
NO PATH
5
0
-95
-100

 

 

 

 

 

 

 

 

 

题意：

给定一个图，有两种边，一种是有向，权值可正可负；另一种无向，权值恒为正；图满足，若a到b有一条有向边，不存在b到a的路。给一个点，求这个点到所有点最短路。

大致思路：

按无向边缩点成一个拓补图，作拓补排序，对于每个拓补点，肯定不会有连续的两段有向边，所以可以用dijkstra。

代码：

#include <cstdio> #include <queue> #include <cstring> #include <iostream> #include <algorithm> using namespace std; #define MAXN 25100 #define MAXM 101000 #define INF 1000000000 struct EDGE { int t, next, v; EDGE( int t = 0, int n = 0, int v = 0 ): t(t), next(n), v(v) {} } edge[5][MAXM]; struct NODE { int t, v; NODE( int t = 0, int v = 0 ): t(t), v(v) {} }; struct cmp { bool operator () ( const NODE &p, const NODE &q ) const { return p.v > q.v; } }; typedef priority_queue< NODE, vector<NODE>, cmp > pq; int n, m1, m2, S; int fa[MAXN], cnt[4], vex[4][MAXN], deg[MAXN], ans[MAXN]; pq topo, dij; bool p1[MAXN], p[MAXN]; int find( int x ) { int y = x; while ( y != fa[y] ) y = fa[y]; while ( x != y ) { int tmp = fa[x]; fa[x] = y; x = tmp; } return y; } void add_edge( int id, int x, int y, int va ) { int cc = ++cnt[id]; edge[id][cc] = EDGE( y, vex[id][x], va ); vex[id][x] = cc; } void dijkstra() { memset( p, 0, sizeof(p) ); while ( !dij.empty() ) { int x = dij.top().t; dij.pop(); if ( p[x] ) continue; p[x] = true; int pos = vex[3][x]; while ( pos != -1 ) { int y = edge[3][pos].t; if ( p[y] == false && ans[x] + edge[3][pos].v < ans[y] ) { ans[y] = ans[x] + edge[3][pos].v; dij.push( NODE( y, ans[y] ) ); } pos = edge[3][pos].next; } } } void build( int t ) { if ( p1[t] == true ) return; p1[t] = true; cnt[3] = 0; memset( vex[3], -1, sizeof(vex[3]) ); while ( !dij.empty() ) dij.pop(); int pos = vex[1][t]; while ( pos != -1 ) { int x = edge[1][pos].t; if ( ans[x] < INF && x != S ) add_edge( 3, S, x, ans[x] ); int po2 = vex[0][x]; while ( po2 != -1 ) { int y = edge[0][po2].t; add_edge( 3, x, y, edge[0][po2].v ); po2 = edge[0][po2].next; } if ( ans[x] < INF ) dij.push( NODE( x, ans[x] ) ); pos = edge[1][pos].next; } dijkstra(); pos = vex[1][t]; while ( pos != -1 ) { int x = edge[1][pos].t; int po2 = vex[2][x]; while ( po2 != -1 ) { int y = edge[2][po2].t; if ( ans[x] < INF && ans[x] + edge[2][po2].v < ans[y] ) ans[y] = ans[x] + edge[2][po2].v; --deg[find(y)]; topo.push( NODE( find(y), deg[find(y)] ) ); po2 = edge[2][po2].next; } pos = edge[1][pos].next; } } int main() { memset( vex, -1, sizeof(vex) ); memset( cnt, 0, sizeof(cnt) ); memset( deg, 0, sizeof(deg) ); memset( p1, 0, sizeof(p1) ); for ( int i = 1; i < MAXN; ++i ) { fa[i] = i; ans[i] = INF; } while ( !topo.empty() ) topo.pop(); scanf( "%d%d%d%d", &n, &m1, &m2, &S ); for ( int i = 1; i <= m1; ++i ) { int a, b, c; scanf( "%d%d%d", &a, &b, &c ); add_edge( 0, a, b, c ); add_edge( 0, b, a, c ); int fa1 = find(a); int fa2 = find(b); fa[fa2] = fa1; } for ( int i = 1; i <= n; ++i ) add_edge( 1, find(i), i, 0 ); for ( int i = 1; i <= m2; ++i ) { int a, b, c; scanf( "%d%d%d", &a, &b, &c ); add_edge( 2, a, b, c ); ++deg[find(b)]; } ans[S] = 0; for ( int i = 1; i <= n; ++i ) if ( fa[i] == i ) topo.push( NODE( i, deg[i] ) ); while ( !topo.empty() ) { build( topo.top().t ); topo.pop(); } for ( int i = 1; i <= n; ++i ) { if ( ans[i] == INF ) puts( "NO PATH" ); else printf( "%d/n", ans[i] ); } return 0; }