最小树形图（tju 2248 UVA 11183 poj 3164）

于 2015-05-14 19:55:29 发布
阅读量567 收藏
点赞数
分类专栏： 图论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u011699990/article/details/45725383
版权

求最小树形图的总权值  

即以固定跟为起点 延给定有向边 可以访问所有的点 并所构成的边权值之和最小 求出这个最小总权值


算法步骤：
① 清除自环，输入的时候判断即可
② 先判断从固定根开始是否可达所有原图中的点。简单搜索加标记位就可以。如果不可就不用说了，肯定没戏。
③ 为除根之外的每个点选定一条最小入边。
（记pre [vi]为该边的起点）
④ 判断这个入边集是否存在有向环，如果不存在，我们很容易证明这个集合就是该图的最小树形图，转⑥，否则接⑤（利用prev数组，枚举为检查过的点作为搜索的起点，做类似DFS的操作）
⑤ 消环。设(u,i,w)表示从u到i的权为w的边。设刚才的有向环缩为新结点new。若u位于环上，并设环中指向u的边权是in[u]。那么对于每条从u出发的边(u, i, w)，在新图中连接(new, i, w)的边，其中new为新加的人工顶点; 对于每条进入u的边(i, u, w)，在新图中建立边(i, new, w-in[u])的边。新图中最小树形图的权加上旧图中被收缩的那个环的权和，就是原图中最小树形图的权。重复③④⑤
⑥ 成功，返回DMST总权值。
补充1：如果无固定根，增加一个节点，连接到所有节点，并且距离一样，即可转化为有固定根。
补充2：本算法只能求最小总权值，但不能求路径。


tju 2248 http://acm.tju.edu.cn/toj/showp2248.html

裸的最小树形图

<span style="font-family:KaiTi_GB2312;font-size:18px;">#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string.h>
#include <string>
#include <vector>
#include <queue>

#define MEM(a,x) memset(a,x,sizeof a)
#define eps 1e-8
#define MOD 10009
#define MAXN 1010
#define MAXM 42010
#define INF 99999999
#define ll __int64
#define bug cout<<"here"<<endl
#define fread freopen("ceshi.txt","r",stdin)
#define fwrite freopen("out.txt","w",stdout)

using namespace std;

int Read()
{
    char c = getchar();
    while (c < '0' || c > '9') c = getchar();
    int x = 0;
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x;
}

void Print(int a)
{
     if(a>9)
         Print(a/10);
     putchar(a%10+'0');
}

struct Edge
{
    int u,v;
    int cost;
}edge[MAXM];
int pre[MAXN],vis[MAXN],id[MAXN];
int in[MAXN];

int zhuliu(int root,int n,int m,Edge edge[])
{
    int res=0;
    int u,v;
    while(1)
    {
        for(int i=0;i<n;i++)
            in[i]=INF;
        for(int i=0;i<m;i++)
            if(edge[i].u!=edge[i].v&&(edge[i].cost-in[edge[i].v]<eps))
        {
            pre[edge[i].v]=edge[i].u;
            in[edge[i].v]=edge[i].cost;
        }
        for(int i=0;i<n;i++)
            if(i!=root&&in[i]==INF)
                return -1;
        int tn=0;
        MEM(id,-1); MEM(vis,-1);
        in[root]=0;
        for(int i=0;i<n;i++)
        {
            res+=in[i];
            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]=tn;
                id[v]=tn++;
            }
        }
        if(tn==0) break;//没有有向环
        for(int i=0;i<n;i++)
            if(id[i]==-1)
                id[i]=tn++;
        for(int i=0;i<m;)
        {
            v=edge[i].v;
            edge[i].u=id[edge[i].u];
            edge[i].v=id[edge[i].v];
            if(edge[i].u!=edge[i].v)
                edge[i++].cost-=in[v];
            else
                swap(edge[i],edge[--m]);
        }
        n=tn;
        root=id[root];
    }
    return res;
}
int g[MAXN][MAXN];
int x[MAXN],y[MAXN];



int main()
{
//    fread;
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(n==0&&m==0) break;
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                g[i][j]=INF;
        int u,v,cost;
        while(m--)
        {
            scanf("%d%d%d",&u,&v,&cost);
            u--; v--;
            if(u==v) continue;
            g[u][v]=min(g[u][v],cost);
        }
        int num=0;
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                if(g[i][j]<INF)
                {
                    edge[num].u=i;
                    edge[num].v=j;
                    edge[num++].cost=g[i][j];
                }
        int ans=zhuliu(0,n,num,edge);
        if(ans==-1) printf("impossible\n");
        else printf("%d\n",ans);
    }
    return 0;
}
</span>

UVA 11183 http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=18548 

