【HDU】5002 Tree 动态树模板题

最新推荐文章于 2015-11-21 23:09:33 发布
原创 最新推荐文章于 2015-11-21 23:09:33 发布 · 980 阅读 · 0 ·
本内容遵循CC 4.0 BY-SA版权协议
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
标签
#HDU
本文详细解析了一道LCT模板题——HDU 5002Tree的解题过程。通过维护set标记、add标记、flip标记，以及最大值、第二大值及其出现次数，作者成功解决该问题并优化了LCT的效率。文章还附带了解题代码，供读者参考和学习。

传送门：【HDU】5002 Tree


题目分析：终于做出来啦~~~~

LCT模板题啊，不过还是很开心~


维护好set标记，add标记，flip标记，维护好最大值，第二大值，以及相应的出现次数就好了~


不过我的LCT效率不高，应该是更新的时候姿势不对，我再去优化优化~


代码如下：


#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std ;

typedef long long LL ;

#define rep( i , a , b ) for ( int i = a ; i < b ; ++ i )
#define For( i , a , b ) for ( int i = a ; i <= b ; ++ i )
#define rev( i , a , b ) for ( int i = a ; i >= b ; -- i )
#define travel( e , H , u ) for ( Edge* e = H[u] ; e ; e = e -> next )
#define clr( a , x ) memset ( a , x , sizeof a )
#define cpy( a , x ) memcpy ( a , x , sizeof a )
#define ls ( o << 1 )
#define rs ( o << 1 | 1 )
#define lson ls , l , m
#define rson rs , m + 1 , r
#define mid ( ( l + r ) >> 1 )
#define root 1 , 1 , n

const int MAXN = 100005 ;
const int INF = 0x3f3f3f3f ;

struct Node* null ;
struct Node {
	Node* c[2] ;
	Node* f ;
	int size ;
	int mmax ;
	int smax ;
	int mmax_num ;
	int smax_num ;
	int v ;
	int flip ;
	int addv ;
	int setv ;
	
	void newnode ( int val ) {
		v = val ;
		mmax = val ;
		smax = -INF ;
		mmax_num = 1 ;
		smax_num = 0 ;
		size = 1 ;
		flip = 0 ;
		addv = 0 ;
		setv = -INF ;
		f = c[0] = c[1] = null ;
	}
	
	void add_val ( int val ) {
		if ( this == null ) return ;
		addv += val ;
		v += val ;
		mmax += val ;
		if ( smax != -INF ) smax += val ;
	}
	
	void set_val ( int val ) {
		if ( this == null ) return ;
		setv = val ;
		addv = 0 ;
		v = val ;
		mmax = val ;
		smax = -INF ;
		mmax_num = size ;
		smax_num = 0 ;
	}
	
	void reverse () {
		if ( this == null ) return ;
		swap ( c[0] , c[1] ) ;
		flip ^= 1 ;
	}
	
	void push_up () {
		size = c[0]->size + 1 + c[1]->size ;
		
		int a[5] = { v , c[0]->mmax , c[1]->mmax , c[0]->smax , c[1]->smax } ;
		sort ( a , a + 5 ) ;
		int cnt = 0 ;
		rep ( i , 1 , 5 ) if ( a[i] != a[cnt] ) a[++ cnt] = a[i] ;
		mmax = a[cnt] ;
		smax = a[cnt - 1] ;
		mmax_num = smax_num = 0 ;
		
		if ( v == mmax ) ++ mmax_num ;
		if ( c[0]->mmax == mmax ) mmax_num += c[0]->mmax_num ;
		if ( c[1]->mmax == mmax ) mmax_num += c[1]->mmax_num ;
		
		if ( smax != -INF ) {
			if ( v == smax ) ++ smax_num ;
			if ( c[0]->mmax == smax ) smax_num += c[0]->mmax_num ;
			if ( c[0]->smax == smax ) smax_num += c[0]->smax_num ;
			if ( c[1]->mmax == smax ) smax_num += c[1]->mmax_num ;
			if ( c[1]->smax == smax ) smax_num += c[1]->smax_num ;
		}
	}
	
