分析发现:连边前后 两个联通块本身的直径不被影响
题意所求可以转化为:第一个连通块最大直径 加连边距离 加第二个连通块最大直径
思路如下:
1.Floyd算法分别求出两个牧场内任意两点的距离
2.对两个牧场分别预处理maxd[i] 意为在该牧场内 与i点联通且距离i点最远的距离
3.枚举所有情况求(maxd1[i]+maxd[2]+两个牧场连边权值)最小值
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
#include<string>
#include<cmath>
#define x first
#define y second
using namespace std;
typedef pair<double, double>PII;
const int N = 200;
const double INF = 1e30;
double dis[N][N], maxd[N];
PII q[N];//存坐标
char g[N][N];
int n;//牧区数目(点的个数)
//1.Floyd算法分别求出两个牧场内任意两点的距离
//2.对两个牧场分别预处理maxd[i] 意为在该牧场内 与i点联通且距离i点最远的距离
//3.枚举所有情况求(maxd1[i]+maxd[2]+两个牧场连边权值)最小值
double get_dist(PII a, PII b)
{
double dx = a.x - b.x;
double dy = a.y - b.y;
return sqrt(dx * dx + dy * dy);
}
int main()
{
cin >> n;
//读入所有点的坐标
for (int i = 0; i < n; i++)cin >> q[i].x >> q[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)//没有自环
{
//从邻接矩阵分析
//如果i,j点之间有边 就计算出来
//如果没有边 就设置为正无穷
if (g[i][j] == '0')dis[i][j] = INF;
else dis[i][j] = get_dist(q[i], q[j]);
}
//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]);
//预处理在原有牧场中 每一个点可以联通到的最远点距离
for (int i = 0; i < n; i++)
for (int j = 0; j < n; j++)
if (dis[i][j] < INF)
maxd[i] = max(maxd[i], dis[i][j]);
//求原牧场最大直径
double res1 = 0;
for (int i = 0; i < n; i++)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, get_dist(q[i], q[j]) + maxd[i] + maxd[j]);
}
printf("%lf", max(res1, res2));
return 0;
}