[BZOJ2631]tree 动态树lct
2631: tree
Time Limit: 30 Sec Memory Limit: 128 MB
Submit: 5171 Solved: 1754
[Submit][Status][Discuss]
Description
一棵n个点的树,每个点的初始权值为1。对于这棵树有q个操作,每个操作为以下四种操作之一:
+ u v c:将u到v的路径上的点的权值都加上自然数c;
- u1 v1 u2 v2:将树中原有的边(u1,v1)删除,加入一条新边(u2,v2),保证操作完之后仍然是一棵树;
* u v c:将u到v的路径上的点的权值都乘上自然数c;
/ u v:询问u到v的路径上的点的权值和,求出答案对于51061的余数。
Input
接下来n-1行每行两个正整数u,v,描述这棵树
接下来q行,每行描述一个操作
Output
Sample Input
1 2
2 3
* 1 3 4
/ 1 1
Sample Output
HINT
数据规模和约定
10%的数据保证,1<=n,q<=2000
另外15%的数据保证,1<=n,q<=5*10^4,没有-操作,并且初始树为一条链
另外35%的数据保证,1<=n,q<=5*10^4,没有-操作
100%的数据保证,1<=n,q<=10^5,0<=c<=10^4
Source
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#define maxn 500005
#define ls(x) t[x].s[0]
#define rs(x) t[x].s[1]
#define ll long long
#define mod 51061
using namespace std;
inline ll read() {
ll x=,f=;char ch=getchar();
for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-;
for(;isdigit(ch);ch=getchar()) x=x*+ch-'';
return x*f;
}
struct data {
ll s[],fa,c,a,v,sum,size;
data(){c=;v=,a=,sum=,size=;}
bool rev;
}t[maxn];
bool isroot(int x) {return ls(t[x].fa)!=x&&rs(t[x].fa)!=x;}
void pushup(int x) {t[x].sum=t[ls(x)].sum+t[rs(x)].sum+t[x].v;t[x].sum%=mod;t[x].size=t[ls(x)].size+t[rs(x)].size+;}
void cal(int x,int ce,int add) {
t[x].sum=t[x].sum*ce+add*t[x].size;t[x].sum%=mod;
t[x].v=t[x].v*ce+add;t[x].v%=mod;
t[x].c*=ce;t[x].c%=mod;
t[x].a=t[x].a*ce+add;t[x].a%=mod;
}
void pushdown(int x) {
if(t[x].rev) {
swap(ls(x),rs(x));
t[ls(x)].rev^=;t[rs(x)].rev^=;t[x].rev^=;
}
int ce=t[x].c,add=t[x].a;
t[x].c=;t[x].a=;
cal(ls(x),ce,add);cal(rs(x),ce,add);
}
void rotate(int x) {
int y=t[x].fa,z=t[y].fa;
bool l=ls(y)!=x,r=l^;
if(!isroot(y)) t[z].s[t[z].s[]==y]=x;
t[x].fa=z;t[y].fa=x;t[y].s[l]=t[x].s[r];
t[t[x].s[r]].fa=y;t[x].s[r]=y;
pushup(y);pushup(x);
}
void pre(int x) {
if(!isroot(x)) pre(t[x].fa);
pushdown(x);
}
void splay(int x) {
pre(x);
while(!isroot(x)) {
int y=t[x].fa,z=t[y].fa;
if(!isroot(y)){
if(ls(y)==x^ls(z)==y) rotate(x);
else rotate(y);
}
rotate(x);
}
}
void access(int x) {for(int y=;x;y=x,x=t[x].fa) {splay(x);t[x].s[]=y;pushup(x);}}
void mroot(int x) {access(x);splay(x);t[x].rev^=;}
void link(int x,int y) {mroot(x);t[x].fa=y;}
void cut(int x,int y) {mroot(x);access(y);splay(y);t[y].s[]=t[x].fa=;pushup(y);}
int n,q;
int main() {
n=read(),q=read();
t[].a=t[].c=t[].sum=t[].size=;
for(int i=;i<n;i++) {int u=read(),v=read();link(u,v);}
for(int i=;i<=q;i++) {
char ch[];scanf("%s",ch);
if(ch[]=='+') {int u=read(),v=read();mroot(u);access(v);splay(v);ll tmp=read();cal(v,,tmp);}
if(ch[]=='-') {int u=read(),v=read();cut(u,v);u=read();v=read();link(u,v);}
if(ch[]=='*') {int u=read(),v=read();mroot(u);access(v);splay(v);ll tmp=read();cal(v,tmp,);}
if(ch[]=='/') {int u=read(),v=read();mroot(u);access(v);splay(v);printf("%lld\n",t[v].sum);}
}
return ;
}
/*
5 100
1 2 1 3 3 4 3 5
+ 4 5 3
/ 4 5
* 4 5 3
/ 4 5
- 3 4 4 5
/ 3 4
*/
[BZOJ2631]tree 动态树lct的更多相关文章
- hdu 5398 动态树LCT
GCD Tree Time Limit: 5000/2500 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Su ...
- hdu 5002 (动态树lct)
Tree Time Limit: 16000/8000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submi ...
- LCT总结——概念篇+洛谷P3690[模板]Link Cut Tree(动态树)(LCT,Splay)
为了优化体验(其实是强迫症),蒟蒻把总结拆成了两篇,方便不同学习阶段的Dalao们切换. LCT总结--应用篇戳这里 概念.性质简述 首先介绍一下链剖分的概念(感谢laofu的讲课) 链剖分,是指一类 ...
- HDU 4718 The LCIS on the Tree (动态树LCT)
The LCIS on the Tree Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Oth ...
- LCT(link cut tree) 动态树
模板参考:https://blog.csdn.net/saramanda/article/details/55253627 综合各位大大博客后整理的模板: #include<iostream&g ...
- 动态树LCT小结
最开始看动态树不知道找了多少资料,总感觉不能完全理解.但其实理解了就是那么一回事...动态树在某种意思上来说跟树链剖分很相似,都是为了解决序列问题,树链剖分由于树的形态是不变的,所以可以通过预处理节点 ...
- [模板] 动态树/LCT
简介 LCT是一种数据结构, 可以维护树的动态加边, 删边, 维护链上信息(满足结合律), 单次操作时间复杂度 \(O(\log n)\).(不会证) 思想类似树链剖分, 因为splay可以换根, 用 ...
- 动态树LCT(Link-cut-tree)总结+模板题+各种题目
一.理解LCT的工作原理 先看一道例题: 让你维护一棵给定的树,需要支持下面两种操作: Change x val: 令x点的点权变为val Query x y: 计算x,y之间的唯一的最短路径的点 ...
- BZOJ 2631 tree 动态树(Link-Cut-Tree)
题目大意:维护一种树形数据结构.支持下面操作: 1.树上两点之间的点权值+k. 2.删除一条边.添加一条边,保证加边之后还是一棵树. 3.树上两点之间点权值*k. 4.询问树上两点时间点的权值和. 思 ...
随机推荐
- Promise使用手册
导读 Promise问世已久, 其科普类文章亦不计其数. 遂本篇初衷不为科普, 只为能够温故而知新. 比如说, catch能捕获所有的错误吗? 为什么有些时候会抛出"Uncaught (in ...
- 如何使用 window.open() 处理ajax请求返回的url: 在本页面打开并防止浏览器拦截
ajax请求中用window.open()打开请求返回url(例如实现下载功能时),可能会因为跨域问题导致浏览器拦截 解决办法是:在请求前,打开一个窗口,请求成功后将返回的url直接赋值给该窗口的hr ...
- vue2.0 vue-cli项目中路由之间的参数传递
1.首先配置路由, import Vue from 'vue' import Router from 'vue-router' Vue.use(Router) export default new R ...
- Java进阶
Java进阶(一)Annotation(注解) Java进阶(二)当我们说线程安全时,到底在说什么 Java进阶(三)多线程开发关键技术 Java进阶(四)线程间通信方式对比 Java进阶(五)NIO ...
- 更新协同开发工具SVN的链接的服务器地址
公司内的协同开发工具使用的SVN,因为换了个服务器需要重置SVN地址,一下子有点措手不及. 研究了下SVN的操作菜单,发现有一个功能“重新定位”,应该就是我要找的了,试了一下果真没错,记录下 第一步: ...
- 2017 Multi-University Training Contest - Team 3 RXD and functions(NTT)
题解: 我是参考的 http://blog.csdn.net/qq_32570675/article/details/76571666 这一篇 orz 原来可以这么变换,涨姿势 代码: #includ ...
- NS10.1 产品技术规范
NS10.1 产品技术规范 产品技术规范==================4层-7层流量管理 4层负载均衡(LB) 支持的协议TCP,UDP,FTP,HTTP,HTTPS,DNS(TC ...
- 左侧导航条+中间显示内容+右侧菜单配置,Bootstrap+AdminLTE+Jquery
1.最近做个导航页面,找了一大堆UI,最终选了AdminLTE,这个UI也是以bootstrap为基础,简单实用,中间内容用jquery的load加载,简单暴力,非常适合快速开发. 2.效果图如下: ...
- 【Linux】线程池
首先,线程池是什么?顾名思义,就是把一堆开辟好的线程放在一个池子里统一管理,就是一个线程池. 其次,为什么要用线程池,难道来一个请求给它申请一个线程,请求处理完了释放线程不行么?也行,但是如果创建线程 ...
- The XOR Largest Pair [Trie]
描述 在给定的N个整数A1,A2--AN中选出两个进行xor运算,得到的结果最大是多少? 输入格式 第一行一个整数N,第二行N个整数A1-AN. 输出格式 一个整数表示答案. 样例输入 3 1 2 3 ...