杭电oj1596--find the safest road(Spfa)

find the safest road

Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 9271 Accepted Submission(s): 3279

Problem Description

XX星球有很多城市，每个城市之间有一条或多条飞行通道，但是并不是所有的路都是很安全的，每一条路有一个安全系数s,s是在 0 和 1 间的实数(包括0，1)，一条从u 到 v 的通道P 的安全度为Safe(P) = s(e1)*s(e2)…*s(ek) e1,e2,ek是P 上的边，现在8600 想出去旅游，面对这这么多的路，他想找一条最安全的路。但是8600 的数学不好，想请你帮忙 ^_^

Input

输入包括多个测试实例，每个实例包括：
第一行：n。n表示城市的个数n<=1000;
接着是一个n*n的矩阵表示两个城市之间的安全系数，(0可以理解为那两个城市之间没有直接的通道)
接着是Q个8600要旅游的路线,每行有两个数字，表示8600所在的城市和要去的城市

Output

如果86无法达到他的目的地，输出"What a pity!",
其他的输出这两个城市之间的最安全道路的安全系数,保留三位小数。

Sample Input

3

1 0.5 0.5 

0.5 1 0.4

0.5 0.4 1 

3

1 2

2 3

1 3

Sample Output

0.500

0.400

0.500

Author

ailyanlu

Source

HDU 2007-Spring Programming Contest - Warm Up （1）

Recommend

8600 | We have carefully selected several similar problems for you: 1598 1142 1690 1385 1548

Spfa 可处理存在负权的问题。

 1 #include <queue>
 2 #include <cstdio>
 3 #include <cstring>
 4 #include <iostream>
 5 using namespace std;
 6 const int INF = 0x3f3f3f3f;
 7 double map[1010][1010], dis[1010]; int vis[1010]; 
 8 int n;
 9 bool Spfa(int a, int b)
10 {
11     queue<int> q;
12     memset(vis, 0, sizeof(vis));
13     for(int i = 1; i <= n; i++)
14         dis[i] = 0.0;
15     vis[a] = 1; dis[a] = 1.0;
16     q.push(a);
17     while(!q.empty())
18     {
19     //    printf("1\n");
20         int temp = q.front();
21         q.pop();    
22         vis[temp] = 0; 
23         for(int i = 1; i <= n; i++)
24         {
25             if(dis[i] < dis[temp] * map[temp][i])
26             {
27                 dis[i] = dis[temp] * map[temp][i];
28                 if(!vis[i])
29                 {
30                     vis[i] = 1;
31                     q.push(i);
32                 }
33             }
34         }
35     } 
36     if(dis[b] == 0.0)
37         printf("What a pity!\n");
38     else
39         printf("%.3lf\n", dis[b]);
40 }
41 int main()
42 {
43     while(~scanf("%d", &n))
44     {
45         for(int i = 1; i <= n; i++)
46             for(int j = 1; j <= n; j++)
47                 scanf("%lf", &map[i][j]);
48         int a, b; double  max = 0.0;
49         int Q;
50         scanf("%d", &Q);
51         for(int i = 0; i < Q; i++)
52         {
53             int a, b;
54             scanf("%d %d", &a, &b);
55             Spfa(a, b);
56         } 
57     }
58     return 0;    
59 }

Floyd:

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 using namespace std;
 5 double map[1010][1010]; 
 6 int a[1010], b[1010];
 7 int n;
 8 void Floyd()
 9 {
10     for(int k = 1; k <= n; k++)
11         for(int i = 1; i <= n; i++)
12             for(int j = 1; j <= n; j++)
13                 if(map[i][j] < map[i][k]*map[k][j])
14                     map[i][j] = map[i][k]*map[k][j];    
15 } 
16 int main()
17 {
18     while(~scanf("%d", &n))
19     {
20         for(int i = 1; i <= n; i++)
21             for(int j = 1; j <= n; j++)
22                 scanf("%lf", &map[i][j]);
23         int Q;
24         scanf("%d", &Q);
25         for(int i = 0; i < Q; i++)
26             scanf("%d %d", &a[i], &b[i]);    
27     //    printf("1\n");
28         Floyd();
29         for(int i = 0; i < Q; i++)
30         {
31             if(!map[a[i]][b[i]])
32                 printf("What a pity!\n");
33             else
34                 printf("%.3lf\n", map[a[i]][b[i]]);  
35         }
36     }
37     return 0;    
38 }