find the safest road
Time Limit : 10000/5000ms (Java/Other) Memory Limit : 32768/32768K (Java/Other)
Total Submission(s) : 23 Accepted Submission(s) : 12
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
代码：
1 #include<stdio.h> 2 #include<string.h> 3 #define MAX(x,y)(x>y?x:y) 4 const int MAXN=1010; 5 double map[MAXN][MAXN],d[MAXN]; 6 int vis[MAXN]; 7 int n; 8 void initial(){ 9 memset(d,0,sizeof(d)); 10 memset(vis,0,sizeof(vis)); 11 } 12 void dijskra(int s){ 13 initial();//写错位置了，应该每次都初始化 14 d[s]=1; 15 while(true){ 16 int k=-1; 17 for(int i=1;i<=n;i++) 18 if(!vis[i]&&(k==-1||d[i]>d[k]))k=i; 19 if(k==-1)break; 20 vis[k]=1; 21 for(int i=1;i<=n;i++) 22 d[i]=MAX(d[i],d[k]*map[k][i]); 23 } 24 } 25 int main(){int Q,a,b; 26 while(~scanf("%d",&n)){ 27 memset(map,0,sizeof(map)); 28 for(int i=1;i<=n;i++) 29 for(int j=1;j<=n;j++) 30 scanf("%lf",&map[i][j]); 31 scanf("%d",&Q); 32 while(Q--){ 33 scanf("%d%d",&a,&b); 34 dijskra(a); 35 if(d[b]==0.0)puts("What a pity!"); 36 else printf("%.3lf ",d[b]); 37 } 38 } 39 return 0; 40 }
SPFA:
1 #include<stdio.h> 2 #include<string.h> 3 #include<queue> 4 using namespace std; 5 const int MAXN=1010; 6 const int MAXM=2000010; 7 int head[MAXM],vis[MAXN]; 8 double dis[MAXN]; 9 struct Edge{ 10 int from,to,next; 11 double value; 12 }; 13 Edge edg[MAXM]; 14 int n,edgnum; 15 queue<int>dl; 16 void initial(){ 17 memset(dis,0,sizeof(dis)); 18 memset(vis,0,sizeof(vis)); 19 while(!dl.empty())dl.pop(); 20 } 21 void add(int u,int v,double value){ 22 Edge E; 23 E.from=u;E.to=v;E.value=value;E.next=head[u]; 24 edg[edgnum]=E; 25 head[u]=edgnum++; 26 } 27 void SPFA(int sx){ 28 initial(); 29 dis[sx]=1; 30 vis[sx]=1; 31 dl.push(sx); 32 while(!dl.empty()){ 33 int k=dl.front(); 34 dl.pop(); 35 vis[k]=0; 36 for(int i=head[k];i!=-1;i=edg[i].next){ 37 int v=edg[i].to; 38 if(dis[k]*edg[i].value>dis[v]){ 39 dis[v]=dis[k]*edg[i].value; 40 if(!vis[v]) 41 dl.push(v),vis[v]=1; 42 } 43 } 44 } 45 } 46 int main(){ 47 int a,b,Q; 48 double v; 49 while(~scanf("%d",&n)){ 50 memset(head,-1,sizeof(head)); 51 edgnum=0; 52 for(int i=1;i<=n;i++) 53 for(int j=1;j<=n;j++){ 54 scanf("%lf",&v); 55 add(i,j,v); 56 } 57 scanf("%d",&Q); 58 while(Q--){ 59 scanf("%d%d",&a,&b); 60 SPFA(a); 61 if(dis[b]==0.0)puts("What a pity!"); 62 else printf("%.3lf ",dis[b]); 63 } 64 } 65 return 0; 66 }