设为首页 加入收藏

TOP

UVA 6199 不定根最小树形图
2014-11-23 20:10:35 来源: 作者: 【大 中 小】 浏览:4次
Tags：UVA 6199 不定 最小 树形

首先是最小树形图的介绍。

看这个博客。最小树形图

上面介绍的很详细了，我就讲一下这道题的题意。

首先给出一些二维点坐标，这些坐标之间构成一些有向图，根据题意，假设两个点a(x1 ,y1) ,b(x2 ,y2) .当y1 <= y2时，他们之间可以连一条有向边，即a -> b。

就是每个点只能连y坐标大于他的点，然后就构成了一张有向图。

最后求出最少的距离可以使得所有的点都连起来。

刚开始以为直接求出两两之间的距离，然后用kruskal求一遍MST就可以了。但是仔细想了一下，这里有向边的限制就使得一些连边的情况是不可行的。

这道题的正解是最小树形图，而且是最裸的。

因为这道题他的根是不确定的，那么我们可以用一个超级源点，作为他的根，将他和所有的点都连起来，边是inf。

#include    
#include    
#include    
#include    
#include    
#include    
#include    
#include    
#include    
#include    
#include    
#include    
#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   
#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 ;  
        //找到每个点的最小入边   
        for (int i = 0 ; i < NE ; i ++ ){  
            int s = ed[i].s ;  
            int e = ed[i].e ;  
            if(ed[i].l < in[e] && s != e){  
                pre[e] = s ;  
                in[e] = ed[i].l ;  
            }  
        }  
        //最小入边   
//        for (int i = 1 ; i < NV ; i ++ ){   
//            cout << i << " : " << in[i] << endl;   
//        }   
        for (int i = 0 ; i < NV ; i ++ ){//除根节点外所有点都找到一条入边   
            if(i == root)continue ;  
            if(in[i] == inf)return -1 ;  
        }  
        //找环   
        int cntnode = 0 ;  
        mem(vis ,-1) ;  
        mem(id ,-1) ;  
        in[root] = 0 ;  
        for (int i = 0 ; i < NV ; i ++ ){  
            ret += in[i] ;  
            int v = i ;  
            while(vis[v] != i && id[v] == -1 && v != root){  
                vis[v] = i ;  
                v = pre[v] ;  
            }  
            if(v != root && id[v] == -1){  
                for (int u = pre[v] ; u != v ; u = pre[u]){  
                    id[u] = cntnode ;  
                }  
                id[v] = cntnode ++ ;  
            }  
        }  
  
        if(cntnode == 0)break ;//无环   
        for (int i = 0 ; i < NV ; i ++ ){  
            if(id[i] == -1)id[i] = cntnode ++ ;  
        }  
        //缩点   
        for (int i = 0 ; i < NE ; i ++ ){  
            int s = ed[i].s ;  
            int e = ed[i].e ;  
            ed[i].s = id[s] ;  
            ed[i].e = id[e] ;  
            if(ed[i].s != ed[i].e){  
                ed[i].l -= in[e] ;  
            }  
        }  
        NV = cntnode ;  
        root = id[root] ;  
    }  
    return ret ;  
}  
int main() {  
    while(cin >> n , n ){  
        int a , b ;  
        for (int i = 1 ; i <= n ; i ++ ){  
            scanf("%lf %lf",&P[i].x ,&P[i].y) ;  
        }  
        init() ;  
        double dis_sum = 0 ;  
        for (int i = 1 ; i <= n ;i ++ ){  
            for (int j = 1 ; j <= n ;j ++ ){  
                if(i == j)continue ;  
                double dis = getdis(i , j) ;  
                if(P[i].y <= P[j].y){  
                    add(i , j , dis) ;  
//                    cout << dis << endl;   
                    dis_sum += dis ;  
                }  
            }  
        }  
        S = 0 ;  
        for (int i = 1 ; i <= n ; i ++ ){  
            add(S , i , inf - 1 ) ;  
        }  
        printf("%.2f\n",Directed_MST(0 , n + 1 , num ) - inf + 1) ;  
    }  
    return 0 ;  
}  

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#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
#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 ;
        //找到每个点的最小入边
        for (int i = 0 ; i < NE ; i ++ ){
            int s = ed[i].s ;
            int e = ed[i].e ;
            if(ed[i].l < in[e] && s != e){
                pre[e] = s ;
                in[e] = ed[i].l ;
            }
        }
        //最小入边
//        for (int i = 1 ; i < NV ; i ++ ){
//            cout << i << " : " << in[i] << endl;
//        }
        for (int i = 0 ; i < NV ; i ++ ){//除根节点外所有点都找到一条入边
            if(i == root)continue ;
            if(in[i] == inf)return -1 ;
        }
        //找环
        int cntnode = 0 ;
        mem(vis ,-1) ;
        mem(id ,-1) ;
        in[root] = 0 ;
        for (int i = 0 ; i < NV ; i ++ ){
            ret += in[i] ;
            int v = i ;
            while(vis[v] != i && id[v] == -1 && v != root){
                vis[v] = i ;
                v = pre[v] ;
            }
            if(v != root && id[v] == -1){
                for (int u = pre[v] ; u != v ; u = pre[u]){
                    id[u] = cntnode ;
                }
                id[v] = cntnode ++ ;
            }
        }

        if(cntnode == 0)break ;//无环
        for (int i = 0 ; i < NV ; i ++ ){
            if(id[i] == -1)id[i] = cntnode ++ ;
        }
        //缩点
        for (int i = 0 ; i < NE ; i ++ ){
            int s = ed[i].s ;
            int e = ed[i].e ;
            ed[i].s = id[s] ;
            ed[i].e = id[e] ;
            if(ed[i].s != ed[i].e){
                ed[i].l -= in[e] ;
            }
        }
        NV = cntnode ;
        root = id[root] ;
    }
    return ret ;
}
int main() {
    while(cin >> n , n ){
        int a , b ;
        for (int i = 1 ; i <= n ; i ++ ){
            scanf("%lf %lf",&P[i].x ,&P[i].y) ;
        }
        init() ;
        double dis_sum = 0 ;
        for (int i = 1 ; i <= n ;i ++ ){
            for (int j = 1 ; j <= n ;j ++ ){
                if(i == j)continue ;
                double dis = getdis(i , j) ;
                if(P[i].y <= P[j].y){
                    add(i , j , dis) ;
//                    cout << dis << endl;
                    dis_sum += dis ;
                }
            }
        }
        S = 0 ;
        for (int i = 1 ; i <= n ; i ++ ){
            add(S , i , inf - 1 ) ;
        }
        printf("%.2f\n",Directed_MST(0 , n + 1 , num ) - inf + 1) ;
    }
    return 0 ;
}

【大 中 小】【打印】 【繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部】
分享到: 
上一篇：HDU 1829 A Bug's Life (并查.. 下一篇：UVA 10602 Editor Nottoobad(贪心..

评论

帐　　号: 密码: (新用户注册)
验 证 码:
表　　情:
内　　容:

·微服务 Spring Boot (2025-12-26 18:20:10)
·如何调整 Redis 内存 (2025-12-26 18:20:07)
·MySQL 数据类型：从 (2025-12-26 18:20:03)
·Linux Shell脚本教程 (2025-12-26 17:51:10)
·Qt教程，Qt5编程入门 (2025-12-26 17:51:07)