	void push_down () {
		if ( flip ) {
			c[0]->reverse () ;
			c[1]->reverse () ;
			flip = 0 ;
		}
		if ( setv != -INF ) {
			c[0]->set_val ( setv ) ;
			c[1]->set_val ( setv ) ;
			setv = -INF ;
		}
		if ( addv ) {
			c[0]->add_val ( addv ) ;
			c[1]->add_val ( addv ) ;
			addv = 0 ;
		}
	}
	
	void link_child ( Node* o , int d ) {
		c[d] = o ;
		o->f = this ;
	}
	
	int is_root () {
		return f == null || this != f->c[0] && this != f->c[1] ;
	}
	
	void sign_down () {
		if ( !is_root () ) f->sign_down () ;
		push_down () ;
	}
	
	void rotate ( int d ) {
		Node* p = f ;
		Node* gp = f->f ;
		p->link_child ( c[d] , !d ) ;
		if ( !p->is_root () ) {
			if ( gp->c[0] == p ) gp->link_child ( this , 0 ) ;
			else gp->link_child ( this , 1 ) ;
		} else f = gp ;
		link_child ( p , d ) ;
		p->push_up () ;
	}
	
	void splay () {
		sign_down () ;
		while ( !is_root () ) {
			if ( f->is_root () ) {
				this == f->c[0] ? rotate ( 1 ) : rotate ( 0 ) ;
			} else {
				if ( f == f->f->c[0] ) {
					this == f->c[0] ? f->rotate ( 1 ) : rotate ( 0 ) ;
					rotate ( 1 ) ;
				} else {
					this == f->c[1] ? f->rotate ( 0 ) : rotate ( 1 ) ;
					rotate ( 0 ) ;
				}
			}
		}
		push_up () ;
	}
	
	void access () {
		Node* o = this ;
		Node* x = null ;
		while ( o != null ) {
			o->splay () ;
			o->link_child ( x , 1 ) ;
			o->push_up () ;
			x = o ;
			o = o->f ;
		}
		splay () ;
	}
	
	void make_root () {
		access () ;
		reverse () ;
	}
	
	Node* find_root () {
		access () ;
		Node* o = this ;
		while ( o->c[0] != null ) {
			o->push_down () ;
			o = o->c[0] ;
		}
		return o ;
	}
	
	void cut () {
		access () ;
		c[0]->f = null ;
		c[0] = null ;
		push_up () ;
	}
	
	void cut ( Node* o ) {
		//if ( o == this || find_root () != o->find_root () ) return ;
		o->make_root () ;
		cut () ;
	}
	
	void link ( Node* o ) {
		//if ( o == this || find_root () == o->find_root () ) return ;
		o->make_root () ;
		o->f = this ;
	}
	
	void add ( Node* o , int val ) {
		//if ( find_root () != o->find_root () ) return ;
		o->make_root () ;
		access () ;
		add_val ( val ) ;
	}
	
	void set ( Node* o , int val ) {
		//if ( find_root () != o->find_root () ) return ;
		o->make_root () ;
		access () ;
		set_val ( val ) ;
	}
	
	void query ( Node* o ) {
		//if ( find_root () != o->find_root () ) return ;
		o->make_root () ;
		access () ;
		if ( !smax_num ) printf ( "ALL SAME\n" ) ;
		else printf ( "%d %d\n" , smax , smax_num ) ;
		//printf ( "%d %d %d %d\n" , mmax , mmax_num , smax , smax_num ) ;
	}
} ;

Node memory_pool[MAXN] ;
Node* cur ;
Node* node[MAXN] ;

int n , q ;

void clear () {
	cur = memory_pool ;
	cur->newnode ( -INF ) ;
	null = cur ++ ;
	null->size = 0 ;
}

