SPOJ QTREE 1-3题解

最新推荐文章于 2019-08-12 16:47:00 发布
最新推荐文章于 2019-08-12 16:47:00 发布
阅读量1.9k 收藏
点赞数
分类专栏： 其他OJ 数据结构 OIer刷题记录 文章标签： SPOJ QTREE 树链剖分 倍增LCA DFS序
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ww140142/article/details/49332595
版权

昨天刷了几道QTREE，感觉码长萌萌哒；

然而因为本人太弱刷不动QTREE4，动态点分治并没有理解上去的能力；

于是暂且弃疗啦，在这里写点题解扔点代码吧；


QTREE1

题意：

给出一个有边权的树；

操作一：改变某条边权；

操作二：查询两点之间路径上最大边权；


题解：

树链剖分的姿势还是挺裸的，想了想没有什么好办法码了一发；

普通的树链剖分维护树上路径，加一个线段树，时间复杂度O(nlog^2n)；


QTREE2

题意：

给出一个有边权的树；

操作一：查询两点之间路径长度；

操作二：查询两点之间路径上经过的第K个点；


题解：

因为是无修的，操作一可以方便的用倍增LCA搞；

操作二在跑过一次LCA之后，YY一下也就是从某个点向上K步之类的东西；

时间复杂度O(nlogn)，手滑打错一个字母调了好久；


QTREE3

题意：

给出一个有点权的树；

每次查询某个点子树中第K小点权是哪个点；


题解：

操作只有一种，子树维护直接在DFS序上就可以了；

第K大用主席树来维护，时间复杂度O(nlogn)，空间复杂度O(nlogn)；

SPOJ的1.5G内存真是太爽啦 (然后到BZOJ上就MLE了)；


代码：


QTREE1：


#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 11000
#define lson l,mid,no<<1
#define rson mid+1,r,no<<1|1
using namespace std;
int next[N<<1],to[N<<1],len[N<<1],head[N],ce;
int n,X[N],Y[N],V[N];
int val[N],deep[N],fa[N],size[N],son[N],top[N],p[N],tot;
int ma[N<<2];
char op[20];
void Pushup(int no)
{
	ma[no]=max(ma[no<<1],ma[no<<1|1]);
}
void update(int l,int r,int no,int k,int val)
{
	if(l==r)
		ma[no]=val;
	else
	{
		int mid=l+r>>1;
		if(k<=mid)	update(lson,k,val);
		else		update(rson,k,val);
		Pushup(no);
	}
}
int query(int l,int r,int no,int st,int en)
{
	if(st<=l&&r<=en)
		return ma[no];
	else
	{
		int mid=l+r>>1;
		if(en<=mid)		return query(lson,st,en);
		else if(st>mid)	return query(rson,st,en);
		else	return max(query(lson,st,en),query(rson,st,en));
	}
}
void add(int x,int y,int v)
{
	to[++ce]=y;
	len[ce]=v;
	next[ce]=head[x];
	head[x]=ce;
}
void dfs1(int x)
{
	deep[x]=deep[fa[x]]+1;
	size[x]=1;
	for(int i=head[x];i;i=next[i])
	{
		if(to[i]!=fa[x])
		{
			fa[to[i]]=x;
			val[to[i]]=len[i];
			dfs1(to[i]);
			size[x]+=size[to[i]];
			if(size[to[i]]>size[son[x]])
				son[x]=to[i];
		}
	}
}
void dfs2(int x,int t)
{
	top[x]=t;
	p[x]=++tot;
	update(1,n,1,p[x],val[x]);
	if(son[x])
		dfs2(son[x],t);
	for(int i=head[x];i;i=next[i])
		if(to[i]!=fa[x]&&to[i]!=son[x])
			dfs2(to[i],to[i]);
}
void init()
{
	memset(head,0,sizeof(head));
	memset(son,0,sizeof(son));
	ce=tot=0;
}
int find(int x,int y)
{
	int ret=0;
	while(top[x]!=top[y])
	{
		if(deep[top[x]]<deep[top[y]])
			swap(x,y);
		ret=max(ret,query(1,n,1,p[top[x]],p[x]));
		x=fa[top[x]];
	}
	if(deep[x]<deep[y])
		swap(x,y);
	if(x!=y)
		ret=max(ret,query(1,n,1,p[son[y]],p[x]));
	return ret;
}
int main()
{
	int c,T,i,j,k,x,y,v;
	scanf("%d",&T);
	for(c=1;c<=T;c++)
	{
		scanf("%d",&n);
		init();
		for(i=1;i<n;i++)
		{
			scanf("%d%d%d",X+i,Y+i,V+i);
			add(X[i],Y[i],V[i]),add(Y[i],X[i],V[i]);
		}
		dfs1(1);
		dfs2(1,0);
		for(i=1;i<n;i++)
			X[i]=deep[X[i]]>deep[Y[i]]?X[i]:Y[i];
		while(scanf("%s",op)!=EOF&&op[0]!='D')
		{
			scanf("%d%d",&x,&y);
			if(op[0]=='C')
			{
				val[X[x]]=y;
				update(1,n,1,p[X[x]],y);
			}
			else
			{
				printf("%d\n",find(x,y));
			}
		}
	}
	return 0;
}


