find the safest road
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模板题，本题是使用dijkstra寻找最大安全系数，注意对安全系数是进行乘法运算。
用一个double数组safety储存从起点到所有城市的最大安全系数，一个布尔数组记录是否固定对应城市的最大安全系数。初始状态将到所有城市的安全系数视为0，从起点开始标记，记录起点到自己的安全系数为1并将其固定，之后每次找到并固定到达起点安全系数最大的城市，固定一个城市后查看以该城市为媒介是否能使还没有固定的城市到达起点的安全系数增高，如果能优化，将优化后的记入该城市对应的safety中。
1 #include <bits/stdc++.h> 2 using namespace std; 3 const int maxn = 1e3+10; 4 int n, q, s; 5 double G[maxn][maxn]; //G记录邻接矩阵 6 double safety[maxn]; //safety记录起点到各城市的最大安全系数 7 bool vis[maxn] = {false}; //vis记录城市最大安全系数是否已经固定 8 9 void dijkstra(int s){ //传入起点 10 memset(safety, 0, sizeof(safety)); //初始化起点到所有城市的最大安全系数为0 11 memset(vis, false, sizeof(vis)); //所有城市都为未固定状态 12 safety[s] = 1.0; //起点到自己的最大安全系数为1 13 for(int i = 0; i < n; i++){ 14 //开始计算并固定城市最大安全系数，需要将每个城市都固定所以循环n次 15 int maxSafetyCity = -1; 16 double maxSafetyNum = -1; 17 for(int j = 0; j < n; j++){ //找到还没有固定的城市中到起点安全系数最大的城市 18 if(!vis[j] && safety[j] > maxSafetyNum){ 19 maxSafetyCity = j; //maxSafeCity记录还没有固定的城市中到起点安全系数最大的城市 20 maxSafetyNum = safety[j]; //maxSafetyNum记录当前最大安全系数 21 } 22 } 23 if(maxSafetyCity == -1){ //如果找不到符合条件的城市，证明其余城市与起点不连通 24 return; 25 } 26 vis[maxSafetyCity] = true; //固定安全系数最大的城市 27 for(int j = 0; j < n; j++){ //寻找以该城市为媒介可以优化的城市 28 if(!vis[j] && G[maxSafetyCity][j] != 0 && G[maxSafetyCity][j] * safety[maxSafetyCity] > safety[j]){ 29 //可以优化的条件为，该城市没有固定 30 //该城市与媒介城市连通 31 //从媒介城市到该城市的安全系数乘以媒介城市到起点的安全系数比不优化前该城市到达起点的安全系数高 32 safety[j] = safety[maxSafetyCity] * G[maxSafetyCity][j]; 33 //优化 34 } 35 } 36 } 37 } 38 int main() 39 { 40 while(scanf("%d", &n) != EOF){ //输入城市数量 41 for(int i = 0; i < n; i++){ 42 for(int j = 0; j < n; j++){ 43 scanf("%lf", &G[i][j]); //输入邻接矩阵 44 } 45 } 46 scanf("%d", &q); //输入旅行路线 47 int ed; 48 for(int i = 0; i < q; i++){ 49 scanf("%d%d", &s, &ed); //输入起点与终点 50 dijkstra(s - 1); //计算起点到所有城市的最大安全系数 51 if(safety[ed - 1] != 0){ 52 printf("%.3f\n", safety[ed - 1]); //格式化输出最大安全系数 53 }else{ 54 printf("What a pity!\n"); //不连通 55 } 56 } 57 } 58 return 0; 59 }