void solve () {
	int op , x , y , a , b , val ;
	clear () ;
	scanf ( "%d%d" , &n , &q ) ;
	For ( i , 1 , n ) {
		scanf ( "%d" , &val ) ;
		cur->newnode ( val ) ;
		node[i] = cur ++ ;
	}
	rep ( i , 1 , n ) {
		scanf ( "%d%d" , &x , &y ) ;
		node[x]->link ( node[y] ) ;
	}
	rep ( i , 0 , q ) {
		scanf ( "%d%d%d" , &op , &x , &y ) ;
		switch ( op ) {
			case 1 :
				scanf ( "%d%d" , &a , &b ) ;
				node[x]->cut ( node[y] ) ;
				node[a]->link ( node[b] ) ;
				break ;
			case 2 :
				scanf ( "%d" , &val ) ;
				node[x]->set ( node[y] , val ) ;
				break ;
			case 3 :
				scanf ( "%d" , &val ) ;
				node[x]->add ( node[y] , val ) ;
				break ;
			case 4 :
				node[x]->query ( node[y] ) ;
				break ;
			default : break ;
		}
	}
}

int main () {
	int T , cas = 0 ;
	scanf ( "%d" , &T ) ;
	while ( T -- ) {
		printf ( "Case #%d:\n" , ++ cas ) ;
		solve () ;
	}
	return 0 ;
}

-------Update 2014.09.28--------


优化了push_up，效率快了近一倍~


代码如下：


#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std ;

typedef long long LL ;

#define rep( i , a , b ) for ( int i = a ; i < b ; ++ i )
#define For( i , a , b ) for ( int i = a ; i <= b ; ++ i )
#define rev( i , a , b ) for ( int i = a ; i >= b ; -- i )
#define travel( e , H , u ) for ( Edge* e = H[u] ; e ; e = e -> next )
#define clr( a , x ) memset ( a , x , sizeof a )
#define cpy( a , x ) memcpy ( a , x , sizeof a )
#define ls ( o << 1 )
#define rs ( o << 1 | 1 )
#define lson ls , l , m
#define rson rs , m + 1 , r
#define mid ( ( l + r ) >> 1 )
#define root 1 , 1 , n

const int MAXN = 100005 ;
const int INF = 0x3f3f3f3f ;

struct Node* null ;
struct Node {
	Node* c[2] ;
	Node* f ;
	int size ;
	int mmax ;
	int smax ;
	int mmax_num ;
	int smax_num ;
	int v ;
	int flip ;
	int addv ;
	int setv ;
	
	void newnode ( int val ) {
		v = val ;
		mmax = val ;
		smax = -INF ;
		mmax_num = 1 ;
		smax_num = 0 ;
		size = 1 ;
		flip = 0 ;
		addv = 0 ;
		setv = -INF ;
		f = c[0] = c[1] = null ;
	}
	
	void add_val ( int val ) {
		if ( this == null ) return ;
		addv += val ;
		v += val ;
		mmax += val ;
		if ( smax != -INF ) smax += val ;
	}
	
	void set_val ( int val ) {
		if ( this == null ) return ;
		setv = val ;
		addv = 0 ;
		v = val ;
		mmax = val ;
		smax = -INF ;
		mmax_num = size ;
		smax_num = 0 ;
	}
	
	void reverse () {
		if ( this == null ) return ;
		swap ( c[0] , c[1] ) ;
		flip ^= 1 ;
	}
	
	void update ( int val , int num ) {
		if ( val == -INF ) return ;
		if ( val < smax ) return ;
		if ( val > mmax ) {
			smax = mmax ;
			smax_num = mmax_num ;
			mmax = val ;
			mmax_num = num ;
		} else if ( val == mmax ) {
			mmax_num += num ;
		} else if ( val > smax ) {
			smax = val ;
			smax_num = num ;
		} else smax_num += num ;
	}
	