QTREE2：


#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 110000
using namespace std;
int next[N<<1],to[N<<1],val[N<<1],head[N],ce;
int fa[N][20],dis[N][20],deep[N];
char op[20];
void add(int x,int y,int v)
{
	to[++ce]=y;
	val[ce]=v;
	next[ce]=head[x];
	head[x]=ce;
}
void dfs(int x)
{
	int i;
	deep[x]=deep[fa[x][0]]+1;
	for(i=1;fa[x][i-1];i++)
		fa[x][i]=fa[fa[x][i-1]][i-1],
		dis[x][i]=dis[x][i-1]+dis[fa[x][i-1]][i-1];
	for(i=head[x];i;i=next[i])
	{
		if(to[i]!=fa[x][0])
		{
			fa[to[i]][0]=x;
			dis[to[i]][0]=val[i];
			dfs(to[i]);
		}
	}
}
int DIST(int x,int y)
{
	if(deep[x]<deep[y])
		swap(x,y);
	int t=15,ret=0;
	while(t>=0)
	{
		if(deep[fa[x][t]]>=deep[y])
			ret+=dis[x][t],x=fa[x][t];
		t--;
	}
	if(x==y)	return ret;
	t=15;
	while(t>=0)
	{
		if(fa[x][t]!=fa[y][t])
			ret+=dis[x][t],x=fa[x][t],
			ret+=dis[y][t],y=fa[y][t];
		t--;
	}
	return ret+dis[x][0]+dis[y][0];
}
int KTH(int x,int y,int k)
{
	int t=15,retx=0,rety=0,tx=x,ty=y;
	if(deep[tx]>deep[ty])
	{
		while(t>=0)
		{
			if(deep[fa[tx][t]]>=deep[ty])
				retx+=(1<<t),tx=fa[tx][t];
			t--;
		}
	}
	else
	{
		while(t>=0)
		{
			if(deep[fa[ty][t]]>=deep[tx])
				rety+=(1<<t),ty=fa[ty][t];
			t--;
		}
	}
	if(tx!=ty)
	{
		t=15;
		while(t>=0)
		{
			if(fa[tx][t]!=fa[ty][t])
				retx+=(1<<t),tx=fa[tx][t],
				rety+=(1<<t),ty=fa[ty][t];
			t--;
		}
		retx++,rety++;
	}
	if(retx+rety+1<k)	return -1;
	if(retx>=k-1)
	{
		k--;
		t=15;
		while(t>=0)
		{
			if(k&1<<t)
				x=fa[x][t];
			t--;
		}
		return x;
	}
	else
	{
		k=rety-k+retx+1;
		t=15;
		while(t>=0)
		{
			if(k&1<<t)
				y=fa[y][t];
			t--;
		}
		return y;
	}
}
void init()
{
	memset(head,0,sizeof(head));
	memset(dis,0,sizeof(dis));
	memset(fa,0,sizeof(fa));
	ce=0;
}
int main()
{
	int c,T,n,m,i,j,k,x,y,v;
	scanf("%d",&T);
	for(c=1;c<=T;c++)
	{
		init();
		scanf("%d",&n);
		for(i=1;i<n;i++)
		{
			scanf("%d%d%d",&x,&y,&v);
			add(x,y,v),add(y,x,v);
		}
		dfs(1);
		while(scanf("%s",op)!=EOF&&op[1]!='O')
		{
			if(op[1]=='I')
			{
				scanf("%d%d",&x,&y);
				printf("%d\n",DIST(x,y));
			}
			else
			{
				scanf("%d%d%d",&x,&y,&k);
				printf("%d\n",KTH(x,y,k));
			}
		}
	}
	return 0;
}

QTREE3： 


