[练习]:LCA练习2

最新推荐文章于 2024-04-09 23:08:42 发布
最新推荐文章于 2024-04-09 23:08:42 发布
阅读量295 收藏
点赞数 1
分类专栏： LCA zoj vijos bzoj codevs 树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xljer_/article/details/72883431
版权
树 同时被 3 个专栏收录
25 篇文章 0 订阅
订阅专栏
bzoj
24 篇文章 0 订阅
订阅专栏
vijos
15 篇文章 0 订阅
订阅专栏

1.bzoj1787紧急集合

找一点使这一点到三个点有最短距离 枚举三个两两LCA即可

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#define inf 999999999
using namespace std;
int n,q,cnt;
int deep[500001],head[500001],fa[500001][20];
bool vis[500001];
struct data{
	int to,next;
}e[1000001];
void ins(int u,int v)
{	
	e[++cnt].to=v;
	e[cnt].next=head[u];
	head[u]=cnt;
}
void insert(int u,int v)
{
	ins(u,v);
	ins(v,u);
}
void dfs(int x)
{
    vis[x]=1;
    for(int i=1;i<=18;i++)
    {
        if(deep[x]<(1<<i))break;
        fa[x][i]=fa[fa[x][i-1]][i-1];
    }
    for(int i=head[x];i;i=e[i].next)
    {
        if(vis[e[i].to])continue;
        deep[e[i].to]=deep[x]+1;
        fa[e[i].to][0]=x;
        dfs(e[i].to);
    }
}
int lca(int x,int y)
{
    if(deep[x]<deep[y])swap(x,y);
    int t=deep[x]-deep[y];
    for(int i=0;i<=18;i++)
       if(t&(1<<i))x=fa[x][i];
    for(int i=18;i>=0;i--)
       if(fa[x][i]!=fa[y][i])
       {x=fa[x][i];y=fa[y][i];}
    if(x==y)return x;
    return fa[x][0];
}
int val(int x,int y,int z)
{
	int p1=lca(x,y),p2=lca(x,z),p3=lca(y,z);
	int ans=inf,tmp,t;
	int q1=lca(p1,z),q2=lca(p2,y),q3=lca(p3,x);
	tmp=deep[x]+deep[y]-deep[p1]+deep[z]-2*deep[q1];
	if(tmp<ans){ans=tmp;t=p1;}
	
	tmp=deep[x]+deep[z]-deep[p2]+deep[y]-2*deep[q2];
	if(tmp<ans){ans=tmp;t=p2;}
	
	tmp=deep[y]+deep[z]-deep[p3]+deep[x]-2*deep[q3];
	if(tmp<ans){ans=tmp;t=p3;}
	
	printf("%d %d\n",t,ans);
}
int main()
{
    scanf("%d%d",&n,&q);
    for(int i=1;i<n;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        insert(u,v);
    }
    dfs(1);
    for(int i=1;i<=q;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        val(x,y,z);
    }
    return 0;
}

2.codevs1036商务旅行

基础题 就求LCA 即可 到一个城市改变now值

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 60005
using namespace std;
int b[N],p[N],nt[N],deep[N],fa[N][25],vis[N];
int n,m;int num=0;
void insert(){
	scanf("%d",&n);
	int x,y;
	for(int i=1;i<=n-1;i++){
		scanf("%d%d",&x,&y);
		num++;
		b[num]=y;
		nt[num]=p[x];
		p[x]=num;

		num++;
		b[num]=x;
		nt[num]=p[y];
		p[y]=num;
	}
}
void dfs(int x){
    vis[x]=1;
	for(int i=1;i<=20;i++){
		if(deep[x]<(1<<i)) break;
		fa[x][i]=fa[fa[x][i-1]][i-1];
	}
	for(int i=p[x];i;i=nt[i]){
		int v=b[i];
		if(vis[v]) continue;
		fa[v][0]=x;
		deep[v]=deep[x]+1;
		dfs(v);
	}
}
int lca(int x,int y){
	if(deep[x]<deep[y]) swap(x,y);
	int t=deep[x]-deep[y];
	for(int i=0;i<=20;i++){
		if(t&(1<<i)) x=fa[x][i];
	}
	if(x==y) return x;
	for(int i=20;i>=0;i--){
		if(fa[x][i]!=fa[y][i])
		{
			x=fa[x][i];y=fa[y][i];
		}
	}
	return fa[x][0];
}
int main(){
	insert();
	dfs(1);
	scanf("%d",&m);
	int now=1,to,ans=0;
	for(int i=1;i<=m;i++){
		scanf("%d",&to);
		int t=lca(now,to);
		int dis=deep[now]+deep[to]-2*deep[t];
		//printf("t=%d %d to %d is %d\n",t,now,to,dis);
		ans+=dis;
		now=to;
	}
	printf("%d",ans);
	return 0;
}

