HDU 1596 find the safest road【dijkstra & floyed & SPFA 】

find the safest road

Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 9314    Accepted Submission(s): 3296

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
  
  

 

 

 

 

 

 

 

floyd算法：

主要是设置断点再分别判断取符合题意的那个。

首先可以将问题转化为在一个二维数组里的问题，设求如何从A到B城市走才能使得路径的安全系数最高，设总有N个城市，每两个城市之间都有一个安全系数。设置i（从1到N）为断点，比较之前AB两城市之间的安全系数与Ai城市与Bi城市之间安全系数的乘积的大小，在这里取大的那一个。当i为断点是遍历每两个城市，之后再是i从1取到N，既得每两个城市之间的最大安全系数。hdoj1596

//floyd算法 

#include<stdio.h>
#include<string.h>
double city[10010][10010];
double max(double a,double b)
{
	return a>b?a:b;
}
void floyd(int n)
{
	int i,j,k;
	for(k=1;k<=n;++k)<span style="font-size: 18px;">//k为断点</span>
	{
		for(i=1;i<=n;++i)
		{
			for(j=1;j<=n;++j)
			{
				city[i][j]=max(city[i][j],city[i][k]*city[k][j]);<span style="font-size: 18px;">//运用自定义函数</span>
			//	if(city[i][j]<city[i][k]*city[k][j])<span style="font-size: 18px;">也可直接比较</span>
			//	city[i][j]=city[i][k]*city[k][j];
			}
		}
	}
}
int main()
{
	int n,m,i,j,s,e;
	while(~scanf("%d",&n))<span style="font-size: 18px;">//n个城市</span>
	{
		for(i=1;i<=n;++i)
		{
			for(j=1;j<=n;++j)
			{
				scanf("%lf",&city[i][j]);
			}
		}
		floyd(n);<span style="font-size: 18px;">//只运用一次，便可求得所有城市之间的最大安全系数</span>
		scanf("%d",&m);<span style="font-size: 18px;">//m组出发地与目的地</span>
		for(i=0;i<m;++i)
		{
			scanf("%d%d",&s,&e); 
		    if(city[s][e]==0)printf("What a pity!\n");
		    else printf("%.3lf\n",city[s][e]);
		}
	}
	return 0;
}

dijkstra算法：

确定出发地A，从A邻接并未访问过的城市中选择距离最近的i（即最符合题意的，在上面的例子中应取安全系数最高的），再以该城市i为断点判断起点A到各个未访问过的城市X之前的距离AX（或安全系数）与A到i的与i到X的距离之和（或者是Ai与iB的安全系数的乘积）孰大孰小，取符合题意的更新数据。之后再在未访问过的城市中继续寻找i直到将所有城市都遍历过，既得从A到其余各城市的最短距离（或最大安全系数）。hdoj1596（还是上面那道题）

//dijkstra算法 

#include<stdio.h>
#include<string.h>
#define max 10010
double city[max][max],saf[max];//saf为记录从起点到各个城市之间的安全系数</span>
int visit[max];//记录该城市是否访问过</span>
void dijkstra(int s,int e,int n)
{
	int i,j,next;//next标记选出的断点</span>
        double ma;//为判断的标准的初始值</span>
	for(i=0;i<=n;++i)
	{
		visit[i]=0;//开始时所有城市都未访问过</span>
	}
	visit[s]=1;//起点标记为已访问</span>
	for(i=1;i<=n;++i)//将city中的从起点开始的数据转移到saf中</span>
	{
		saf[i]=city[s][i];
	}
	for(i=1;i<=n;++i)
	{
		ma=0;//求最大安全系数，初始值可以置为0</span>
		for(j=1;j<=n;++j)//求断点</span>
		{
			if(!visit[j]&&saf[j]>ma)
			{
				ma=saf[j];
				next=j;
			}
		} 
		visit[next]=1;//在这里已将i为断点即以已取的形式存在</span>
		for(j=1;j<=n;++j)//遍历其他未访问的城市</span>
		{
			if(!visit[j]&&saf[j]<saf[next]*city[next][j])
			{
				saf[j]=saf[next]*city[next][j];
			}
		}
	}
	if(saf[e]==0)printf("What a pity!\n");
    else printf("%.3lf\n",saf[e]);
}
int main()
{
	int i,j,n,m,s,e;
	while(~scanf("%d",&n))
	{
		for(i=1;i<=n;++i)
		{
			for(j=1;j<=n;++j)
			{
				scanf("%lf",&city[i][j]);
			}
		}
		scanf("%d",&m);
		while(m--)
		{
			scanf("%d%d",&s,&e);
			if(s==e)printf("1.000\n"); 
			else dijkstra(s,e,n); //即需要对于每个不同的起点调用一次函数</span>
		} 
	}
	return 0;
}

 

 

 

SPFA

#include <iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define maxn 1100
#define maxm 1000000+10
#define inf 0x3f3f3f3f
using namespace std;
int n,q,cnt,head[maxn];
struct node
{
    int from,to;
    double val;
    int next;
};
node edge[maxm];
void initialize()
{
    cnt=0;
    memset(head,-1,sizeof(head));
}
void add(int a,int b,double c)
{
    node E={a,b,c,head[a]};
    edge[cnt]=E;
    head[a]=cnt++;
}
void spfa(int s,int e)
{
    queue<int>q;
    int vis[maxn];
    double dis[maxn];
    memset(vis,0,sizeof(vis));
    memset(dis,0,sizeof(dis));
    dis[s]=1;
    vis[s]=1;
    q.push(s);
    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])
                {
                    q.push(v);
                    vis[v]=1;
                }
            }
        }
    }
    if(dis[e]==0)
        printf("What a pity!\n");
    else
        printf("%.3lf\n",dis[e]);
}
int main()
{
    int n,s,e;
    double f;
    while(~scanf("%d",&n))
    {
        initialize();
        for(int i=1;i<=n;++i)
        {
            for(int j=1;j<=n;++j)
            {
                scanf("%lf",&f);
                add(i,j,f);
            }
        }
        scanf("%d",&q);
        while(q--)
        {
            scanf("%d%d",&s,&e);
            spfa(s,e);
        }
    }
    return 0;
}