The Fourth Dimension Space

枯叶北风寒,忽然年以残,念往昔,语默心酸。二十光阴无一物,韶光贱,寐难安; 不畏形影单,道途阻且慢,哪曲折,如渡飞湍。斩浪劈波酬壮志,同把酒,共言欢! -如梦令

pku 3164 最小树形图(Zhu-Liu Algorithm)

所谓的最小树形图,就是让你在一个有向图中找一个根和由根扩展出来的一颗有向的生成树,并且使这棵树的权值最小。
令人欣喜的是,这个算法是由中国人提出来的,这也是我正式学习的第一个中国人提出来的算法^_^

最小树形图算法
(Zhu-Liu Algorithm)

1.       设最小树形图的总权值为cost,置cost0

2.       除源点外,为其他所有节点Vi找一条权值最小的入边,加入集合TT就是最短边的集合。加边的方法:遍历所有点到Vi的边中权值最小的加入集合T,记pre[Vi]为该边的起点,mincost[Vi]为该边的权值。

3.       检查集合T中的边是否存在有向环,有则转到步骤4,无则转到步骤5。这里需要利用pre数组,枚举检查过的点作为搜索的起点,类似dfs的操作判断有向环。

4.       将有向环缩成一个点。设环中有点{Vk1,Vk2,…,Vki}i个点,用Vk代替缩成的点。在压缩后的图中,更新所有不在环中的点VVk的距离:

map[V][Vk] = min {map[V][Vkj]-mincost[Vki]} 1<=j<=i

map[Vk][V] = min {map[Vkj][V]}           1<=j<=I

5.       cost加上T中有向边的权值总和就是最小树形图的权值总和。

源代码如下:(参考了网上的代码)
#include<iostream>
#include
<cstring>
#include
<cstdio>
#include
<cmath>
using namespace std;

#define INF 99999999
#define min( a, b ) ( (a)< (b)?(a): (b) )

struct point
{

    
double x;
    
double y;
}
p[200];

int pre[200];//记录该节点的前驱
double graph[200][200], ans;//图数组和结果
bool   visit[110], circle[110];//visit记录该点有没有被访问过,circle记录改点是不是在一个圈里
int    n, m, root;//顶点数+边数+根节点标号

void dfs( int t )//一个深度优先搜索,搜索出一个最大的联通空间
{
    
int i;
    visit[t]
= true;
    
for(i= 1; i<= n; ++i )
    
{
        
if!visit[i] && graph[t][i]!= INF )
            dfs( i );
    }

}


bool check()//这个函数用来检查最小树形图是否存在,即如果存在,那么一遍dfs后,应该可以遍历到所有的节点
{
    memset( visit, 
falsesizeof(visit) );
    dfs( root );
    
    
forint i= 1; i<= n; ++i )
    
{
        
if!visit[i] ) 
            
return false;
    }

    
return true;
}


double dist( 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) );
}


int exist_circle()//判断图中是不是存在有向圈
{
    
int i;
    
int j;
    root
= 1; pre[root]= root;
    
for(i= 1; i<= n; ++i )
    
{
        
if!circle[i] && i!= root )
        
{
            pre[i]
= i; graph[i][i]= INF;
            
            
for(j= 1; j<= n; ++j )
            
{
                
if!circle[j] && graph[j][i]< graph[pre[i]][i] )
                    pre[i]
= j;
            }

        }

    }
//这个for循环负责找出所有非根节点的前驱节点
    for( i= 1; i<= n; ++i )
    
{
        
if( circle[i] ) 
            
continue;
        memset( visit, 
falsesizeof(visit) );
        
int j= i;
        
while!visit[j] ) 
        

            visit[j]
= true
            j
= pre[j]; 
        }

        
if( j== root )
            
continue;
        
        
return j;
    }
//找圈过程,最后返回值是圈中的一个点
    
    
return -1;//如果没有圈,返回-1
}



void  update( int t )//缩圈之后更新数据
{
    
int i;
    
int j;
    ans
+= graph[pre[t]][t];
    
for(i=pre[t]; i!= t; i= pre[i] )
    
{
        ans
+= graph[pre[i]][i];
        circle[i]
= true;
    }
//首先把圈里的边权全部加起来,并且留出t节点,作为外部接口
    
    
for(i= 1; i<= n; ++i )
        
if!circle[i] && graph[i][t]!= INF )
            graph[i][t]
-= graph[pre[t]][t];
    
//上面这个for循环的作用是对t节点做更新操作,为什么要单独做?你可以看看线面这个循环的跳出条件。
    for(j= pre[t]; j!= t; j= pre[j] )
        
forint i= 1; i<= n; ++i )
        
{
            
if( circle[i] ) 
                
continue;
            
if( graph[i][j]!= INF )
            graph[i][t]
= min( graph[i][t], graph[i][j]- graph[pre[j]][j] );
            
//////////////////////////////////////////////////////////////////////////
            graph[t][i]= min( graph[j][i], graph[t][i] );
        }

    
//这个循环对圈中的其他顶点进行更新
}


void solve()
{
    
int j;
    memset( circle, 
falsesizeof(circle) );
    
while( ( j= exist_circle() )!= -1 ) 
        update( j );
    
    
for( j= 1; j<= n; ++j )
        
if( j!= root && !circle[j] )
            ans
+= graph[pre[j]][j];
    
    printf(
"%.2lf\n", ans );
}


int main()
{
    
int i;
    
while( scanf("%d%d",&n,&m)!= EOF )
    
{
        
for(i= 0; i<= n; ++i )
            
forint j= 0; j<= n; ++j )
                graph[i][j]
= INF;
        
        
for(i= 1; i<= n; ++i )
            scanf(
"%lf%lf",&p[i].x, &p[i].y );
        
        
for(i= 0; i< m; ++i )
        
{
            
int a, b;
            scanf(
"%d%d",&a,&b);
            graph[a][b]
= dist( a, b );
        }

        
        root
= 1
        ans
= 0;
        
if!check() ) 
            printf(
"poor snoopy\n");
        
else  
            solve();
    }

    
    
return 0;
}



posted on 2009-10-29 19:38 abilitytao 阅读(2057) 评论(2)  编辑 收藏 引用

评论

# re: pku 3164 最小树形图(Zhu-Liu Algorithm) 2009-10-30 11:58 淘宝网首页

学习!学习~~  回复  更多评论   

# re: pku 3164 最小树形图(Zhu-Liu Algorithm) 2009-11-03 16:31 argmax

这个算法应该叫Chu-Liu-Edmond算法吧!  回复  更多评论   


只有注册用户登录后才能发表评论。
【推荐】超50万行VC++源码: 大型组态工控、电力仿真CAD与GIS源码库
网站导航: 博客园   IT新闻   BlogJava   知识库   博问   管理