poj3237--Tree(树链剖分+线段树)

2024-08-25 00:58

本文主要是介绍poj3237--Tree(树链剖分+线段树),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目链接:poj--3237

题意很简单,给出n个节点的一棵树,有三种操作:

1、C修改第i条边的值为v

2、N改变节点a到b内边的权值的符号(取反)

3、Q询问节点a到b内权值的最大值

首先树链剖分,将边整合到线段树上,线段树数组cl,因为存在取反操作,所以最大值可能是由最小值取反得到,所以记录最大和最小值,cl[i][0]记录第i段的最大值,cl[i][1]记录最小值,lazy做标记,该段是否取反,两次取反==没有取反

 

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;
#define LL __int64
#define INF 0x3f3f3f3f
#define maxn 11000
#define int_now int l,int r,int rt
#define now l,r,rt
#define lson l,(l+r)/2,rt<<1
#define rson (l+r)/2+1,r,rt<<1|1
struct node{
int u , v , w ;
int next ;
} edge[maxn<<1] ;
int head[maxn] , cnt , vis[maxn] ;
int num[maxn] , dep[maxn] , fa[maxn] , son[maxn] , top[maxn] , w[maxn] , step ;
int cl[maxn<<2][2] , lazy[maxn<<2] , n ;//cl[i][0]最大值,cl[i][1]最小值
void add(int u,int v,int w) {
edge[cnt].u = u ; edge[cnt].v = v ; edge[cnt].w = w ;
edge[cnt].next = head[u] ; head[u] = cnt++ ;
return ;
}
void dfs1(int u) {
int i , v ;
num[u] = 1 ;
son[u] = -1 ;
for(i = head[u] ; i != -1 ; i = edge[i].next) {
v = edge[i].v ;
if( vis[v] ) continue ;
vis[v] = 1 ;
dep[v] = dep[u] + 1 ;
fa[v] = u ;
dfs1(v) ;
num[u] += num[v] ;
if( son[u] == -1 || ( num[son[u]] < num[v] ) )
son[u] = v ;
}
return ;
}
void dfs2(int u) {
if( son[u] == -1 ) return ;
top[son[u]] = top[u] ;
w[ son[u] ] = step++ ;
vis[ son[u] ] = 1 ;
dfs2(son[u]) ;
int i , v ;
for(i = head[u] ; i != -1 ; i = edge[i].next) {
v = edge[i].v ;
if( vis[v] ) continue ;
vis[v] = 1 ;
top[v] = v ;
w[v] = step++ ;
dfs2(v) ;
}
return ;
}
void dfs() {
memset(vis,0,sizeof(vis)) ;
vis[1] = 1 ; dep[1] = 1 ; fa[1] = -1 ;
dfs1(1) ;
memset(vis,0,sizeof(vis)) ;
vis[1] = 1 ; top[1] = 1 ; step = 1 ;
dfs2(1) ;
return ;
}
void swap1(int rt) {
cl[rt][0] = -cl[rt][0] ;
cl[rt][1] = -cl[rt][1] ;
swap(cl[rt][0],cl[rt][1]) ;
}
void push_up(int_now) {
cl[rt][0] = max(cl[rt<<1][0],cl[rt<<1|1][0]) ;
cl[rt][1] = min(cl[rt<<1][1],cl[rt<<1|1][1]) ;
if( lazy[rt] )
swap1(rt) ;
return ;
}
void push_down(int_now) {
if( lazy[rt] ) {
lazy[rt] = 0 ;
lazy[rt<<1] = 1 - lazy[rt<<1] ;
lazy[rt<<1|1] = 1 - lazy[rt<<1|1] ;
swap1(rt<<1) ;
swap1(rt<<1|1) ;
}
return ;
}
void update1(int k,int x,int_now) {
if( l == k && r == k ) {
lazy[rt] = 0 ;
cl[rt][0] = cl[rt][1] = x ;
return ;
}
push_down(now) ;
if( k <= (l+r)/2 )
update1(k,x,lson) ;
else
update1(k,x,rson) ;
push_up(now) ;
return ;
}
void update2(int ll,int rr,int_now) {
if( ll > r || rr < l ) return ;
if( ll <= l && rr >= r ) {
lazy[rt] = 1 - lazy[rt] ;
swap1(rt) ;
return ;
}
push_down(now) ;
update2(ll,rr,lson) ;
update2(ll,rr,rson) ;
push_up(now) ;
return ;
}
int query(int ll,int rr,int_now,int sum) {
if( ll > r || rr < l ) return -INF ;
if( ll <= l && rr >= r ) {
if( sum%2 ) return -cl[rt][1] ;
else return cl[rt][0] ;
}
sum += lazy[rt] ;
return max( query(ll,rr,lson,sum),query(ll,rr,rson,sum) ) ;
}
void solve(int u,int v,int k) {
int f1 , f2 , ans = -INF ;
while( u != v ) {
if( dep[u] > dep[v] )
swap(u,v) ;
f1 = top[u] ; f2 = top[v] ;
if( f1 == f2 ) {
if( k )
ans = max(ans,query(w[son[u]],w[v],1,step-1,1,0)) ;
else
update2(w[son[u]],w[v],1,step-1,1) ;
v = u ;
}
else if( dep[f1] > dep[f2] ) {
if( k )
ans = max(ans,query(w[f1],w[u],1,step-1,1,0)) ;
else
update2(w[f1],w[u],1,step-1,1) ;
u = fa[f1] ;
}
else{
if( k )
ans = max(ans,query(w[f2],w[v],1,step-1,1,0)) ;
else
update2(w[f2],w[v],1,step-1,1);
v = fa[f2] ;
}
}
if( ans == -INF ) ans = 0 ;
if( k )
printf("%d\n", ans) ;
return ;
}
char str[100] ;
int main() {
int t , i , j , k ;
int u , v , s ;
scanf("%d", &t) ;
while( t-- ) {
memset(head,-1,sizeof(head)) ;
cnt = 0 ;
scanf("%d", &n) ;
for(i = 1 ; i < n ; i++) {
scanf("%d %d %d", &u, &v, &s) ;
add(u,v,s) ;
add(v,u,s) ;
}
dfs() ;
memset(cl,0,sizeof(cl)) ;
memset(lazy,0,sizeof(lazy)) ;
for(i = 0 ; i < n-1 ; i++) {
if( dep[ edge[i*2].u ] > dep[ edge[i*2].v ] )
swap(edge[i*2].u,edge[i*2].v  ) ;
update1(w[ edge[i*2].v ],edge[i*2].w,1,step-1,1) ;
}
while( scanf("%s", str) ) {
if( str[0] == 'D' ) break ;
if( str[0] == 'C' ) {
scanf("%d %d", &i, &k) ;
update1(w[ edge[(i-1)*2].v ],k,1,step-1,1) ;
}
else if( str[0] == 'N' ) {
scanf("%d %d", &i, &j) ;
solve(i,j,0) ;
}
else{
scanf("%d %d", &i, &j) ;
solve(i,j,1) ;
}
}
}
return 0 ;
}


 

