Floyd - 牛的旅行 - 洛谷 P1522

最新推荐文章于 2023-08-03 20:43:36 发布

njuptACMcxk

最新推荐文章于 2023-08-03 20:43:36 发布

阅读量153

点赞数

分类专栏：最短路文章标签：算法图论 Floyd算法 Floyd算法 ACM 程序设计竞赛题目解读最短路

本文链接：https://blog.csdn.net/njuptACMcxk/article/details/107084178

版权

最短路专栏收录该内容

48 篇文章 0 订阅

订阅专栏

Floyd - 牛的旅行 - 洛谷 P1522

题意：

$有多个牧区 (连通块) ，每个牧区的直径为牧区中相距最远的两点的距离。$

$这里的距离都是指的两点间的最短距离。$

$给定一个 n 个点的图，可能有多个连通块。$

$要在某两个连通块之间添加一条边，使得两个连通块合并为一个连通块。$

$要求合并之后，所有连通块中，最大直径的最小值。$

输入格式

第 1 行：一个整数 N, 表示牧区数；

第 2 到 N+1 行：每行两个整数 X,Y，表示 N 个牧区的坐标。每个牧区的坐标都是不一样的。

第 N+2 行到第 2*N+1 行：每行包括 N 个数字 ( 0或1 ) 表示一个对称邻接矩阵。

例如，题目描述中的两个牧场的矩阵描述如下：

  A B C D E F G H 
A 0 1 0 0 0 0 0 0 
B 1 0 1 1 1 0 0 0 
C 0 1 0 0 1 0 0 0 
D 0 1 0 0 1 0 0 0 
E 0 1 1 1 0 0 0 0 
F 0 0 0 0 0 0 1 0 
G 0 0 0 0 0 1 0 1 
H 0 0 0 0 0 0 1 0

输入数据中至少包括两个不连通的牧区。

输出格式

只有一行，包括一个实数，表示所求答案。

数字保留六位小数。

数据范围

$1≤N≤150, 0≤X,Y≤10^5$

输入样例：

输出样例：

22.071068

分析：

$对于每个连通块，计算直径都得先计算连通块内各点间的最短距离。$

$所以，建图后，在每个连通块内部跑一遍 f l o y d 。$

$d i s [i] [j] 保存点 i 和点 j 之间的最短距离。$

$要求所有连通块的最大直径。有两种情况：$

$① 、合并后的连通块的直径不是最大的：$

$\qquad最大的直径在原先未合并的几个连通块中取。$

$\qquad用res1来更新合并前所有连通块中的最大直径，即max_{1≤i,j≤n}(dis[i][j])。$

$② 、合并后的连通块的直径是最大的：$

$\qquad假设在点i和点j之间连一条边，距离为dis[i][j]，max[i]表示i所在连通块中，到i点的最大距离。$

$\qquad则合并后的连通块的直径为max[i]+dis[i][j]+max[j]。$

$\qquad因此，我们还要预处理一个数组max[i]。$

$\qquad由于我们要使得最大直径最小，因为情况①是固定的，所以我们要使得情况②得到的新的直径尽量最小。$

$\qquad用res2来更新合并后的直径的最小值，即min_{1≤i,j≤n}(max[i]+dis[i][j]+max[j])。$

$最后在 r e s 1 和 r e s 2 中取一个较大值即可。$

注意：

$跑过 f l o y d 后会建立新的连通关系，若 d i s [i] [j] < i n f 表示点 i 和 j 连通，否则表示不连通。$

$\%lf默认小数点后保留6位。$

代码：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>

#define P pair<int,int>
#define x first
#define y second

using namespace std;

const int N=160;
const double inf=1e20;

int n;
P V[N];
char g[N][N];
double dis[N][N],maxd[N];

double get_dist(P a,P b)
{
    double dx=a.x-b.x,dy=a.y-b.y;
    return sqrt(dx*dx+dy*dy);
}

void floyd()
{
    for(int k=0;k<n;k++)
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
}

int main()
{
    cin>>n;
    for(int i=0;i<n;i++) cin>>V[i].x>>V[i].y;
    
    for(int i=0;i<n;i++) cin>>g[i];
    
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            if(i!=j)
            {
                if(g[i][j]=='1') dis[i][j]=get_dist(V[i],V[j]);
                else dis[i][j]=inf;
            }
            
    floyd();
    
    double res1=0;
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
            if(dis[i][j]<inf)   //连通不能通过g矩阵，因为跑过floyd后可能间接连通
                maxd[i]=max(maxd[i],dis[i][j]);
        res1=max(res1,maxd[i]);
    }
            
    double res2=inf;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            if(dis[i][j]>=inf)
                res2=min(res2,maxd[i]+get_dist(V[i],V[j])+maxd[j]);
                
    printf("%lf\n",max(res1,res2));
    
    return 0;
}