树的最长链-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 ...
随机推荐
- PHP APP端微信支付
前面已经写了手机APP支付宝支付,今天再把手机APP微信支付补上,前期的准备工作在这里就不多说了,可以参考微信支付开发文档,一定要仔细阅读开发文档,可以让你少踩点坑:准备工作完成后就是配置参数,调用统 ...
- Hive-1.2.1_05_案例操作
1. 建库建表 # 建库 create database exercise; # 建表 create table student(Sno int,Sname string,Sex string,Sag ...
- 使用Maven分环境打包:dev sit uat prod
使用Maven管理的项目,经常需要根据不同的环境打不同的包,因为环境不同,所需要的配置文件不同,比如database的连接信息,相关属性等等. 在Maven中,我们可以通过P参数和profiles元素 ...
- ES5-ES6-ES7_let关键字声明变量
let命令的介绍 let是ECMAScript6中新增的关键字,用于声明变量.它的用法类似于var var a = 3 let b = 4 let变量的声明 let 命令的特点不允许在同一作用域下声明 ...
- 深入浅出RxJava(二:操作符)
看完这篇blog,我相信你肯定想立即在你的项目中使用RxJava了,这篇blog将介绍许多RxJava中的操作符,RxJava的强大性就来自于它所定义的操作符. 首先先看一个例子: 准备工作 假设我有 ...
- 转://UDEV简介及配置过程
在Linux环境下安装Oracle11g RAC时,OS层面配置好多路径软件后(multipath),下一步就需要配置udev或asmlib来处理共享分区(Lun),以便Orace ASM能够看到这些 ...
- warning: ISO C++ forbids converting a string constant to 'char*'
第1种字符串赋值方式: char * fileName="./2017-09-02-10-34-10.xml";//这一种字符串赋值方式已经被ISO禁止了 第2种字符串赋值方式: ...
- 学习 JS滚轮事件(mousewheel/DOMMouseScroll)
学习 JS滚轮事件(mousewheel/DOMMouseScroll) 1-1 滚轮事件兼容性的差异 IE,chrome,safari 浏览器都使用 onmousewheel, 只有firefo ...
- CVE-2018-14634 - Linux create_elf_tables()中的整型溢出 - 翻译
原文:https://seclists.org/oss-sec/2018/q3/274 摘要 Qualys研究实验室的安全团队发现一个位于Linux内核函数create_elf_tables()中的整 ...
- Java多线程编程模式实战指南之Promise模式
Promise模式简介(转) Promise模式是一种异步编程模式 .它使得我们可以先开始一个任务的执行,并得到一个用于获取该任务执行结果的凭据对象,而不必等待该任务执行完毕就可以继续执行其他操作.等 ...