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[1010], sa[1010][1010];
int vis[1010];
double max(double x, double y)
{
    return x > y ? x : y;
}
void dijkstra(int s, int t)
{
    int u, v;
    for(u = 1; u <= n; u++)
    {
        vis[u] = 0;
        safe[u] = 0.0;
    }
    safe[s] = 1.0;
    while(true)
    {
        v= -1;
        for(u = 1; u <= n; u++)
            if(!vis[u] && (safe[u]>safe[v] || v==-1))
                v = u;
        if(v == -1)
            break;
        vis[v] = 1; 
        for(u = 1; u <= n; u++)
        {
            safe[u] = max(safe[u], safe[v]*sa[v][u]);
        }
    }
    if(safe[t] == 0.0)
        printf("What a pity!
");
    else
        printf("%.3f
", safe[t]);
}
int main()
{
    double s;
    while(~scanf("%d", &n))
    {
        for(int i = 1; i <= n; i++)
            for(int j = 1; 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 0;
}

 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 = 1; i <= n; i++)
    {
        dis[i] = 0.0;
    }
    memset(vis, 0, sizeof(vis));
    q.push(s);
    dis[s] = 1.0;
    vis[s] = 1;
    while(!q.empty())
    {
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for(int i = head[u]; i != -1; 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] = 1;
                    q.push(v);
                }
            }
        } 
    }
    if(dis[e] == 0.0)
        printf("What a pity!
");
    else
        printf("%.3f
", dis[e]);
}
int main()
{
    while(~scanf("%d", &n))
    {
        double a;
        cnt = 0;
        memset(head, -1, sizeof(head));
        for(int i = 1; i <= n; i++)
            for(int j = 1; 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 0;
}