find the safest road
Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 9081 Accepted Submission(s): 3193
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
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#define mx 0x3f3f3f3f
using namespace std;
double map[1010][1010];
double dis[1010];
int vis[1010];
int n;
void djs(int x)
{
int i,j,k;
memset(vis,0,sizeof(vis));
for(i=1;i<=n;i++)
dis[i]=0; //题中0表示不通
dis[x]=1;
for(j=1;j<n;j++)
{
k=-1;
for(i=1;i<=n;i++)
{
if(!vis[i]&&(k==-1||dis[i]>dis[k]))//为了找最安全的也就是找安全系数最大的，所以找最大的
k=i;
}
if(k==-1)
break;
vis[k]=1;
for(i=1;i<=n;i++)
dis[i]=max(dis[i],dis[k]*map[k][i]);//因为安全系数是相乘的关系
}
}
int main(){
int a,b,c,i,j;
int n1,a1,b1;
double m;
while(scanf("%d",&n)!=EOF)
{
memset(map,mx,sizeof(map));
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
{
scanf("%lf",&map[i][j]);
}
scanf("%d",&n1);
while(n1--)
{
scanf("%d%d",&a1,&b1);
djs(a1);
if(dis[b1]==0)//如果不通
printf("What a pity!\n");
else
printf("%.3lf\n",dis[b1]);
}
}
return 0;
}
//SPFA解
#include<stdio.h>
#include<string.h>
#include<queue>
#include<algorithm>
#define INF 0x3f3f3f3f
using namespace std;
int n;
struct Edge
{
int from,to;
double val;
int next;
}edge[2100000];//这块很坑，给的点的数量是1000，变数应该有1000*1000条
int head[1100],edgenum;
void add(int u,int v,double w)
{
Edge E={u,v,w,head[u]};
edge[edgenum]=E;
head[u]=edgenum++;
}
int a,b;
double dis[1100];
int vis[1100];
void SPFA(int x)
{
queue<int>q;
memset(vis,0,sizeof(vis));
memset(dis,0,sizeof(dis));
q.push(x);
vis[x]=1;
dis[x]=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);
}
}
}
}
if(dis[b]==0)
printf("What a pity!\n");
else
printf("%.3lf\n",dis[b]);
}
int main(){
while(scanf("%d",&n)!=EOF)
{
int i,j,n1;
double m;
memset(head,-1,sizeof(head));
edgenum=0;
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
{
scanf("%lf",&m);
if(m)
add(i,j,m);
}
}
scanf("%d",&n1);
while(n1--)
{
scanf("%d%d",&a,&b);
SPFA(a);
}
}
return 0;
}
