航空公司（最小生成树）Prim算法

航空公司（最小生成树）Prim算法

题目描述

A国是一个拥有很多岛屿的国家。岛屿上的风景迷人且各有特点。不但如此，这些岛屿都是圆形的。所有的一切吸引了很多国内外的游客，A国当然不会放弃这个发展经济的好机会了。为了更好发展旅游业，A国决定由S航空公司来设计生产用于往返于各个岛屿之间的旅游小型机。
A国拥有大量的能源资源，但是缺乏航空技术。航空公司的工程师们在设计旅游小型机的时候总是希望能够减小油箱的大小而降低设计难度。由于每个岛屿上都有加油站，即使油箱容积小点，也可以通过到一些中转岛屿加油的办法往返于任意两个岛屿之间。A国拥有大量的能源，因而中转带来的资源浪费它们并不关心。已经知道，飞行1km需要1个单位的汽油。如果油箱过小，可能出导致飞行中由于汽油竭尽而坠入大海的危险。另外如果油箱过大，又会增加设计难度。因此，工程师们希望知道最小的油箱容积是多少？
岛屿的数目是那么多，工程师们一时都很难计算出：旅游小型机最小需要多大的油箱才能够保证可以往返于任意两个岛屿之间？聪明的你能否编写一个程序来帮助繁忙的工程师们呢。

输入

第一行：一个整数N (1<N<=1000)
接下来N行,每行包括3个整数(x,y,r).分别表示N个岛屿的坐标(x,y)及其它们的半径

输出

输出一个整数V （V是保证可以往返于任意两个岛屿之间的旅游小型机的最小油箱容积）。

样例输入 Copy

3
0 0 1
3 0 1
0 4 1

样例输出 Copy

2

提示

输入样例中A岛位于(0,0),半径为1;B岛位于(3,0),半径为1;C岛位于(0,4),半径为1.从A岛到达B岛的最近距离需要1个单位的汽油。从A岛到达C岛的最近距离需要2个单位的汽油。直接从B岛到达C岛需要3个单位的汽油，也就是说需要3个单位大小的油箱。但是，假如先从B岛飞到A岛，然后再从A岛飞到C岛，那么就只需要一个2个单位大小的油箱。所以，样例的正确输出应当是2。

思路：其实就是求一个图里，把它转化位最小生成树，最长边是多大。一棵树可以互达。

#pragma GCC optimize(3 , "Ofast" , "inline")
 
#include <bits/stdc++.h>
 
#define rep(i , a , b) for(register int i=(a);i<=(b);i++)
#define rop(i , a , b) for(register ll i=(a);i<(b);i++)
#define per(i , a , b) for(register ll i=(a);i>=(b);i--)
#define por(i , a , b) for(register ll i=(a);i>(b);i--)
 
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll , ll> pi;
const int maxn = 1e3+10;
const double eps = 1e-10;
int x[maxn],y[maxn],r[maxn];
int e[maxn][maxn];
int d[maxn];
bool v[maxn];
int n;
template<class T>
inline void read (T &x) {
    x = 0;
    int sign = 1;
    char c = getchar ();
    while (c < '0' || c > '9') {
        if ( c == '-' ) sign = - 1;
        c = getchar ();
    }
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar ();
    }
    x = x * sign;
}
void prime(){
    memset (d,0x3f, sizeof (d));
    memset (v,0, sizeof (v));
    d[1]=0;
    rep (i,1,n-1) {
        int x=0;
        rep(j,1,n) {
            if(!v[j]&&(x==0||d[j]<d[x])) x=j;
        }
        v[x]=1;
        rep (y,1,n) {
            if(!v[y]) d[y]=min (d[y],e[x][y]);
        }
    }
}
 
int main () {
    read (n);
    rep (i,1,n) {
        read (x[i]);read (y[i]);read (r[i]);
    }
    double dist = 0;
    rep (i,1,n) {
        rep (j,i+1,n) {
            int d = (x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]);
            double m = sqrt (d*1.0);
            m=m-r[i]-r[j];
            if(m<0) m=0;
            e[i][j]=ceil (m);
            e[j][i]=ceil (m);
        }
    }
    prime();
    int ans = 0;
    rep (i,1,n) {
        ans=max (ans,d[i]);
    }
    cout<<ans<<endl;
    return 0;
}