find the safest road
Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 9162 Accepted Submission(s): 3236
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所在的城市和要去的城市
第一行: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)<span style="font-size:18px;"><span style="font-size:18px;">#include<cstdio> #include<algorithm> using namespace std; double map[1010][1010],d[1010]; bool mark[1010]; int n; void dijkstra(int a) { for(int i=1;i<=n;i++) { d[i]=0; mark[i]=false; } d[a]=1; while(true) { int v=-1; for(int i=1;i<=n;i++) { if(!mark[i]&&(v==-1||d[i]>d[v])) v=i; } if(v==-1) break; mark[v]=true; for(int i=1;i<=n;i++) d[i]=max(d[i],d[v]*map[v][i]); } } int main() { double c; int a,b; while(~scanf("%d",&n)) { int m; for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) { scanf("%lf",&c); map[i][j]=map[j][i]=c; } } scanf("%d",&m); for(int i=1;i<=m;i++) { scanf("%d%d",&a,&b); dijkstra(a); if(d[b]==0) printf("What a pity!\n"); else printf("%.3lf\n",d[b]); } } return 0; } </span></span>
源代码(SPFA)<span style="font-size:18px;">//1005 #include<cstdio> #include<cstring> #include<queue> #include<algorithm> using namespace std; int head[1000*1000],vis[1010];//head数组开的小会超内存 double dis[1010]; int n,m,edgenum; struct Edge { int from,to; double val; int next; }; Edge edge[1000*1000]; void addedge(int u,int v,double w) { Edge E={u,v,w,head[u]}; edge[edgenum]=E; head[u]=edgenum++; } void SPFA(int s) { queue<int> Q; memset(dis,0,sizeof(dis)); memset(vis,0,sizeof(vis)); Q.push(s); dis[s]=1; 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); } } } } } int main() { double k; int a,b; while(~scanf("%d",&n)) { edgenum=0; memset(head,-1,sizeof(head)); for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++)//从1开始 { scanf("%lf",&k); addedge(i,j,k); } } scanf("%d",&m); while(m--) { scanf("%d%d",&a,&b); SPFA(a); if(dis[b]==0) printf("What a pity!\n"); else printf("%.3lf\n",dis[b]); } } return 0; }</span>