	void push_up () {
		size = c[0]->size + 1 + c[1]->size ;
		mmax = smax = -INF ;
		mmax_num = smax_num = 0 ;
		update ( v , 1 ) ;
		update ( c[0]->mmax , c[0]->mmax_num ) ;
		update ( c[0]->smax , c[0]->smax_num ) ;
		update ( c[1]->mmax , c[1]->mmax_num ) ;
		update ( c[1]->smax , c[1]->smax_num ) ;
	}
	
	void push_down () {
		if ( flip ) {
			c[0]->reverse () ;
			c[1]->reverse () ;
			flip = 0 ;
		}
		if ( setv != -INF ) {
			c[0]->set_val ( setv ) ;
			c[1]->set_val ( setv ) ;
			setv = -INF ;
		}
		if ( addv ) {
			c[0]->add_val ( addv ) ;
			c[1]->add_val ( addv ) ;
			addv = 0 ;
		}
	}
	
	void link_child ( Node* o , int d ) {
		c[d] = o ;
		o->f = this ;
	}
	
	int is_root () {
		return f == null || this != f->c[0] && this != f->c[1] ;
	}
	
	void sign_down () {
		if ( !is_root () ) f->sign_down () ;
		push_down () ;
	}
	
	void rotate ( int d ) {
		Node* p = f ;
		Node* gp = f->f ;
		p->link_child ( c[d] , !d ) ;
		if ( !p->is_root () ) {
			if ( gp->c[0] == p ) gp->link_child ( this , 0 ) ;
			else gp->link_child ( this , 1 ) ;
		} else f = gp ;
		link_child ( p , d ) ;
		p->push_up () ;
	}
	
	void splay () {
		sign_down () ;
		while ( !is_root () ) {
			if ( f->is_root () ) {
				this == f->c[0] ? rotate ( 1 ) : rotate ( 0 ) ;
			} else {
				if ( f == f->f->c[0] ) {
					this == f->c[0] ? f->rotate ( 1 ) : rotate ( 0 ) ;
					rotate ( 1 ) ;
				} else {
					this == f->c[1] ? f->rotate ( 0 ) : rotate ( 1 ) ;
					rotate ( 0 ) ;
				}
			}
		}
		push_up () ;
	}
	
	void access () {
		Node* o = this ;
		Node* x = null ;
		while ( o != null ) {
			o->splay () ;
			o->link_child ( x , 1 ) ;
			o->push_up () ;
			x = o ;
			o = o->f ;
		}
		splay () ;
	}
	
	void make_root () {
		access () ;
		reverse () ;
	}
	
	Node* find_root () {
		access () ;
		Node* o = this ;
		while ( o->c[0] != null ) {
			o->push_down () ;
			o = o->c[0] ;
		}
		return o ;
	}
	
	void cut () {
		access () ;
		c[0]->f = null ;
		c[0] = null ;
		push_up () ;
	}
	
	void cut ( Node* o ) {
		//if ( o == this || find_root () != o->find_root () ) return ;
		o->make_root () ;
		cut () ;
	}
	
	void link ( Node* o ) {
		//if ( o == this || find_root () == o->find_root () ) return ;
		o->make_root () ;
		o->f = this ;
	}
	
	void add ( Node* o , int val ) {
		//if ( find_root () != o->find_root () ) return ;
		o->make_root () ;
		access () ;
		add_val ( val ) ;
	}
	
	void set ( Node* o , int val ) {
		//if ( find_root () != o->find_root () ) return ;
		o->make_root () ;
		access () ;
		set_val ( val ) ;
	}
	
	void query ( Node* o ) {
		//if ( find_root () != o->find_root () ) return ;
		o->make_root () ;
		access () ;
		if ( !smax_num ) printf ( "ALL SAME\n" ) ;
		else printf ( "%d %d\n" , smax , smax_num ) ;
		//printf ( "%d %d %d %d\n" , mmax , mmax_num , smax , smax_num ) ;
	}
} ;

