find the safest road HDU杭电1596【Dijkstra || SPFA】

http://acm.hdu.edu.cn/showproblem.php?pid=1596


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
 

SPFA

#include <cstdio>
#include <cstring>
#include <queue>
#define MAXN 1100
using namespace std;

double map[MAXN][MAXN];
int vis[MAXN];//判断是否加入队列了 
int num;
double low[MAXN];//存最短路径 
int s, e;
int M, N;
void SPFA()
{
	int i, j;
	queue<int> Q;
	memset(low, 0, sizeof(low)); 
	memset(vis, 0, sizeof(vis));	
	vis[s] = 1;
	low[s] = 1;
	Q.push(s);
	while(!Q.empty())
	{
		int u = Q.front();
		Q.pop();
		vis[u] = 0;//出队列了,不在队列就变成0 
		for(i = 1; i <= N; ++i)
		{
			
			if(low[i] < low[u] * map[u][i])
			{
				low[i] = low[u] * map[u][i];
				if(!vis[i])				 
			 	{
			 		vis[i]=1;
			 		Q.push(i);
			 	}
			}
		}
	}
	if(low[e] == 0) printf("What a pity!\n");
	else printf("%.3f\n",low[e]);
}
int main()
{
	int u, v;
	double w;

	while(~scanf("%d", &N))
	{
		num=0;
		for(int i=1 ;i <= N; ++i)
		{
			for(int j=1; j <= N; ++j)
			{
				scanf("%lf",&w);
				map[i][j]=w;
			}
		}
		scanf("%d", &M);
		while(M--)
		{
			scanf("%d%d", &s, &e);
			SPFA();
		}	
		
	}
	return 0;
}


Dijkstra

#include<stdio.h>
#include<string.h>
double map[1010][1010];
double dis[1010];
bool used[1010];
int n;
int i,j;
void dijkstra(int u)
{
	memset(used,0,sizeof(used));
	for(i=1;i<=n;++i)
		dis[i]=0;

	int pos=u;
	for(i=1;i<=n;++i)//第一次给dis赋值 
	{
		dis[i]=map[u][i];
	}
	dis[u]=1;
	used[u]=1;
	for(i=1;i<n;++i)//再找n-1个点
	{
		double max=0;
		for(j=1;j<=n;++j)
		{
			if(!used[j]&&max<dis[j])
			{
				max=dis[j];
				pos=j;
			}
		} 		
		used[pos]=1;
		dis[pos]=max;
		for(j=1;j<=n;++j)//把dis数组更新,也叫松弛
		{
			if(!used[j]&&dis[j]<map[pos][j]*dis[pos])
			{
				dis[j]=map[pos][j]*dis[pos];
			}
		}
	}
}
int main()
{
	int m;
	double w;
	int u,v;
	while(~scanf("%d",&n))
	{		
		for(i=1;i<=n;++i)
      		for(j=1;j<=n;++j)
        	{
        		scanf("%lf",&w);
        		map[i][j]=w;
        	}
        	
        scanf("%d",&m);
        while(m--)
        {
        	scanf("%d%d",&u,&v);
           	dijkstra(u);
           	if(dis[v]==0) printf("What a pity!\n");
           	else
        		printf("%.3lf\n",dis[v]);
        }		
	}
	return 0;
}

/*
3
1 0.5 0.5
0.5 1 0.2
0.5 0.2 1
3
1 2
2 3
1 3
*/

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值