#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 110000
#define M 10000000
#define lson l,mid,ls[no]
#define rson mid+1,r,rs[no]
using namespace std;
int next[N<<1],to[N<<1],head[N],ce;
int val[N],dis[N],which[N];
int n,L[N],R[N],tim;
int size[M],ls[M],rs[M],root[N],tot;
void Insert(int l,int r,int &no,int val)
{
	int p=++tot;
	ls[p]=ls[no],rs[p]=rs[no],size[p]=size[no];
	no=p;
	size[no]++;
	if(l==r)	return ;
	int mid=l+r>>1;
	if(val<=mid)	Insert(lson,val);
	else			Insert(rson,val);
}
int query(int l,int r,int nol,int nor,int k)
{
	if(l==r)
		return l;
	else
	{
		int mid=l+r>>1;
		if(k<=size[ls[nor]]-size[ls[nol]])
			return query(l,mid,ls[nol],ls[nor],k);
		else
			return query(mid+1,r,rs[nol],rs[nor],k-size[ls[nor]]+size[ls[nol]]);
	}
}
void add(int x,int y)
{
	to[++ce]=y;
	next[ce]=head[x];
	head[x]=ce;
}
void dfs(int x,int pre)
{
	L[x]=++tim;
	root[tim]=root[tim-1];
	Insert(1,n,root[tim],val[x]);
	for(int i=head[x];i;i=next[i])
		if(to[i]!=pre)
			dfs(to[i],x);
	R[x]=tim;
}
int main()
{
	int m,i,j,k,x,y;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
	{
		scanf("%d",val+i);
		dis[i]=val[i];
	}
	sort(dis+1,dis+n+1);
	for(i=1;i<=n;i++)
	{
		val[i]=lower_bound(dis+1,dis+n+1,val[i])-dis;
		which[val[i]]=i;
	}
	for(i=1;i<n;i++)
	{
		scanf("%d%d",&x,&y);
		add(x,y),add(y,x);
	}
	dfs(1,0);
	scanf("%d",&m);
	for(i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&k);
		printf("%d\n",which[query(1,n,root[L[x]-1],root[R[x]],k)]);
	}
	return 0;
}