3.codevs2370小机房的树

同样的基础题 比上面的还简单

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define M 100005
#define N 100005
using namespace std;
int n,m;int num=0;
int b[M],w[M],nt[M],p[N];
int fa[M][20],d[M],val[N];
void insert(){
	scanf("%d",&n);
	for(int i=1;i<=n-1;i++){
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		num++;
		b[num]=y;
		w[num]=z;
		nt[num]=p[x];
		p[x]=num;

		num++;
		b[num]=x;
		w[num]=z;
		nt[num]=p[y];
		p[y]=num;
	}
}
void dfs(int x){
	for(int i=1;i<=20;i++){
		if(d[x]<(1<<i))break;
		fa[x][i]=fa[fa[x][i-1]][i-1];
	}
	for(int i=p[x];i>0;i=nt[i]){
		int v=b[i];
		if(v==fa[x][0]) continue;//无向边 不走已经走过的父节点 
		fa[v][0]=x;
		d[v]=d[x]+1;//求深度=父节点深度+1 
		val[v]=val[x]+w[i];//求与根结点距离 和深度类似
        dfs(v);
	}
}
int lca(int x,int y){//求LCA 
	int h;
	if(d[x]<d[y]) swap(x,y);
	for(h=0;(1<<h)<=d[x];h++);
	h--;
	for(int i=h;i>=0;i--){
		if(d[x]-(1<<i)>=d[y])
			x=fa[x][i];
	}
	if(x==y) return x;
	for(int i=h;i>=0;i--){
		if(fa[x][i]!=fa[y][i]){
			x=fa[x][i];
			y=fa[y][i];
		}
	}
	return fa[x][0]; 
}	 
int main(){
	insert();
	dfs(0);
	scanf("%d",&m);
	int x,y;
	for(int i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		printf("%d\n",val[x]+val[y]-2*val[lca(x,y)]);
	}
	
	return 0;
}

4.vijos1460拉力赛

记录深度数组和求LCA即可

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 10005
#define ll long long
using namespace std;
int n,m;
int b[N],p[N],nt[N],fa[N][21],d[N];
ll w[N];
ll val[N];
ll cnt,ans;
void insert(){
	scanf("%d%d",&n,&m);
	int x,y,z;
	for(int i=1;i<=n-1;i++){
		scanf("%d%d%d",&x,&y,&z);
		b[i]=y;
		w[i]=z;
		nt[i]=p[x];
		p[x]=i;
	}
}
void dfs(int x){
	for(int i=1;i<=20;i++){
		if(d[x]<(1<<i)) break;
		fa[x][i]=fa[fa[x][i-1]][i-1];
	}
	for(int i=p[x];i;i=nt[i]){
		int v=b[i];
		if(v==fa[x][0]) continue;
		fa[v][0]=x;
		d[v]=d[x]+1;
		val[v]=val[x]+w[i];
		dfs(v);
	}
}
bool pdlca(int x,int y){
	int t=d[x]-d[y];
	for(int i=0;i<=20;i++)
		if(t&(1<<i)) x=fa[x][i];
	if(x==y) return 1;
	else return 0;
}
int main(){
	insert();
	dfs(1);
	int u,v;
	for(int i=1;i<=m;i++){
		scanf("%d%d",&u,&v);
		if(d[v]<=d[u]) continue;
		if(!pdlca(v,u)) continue;
		cnt++;
		ans+=val[v]-val[u];
	}
	printf("%lld\n",cnt);
	printf("%lld",ans);
	return 0;
}

5.zoj3195 design the city

求树上三个点到某点 

有最小距离和输出这个距离和方法是 

会求两点最小距离 就会求三点距离=三个两两之间的最短距离相加再除2 

