题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1596
解题思路：Floyd 算法（邻接矩阵）
1 /************************************************************************** 2 user_id: SCNU20102200088 3 problem_id: hdoj 1596 4 problem_name: find the safest road 5 **************************************************************************/ 6 7 #include <algorithm> 8 #include <iostream> 9 #include <iterator> 10 #include <iomanip> 11 #include <sstream> 12 #include <fstream> 13 #include <cstring> 14 #include <cstdlib> 15 #include <climits> 16 #include <bitset> 17 #include <string> 18 #include <vector> 19 #include <cstdio> 20 #include <cctype> 21 #include <ctime> 22 #include <cmath> 23 #include <queue> 24 #include <stack> 25 #include <list> 26 #include <set> 27 #include <map> 28 using namespace std; 29 30 //线段树 31 #define lson l,m,rt<<1 32 #define rson m+1,r,rt<<1|1 33 34 //手工扩展栈 35 #pragma comment(linker,"/STACK:102400000,102400000") 36 37 const double EPS=1e-9; 38 const double PI=acos(-1.0); 39 const double E=2.7182818284590452353602874713526; //自然对数底数 40 const double R=0.5772156649015328606065120900824; //欧拉常数:(1+1/2+...+1/n)-ln(n) 41 42 const int x4[]={-1,0,1,0}; 43 const int y4[]={0,1,0,-1}; 44 const int x8[]={-1,-1,0,1,1,1,0,-1}; 45 const int y8[]={0,1,1,1,0,-1,-1,-1}; 46 47 typedef long long LL; 48 49 typedef int T; 50 T max(T a,T b){ return a>b? a:b; } 51 T min(T a,T b){ return a<b? a:b; } 52 T gcd(T a,T b){ return b==0? a:gcd(b,a%b); } 53 T lcm(T a,T b){ return a/gcd(a,b)*b; } 54 55 /////////////////////////////////////////////////////////////////////////// 56 //Add Code: 57 const int maxn=1005; 58 int n; 59 double dist[maxn][maxn],f[maxn][maxn]; 60 61 void Floyd(){ 62 int i,j,k; 63 for(i=1;i<=n;i++){ 64 for(j=1;j<=n;j++) dist[i][j]=f[i][j]; 65 } 66 for(k=1;k<=n;k++){ 67 for(i=1;i<=n;i++){ 68 for(j=1;j<=n;j++){ 69 if(dist[i][k]*dist[k][j]>dist[i][j]) dist[i][j]=dist[i][k]*dist[k][j]; 70 } 71 } 72 } 73 } 74 /////////////////////////////////////////////////////////////////////////// 75 76 int main(){ 77 std::ios::sync_with_stdio(false); 78 //freopen("in.txt","r",stdin); 79 //freopen("out.txt","w",stdout); 80 /////////////////////////////////////////////////////////////////////// 81 //Add Code: 82 int q,i,j,start,end; 83 while(scanf("%d",&n)!=EOF){ 84 for(i=1;i<=n;i++){ 85 for(j=1;j<=n;j++) scanf("%lf",&f[i][j]); 86 } 87 Floyd(); 88 scanf("%d",&q); 89 while(q--){ 90 scanf("%d%d",&start,&end); 91 if(dist[start][end]==0) printf("What a pity! "); 92 else printf("%.3lf ",dist[start][end]); 93 } 94 } 95 /////////////////////////////////////////////////////////////////////// 96 return 0; 97 } 98 99 /************************************************************************** 100 Testcase: 101 Input: 102 3 103 1 0.5 0.5 104 0.5 1 0.4 105 0.5 0.4 1 106 3 107 1 2 108 2 3 109 1 3 110 Output: 111 0.500 112 0.400 113 0.500 114 **************************************************************************/