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;
}