<span style="font-family:KaiTi_GB2312;font-size:18px;">#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string.h>
#include <string>
#include <vector>
#include <queue>

#define MEM(a,x) memset(a,x,sizeof a)
#define eps 1e-8
#define MOD 10009
#define MAXN 1010
#define MAXM 42010
#define INF 99999999
#define ll __int64
#define bug cout<<"here"<<endl
#define fread freopen("ceshi.txt","r",stdin)
#define fwrite freopen("out.txt","w",stdout)

using namespace std;

int Read()
{
    char c = getchar();
    while (c < '0' || c > '9') c = getchar();
    int x = 0;
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x;
}

void Print(int a)
{
     if(a>9)
         Print(a/10);
     putchar(a%10+'0');
}

struct Edge
{
    int u,v;
    int cost;
}edge[MAXM];
int pre[MAXN],vis[MAXN],id[MAXN];
int in[MAXN];

int zhuliu(int root,int n,int m,Edge edge[])
{
    int res=0;
    int u,v;
    while(1)
    {
        for(int i=0;i<n;i++)
            in[i]=INF;
        for(int i=0;i<m;i++)
            if(edge[i].u!=edge[i].v&&(edge[i].cost-in[edge[i].v]<eps))
        {
            pre[edge[i].v]=edge[i].u;
            in[edge[i].v]=edge[i].cost;
        }
        for(int i=0;i<n;i++)
            if(i!=root&&in[i]==INF)
                return -1;
        int tn=0;
        MEM(id,-1); MEM(vis,-1);
        in[root]=0;
        for(int i=0;i<n;i++)
        {
            res+=in[i];
            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]=tn;
                id[v]=tn++;
            }
        }
        if(tn==0) break;//没有有向环
        for(int i=0;i<n;i++)
            if(id[i]==-1)
                id[i]=tn++;
        for(int i=0;i<m;)
        {
            v=edge[i].v;
            edge[i].u=id[edge[i].u];
            edge[i].v=id[edge[i].v];
            if(edge[i].u!=edge[i].v)
                edge[i++].cost-=in[v];
            else
                swap(edge[i],edge[--m]);
        }
        n=tn;
        root=id[root];
    }
    return res;
}
int g[MAXN][MAXN];
int x[MAXN],y[MAXN];



int main()
{
//    fread;
    int tc;
    int cs=1;
    scanf("%d",&tc);
    while(tc--)
    {
        int n,m;
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                g[i][j]=INF;
        int u,v,cost;
        while(m--)
        {
            scanf("%d%d%d",&u,&v,&cost);
            if(u==v) continue;
            g[u][v]=min(g[u][v],cost);
        }
        int num=0;
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                if(g[i][j]<INF)
                {
                    edge[num].u=i;
                    edge[num].v=j;
                    edge[num++].cost=g[i][j];
                }
        int ans=zhuliu(0,n,num,edge);
        printf("Case #%d: ",cs++);
        if(ans==-1) printf("Possums!\n");
        else printf("%d\n",ans);
    }
    return 0;
}
</span>

poj 3164 http://poj.org/problem?id=3164

此题中cost是两个点之间的距离 

<span style="font-family:KaiTi_GB2312;font-size:18px;">#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string.h>
#include <string>
#include <vector>
#include <queue>

#define MEM(a,x) memset(a,x,sizeof a)
#define eps 1e-8
#define MOD 10009
#define MAXN 1010
#define MAXM 12010
#define INF 99999999
#define ll __int64
#define bug cout<<"here"<<endl
#define fread freopen("ceshi.txt","r",stdin)
#define fwrite freopen("out.txt","w",stdout)

using namespace std;

int Read()
{
    char c = getchar();
    while (c < '0' || c > '9') c = getchar();
    int x = 0;
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x;
}

void Print(int a)
{
     if(a>9)
         Print(a/10);
     putchar(a%10+'0');
}

struct Edge
{
    int u,v;
    double cost;
}edge[MAXM];
int pre[MAXN],vis[MAXN],id[MAXN];
double in[MAXN];

double zhuliu(int root,int n,int m,Edge edge[])
{
    double res=0;
    int u,v;
    while(1)
    {
        for(int i=0;i<n;i++)
            in[i]=INF;
        for(int i=0;i<m;i++)
            if(edge[i].u!=edge[i].v&&(edge[i].cost-in[edge[i].v]<eps))
        {
            pre[edge[i].v]=edge[i].u;
            in[edge[i].v]=edge[i].cost;
        }
        for(int i=0;i<n;i++)
            if(i!=root&&in[i]==INF)
                return -1;
        int tn=0;
        MEM(id,-1); MEM(vis,-1);
        in[root]=0;
        for(int i=0;i<n;i++)
        {
            res+=in[i];
            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]=tn;
                id[v]=tn++;
            }
        }
        if(tn==0) break;//没有有向环
        for(int i=0;i<n;i++)
            if(id[i]==-1)
                id[i]=tn++;
        for(int i=0;i<m;)
        {
            v=edge[i].v;
            edge[i].u=id[edge[i].u];
            edge[i].v=id[edge[i].v];
            if(edge[i].u!=edge[i].v)
                edge[i++].cost-=in[v];
            else
                swap(edge[i],edge[--m]);
        }
        n=tn;
        root=id[root];
    }
    return res;
}
double g[MAXN][MAXN];
int x[MAXN],y[MAXN];

