poj 2253 求可达路上的最大边，dijkstra，spfa

最新推荐文章于 2019-09-05 20:11:00 发布

原创最新推荐文章于 2019-09-05 20:11:00 发布 · 267 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#path #distance #struct #优化

图专栏收录该内容

12 篇文章

订阅专栏

本文详细介绍了两种经典的最短路径算法：Dijkstra算法和SPFA算法。通过具体实例讲解了如何利用这两种算法寻找两点间最大权重边的路径，并提供了完整的C++代码实现。

/*

http://acm.pku.edu.cn/JudgeOnline/problem?id=2253

求可达路上的权值最大的边

在dijkstra中，d[i]表示源点s到i的最短路径，这里根据题目要求修改d[i]为“未访问”节点和“已访问”节点，两个点集之间的最短路

则松弛条件为：d[i]=min(d[i],path[mark][i]);

在算浮点数时，可将所有的路径都用int来表示，最后才开方。

spfa是bellman_ford的优化，省去很多冗余的操作

期望的时间复杂度O(ke)， 其中k为所有顶点进队的平均次数，可以证明k一般小于等于2。

实现方法：建立一个队列，初始时队列里只有起始点

建立一个path记录起始点到所有点的最短路径（该表格的初始值要赋为极大值，该点到他本身的路径赋为0）。

然后执行松弛操作，用队列里有的点去刷新起始点到所有点的最短路，如果刷新成功且被刷新点不在队列中则把该点加入到队列最后。

重复执行直到队列为空

判断有无负环：如果某个点进入队列的次数超过N次则存在负环（SPFA无法处理带负环的图）

*/

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<vector>
#include<set>
#include<queue>
#include<map>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<functional>
using namespace std;
const int INF=(1<<31)-1;
const double INF_DOUBLE=1e8;
const double EPS=1e-8;
const int N=210; 
int n;
int path[N][N];
struct point
{
  int x,y;     
}p[N];

int calDis(int i,int j)
{
   return (p[i].x-p[j].x)*(p[i].x-p[j].x)+(p[i].y-p[j].y)*(p[i].y-p[j].y);      
}

double dijkstra(int s,int t)
{
       int d[N];
       bool v[N];
       memset(v,false,sizeof(v));
       for(int i=1;i<=n;i++)d[i]=path[s][i];
       v[s]=true;
       int ans=-1;
       while(true)
       {
         int Min=INF;
         int mark=-1;
         for(int i=1;i<=n;i++)
          if(!v[i]&&d[i]<Min){Min=d[i];mark=i;}
         if(mark==-1)return ans;
         v[mark]=true;
         ans=max(ans,Min);
         if(mark==t)return ans; 
         for(int i=1;i<=n;i++)
          if(!v[i]&&d[i]>path[mark][i])
           d[i]=path[mark][i];                   
       }      
}

double spfa(int s,int t)
{
    int d[N];
    bool in_q[N];
    for(int i=1;i<=n;i++)d[i]=INF;
    memset(in_q,false,sizeof(in_q));
    d[s]=0;
    queue<int>q;
    q.push(s);    
    in_q[s]=true;
    while(!q.empty())
    {
       int temp=q.front(); 
       q.pop();
       in_q[temp]=false;
       for(int i=1;i<=n;i++)
        if(d[i]>max(d[temp],path[temp][i]))
        {
          d[i]=max(d[temp],path[temp][i]);
          if(!in_q[i])
          {
             q.push(i);
             in_q[i]=true;
          }
        }                
    }
    return d[t];
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int T=1;
    while(cin>>n,n)
    {
       int a,b;
       for(int i=1;i<=n;i++){scanf("%d %d",&p[i].x,&p[i].y);}
       for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
         path[i][j]=calDis(i,j);
       printf("Scenario #%d/nFrog Distance = %.3lf/n/n",T++,sqrt(dijkstra(1,2)));                       
    }
    //system("pause");
    return 0;
}