Node memory_pool[MAXN] ;
Node* cur ;
Node* node[MAXN] ;

int n , q ;

void clear () {
	cur = memory_pool ;
	cur->newnode ( -INF ) ;
	null = cur ++ ;
	null->size = 0 ;
}

void solve () {
	int op , x , y , a , b , val ;
	clear () ;
	scanf ( "%d%d" , &n , &q ) ;
	For ( i , 1 , n ) {
		scanf ( "%d" , &val ) ;
		cur->newnode ( val ) ;
		node[i] = cur ++ ;
	}
	rep ( i , 1 , n ) {
		scanf ( "%d%d" , &x , &y ) ;
		node[x]->link ( node[y] ) ;
	}
	rep ( i , 0 , q ) {
		scanf ( "%d%d%d" , &op , &x , &y ) ;
		switch ( op ) {
			case 1 :
				scanf ( "%d%d" , &a , &b ) ;
				node[x]->cut ( node[y] ) ;
				node[a]->link ( node[b] ) ;
				break ;
			case 2 :
				scanf ( "%d" , &val ) ;
				node[x]->set ( node[y] , val ) ;
				break ;
			case 3 :
				scanf ( "%d" , &val ) ;
				node[x]->add ( node[y] , val ) ;
				break ;
			case 4 :
				node[x]->query ( node[y] ) ;
				break ;
			default : break ;
		}
	}
}

int main () {
	int T , cas = 0 ;
	scanf ( "%d" , &T ) ;
	while ( T -- ) {
		printf ( "Case #%d:\n" , ++ cas ) ;
		solve () ;
	}
	return 0 ;
}