这篇关于poj3237--Tree(树链剖分+线段树)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1104109

相关文章

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in

hdu1689(线段树成段更新)

两种操作:1、set区间[a,b]上数字为v;2、查询[ 1 , n ]上的sum 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdl

hdu 1754 I Hate It(线段树,单点更新,区间最值)

题意是求一个线段中的最大数。 线段树的模板题,试用了一下交大的模板。效率有点略低。 代码: #include <stdio.h>#include <string.h>#define TREE_SIZE (1 << (20))//const int TREE_SIZE = 200000 + 10;int max(int a, int b){return a > b ? a :

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while

poj 1127 线段相交的判定

题意: 有n根木棍,每根的端点坐标分别是 px, py, qx, qy。 判断每对木棍是否相连,当他们之间有公共点时,就认为他们相连。 并且通过相连的木棍相连的木棍也是相连的。 解析: 线段相交的判定。 首先,模板中的线段相交是不判端点的,所以要加一个端点在直线上的判定; 然后,端点在直线上的判定这个函数是不判定两个端点是同一个端点的情况的,所以要加是否端点相等的判断。 最后

HDU4737线段树

题目大意:给定一系列数,F(i,j)表示对从ai到aj连续求或运算,(i<=j)求F(i,j)<=m的总数。 const int Max_N = 100008 ;int sum[Max_N<<2] , x[Max_N] ;int n , m ;void push_up(int t){sum[t] = sum[t<<1] | sum[t<<1|1] ;}void upd

zoj 1721 判断2条线段(完全)相交

给出起点,终点,与一些障碍线段。 求起点到终点的最短路。 枚举2点的距离,然后最短路。 2点可达条件:没有线段与这2点所构成的线段(完全)相交。 const double eps = 1e-8 ;double add(double x , double y){if(fabs(x+y) < eps*(fabs(x) + fabs(y))) return 0 ;return x + y ;

圆与线段的交点

poj 3819  给出一条线段的两个端点,再给出n个圆,求出这条线段被所有圆覆盖的部分占了整条线段的百分比。 圆与线段的交点 : 向量AB 的参数方程  P = A + t * (B - A)      0<=t<=1 ; 将点带入圆的方程即可。  注意: 有交点 0 <= t <= 1 ; 此题求覆盖的部分。 则 若求得 t  满足 ; double ask(d

Codeforces 482B 线段树

求是否存在这样的n个数; m次操作,每次操作就是三个数 l ,r,val          a[l] & a[l+1] &......&a[r] = val 就是区间l---r上的与的值为val 。 也就是意味着区间[L , R] 每个数要执行 | val 操作  最后判断  a[l] & a[l+1] &......&a[r] 是否= val import ja