double dis(int u,int v)
{
    return sqrt((double)((x[u]-x[v])*(x[u]-x[v])+(y[u]-y[v])*(y[u]-y[v])));
}

int main()
{
//    fread;
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                g[i][j]=INF;
        for(int i=0;i<n;i++)
            scanf("%d%d",&x[i],&y[i]);
        while(m--)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            u--; v--;
            if(u==v) continue;
            double len=dis(u,v);
            if(g[u][v]-len>eps)
                g[u][v]=len;
        }
        int num=0;
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                if(g[i][j]<INF)
        {
            edge[num].u=i; edge[num].v=j; edge[num++].cost=g[i][j];
        }
        double ans=zhuliu(0,n,num,edge);
        if(ans==-1) printf("poor snoopy\n");
        else printf("%.2lf\n",ans);
    }
    return 0;
}
</span>


理想主义文艺青年
关注
专栏目录
poj 3164 && tju 2248 最小树形图 朱刘算法
Join Hands
09-04 7727
个人觉得这个博客把这个算法说的比较详细了，直接搬过来吧，我再阐述一遍的话没有人家说的好，还容易说错。 ========================== 分割线之下摘自Sasuke_SCUT的blog======================================
最小树形图
_Griefs
08-24 268
日常%一发mmh学长。 时间复杂度为O（n * m）。 据说tarjan大神有一个（m+nlogn）的算法，但是找不到资料。。。 ①P4716 【模板】最小树形图： 题目背景 这是一道模板题。 题目描述 给定包含 n 个结点， m 条有向边的一个图。试求一棵以结点 r 为根的最小树形图，并输出最小树形图每条边的权值之和，如果没有以 r 为根的最小树形图，输出 −1。 输入格式 第一行包含三个整数 ...
poj 3164 Command Network(最小树形图模板题)
D_Double's Journey
10-06 5446
链接： http://poj.org/problem?id=3164 题目： Command Network Time Limit: 1000MS   Memory Limit: 131072K Total Submissions: 8922   Accepted: 2609 Description
【最小树形图】最小树形图复习
void* ptr = &OI;
06-25 2275
模板题： poj3164 Command Network /***********************************\ * @prob: poj3164 Command Network * * @auth: Wang Junji * * @stat: Accepted. * * @date: June.
TJU计科数字逻辑实验
09-17
《TJU计科数字逻辑实验详解》 天津大学计算机科学与技术专业的一门核心课程——数字逻辑与数字系统，其实践环节尤为重要。本实验旨在深入理解和应用数字电路的基本原理，通过一系列精心设计的项目，让学生掌握数字...
TJU-智算-算法设计与分析作业
09-17
第七次作业可能涉及图论问题，如最小生成树、最短路径等；第一次作业可能从基础的递归和分治策略入手。 "3020244271-彭筱涵-分治算法作业.pdf"则直接点明了主题，分治算法(Divide and Conquer)是一种将大问题划分为...
TJU-计科数据库实验（物流数据库ER模型与关系模型转换）
09-17
- 使用ER图表示实体及其之间的关系，比如实体间的多对多、一对多和一对一关系。 3. **关系模型转换**： - 根据ER模型，转换为关系表结构。例如，货主和承运方表可以直接对应，运单表需要关联货主和承运方的id，...
AR_tju:北洋AR，建立在天津大学录取通知书上，然后可以将学校的图像显示出来
05-11
AR_tju 北洋AR，建立在天津大学录取通知书上，然后可以将学校的图像显示出来 ### 1。 简介增强现实（AR）和虚拟现实（VR）是计算机图形学中的两个基本领域。 AR是物理，真实世界环境的实时直接或间接视图，其元素...
考研复习-英语二真题考试题集-带答案
09-14
英语二考研真题复习资料，带答案版
2024中美独角兽公司发展分析报告.pdf
09-14
全球各大洲独角兽企业分布、中美独角兽企业对比（数量、估值、新增及退榜情况、行业分布、所在城市）、
C++ 中的异步编程模型是什么
09-14
在C++中，异步编程模型是处理并发任务、提高程序性能和响应性的关键技术。以下是C++中实现异步编程的几种主要方式： 每种异步编程模型都有其适用场景和优缺点。选择合适的模型可以提高代码的可读性、可维护性和性能。随着C++标准的不断发展，异步编程模型也在不断进化，为开发者提供了更多的工具和选择。 在实际开发中，应根据具体需求选择合适的异步编程模型。例如，对于简单的异步任务，回调函数可能是最直接的选择；而对于需要结构化错误处理和结果获取的复杂异步任务，std::async和std::future可能更合适；在需要高效资源管理的场景下，线程池是一个不错的选择；而对于需要编写大量异步代码的现代应用程序，协程提供了一种更简洁、更直观的解决方案。 总之，C++中的异步编程模型是多核和高并发环境下提高程序性能的重要工具。通过合理使用这些模型，开发者可以构建出更高效、更可靠的软件系统。
正则表达式Regex是一种文本模式.docx
09-14
正则表达式（Regular Expression，简称Regex或Regexp）是一种文本模式，包括普通字符（例如，a 到 z 之间的字母）和特殊字符（称为"元字符"）。正则表达式使用单个字符串来描述、匹配一系列符合某个句法规则的字符串。正则表达式是强大的文本处理工具，广泛用于搜索、编辑或操作文本和数据。 基本组成 普通字符：大多数字符，包括所有大写和小写字母、所有数字、所有标点符号和一些其他符号，都是普通字符。正则表达式中的普通字符表示它们自身。例如，正则表达式test会匹配字符串"test"。 特殊字符（元字符）：一些字符在正则表达式中具有特殊的意义，如^、$、.、*、+、?、{、}、[、]、|、\等。这些特殊字符用于表示在搜索文本时要匹配的一个或多个字符。例如，.匹配除换行符之外的任何单个字符。 字符类：字符类允许你指定一组字符中的任何一个字符。例如，[abc]匹配"a"、"b"或"c"中的任意一个字符。你也可以使用范围，如[a-z]匹配任何小写字母。 预定义字符类：正则表达式提供了一些预定义字符类，用于匹配常见的字符集合。例如，\d匹配任何数字（等价于[0-9]），\s匹配任
基于struts+sqlserver网络购物系统毕业课程源码设计+论文资料
09-14
编号：150 系统功能： 展示网站最新的商品信息。 展示网站特价的商品信息。 为用户提供修改个人资料和查看在网站操作情况的平台。 提供用户在网站上购物的平台。 展示网站发布的公告信息。 展示商品的销量排行。 展示网站的友情链接信息。 对商品详细信息以及分类信息进行管理。 对用户基本资料、交易制度、消费情况及留言信息进行管理。 对用户提交的订单进行管理。 对管理员信息、网站公告信息、商业资讯信息及友情链接信息进行管理。 系统运行稳定，具有强大的数据处理能力。 操作注意事项 （1）本系统的用户名为：admin，密码为：admin （2）admin管理员的信息不能删除。 （3）用户注册登录后，可进行商品购买、商品信息查看以及订单查询操作。 （4）后台登陆地址：http://localhost:8080/WebShopping/bg-land.jsp。
基于ssm的高校毕业生就业管理系统设计与实现.docx
09-14
基于ssm的高校毕业生就业管理系统设计与实现.docx
退工、裁员培训讲义.ppt
09-14
退工、裁员培训讲义.ppt
认知实习decesion tree代码
最新发布
09-14
认知实习decesion tree代码
Java 算法：冒泡，选择，插入排序算法
09-14
Java 算法：冒泡，选择，插入排序算法
网约车出行预约系统数据库设计与数据操作
09-14
本课程设计旨在通过对网约车出行预约系统的数据库设计和数据操作的研究，提高对数据库系统的理解和应用能力。课程内容包括需求分析、数据库设计、数据操作等方面，通过实践操作，使掌握数据库设计的基本方法和技巧，提高实际操作能力。课程将从需求分析开始，明确网约车出行预约系统的功能需求和性能需求，为后续的数据库设计提供依据。接下来进行数据库设计，包括概念结构设计、逻辑结构设计和物理结构设计，确保数据库结构的合理性和高效性。在数据操作方面，将涵盖数据的增删改查、事务处理等内容，使学生能够熟练地进行数据操作和管理。通过理论教学和实践操作相结合的方式，使全面掌握网约车出行预约系统数据库设计和数据操作的方法和技巧，为今后从事相关领域的工作打下坚实的基础。
tju图卷积神经网络(GCN) 实验
06-08
好的，您需要我为您讲解如何在图卷积神经网络(GCN)上进行实验，是吗？那么，首先我们需要明确一些基本概念。 图卷积神经网络(GCN)是一种适用于图像数据的深度学习模型。与传统的卷积神经网络(CNN)不同，GCN用于处理...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值