Program4_P

 现在做的是第四专题编号为1016的试题，具体内容如下所示：

Problem P

Time Limit : 10000/5000ms (Java/Other)   Memory Limit : 32768/32768K (Java/Other)

Total Submission(s) : 27   Accepted Submission(s) : 6

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

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

 

Output

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

 

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
  
  

简单题意：

找一条最安全的路

解题思路：

找一条最安全的路则需要安全度最大，因为ek是P上的边，而且s是在0~1之间的数，相乘的越少则安全度越大，所以还是一个最短路径的问题，可以用dijstral算法求得。

编写代码：

           #include <stdio.h>  
    #include <iostream>  
    #include <string.h>  
    #define maxn 10000  
    using namespace std;  
    int vis[10100],n,mm;  
    double map[1010][1010],d[10100];  
    void dijkstra(int start)  
    {  
        int i,y,x;  
        double m;  
        memset(vis,0,sizeof(vis));  
        for(i=1;i<=n;i++)  
            d[i]=map[start][i];  
        for(i=1;i<=n;i++)  
        {  
            m=-maxn;  
            for(y=1;y<=n;y++)  
                if(!vis[y]&&d[y]>m)  
                {  
                    m=d[y];  
                    x=y;  
                }  
            vis[x]=1;  
            for(y=1;y<=n;y++)  
                d[y]=max(d[y],d[x]*map[x][y]);  
        }  
    }  
    int main ()  
    {  
        int m,i,j;  
        while (scanf("%d",&n)!=EOF)  
        {  
            for(i=1;i<=n;i++)  
                for(j=1;j<=n;j++)  
                    scanf("%lf",&map[i][j]);  
            scanf("%d",&m);  
            while (m--)  
            {  
                scanf("%d %d",&i,&j);  
                dijkstra(i);  
                if (d[j]>0)   
                    printf ("%.3lf\n",d[j]);  
                else  
                    printf ("What a pity!\n");  
            }  
        }  
        return 0;  
    }