树的最长链-POJ 1985 树的直径(最长链)+牛客小白月赛6-桃花
求树直径的方法在此转载一下大佬们的分析:
可以随便选择一个点开始进行bfs或者dfs,从而找到离该点最远的那个点(可以证明,离树上任意一点最远的点一定是树的某条直径的两端点之一;树的直径:树上的最长简单路径)。再从找到的点出发,找到据该点的最远点,那么这两点就确定了树的一条直径,两点间距即为所求距离。
无意中看到一道水题,也就是POJ 1383
题目中给出了一个无环的迷宫,求出其中最长的一条路
我们知道无环图本质上可以认为就是树,所以此题完全可以使用树的最长链算法
即:随便从某个节点C开始DFS或BFS找到最远的一个点A,再从点A开始DFS或BFS找到最远的一点B,那么路径A->B必然是树上的最长路径。
这个算法很多人都知道并且当做结论使用,但很少看到有人给出正确性证明
所以本人简要分析后,给出一个简单的文字性证明:
(最好拿出纸笔画个图,有助于理解下面这段抽象的语死早描述)
首先我们知道对于一个图我们可以任取一个起始点C作为根来构造生成树(其实本来就是树,就是指定一个根C)
那么我们接下来可以分两种情况来讨论:
1. 若最长链包含树根(起始点C),那么DFS/BFS一次必然能够找到最长链的一端,再DFS/BFS一次必然找到最长链另一端。这种情况很容易理解。
2. 若最长链不包含树根(起始点C),在这里我们可以先做出【断言1】:从起始点C开始DFS/BFS查找最远节点A,一定会在中途遇到最长链上的点(暂且把首次遇到的最长链上的那个点称作TR);只要遇到了最长链上的点之后,子问题变为从TR作为根寻找最长链,然后就和上面的情况1相同了。
为什么【断言1】一定成立?首先如果最长链不包含起始点C,则最长链必然完全包含于C手下的某个子树T中(而且这个子树的根就是上面所说的TR)。然后用反证法,假设这个子树T的最深节点B的对于起始点C的深度为h(B),那么如果说从起始点C一直走到某个尽头A而没有经过TR的话,那么我们可以做出【断言2】:点A的深度h(A)必定小于等于h(B);与前提“A是距离起始点C最远的节点”矛盾。
所以只要【断言2】成立则上述反证成功。
而说明【断言2】仍然是反证法,如果点A的深度h(A)大于h(B),那么路径(A->C->TR)的长度至少为h(A)+1且大于h(B),所以路径(A->C->TR->B)的长度一定比子树T中的最长链更长!与前提“最长链一定完全包含于C手下的子树T中”矛盾。
至此得证。原文出处https://blog.csdn.net/macuilxochitl/article/details/19157579
这里讲POJ-1985
题自己看
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<queue>
#define INF 1e9+7
using namespace std;
const int MAXN = ;
const int MAXM = ;
int n,m;
struct EDGE
{
int v,next,w;
}edge[MAXN];
int head[MAXN],e;
int q[MAXN],vis[MAXN],d[MAXN];
void init()
{
e=;
memset(head,-,sizeof(head));
}
void add(int u,int v,int w)//链式前向星部分
{
edge[e].v=v;
edge[e].w=w;//权重
edge[e].next=head[u];//下一条边
head[u]=e++;//u为节点边的头节点
}
void bfs(int src)
{
for(int i=;i<=n;i++)vis[i]=,d[i]=INF;
int h=,t=;
queue<int>q;
q.push(src);
vis[src]=;
d[src]=;
while(!q.empty())
{
int u=q.front();
q.pop();
for (int i=head[u];i!=-;i=edge[i].next)
{
int v = edge[i].v;
int w = edge[i].w;
if (d[u]+w<d[v]){
d[v]=d[u]+w;
if (!vis[v])
{
q.push(v);//把它存下来
vis[v]=;
}
}
}
}
}
int main(){
int u,v,w;
char k;
scanf("%d%d",&n,&m);
init();
for(int i=;i<=m;i++)
{
scanf("%d%d%d%*s",&u,&v,&w);
add(u,v,w);
add(v,u,w);
}
bfs();
int pos=-,mx=-;
for (int i=;i<=n;i++)//找到距离这个点最大的并记录下距离它最远的点
if (d[i]>mx)
{
mx=d[i];
pos=i;
}
bfs(pos);//从距离它最远的点出发,再次bfs找这样就能找到直径
mx=-;
for (int i=;i<=n;i++)
if (d[i]>mx)mx=d[i];
printf("%d\n",mx);
return ;
}
类似题-牛客小白月赛6-桃花,这题没有边权值,可以给边一个权值1
这样答案就是最长链长度+1
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define INF 1e9+7
using namespace std;
const int MAXN = ;
const int MAXM = ;
int n,m;
struct EDGE
{
int v,next,w;
}edge[MAXN];
int head[MAXN],e;
int q[MAXN],vis[MAXN],d[MAXN];
void init()
{
e=;
memset(head,-,sizeof(head));
}
void add(int u,int v,int w)//链式前向星部分
{
edge[e].v=v;
edge[e].w=w;//权重
edge[e].next=head[u];//下一条边
head[u]=e++;//u为节点边的头节点
}
void bfs(int src)
{
for(int i=;i<=n;i++)vis[i]=,d[i]=INF;
int h=,t=;
q[t++]=src;
vis[src]=;
d[src]=;
while(h<t)
{
int u=q[h++];
for (int i=head[u];i!=-;i=edge[i].next)
{
int v = edge[i].v;
int w = edge[i].w;
if (d[u]+w<d[v]){
d[v]=d[u]+w;
if (!vis[v])
{
q[t++]=v;//把它存下来
vis[v]=;
}
}
}
}
}
int main(){
int u,v,w;
char k;
scanf("%d",&n);
init();
for(int i=;i<=n-;i++)
{
scanf("%d%d",&u,&v);
add(u,v,);
add(v,u,);
}
bfs();
int pos=-,mx=-;
for (int i=;i<=n;i++)
if (d[i]>mx)
{
mx=d[i];
pos=i;
}
bfs(pos);
mx=-;
for (int i=;i<=n;i++)
if (d[i]>mx)mx=d[i];
printf("%d\n",mx+);
return ;
}
树的最长链-POJ 1985 树的直径(最长链)+牛客小白月赛6-桃花的更多相关文章
- 树状数组求区间和模板 区间可修改 参考题目:牛客小白月赛 I 区间
从前有个东西叫树状数组,它可以轻易实现一些简单的序列操作,比如单点修改,区间求和;区间修改,单点求值等. 但是我们经常需要更高级的操作,比如区间修改区间查询.这时候树状数组就不起作用了,只能选择写一个 ...
- 牛客网 牛客小白月赛5 I.区间 (interval)-线段树 or 差分数组?
牛客小白月赛5 I.区间 (interval) 休闲的时候写的,但是写的心情有点挫,都是完全版线段树,我的一个队友直接就水过去了,为啥我的就超内存呢??? 试了一晚上,找出来了,多初始化了add标记数 ...
- 牛客小白月赛19 E 「火」烈火燎原 (思维,树)
牛客小白月赛19 E 「火」烈火燎原 (思维,树) 链接:https://ac.nowcoder.com/acm/contest/2272/E来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空 ...
- 牛客小白月赛16 小石的妹子 二分 or 线段树
牛客小白月赛16 这个题目我AC之后看了一下别人的题解,基本上都是线段树,不过二分也可以. 这个题目很自然就肯定要对其中一个进行排序,排完序之后再处理另外一边,另一边记得离散化. 怎么处理呢,你仔细想 ...
- 牛客小白月赛6C-桃花(DFS/BFS求树的直径)
链接:https://www.nowcoder.com/acm/contest/136/C 来源:牛客网 桃花 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言 ...
- 牛客小白月赛12 F 华华开始学信息学 (分块+树状数组)
链接:https://ac.nowcoder.com/acm/contest/392/F来源:牛客网 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/C++ 32768K,其他语言65536K ...
- 牛客小白月赛6 F 发电 树状数组单点更新 求区间乘积 模板
链接:https://www.nowcoder.com/acm/contest/136/F来源:牛客网 HA实验是一个生产.提炼“神力水晶”的秘密军事基地,神力水晶可以让机器的工作效率成倍提升. ...
- 牛客小白月赛12 H 华华和月月种树 (离线dfs序+线段树)
链接:https://ac.nowcoder.com/acm/contest/392/H 来源:牛客网 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/C++ 131072K,其他语言2621 ...
- 牛客小白月赛16 H小阳的贝壳 (线段树+差分数组)
链接:https://ac.nowcoder.com/acm/contest/949/H来源:牛客网 题目描述 小阳手中一共有 n 个贝壳,每个贝壳都有颜色,且初始第 i 个贝壳的颜色为 colico ...
随机推荐
- Hive-1.2.1_03_DDL操作
Hive官方文档:Home-UserDocumentation Hive DDL官方文档:LanguageManual DDL 参考文章:Hive 用户指南 注意:各个语句的版本时间,有的是在 hiv ...
- Arcgis10.3在添加XY数据时出现问题
准备通过excel表格(xls格式)中的经纬度生成点数据,但是选择数据的时候报错:连接到数据库失败,常规功能故障,外部表不是预期的格式.如下图所示: 解决方法: 将xls格式的表格另存为csv格式,重 ...
- 一个CSS值转REM的Sublime Text插件
CSSREM 一个CSS的px值转rem值的Sublime Text 3自动完成插件. 插件效果如下: 安装 下载本项目,比如:git clone https://github.com/flashli ...
- Usaco 2019 Jan Platinum
Usaco 2019 Jan Platinum 要不是昨天老师给我们考了这套题,我都不知道usaco还有铂金这么一级. 插播一则新闻:杨神坚持认为铂金比黄金简单,原因竟是:铜 汞 银 铂 金(金属活动 ...
- P1754 球迷购票问题
题目背景 盛况空前的足球赛即将举行.球赛门票售票处排起了球迷购票长龙. 按售票处规定,每位购票者限购一张门票,且每张票售价为50元.在排成长龙的球迷中有N个人手持面值50元的钱币,另有N个人手持面值1 ...
- php 怎么在foreach中循环数组的时候添加元素的属性
foreach($arr as $key => &$vo){ //注意,由于上面遍历的时候写了地址传值符&, //所以下面可以直接给$vo 赋值;如果不写&符号,下面这样 ...
- php单例模式实现对象只被创建一次
这是我在php面试题中遇到的一道试题,单例模式按字面来看就是某一个类只有一个实例,这样做的好处还是很大的,比如说数据库的连接,我们只需要实例化一次,不需要每次都去new了,这样极大的降低了资源的耗费. ...
- 天融信资料下载官方FTP服务器
2017年更新:最近发现天融信的FTP服务器变更为了ftp://ftp.topsec.com.cn/,用户名:topsec,密码:topsec2016(不知道后期密码会不会改成topsec2017,如 ...
- LeetCode算法题详解之两个数组的交集
题目背景: 这个与我们高中时期学习的交集是一样的,顺便复习一下相关的数学知识有助于更好的理解. 交集的定义: 对于两个集合A和B,定义A和B的交集为C,其中C={x|x属于A且X属于B},记作A∩B. ...
- MyBatis之反射技术+JDK动态代理+cglib代理
一.反射 引用百度百科说明: JAVA反射机制是在运行状态中,对于任意一个类,都能够知道这个类的所有属性和方法:对于任意一个对象,都能够调用它的任意方法和属性:这种动态获取信息以及动态调用对象方法的功 ...