首先是最小树形图的介绍。
    看这个博客。最小树形图
    上面介绍的很详细了，我就讲一下这道题的题意。
    首先给出一些二维点坐标，这些坐标之间构成一些有向图，根据题意，假设两个点a（x1 ,y1） ,b（x2 ,y2） .当y1 <= y2时，他们之间可以连一条有向边，即a -> b.
    就是每个点只能连y坐标大于他的点，然后就构成了一张有向图。
    最后求出最少的距离可以使得所有的点都连起来。
    刚开始以为直接求出两两之间的距离，然后用kruskal求一遍MST就可以了。但是仔细想了一下，这里有向边的限制就使得一些连边的情况是不可行的。
    这道题的正解是最小树形图，而且是最裸的。
    因为这道题他的根是不确定的，那么我们可以用一个超级源点，作为他的根，将他和所有的点都连起来，边是inf.
    #include <set>
    #include <map>
    #include <stack>
    #include <cmath>
    #include <queue>
    #include <cstdio>
    #include <string>
    #include <vector>
    #include <iomanip>
    #include <cstring>
    #include <iostream>
    #include <algorithm>
    #define Max 2505
    #define FI first
    #define SE second
    #define ll long long
    #define PI acos（-1.0）
    #define inf 0x3fffffff
    #define LL（x） （ x 《 1 ）
    #define bug puts（"here"）
    #define PII pair<int,int>
    #define RR（x） （ x 《 1 | 1 ）
    #define mp（a,b） make_pair（a,b）
    #define mem（a,b） memset（a,b,sizeof（a））
    #define REP（i,s,t） for（ int i = （ s ） ; i <= （ t ） ; ++ i ）
    using namespace std;
    inline void RD（int &ret） {
    char c;
    int flag = 1 ;
    do {
    c = getchar（）；
    if（c == '-'）flag = -1 ;
    } while（c < '0' || c > '9'） ;
    ret = c - '0';
    while（（c=getchar（）） >= '0' && c <= '9'）
    ret = ret * 10 + （ c - '0' ）；
    ret *= flag ;
    }
    inline void OT（int a） {
    if（a >= 10）OT（a / 10） ;
    putchar（a % 10 + '0'） ;
    }
    inline void RD（double &ret） {
    char c ;
    int flag = 1 ;
    do {
    c = getchar（） ;
    if（c == '-'）flag = -1 ;
    } while（c < '0' || c > '9'） ;
    ll n1 = c - '0' ;
    while（（c = getchar（）） >= '0' && c <= '9'） {
    n1 = n1 * 10 + c - '0' ;
    }
    ll n2 = 1 ;
    while（（c = getchar（）） >= '0' && c <= '9'） {
    n1 = n1 * 10 + c - '0' ;
    n2 *= 10 ;
    }
    ret = flag * （double）n1 / （double）（n2） ;
    }
    /*********************************************/
    #define N 1005
    struct PP{
    double x , y ;
    }P[N] ;
    double getdis（int i ,int j）{
    return sqrt（（P[i].x - P[j].x） * （P[i].x - P[j].x） + （P[i].y - P[j].y） * （P[i].y - P[j].y）） ;
    }
    struct kdq{
    int s , e ;
    double l ;
    }ed[N * N] ;
    int num ;
    void add（int s ,int e ,double l）{
    ed[num].s = s ;
    ed[num].e = e ;
    ed[num].l = l ;
    num ++ ;
    }
    void init（）{
    num = 0 ;
    }
    int n ;
    int S ;
    int pre[N] , id[N] , vis[N] ;
    double in[N] ;
    double Directed_MST（int root ,int NV , int NE）{
    double ret = 0 ;
    while（1）{
    for （int i = 0 ; i < NV ; i ++ ）in[i] = inf ;