别再为差分时钟头疼了！手把手教你用Xilinx FPGA的IBUFGDS和PLL搞定它
最新发布
weixin_30521161的博客
04-17 343
本文详细解析了Xilinx FPGA中差分时钟的处理方案，重点介绍了IBUFGDS原语和PLL硬核的使用方法。通过实战代码示例和配置技巧，帮助工程师高效实现差分到单端时钟的转换，解决高速设计中的信号完整性问题，并优化时钟管理策略。
动态树模板(HDU4010题)
ACdreamer
05-13 1695
题目：Query on The Trees   #include #include #include #include #include const int maxn=100011; const int INF=0x7fffffff; using namespace std; struct SplayTree { int val,mn,lazy;
BZOJ 2049 动态树-模板题
一花一世界 一叶一菩提
04-05 528
BZOJ 2049:[Sdoi2008]Cave 洞穴勘测   辉辉有一台监测仪器可以实时将通道的每一次改变状况在辉辉手边的终端机上显示：如果监测到洞穴u和洞穴v之间出现了一条通道，终端机上会显示一条指令 Connect u v 。如果监测到洞穴u和洞穴v之间的通道被毁，终端机上会显示一条指令 Destroy u v 。经过长期的艰苦卓绝的手工推算，辉辉发现一个奇怪的现象：无论通道怎么改变，任
Link-Cut Tree(LCT) 模板总结 & 水题/模板题 动态树
とある黄鱼のOIブログ
11-21 2924
1036: [ZJOI2008]树的统计Count Time Limit: 10 Sec  Memory Limit: 162 MB Submit: 9352  Solved: 3783 [Submit][Status][Discuss] Description 一棵树上有n个节点，编号分别为1到n，每个节点都有一个权值w。我们将以下面的形式来要求你对这棵树完成一些操作： I. CH
2014鞍山网络预选赛1006（LCT模板题）hdu5002
cq_phqg
09-13 2015
Tree Time Limit: 16000/8000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 41    Accepted Submission(s): 10 Problem Description You are given a tree wit
【HDU】4677 Query on Graph 动态树
poursoul
12-04 1040
传送门：【HDU】4677 Query on Graph 题目分析： 首先将询问离线。 将所有的边固定形式：左端点小，右端点大。 将边按照右端点从小到大排序，将询问按照右端点从小到大排序。 然后从左到右遍历询问，将边的右端点小于等于询问右端点的都插入到动态树上。 因为动态树要保证是一棵树，于是我们要在加一条边成环的时候去掉环上的一条边来保证动态树依旧是树结构。注意到加边成环时再
hdu 5002 Tree
pkwv
10-03 712
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5002
HDU 5002 Tree（动态树LCT）（2014 ACM/ICPC Asia Regional Anshan Online）
weixin_30699831的博客
09-18 139
Problem Description You are given a tree with N nodes which are numbered by integers 1..N. Each node is associated with an integer as the weight.Your task is to deal with M operations of 4 types:1.De...
LCT - hdu4010 Query on The Trees
Keith_windniw的专栏
09-09 445
题目： http://acm.hdu.edu.cn/showproblem.php?pid=4010 题意： 给一棵树，四种操作 1：添加(x,y)边，若x,y在同一棵树上，则操作非法 2：删除(x,y)边，若x = y 或者 x,y不在同一棵树上，则操作非法 3：x--y的路径上每一个结点的权值增加w，若x,y不在同一棵树上，则操作非法（注意这里给的顺序是w x y 不是
【codeforces】482E. ELCA 动态树
poursoul
11-11 2212
传送门：【codeforces】482E. ELCA
【HDU】5314 Happy King【动态树（点分治）】
poursoul
07-27 1981
传送门：【HDU】5314 Happy King算法分析： 点分治，子树去重类可做。对于一个重心，取出所有路径按照最大值排序，扫下最小值，插入到线段树里维护个数，双指针扫扫解决。然后对所有重心的儿子对应的子树再做一次去重。这个方法我是今早在床上想了下想到的。 我是动态树，动态树比较裸，对点权排序，然后双指针扫扫，接下来就是链接两棵子树或者拆分两棵子树时的子树的siz。子树siz可以用两个变量维护
【codeforces】2014-2015ACM-ICPC CERC 14 Problem J: Pork barrel 【LCT+主席树】
poursoul
12-24 1956
传送门：【codeforces】2014-2015ACM-ICPC CERC 14 Problem J: Pork barrel 题目分析：pushup内写错了一直没发现。。。把maxidx写成idx然后查了一下午没查出来，到晚饭后才发现。。然后再纠正了数组大小，终于AC了。。。写了我一下午，就因为这么一个小错误T  T 今天是平安夜，祝大家平安夜快乐！ 这题是要维护
【HDU】5333 Undirected Graph【LCT+BIT】
poursoul
07-31 1446
传送门：【HDU】5333 Undirected Graphmy  code:my~~code:#pragma comment(linker, "/STACK:1024000000") #include <stdio.h> #include <string.h> #include <map> #include <algorithm> using namespace std ;typedef long
【BZOJ】2049 [Sdoi2008]Cave 洞穴勘测 LCT入门题
poursoul
10-18 1387
传送门：【BZOJ】2049 [Sdoi2008]Cave 洞穴勘测 题目分析：模板题 dai m r
【BZOJ】2631 tree LCT入门题
poursoul
09-27 1120
传送门：【BZOJ】2631 tree 题目分析：LCT模板题~
【SPOJ】 4155. OTOCI LCT
poursoul
11-01 1072
传送门：【SPOJ】 4155. OTOCI 题目分析：dong
【HDU】4010 Query on The Trees 动态树之Link Cut Tree（LCT）
poursoul
09-26 980
传送门：【HDU】4010 Query on The Trees 【HDU】4010 Query on The Trees
【BZOJ】2002 [Hnoi2010]Bounce 弹飞绵羊 LCT入门题
poursoul
09-30 939
传送门：【BZOJ】2002 [Hnoi2010]Bounce 弹飞绵羊
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值