poj 3237 树链剖分模板(用到线段树lazy操作)
/*
本体在spoj375的基础上加了一些操作,用到线段树的lazy操作模板类型
*/
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<map>
#include<string.h>
#include<stdlib.h>
#include<math.h>
using namespace std;
#define N 11000
#define inf 0x3fffffff
int head[N];
int son[N];//记录与当前点相连的数目最多的子节点的下标
int fa[N];//记录上一个父节点
int siz[N];//记录当前节点的子节点的数目
int top[N];//当前链的最顶层
int f[N];//重新标记
int fp[N];//记录重新标记前的点
int deep[N];//深度
int w[N];// 记录当前点与其父节点的关系
int nu,yong;
int Max;
struct nodee
{
int u,v,w,next;
} bian[N*4],ff[N];
void init()
{
yong=nu=0;
memset(head,-1,sizeof(head));
memset(son,-1,sizeof(son));
}
void addedge(int u,int v,int w)
{
bian[yong].u=u;
bian[yong].v=v;
bian[yong].next=head[u];
head[u]=yong++;
}
void dfs(int u,int father,int d)
{
deep[u]=d;//记录深度
siz[u]=1;//初始化
fa[u]=father;//记录父节点
int i;
for(i=head[u]; i!=-1; i=bian[i].next)
{
int v=bian[i].v;
if(v!=father)
{
dfs(v,u,d+1);//
siz[u]+=siz[v];//回溯时累加数目
if(son[u]==-1||siz[son[u]]<siz[v])//son[u]记录当前点相连的点的子节点数目最多的点,即重链
son[u]=v;//u-v
}
}
}
void getnu(int u,int cnt)
{
top[u]=cnt;//记录当前链的顶端
f[u]=nu++;//重新标记
fp[f[u]]=u;//记录标记前的编号
if(son[u]==-1)return ;//如果他没有儿子节点
getnu(son[u],cnt);//重链
int i;
for(i=head[u]; i!=-1; i=bian[i].next)
{
int v=bian[i].v;
if(v!=fa[u]&&v!=son[u])//排除与u相连的重链
getnu(v,v);//轻链
}
return;
}
//******以上求重链和轻链以及各部分的初始化,下面是线段树部分和查询询问**********//
struct node
{
int l,r,maxx,yanchi,minn;
} tree[N*4];
int Ma(int v,int vv)
{
return v>vv?v:vv;
}
int Min(int v,int vv) {
return v>vv?vv:v;
}
void pushdown(int t) {
while(tree[t].yanchi) {
tree[t*2].maxx=-tree[t*2].maxx;
tree[t*2].minn=-tree[t*2].minn;
swap(tree[t*2].maxx,tree[t*2].minn);
tree[t*2].yanchi^=1;
tree[t*2+1].maxx=-tree[t*2+1].maxx;
tree[t*2+1].minn=-tree[t*2+1].minn;
swap(tree[t*2+1].maxx,tree[t*2+1].minn);
tree[t*2+1].yanchi^=1;
tree[t].yanchi=0;
}
return ;
}
void pushup(int t) //回溯时更新最大值和最小值
{
tree[t].maxx=Ma(tree[t*2].maxx,tree[t*2+1].maxx);
tree[t].minn=Min(tree[t*2].minn,tree[t*2+1].minn);
}
void build(int t,int l,int r)//建树
{
tree[t].l=l;
tree[t].r=r;
tree[t].yanchi=0;
if(tree[t].l==tree[t].r)
{
tree[t].maxx=tree[t].minn=w[tree[t].l];//记录边权值
return ;
}
int mid=(l+r)>>1;
build(t*2,l,mid);
build(t*2+1,mid+1,r);
pushup(t);
}
void qury(int t,int l,int r)//询问区间的最大值
{
if(tree[t].l==l&&tree[t].r==r)//如果查到
{
Max=Ma(Max,tree[t].maxx);//
return ;
}
int mid=(tree[t].l+tree[t].r)>>1;
pushdown(t);
if(r<=mid)
qury(t*2,l,r);
else if(l>mid)
qury(t*2+1,l,r);
else
{
qury(t*2,l,mid);
qury(t*2+1,mid+1,r);
}
pushup(t);
}
int findmax(int u,int v)//查找最大值
{
int f1=top[u];//得到顶端编号值
int f2=top[v];
int ans=-inf;//初始化最小值
while(f1!=f2)//结束条件,再通一个重链上
{
if(deep[f1]<deep[f2])//从最深层开始网上
{
swap(f1,f2);//交换
swap(u,v);
}
Max=-inf;
qury(1,f[f1],f[u]);//询问重新编号后的f[u]和其顶端节点之间的最大值,从而使其从f[u]跳到顶端
ans=Ma(ans,Max);//ans储存最大值
u=fa[f1];//从f1向上跳一步,不管当前链是轻链还是重链
f1=top[u];//得到跳一步后的顶端节点,继续比较
}//
if(v==u)return ans;//如果在同一点就直接返回
if(deep[u]>deep[v])swap(u,v);//得到u,v之间的最小值
Max=-inf;
qury(1,f[son[u]],f[v]);//求出u的子节点的f[v]---f[u的最大数目子节点];,因为此时他们在一个重链上
ans=Ma(ans,Max);//求出最大值
return ans;
}
void update(int t,int x,int y)//更新
{
if(tree[t].l==x&&tree[t].r==x)//
{
tree[t].maxx=tree[t].minn=y;
return ;
}
pushdown(t);
int mid=(tree[t].l+tree[t].r)/2;
if(x<=mid)
update(t*2,x,y);
else
update(t*2+1,x,y);
pushup(t);
}
void chant(int t,int l,int r) {
if(tree[t].l==l&&tree[t].r==r) {
tree[t].maxx=-tree[t].maxx;
tree[t].minn=-tree[t].minn;
swap(tree[t].maxx,tree[t].minn);
tree[t].yanchi^=1;
return ;
}
pushdown(t);
int mid=(tree[t].l+tree[t].r)/2;
if(r<=mid)chant(t*2,l,r);
else if(l>mid)chant(t*2+1,l,r);
else {
chant(t*2,l,mid);
chant(t*2+1,mid+1,r);
}
pushup(t);
}
void change(int u,int v) {
int f1=top[u];
int f2=top[v];
while(f1!=f2) {
if(deep[f1]<deep[f2]) {
swap(f1,f2);
swap(u,v);
}
chant(1,f[f1],f[u]);
u=fa[f1];
f1=top[u];
}
if(u==v)return ;
if(deep[u]>deep[v])swap(u,v);
chant(1,f[son[u]],f[v]);
return ;
}
int main()
{
int T,i,n;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
init();
for(i=1; i<n; i++)
{
scanf("%d%d%d",&ff[i].u,&ff[i].v,&ff[i].w);
addedge(ff[i].u,ff[i].v,ff[i].w);
addedge(ff[i].v,ff[i].u,ff[i].w);
}
dfs(1,1,0);//深搜的到 每个节点的深度,父节点和其子节点的数目(包括本身),还有最大数目的子节点的编号
getnu(1,1);//得到重链或者轻链的顶端和重新编号并记录重新编号前的值。如果是轻链的话
for(i=1; i<n; i++)
{
if(deep[ff[i].u]<deep[ff[i].v])
swap(ff[i].u,ff[i].v);//得到深度最大的节点
w[f[ff[i].u]]=ff[i].w;//记录重新编号后的当前点与上一个点的权值
}
build(1,1,nu-1);//建树
char s[222];
int x,y;
while(scanf("%s",s),strcmp(s,"DONE"))
{
scanf("%d%d",&x,&y);
if(s[0]=='Q')
printf("%d\n",findmax(x,y));//找区间最大值
else
if(s[0]=='C')
update(1,f[ff[x].u],y);//更换区间中某个值,ff[x].u是深度较大的数,所以不会出现越界情况,即f[ff[x].u]不为0
else
change(x,y);
}
}
return 0;
}
poj 3237 树链剖分模板(用到线段树lazy操作)的更多相关文章
- Qtree3题解(树链剖分(伪)+线段树+set)
外话:最近洛谷加了好多好题啊...原题入口 这题好像是SPOJ的题,挺不错的.看没有题解还是来一篇... 题意: 很明显吧.. 题解: 我的做法十分的暴力:树链剖分(伪)+线段树+\(set\)... ...
- Luogu 2590 [ZJOI2008]树的统计 / HYSBZ 1036 [ZJOI2008]树的统计Count (树链剖分,LCA,线段树)
Luogu 2590 [ZJOI2008]树的统计 / HYSBZ 1036 [ZJOI2008]树的统计Count (树链剖分,LCA,线段树) Description 一棵树上有n个节点,编号分别 ...
- 【bzoj4999】This Problem Is Too Simple! 树链剖分+动态开点线段树
题目描述 给您一颗树,每个节点有个初始值. 现在支持以下两种操作: 1. C i x(0<=x<2^31) 表示将i节点的值改为x. 2. Q i j x(0<=x<2^31) ...
- BZOJ 3531 [Sdoi2014]旅行 树链剖分+动态开点线段树
题意 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天面条神教.隐形独角兽教.绝地教都是常见的信仰. 为了方便,我们用 ...
- [bzoj 3531][SDOI2014]旅行(树链剖分+动态开点线段树)
题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3531 分析: 对于每个颜色(颜色<=10^5)都建立一颗线段树 什么!那么不是M ...
- 洛谷P3313 [SDOI2014]旅行(树链剖分 动态开节点线段树)
题意 题目链接 Sol 树链剖分板子 + 动态开节点线段树板子 #include<bits/stdc++.h> #define Pair pair<int, int> #def ...
- 刷题总结——骑士的旅行(bzoj4336 树链剖分套权值线段树)
题目: Description 在一片古老的土地上,有一个繁荣的文明. 这片大地几乎被森林覆盖,有N座城坐落其中.巧合的是,这N座城由恰好N-1条双 向道路连接起来,使得任意两座城都是连通的.也就是说 ...
- 【BZOJ3531】[Sdoi2014]旅行 树链剖分+动态开点线段树
[BZOJ3531][Sdoi2014]旅行 Description S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天 ...
- bzoj3531——树链剖分+动态开点线段树
3531: [Sdoi2014]旅行 Time Limit: 20 Sec Memory Limit: 512 MB Description S国有N个城市,编号从1到N.城市间用N-1条双向道路连 ...
- BZOJ4999: This Problem Is Too Simple!树链剖分+动态开点线段树
题目大意:将某个节点的颜色变为x,查询i,j路径上多少个颜色为x的点... 其实最开始一看就是主席树+树状数组+DFS序...但是过不去...MLE+TLE BY FCWWW 其实树剖裸的一批...只 ...
随机推荐
- (前缀和 内存分配)51NOD 1081 子段求和
给出一个长度为N的数组,进行Q次查询,查询从第i个元素开始长度为l的子段所有元素之和. 例如,1 3 7 9 -1,查询第2个元素开始长度为3的子段和,1 {3 7 9} -1.3 + 7 + 9 ...
- RT-Thread 设备驱动-硬件定时器浅析与使用
RT-Thread 4.0.0 访问硬件定时器设备 应用程序通过 RT-Thread 提供的 I/O 设备管理接口来访问硬件定时器设备,相关接口如下所示: 函数 描述 rt_device_find() ...
- $割点割顶tarjan$
原题 #include <bits/stdc++.h> using namespace std; typedef long long LL; inline LL read () { LL ...
- C语言小项目-火车票订票系统
list.h #ifndef __LIST_H__ #define __LIST_H__ #include "stdafx.h" #include <stdio.h> ...
- linux学习之路7 linux下获取帮助
help 帮助 ls -h或者ls - -help man 最常用的帮助命令 man (+数字 )+命令 (数字代表文档帮助类型) man -k 关键字 可以用来查询包含该关键字的文档 info 与m ...
- MySQL的两种存储引擎storage engine特点和对比
MyISAM 优点:快速读取数据, 占用空间小 缺点:不支持事务,外键 (表级别锁) InnoDB 优点:支持事务,外键; 高性能(CPU效率高) 缺点: 慢,占空间 (行级别锁)
- [转]iOS WebKit browsers and auto-zooming form controls
问题描述:https://github.com/jquery/jquery-mobile/issues/2581 本文转自:http://www.456bereastreet.com/archive/ ...
- webpack学习汇总
一. 安装 window : 附件 --- 命令提示符 1:node -------- http://pan.baidu.com/s/1boFor3D node -v : 查看版本: npm conf ...
- Spring.Net学习笔记(0)-错误汇总
1.错误一:ObjectDefinitionStoreException "Spring.Objects.Factory.ObjectDefinitionStoreException&quo ...
- Scala-基础-函数(1)
import junit.framework.TestCase //函数(1) class Demo5 extends TestCase { def testDemo(){ println(" ...