ww140142
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
SPOJ QTREE4 Query on a tree IV
Occult
05-13 1178
Description You are given a tree (an acyclic undirected connected graph) with N nodes, and nodes numbered 1,2,3...,N. Each edge has an integer value assigned to it(note that the value can be nega
SPOJ QTREE3 - Query on a tree again!
tangzhide123yy的博客
08-17 257
题目链接：https://www.spoj.com/problems/QTREE3/ 感谢Konjak谷弱博客的启发 题目分析 0 i : change the color of the i-th node (from white to black, or from black to white); 1 v : ask for the id of the first black node ...
SPOJ Qtree系列
weixin_30657541的博客
08-12 82
Qtree1 将边权变为这条边连接的两个点中深度更深的点的点权，这样就可以变为带修改链上最大点权。直接树链剖分即可。 下面是一份C语言代码 #include<stdio.h> #include<string.h> #define MAXN 10001 inline int read(){ int a = 0; int f = 0; char...
spoj Qtree1
ehdhg13455的博客
06-25 147
说几个需要注意的问题： 1.多组样例时初始化问题 首先 lca 的 f 数组和 根的 d 值需要初始化 其次，树剖的son值和num需要初始化 2.把边赋给点时需要注意的问题 首先，根节点对应的边应该注意赋初值 其次，线段树是对编号之后的树进行建树，应该注意。 最后，因为lca的边不应该被算上，所以处理方式和点的不太一样，注意一下。 #include<bits/stdc++...
SPOJ QTree 1
lemonoil的博客
07-07 505
树链剖分+线段树裸题，注意边权维护的细节。#include<algorithm> #include<cstring> #include<cstdio> #include<cmath> #define clr(x) memset((x),0,sizeof(x)) const int N=200020,inf=~0U>>1; using namespace std; char cmd[10]; stru
SPOJ-2.zip_SPOJ_SPOJ-TRIKA.cpp_online judge_sPOJ-Solution_spoj2
09-24
3. **SPOJ-GCPC11A.cpp**：可能对应的是某个编程竞赛的问题，可能涉及到动态规划、贪心算法等，要求对算法的选取和实现有较高要求。解决这类问题需要对问题背后的数学模型有深入理解。 4. **SPOJ-SUMSUMS.cpp**：从...
SPOJ-BACKUP-TOOL:在 SPOJ 上下载已提交代码的脚本
05-30
SPOJ-备份工具 介绍 在 Sphere Online Judge ( ) 中，您可以尝试所给的具有挑战性的问题。 它还使您能够查看和下载您自己的解决方案。 工具 SPOJ_BACKUP 备份所有已接受的提交并将它们保存在脚本所在的计算机位置。...
Ultimate SPOJ Rank Tracker-crx插件
04-03
语言:English (UK) 分机检查spoj用户的排名。 此扩展名有2个选项： - 在选择的比赛中显示用户排名 - 在选择比赛中最多... 15这意味着将有3个用户比较和3场比赛。第一个比赛中第一个用户的排名将显示在扩展名图标中，并
SPOJ Problems Score-crx插件
04-06
语言:English (United States) 在问题页面本身上显示spoj... spoj有超过20,000个问题。 通常，要了解问题的分数，您需要在问题页面中找到该问题并在“用户”列中悬停。 这是一个Chrome扩展，在问题页面本身上显示得分。
SPOJ QTREE 系列
yamiedie_
08-04 1301
QTREE Query on a tree 树链剖分：QTREE LCT:QTREE QTREE2 Query on a tree II 倍增LCA:QTREE2 PTO7J Query on a tree III dfs序+主席树：QTREE3 QTREE4 Query on a tree IV 边分治+堆：QTREE4
Spoj2798 Qtree3
Canopus
10-13 196
Description You are given a tree (an acyclic undirected connected graph) with NN nodes. The tree nodes are numbered from 1 to NN. In the start, the color of any node in the tree is white. We will ask
SPOJ - QTREE
常胜将军的博客
11-19 237
题目：点击打开链接 题意：给你一颗树，现在有两个操作，一种是改变某条边的权值，一种是查询点u到v之间的路径的最大边权。 分析：树链剖分模板题，将点权转化为边权。（树链剖分的模板题了。唯一不同的是这次权值是在边上，这样我们把每一条边的边权转移记录到深度较大的点上，因为每个点向上只有一条边，向下可能有多条边。这样我们就还是对点操作了。注意最后tp1 == tp2时，可能x != y,这样我们还需要查询...
spoj 2398 Qtree3
mlzmlz95的专栏
02-23 838
Description给出一棵树，树节点的颜色初始时为白色，有两种操作：0 x:把x号节点颜色取反1 x:询问1到x路径上第一个黑点编号Solution最近想练练剖分和线段树，于是想到来做做Qtree系列，不会Lct嘤嘤嘤>_<这题比较裸，直接剖分然后维护区间最浅的黑色点标号就可以了注意到把1当做根节点，其实剖分后每个区间左区间显然是更浅的，如果有黑点直接统计答案就可以了Code#include <
[SPOJ2798]QTREE3——[树链剖分]
sleepyNick
09-14 643
【题目描述】 给出N个点的一棵树（N-1条边），节点有白有黑，初始全为白。有两种操作： 0 i : 改变某点的颜色（原来是黑的变白，原来是白的变黑） 1 v : 询问1到v的路径上的第一个黑点，若无，输出-1 【输入格式】 第一行 N，Q，表示N个点和Q个操作 第二行到第N行N-1条无向边 再之后Q行，每行一个操作”0 i” 或者”1 v” (1 ≤ i, v ≤ N). 【输出格...
SPOJ QTREE(树链剖分)
俯瞰风景
10-15 571
题意：有一棵树，n个点，n - 1条边，边有权值，两种操作，query a b，询问a到b所有路径中的最大值，change i a，把第i条边权值设置a。 题解：线段树维护最大值。#include <cstdio> #include <cstring> #include <algorithm> using namespace std; const int N = 10005; const int
2018.10.06 spoj QTREE3（树链剖分）
苟为蒟蒻又何妨
10-06 216
传送门 发现把树投射到序列上。 发现每个黑点都只会对其子树产生影响。 再继续想想好像在投射到序列上之后好像就是查询从点uuu到根路径上深度最浅的黑点。 也就是logloglog段线段树上最靠左的黑点。 这个直接在线段树上二分到叶子结点返回下标就行了。 代码： #include&lt;bits/stdc++.h&gt; #define N 100005 #define lc (p&lt;&lt;1)...
SPOJ QTREE3 lct
九野的博客
07-05 1315
题目链接 题意： 给定n个点 q个询问 下面n-1行给出树边，点有黑或白色，初始化为白色 下面q行： 询问有2种： 1、 0 x 把x点黑变白，白变黑 2、1 x 询问Path(1,x)路径上第一个黑点的点标， 若不存在黑点则输出-1 思路： lct裸题 #include #include #include #include #include #include #i
SPOJ QTREE5
Cyhlnj
03-24 358
题意 一棵nnn个点的树，点从111到nnn编号。每个点可能有两种颜色：黑或白。 我们定义dist(a,b)dist(a,b)dist(a,b)为点aaa至点bbb路径上的边个数。 一开始所有的点都是黑色的。 要求作以下操作： 0i0i0 i 将点iii的颜色反转（黑变白，白变黑） 1v1v1 v 询问dist(u,v)dist(u,v)dist(u,v)的最小值，uuu与vvv可以相...
洛谷SPOJ如何注册
最新发布
07-27
chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [luogu p...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值