hdoj 1596 find the safest rode

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

 

思路：算法还是利用dijkstra算法，只要改成求的是最大距离就好了，要注意数的类型

 

 #include <stdio.h>
 #include <algorithm>
 int n, q;
 double safe[], sa[][];
 int vis[];
 double max(double x, double y)
 {
     return x > y ? x : y;
 }
 void dijkstra(int s, int t)
 {
     int u, v;
     for(u = ; u <= n; u++)
     {
         vis[u] = ;
         safe[u] = 0.0;
     }
     safe[s] = 1.0;
     while(true)
     {
         v= -;
         for(u = ; u <= n; u++)
             if(!vis[u] && (safe[u]>safe[v] || v==-))
                 v = u;
         if(v == -)
             break;
         vis[v] = ;
         for(u = ; u <= n; u++)
         {
             safe[u] = max(safe[u], safe[v]*sa[v][u]);
         }
     }
     if(safe[t] == 0.0)
         printf("What a pity!\n");
     else
         printf("%.3f\n", safe[t]);
 }
 int main()
 {
     double s;
     while(~scanf("%d", &n))
     {
         for(int i = ; i <= n; i++)
             for(int j = ; j <= n; j++)
             {
                 scanf("%lf", &s);
                 if(i <= j)
                     sa[i][j] = sa[j][i] = s;
             }
         scanf("%d", &q);
         while(q--)
         {
             int a, b;
             scanf("%d%d", &a, &b);
             dijkstra(a, b);
         }
     }
     return ;
 }

spfa算法代码：

 #include <stdio.h>
 #include <string.h>
 #include <queue>
 #define INF 0x3f3f3f3f
 #define N 1010
 #define M 10000010
 using namespace std;
 int n, cnt;
 int vis[N], head[N];
 double dis[N];
 struct node
 {
     int from, to, next;
     double val;
 }edge[M];
 void add(int x, int y, double z)
 {
     node e = {x, y, head[x], z};
     edge[cnt] = e;
     head[x] = cnt++;
 }
 void spfa(int s, int e)
 {
     queue<int>q;
     for(int i = ; i <= n; i++)
     {
         dis[i] = 0.0;
     }
     memset(vis, , sizeof(vis));
     q.push(s);
     dis[s] = 1.0;
     vis[s] = ;
     while(!q.empty())
     {
         int u = q.front();
         q.pop();
         vis[u] = ;
         for(int i = head[u]; i != -; i = edge[i].next)
         {
             int v = edge[i].to;
             if(dis[v] < dis[u]*(edge[i].val))
             {
                 dis[v] = dis[u]*(edge[i].val);
                 if(!vis[v])
                 {
                     vis[v] = ;
                     q.push(v);
                 }
             }
         }
     }
     if(dis[e] == 0.0)
         printf("What a pity!\n");
     else
         printf("%.3f\n", dis[e]);
 }
 int main()
 {
     while(~scanf("%d", &n))
     {
         double a;
         cnt = ;
         memset(head, -, sizeof(head));
         for(int i = ; i <= n; i++)
             for(int j = ; j <= n; j++)
             {
                 scanf("%lf", &a);
                 add(i, j, a);
             }

         int b;
         scanf("%d", &b);
         while(b--)
         {
             int start, end;
             scanf("%d%d", &start, &end);
             spfa(start, end);
         }
     }
     return ;
 }