(不明觉厉。。。)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#define N 100005
using namespace std;
int b[N],p[N],nt[N],fa[N][30],val[N],w[N],deep[N];
int n,q,num=0;
void insert(){
    num=0;
	for(int i=1;i<=n-1;i++)
	{
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		num++;
		b[num]=y;
		w[num]=z;
		nt[num]=p[x];
		p[x]=num;
		num++;
		b[num]=x;
		w[num]=z;
		nt[num]=p[y];
		p[y]=num;
	}
}
void dfs(int x,int h){
	deep[x]=h;
	for(int i=1;i<=26;i++){
		if(deep[x]<(1<<i)) break;
		fa[x][i]=fa[fa[x][i-1]][i-1];
	}
	for(int i=p[x];i;i=nt[i]){
		int v=b[i];
		if(v==fa[x][0]) continue;
		fa[v][0]=x;
		val[v]=val[x]+w[i];
		dfs(v,h+1);
	}
}
int lca(int x,int y){
	if(deep[x]<deep[y]) swap(x,y);
	int t=deep[x]-deep[y];
	for(int i=0;i<=26;i++)
		if(t&(1<<i)) x=fa[x][i];
	if(x==y) return x;
	for(int i=26;i>=0;i--){
		if(fa[x][i]!=fa[y][i])
		{
			x=fa[x][i];y=fa[y][i];
		}
	}
	return fa[x][0];
}
void clear(){
	memset(fa,0,sizeof(fa));
	memset(p,0,sizeof(p));
	memset(val,0,sizeof(val));
}
int query(int x,int y){
	int t=lca(x,y);
	return val[x]+val[y]-2*val[t];
}
int deal(int x,int y,int z){
	int ans=0;
	ans=query(x,y)+query(y,z)+query(x,z);
	ans=ans>>1;
	return ans;
}
int main(){
	int x,y,z;
	int tot=0;
	while(scanf("%d",&n)!=EOF){
		if(tot++)  {printf("\n");}//如果最后多了一个空行就PE了  
		clear();
		insert();
		dfs(0,0);
		scanf("%d",&q);
		for(int i=1;i<=q;i++){
			scanf("%d%d%d",&x,&y,&z);
			printf("%d\n",deal(x,y,z));
		}
	}
	return 0;
}



