題目描述
平面上有 \(n\) 個點,第 \(i\) 個點座標為 \((x_i, y_i)\)。連線 \(i, j\) 兩點的邊權為 \(\sqrt{(x_i - x_j) ^ 2 + (y_i - y_j) ^ 2}\)。求最小生成樹的邊權之和。
輸入格式
第一行一個整數 \(n\)。
接下來 \(n\) 行,每行輸入兩個整數 \(x_i, y_i\)。
輸出格式
輸出一行一個實數,表示答案。
當你的答案與標準輸出的絕對誤差或相對誤差在 \(10^{-6}\) 內時,就會被視為正確。
輸入輸出樣例
輸入 #1
4
0 0
1 2
-1 2
0 4
輸出 #1
6.472136
說明/提示
樣例解釋 1
該樣例中,最小生成樹如下圖所示:
邊權之和為 \(2 \sqrt{5} + 2 \approx 6.47213595500\)
資料規模與約定
對於 \(50\%\) 的資料,\(n \le 5000\)。
對於 \(100\%\) 的資料,\(3 \le n \le 10 ^ 5\),\(\lvert x_i \rvert, \lvert y_i \rvert \le 10 ^ 5\)。
分析
邊數太多,肯定不能用 \(Kruskal\)
\(n^2\) 的 \(Prim\) 也過不去
所以可以用 \(Boruvka\) 演算法
找到某一個點距離最近的點用 \(kdtree\) 查詢就行了
查詢的時候加一個剪枝,初始的答案不要置為無窮大,要設為當前聯通塊內的最優解
程式碼
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define rg register
inline int read(){
rg int x=0,fh=1;
rg char ch=getchar();
while(ch<'0' || ch>'9'){
if(ch=='-') fh=-1;
ch=getchar();
}
while(ch>='0' && ch<='9'){
x=(x<<1)+(x<<3)+(ch^48);
ch=getchar();
}
return x*fh;
}
const int maxn=1e5+5;
int cnt,orz,n,rt,x[maxn],y[maxn];
struct KDT{
int mn[2],mx[2],d[2],lc,rc,id,col;
friend bool operator <(const KDT& A,const KDT& B){
return A.d[orz]<B.d[orz];
}
}tr[maxn],jl[maxn];
void push_up(rg int da){
rg int lc=tr[da].lc,rc=tr[da].rc;
for(rg int i=0;i<2;i++){
tr[da].mn[i]=tr[da].mx[i]=tr[da].d[i];
if(lc){
tr[da].mx[i]=std::max(tr[da].mx[i],tr[lc].mx[i]);
tr[da].mn[i]=std::min(tr[da].mn[i],tr[lc].mn[i]);
}
if(rc){
tr[da].mx[i]=std::max(tr[da].mx[i],tr[rc].mx[i]);
tr[da].mn[i]=std::min(tr[da].mn[i],tr[rc].mn[i]);
}
}
}
int build(rg int l,rg int r,rg int pl){
orz=pl;
rg int mids=(l+r)>>1;
std::nth_element(jl+l,jl+mids,jl+r);
tr[mids]=jl[mids];
if(l<mids) tr[mids].lc=build(l,mids-1,!pl);
if(mids<r) tr[mids].rc=build(mids+1,r,!pl);
push_up(mids);
return mids;
}
long long getdis(rg int ax,rg int ay,rg int bx,rg int by){
return 1LL*(ax-bx)*(ax-bx)+1LL*(ay-by)*(ay-by);
}
long long sqr(rg int aa){
return 1LL*aa*aa;
}
long long mindis(rg int da,rg int xx,rg int yy){
rg long long mans=0;
mans+=sqr(std::max(0,xx-tr[da].mx[0])+std::max(0,tr[da].mn[0]-xx));
mans+=sqr(std::max(0,yy-tr[da].mx[1])+std::max(0,tr[da].mn[1]-yy));
return mans;
}
long long nans=0x3f3f3f3f3f3f3f3f;
int haha=0;
void cx(rg int da,rg int xx,rg int yy,rg int zz){
if(!da) return;
rg long long tmp=getdis(xx,yy,tr[da].d[0],tr[da].d[1]);
if(tr[da].col!=zz && nans>tmp){
nans=tmp,haha=tr[da].id;
}
tmp=mindis(da,xx,yy);
if(tmp>=nans) return;
rg long long tmp1=0x3f3f3f3f3f3f3f3f,tmp2=0x3f3f3f3f3f3f3f3f;
if(tr[da].lc) tmp1=mindis(tr[da].lc,xx,yy);
if(tr[da].rc) tmp2=mindis(tr[da].rc,xx,yy);
if(tmp1<tmp2){
if(tmp1<=nans) cx(tr[da].lc,xx,yy,zz);
if(tmp2<=nans) cx(tr[da].rc,xx,yy,zz);
} else {
if(tmp2<=nans) cx(tr[da].rc,xx,yy,zz);
if(tmp1<=nans) cx(tr[da].lc,xx,yy,zz);
}
}
double ans=0;
int bes[maxn],fa[maxn],tot;
long long bes2[maxn];
int zhao(rg int xx){
if(xx==fa[xx]) return xx;
return fa[xx]=zhao(fa[xx]);
}
double getdis2(rg int ax,rg int ay,rg int bx,rg int by){
return sqrt(1.0*(ax-bx)*(ax-bx)+1.0*(ay-by)*(ay-by));
}
int main(){
n=read();
for(rg int i=1;i<=n;i++) jl[i].d[0]=read(),jl[i].d[1]=read(),jl[i].id=jl[i].col=i,x[i]=jl[i].d[0],y[i]=jl[i].d[1];
for(rg int i=1;i<=n;i++) fa[i]=i;
rt=build(1,n,0);
rg int tmp=0;
while(tot<n-1){
for(rg int i=1;i<=n;i++) bes[i]=0,bes2[i]=0x3f3f3f3f3f3f3f3f;
for(rg int i=1;i<=n;i++){
tmp=zhao(i);
nans=bes2[tmp],haha=-1;
cx(rt,x[i],y[i],tmp);
if(haha==-1) continue;
if(bes2[tmp]>nans){
bes2[tmp]=nans;
bes[tmp]=haha;
} else if(bes2[tmp]==nans){
if(bes[tmp]<haha) bes[tmp]=haha;
}
}
for(rg int i=1;i<=n;i++){
tmp=zhao(i);
if(bes[tmp] && tmp!=zhao(bes[tmp])){
fa[tmp]=zhao(bes[tmp]);
ans+=sqrt((double)bes2[tmp]);
tot++;
}
}
for(rg int i=1;i<=n;i++){
tr[i].col=zhao(tr[i].id);
}
}
printf("%.6f\n",ans);
return 0;
}