POJ1097 Roads Scholar

1. #include<vector>
2. #include<iostream>
3. #include<iomanip>
4. #include<string>
6. using namespace std;
8. const int MAXN = 30;
9. const int MAXDIST = 100000000;
11. int distanc[MAXN][MAXN];
13. //从交点a到b的最短距离
14. int dist[MAXN][MAXN];
15. //从交点a到b最段路径要经过的下一个交点
16. int next[MAXN][MAXN];
17. //在交点a处的城市的名称
18. string name[MAXN];
19. //all the Cs that the shortest path from A to C will pass B.
20. vector<int> passby[MAXN][MAXN];
22. //调试用方法， 打印二维数组的内容
23. void ShowMatrix(int matrix[MAXN][MAXN], int rows, int cols)
24. {
25.     int i, j;
26.     for ( i = 0; i < rows; i ++ )
27.     {
28.         for ( j = 0; j < cols; j ++ )
29.             cout << setw(11) <<  matrix[i][j] << " ";
30.         cout << endl << endl << endl;
31.     }
32.     cout << " ------------------------ " << endl;
33. }
34. //调试方法， 打印出从节点a 到 b的最段路径
35. void FindPath(int a, int b)
36. {
37.     while ( a != b )
38.     {
39.         cout << a << "->" ;
40.         a = next[a][b];
41.     }
42.     cout << b << endl;
43. }
44. //对道路通往的城市排序
45. //排序方法， 按 （距离， 城市名）排序
46. void sortD(vector<int> &verc, int a, int b, int dis)
47. {
48.     int i, j, di, dj, t;
49.     for ( i = 0; i < verc.size(); i ++)
50.     {       
51.         for ( j = i + 1; j < verc.size(); j ++)
52.         {           
53.             di = ( dist[a][verc[i]] - dis + 50 ) / 100;
54.             dj = ( dist[a][verc[j]] - dis + 50 ) / 100;
55.             if ( di > dj || di == dj && name[verc[i]]  > name[verc[j]])
56.             {
57.                 t = verc[i];
58.                 verc[i] = verc[j];
59.                 verc[j] = t;
60.             }
61.         }
62.     }
63. }
65. int main()
66. {
67.     int n,/*Number of intersections*/
68.         m,/*Number of roads*/
69.         k,/*Number of cities*/ 
70.         s,/*Number of Signs*/
71.         a, b, dis, i, j;
72.     float d;
74.     cin >> n >> m >> k;
75.     for ( a = 0; a < n; a ++ )
76.         for ( b = 0; b < n ; b ++ )
77.             dist[a][b] = distanc[a][b] = MAXDIST;   
78.     for ( a = 0; a < n; a++)
79.         dist[a][a] = distanc[a][a] = 0;
80.     
81.     //读入铁路图信息
82.     for ( i = 0; i < m; i ++ )
83.     {
84.         cin >> a >> b >> d;
85.         dis = (int) (d * 100 + 0.5f); //加0.5避免浮点数精度错误
86.         distanc[a][b] = distanc[b][a] 
87.             = dist[a][b] = dist[b][a] = dis;
88.     }
89.     //读入城市名称
90.     for ( i = 0; i < k; i ++ )
91.     {
92.         cin >> a;
93.         cin >> name[a];
94.     }
96.     /************************************************************************/
97.     /* 以下程序段计算交点间最短路径并记录路径经由的下一个交点               */
98.     /************************************************************************/
99.     for ( a = 0; a < n; a ++ )
100.         for ( b = 0; b < n; b ++ )
101.             if ( dist[a][b] < MAXDIST )
102.                 next[a][b] = b;
103.     //Dijstra's algorithm.
104.     for ( s = 0; s < n; s ++ )
105.         for ( a = 0; a < n; a ++ )
106.             if ( dist[a][s] < MAXDIST )
107.                for ( b = 0; b < n; b ++ )
108.                    if ( dist[a][b] > dist [a][s] + dist[s][b] )
109.                    {
110.                        dist[a][b] = dist[a][s] + dist[s][b];
111.                    }
112.     //最短路径下一步节点
113.     for ( a = 0; a < n; a ++ )
114.         for ( b = 0; b < n; b ++)
115.             if ( a != b)
116.                 for ( s = 0; s < n; s ++)
117.                     if ( s != a && s != b &&  distanc[a][s] < MAXDIST 
118.                         && dist[a][b] == dist[a][s] + dist[s][b] )
119.                     {
120.                         next[a][b] = s;
121.                         break;
122.                     }
123.     
124.     for ( a = 0; a < n; a ++ )
125.         for ( b = 0; b < n; b ++ )
126.             if ( a != b )
127.             {
128.                 passby[a][next[a][b]].push_back(b);
129.             }
130.     
131.     cin >> s;
132.     for ( i = 0 ; i < s ; i ++ )
133.     {
134.         cin >> a >> b >> d;
135.         dis = (int) ( d * 100 + 0.5f );
136.         sortD( passby[a][b], a, b, dis );
137.         for (j = 0; j < passby[a][b].size(); j++)
138.         {
139.             k = passby[a][b][j];
140.             if ( name[k] != "" )
141.                 cout << setw(20) << setiosflags(ios_base::left) << name[k] 
142.                      << ( dist[a][k] - dis + 50 ) / 100 << endl;            
143.         }
144.         cout << endl;
145.     }
146.     return 0;
147. }