Lndulge.
关注
专栏目录
LCA练习题
卡西莫多
11-11 1166
poj1330题意：求祖先，解释了啥是LCA。 想法：模板。#include<cstdio> #include<algorithm> #include<cstring> #include<vector> using namespace std; #define M 100009 int f[M],in[M],ancestor[M]; vector<int> tree[M],query[M]; in
LCA四种写法其一：倍增
m0_62021646的博客
10-08 774
文章目录众所周知一、倍增区间倍增二、LCA2.读入数据总结 众所周知 LCA有四种写法 一、倍增 fa[n][logn]，fa[i][j]表示i节点的第2^j个父亲。 显然：fa[i][j]=fa[ fa[i][j-1] ][j-1] 即i的第2j个父亲是i的第2(j-1)个父亲的第2(j-1)个父亲。 于是，求i的第k个父亲的复杂度O(k)—>O(logk) 求i的第k个父亲： 1. for(int j=0; j<=log2(k); j++) 2. if((1<<j)&a.
【蓝桥杯练习】tarjan算法求解LCA
m0_59385870的博客
04-03 463
函数中就能得出每个询问的答案，当然对于起始和结束都在同一个节点的情况下,它的答案就是当前节点的出度，我们可以进行特判一下。(最近公共祖先)模型的题目,我们可以使用多种方法来解决该问题，这里我们使用更好写的离线的。：每个点遍历一次,复杂度级别。，所以我们需要减去两倍的。的时间，所以还得加上一个。但很明显，我们多加了两。还是一道比较明显的求。但延迟时间还包括经过。
[练习] LCA练习1(最基础)
Edt!
06-02 392
hdu 2586 How Far Away? lca求最短距离（模板题） #include #include #include #include #include #include #define N 40005 using namespace std; int n,m;int num=0; int b[2*N],w[2*N],nt[2*N],p[2*N]; int fa[N][20],d[N
剑指offer所有编程练习总结分析
热门推荐
张彦峰的博客
07-25 168万+
给出剑指offer所有题目描述和对应的解题思路，并且在最有给出对应的代码展示和测试内容 由于牛客剑指Offer题目经过整理和筛选，是面试备考的经典资源之一，因此它的重要性在面试准备过程中是很高的。做这些题目可以帮助面试者熟悉常见的面试题型、巩固基础知识、提升编程能力，为应对面试提供了很好的准备。然而，还需注意不仅仅局限于做题，理解背后的原理和思想同样重要，以便能够在面试时灵活运用和展示自己的能力。
2016夏季练习——LCA
kiss forever
07-20 226
LCA，在线
【离线Tarjan 模板】LCA tarjan 算法 练习： hdu 2586 + poj 1986
qq_24664053的博客
05-16 494
个人一直觉得kuangbin大牛的LCA板子 感觉有点不太好用，于是在这个大牛的博客帮助下，更换了自己的板子； http://www.cnblogs.com/scau20110726/archive/2013/05/26/3100265.htmlhdu 2586， 给一颗树 ，每一条边有权值，然后m次询问，每一次询问两个点：i，j。 输出 i，j的最短距离。#include <iostream>
牛客练习赛55 E 树
Ajwallet
12-14 323
树形dp，二次换根(算是吧
P3379 【模板】最近公共祖先（LCA）|倍增法练习
最新发布
qq_61814350的博客
04-09 507
定义dep数组记录深度，edge二维数组存结点u的邻接点，fa[u][i]存结点u的i级祖先结点（i从0到19）dfs建表：记录深度dep[u]，为父节点深度+1。每个结点的i级祖先（通过递推获得，u的i级祖先等于u的i-1级祖先的 i-1级祖先，相当于我一半一半的跳），所以级数从小到大枚举。两者计算完毕，向下遍历所有孩子结点，注意，要判断孩子是否为父节点，防止往回走。
编程题练习2
qq_37453637的博客
05-03 318
1. 正整数A和正整数B 的最小公倍数是指 能被A和B整除的最小的正整数值，设计一个算法，求输入A和B的最小公倍数。 输入描述：输入两个正整数A和B。 输出描述：输出A和B的最小公倍数。 示例： 输入：5 7 输出：35 参考代码： import java.util.Scanner; //最小公倍数 = 两数之积除以最大公约数，这里使用碾转相除法进行最大公约数的求解：即a与b的最大公约数可以转化为a、b之间的余数为两者之间最小的数之间的公约数。 // 所以对于输入的两个数进行连续求余，直到余数为0，求余
RMQ与LCA问题详解：从提出到解决
2. Tarjan算法：用于解决LCA问题，预处理时间为O(N)，在线查询时间为O(Nα(N)+Q)，其中α(N)是逆阿克曼函数，通常增长极慢，因此查询效率较高，空间消耗为O(N)。 3. ±1RMQ算法：专门处理±1RMQ问题，预处理时间为O...
【LCA练习 倍增 tarjin】 洛谷P3379【模板】最近公共祖先 P1967货车运输
wldcmzy的博客
04-17 161
洛谷P3379【模板】最近公共祖先 构造fa数组的正确写法是这样的 for(int i=1; i<=20; i++){ fa[x][i] = fa[fa[x][i-1]][i-1]; } 然而第一遍错写成了这样 for(int i=1; i<=20; i++){ fa[x][i] = fa[father][i-1]; } 这显然是错误的，不是方法他的父亲的上2^i-1个节点，而是访问它的上2^i-1个节点的上2^i-1个节点 引以为戒 #include <iostream> #i
算法面试攻略：从RMQ到LCA，数据结构与线段树解析
2. 复习基础面试题：通过研究如微软面试100题这样的题库，了解常见题型和考察点。 3. 数据结构基础：深入学习数据结构，包括链表、树、图、排序算法等，这对于解决实际问题至关重要。 4. 学习《算法导论》：这本书...
[poj 3321]:Apple Tree(树状数组/线段树 和dfs序)
Edt!
05-24 1155
Apple Tree Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 29084 Accepted: 8633 Description There is an apple tree outside of kaka’s house. Every autumn, a lot of apples wi
[hdu&poj&洛谷] 经典线段树练习题
Edt!
06-01 763
1.hdu1166敌兵布阵 #include #include #define maxn 50000 int ans; struct node { int left,right,sum; int mid() { return (left+right)>>1; } }tree[maxn*4]; void btree(int left,int r
[vijos 1448]: 校门外的树(树状数组/线段树)
Edt!
05-23 629
描述校门外有很多树，有苹果树，香蕉树，有会扔石头的，有可以吃掉补充体力的…… 如今学校决定在某个时刻在某一段种上一种树，保证任一时刻不会出现两段相同种类的树，现有两个操作： K=1，K=1，读入l、r表示在区间[l,r]中种上一种树，每次操作种的树的种类都不同 K=2，读入l,r表示询问l~r之间能见到多少种树 （l,r>0） 格式输入格式第一行n,m表示道路总长为n，共有m个操作 接下
[笔记]: 记忆化搜索+hash优化&树形dp
Edt!
06-07 559
1.记忆化搜索 记忆化搜索就是搜索过程中记下搜索的值 很简单 没什么好说的 下面是一个简单的斐波拉契数列的记忆化搜索代码 #include #include #include #include #include #include #define N 100000 using namespace std; int f[N]; int dfs(int x){ if(f[x]) return f
[bzoj2144]: 跳跳棋
Edt!
06-02 543
2144: 跳跳棋 Time Limit: 10 Sec  Memory Limit: 259 MB Submit: 689  Solved: 326 [Submit][Status][Discuss] Description 跳跳棋是在一条数轴上进行的。棋子只能摆在整点上。每个点不能摆超过一个棋子。我们用跳跳棋来做一个简单的游戏：棋盘上有3颗棋子，分别在a，b，c这三个位置。我们要
[noip2013] 货车运输(最大生成树+并查集+LCA)
Edt!
06-02 521
货车运输 描述 A 国有 n 座城市，编号从 1 到 n，城市之间有 m 条双向道路。每一条道路对车辆都有重量限制，简称限重。现在有 q 辆货车在运输货物，司机们想知道每辆车在不超过车辆限重的情况下，最多能运多重的货物。 格式 输入格式 第一行有两个用一个空格隔开的整数 n，m，表示 A 国有 n 座城市和 m 条道路。 接下来 m 行每行 3 个整数 x、
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值