find the safest road

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

 

注意：这个题是最小路径的变形，只是把求起点到终点最小长度变成了安全系数的乘积的最大，直接套用模板注意题中所说的，0可以理解为那两个城市之间没有直接的通道，为最后结果的判断提供了依据。

 

 

#include<stdio.h>#include<string.h>#define MAX(a,b) a>b?a:b#define INF 0x3f3f3fint n;//n是城市的个数 double map[1100][1100];double cost[1100][1100];double d[1100];int mark[1100];void dijkstra(int s,int e){memset(d,0,sizeof(d)); int i,j;for(i=1;i<=n;i++){d[i]=-100;mark[i]=0;}d[s]=1<span style="color:#ff0000;">;//对起点的处理 </span>int u;while(1){int v=-1;for(u=1;u<=n;u++)  {if(mark[u]==0&&(v==-1||d[u]>d[v]))   v=u;  }  if(v==-1)    {    if(d[e])      printf("%.3lf\n",d[e]);//若安全度为0，则不能到达     else      printf("What a pity!\n");    break;    }    mark[v]=1;    for(u=1;u<=n;u++)          {    //if(mark[u]==0&&d[u]>d[v]*map[v][u])      // d[u]=d[v]*map[v][u];   d[u]=MAX(d[u],d[v]*map[v][u]);//<span style="color:#ff0000;">需要注意的地方 </span>    }}}int main(){while(scanf("%d",&n)!=EOF){        //memset(map,0,sizeof(map));int i,j;for(i=1;i<=n;i++)  for(j=1;j<=n;j++)     scanf("%lf",&map[i][j]);//构建了一个邻接矩阵int Q,be,end;scanf("%d",&Q);//WA了无数次原来少写了这一句 ,无语了 while(Q--){<span style="color:#33ff33;">//for(i=1;i<=n;i++)  //for(j=1;j<=n;j++)原来以为这儿错了，原来不是，只是少了scanf("%d",&Q);     //map[i][j]=cost[i][j];//每一次调用dijkstra函数 的时候 map的值都会改变 ，会导致结果错误，所以每用一次，复制一次 </span>scanf("%d%d",&be,&end);//s是起点城市，e是终点城市dijkstra(be,end); }   }    return 0;} 

 

 SPFA  不知道为什么超时？？？

 

 

//超时 #include<stdio.h>#include<string.h>#include<queue>#include<algorithm>#define MAX 1000 #define INF 0x3f3f3fusing namespace std;double d[MAX];int head[MAX];int mark[MAX];double map[MAX];int edgenum;int n;struct Edge{int from,to,next;double val;}edge[MAX*10];//void addedge(int u,int v,int w)//w是权值，是double型 void addedge(int u,int v,double w){  <span style="color:#009900;"> /* Edge E={u,v,w,head[u]};//</span>

<span style="color:#009900;">edge[edgenum]=E;       head[u]=edgenum++;*/</span>    Edge E;E.from=u;E.to=v;E.val=w;E.next=head[u]; edge[edgenum]=E;//错误在这里！少了这一句 head[u]=edgenum++;}void SPFA(int sx,int ex){memset(mark,0,sizeof(mark));<span style="color:#009900;">//memset(d,INF,sizeof(d));</span>   <span style="color:#ff0000;"> memset(d,0,sizeof(d));//这一</span><span style="color:#cc0000;">点要全部赋值为0，因为刚开始假设所有的点都不能到 。</span> queue<int>q;q.push(sx);     <span style="color:#33cc00;">  </span><span style="color:#009900;"> //d[sx]=0;</span>    d[sx]=1;//<span style="color:#cc0000;">与之前不一样的地方 ，若自己与自己的权值是0的话，后面是相乘的关系，结果都为0。 </span>mark[sx]=1;while(!q.empty()){int u=q.front();   q.pop();   mark[u]=0;   for(int i=head[u];i!=-1;i=edge[i].next)   {     int v=edge[i].to;     if(d[v]<d[u]*edge[i].val)     {      d[v]=d[u]*edge[i].val;      if(mark[v]==0)      {      q.push(v);      mark[v]=1;      }     }   }}if(d[ex]!=0)   printf("%.3lf\n",d[ex]);else   printf("What a pity!\n");}int main(){while(scanf("%d",&n)!=EOF){memset(head,-1,sizeof(head));edgenum=0;int i,j;double a;for(i=1;i<=n;i++){for(j=1;j<=n;j++){scanf("%lf",&a);addedge(i,j,a);<span style="color:#33cc00;">//addedge(j,i,a);</span>}}int Q,b,c;scanf("%d",&Q);while(Q--){scanf("%d%d",&b,&c);SPFA(b,c);}}return 0;}