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 0x3f3f3f
int 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 0x3f3f3f
using 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;
}