裸LCT。。QAQ写了三遍没写对

真是老了。。QAQ

主要错的地方是

init:

size[i] = sum[i] = val[i] = mul[i] = 1;

pushdown:

注意判断左右儿子是否为空

splay:

前面有pushdown, stack..

while(!isroot(p)){

int x = fa[p], y = fa[x];

if(!isroot(x))...

}

Access:

while(u){

...

c[u][] = t;

u = fa[u];
}

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 100010
using namespace std; typedef long long ll; const int md = 51061; int n, q;
int c[maxn][2], fa[maxn], val[maxn], mul[maxn], add[maxn];
int sum[maxn], top, st[maxn], size[maxn];
bool rev[maxn];
namespace splay{
inline void Mul(int& x, ll y){
x = x * y % md;
}
inline void Add(int& x, ll y){
x = (x + y) % md;
} void init(){
for(int i = 1; i <= n; i ++)
val[i] = mul[i] = 1, add[i] = 0;
} inline void pushup(int x){
sum[x] = (sum[c[x][0]] + sum[c[x][1]] + val[x]) % md;
size[x] = size[c[x][0]] + size[c[x][1]] + 1;
} inline void pushdown(int x){
if(rev[x]){
if(c[x][0])rev[c[x][0]] ^= 1;
if(c[x][1])rev[c[x][1]] ^= 1;
swap(c[x][0], c[x][1]);
rev[x] = 0;
} int lc = c[x][0], rc = c[x][1];
if(lc != 0){
if(mul[x] != 1){
Mul(mul[lc], mul[x]);
Mul(add[lc], mul[x]);
Mul(val[lc], mul[x]);
Mul(sum[lc], mul[x]);
}
if(add[x] != 0){
Add(add[lc], add[x]);
Add(val[lc], add[x]);
Add(sum[lc], (ll)size[lc] * add[x]);
}
} if(rc != 0){
if(mul[x] != 1){
Mul(mul[rc], mul[x]);
Mul(add[rc], mul[x]);
Mul(val[rc], mul[x]);
Mul(sum[rc], mul[x]);
}
if(add[x] != 0){
Add(add[rc], add[x]);
Add(val[rc], add[x]);
Add(sum[rc], (ll)size[rc] * add[x]);
}
}
mul[x] = 1, add[x] = 0;
} inline bool isroot(int x){
return c[fa[x]][0] != x && c[fa[x]][1] != x;
} inline void rotate(int p, int x){
int mark = p == c[x][1], y = c[p][mark^1], z = fa[x];
if(c[z][0] == x) c[z][0] = p;
if(c[z][1] == x) c[z][1] = p;
if(y) fa[y] = x;
c[p][mark^1] = x, fa[p] = z;
c[x][mark] = y, fa[x] = p;
pushup(x);
} inline void splay(int p){
st[top = 1] = p;
for(int i = p; !isroot(i); i = fa[i])
st[++ top] = fa[i];
while(top) pushdown(st[top --]);
while(!isroot(p)){
int x = fa[p], y = fa[x];
if(isroot(x))rotate(p, x);
else if(p == c[x][0] ^ x == c[y][0])
rotate(p, x), rotate(p, y);
else rotate(x, y), rotate(p, x);
}
pushup(p);
} inline void update(int x, int mu, int ad){
if(mu != 1){
Mul(val[x], mu);
Mul(mul[x], mu);
Mul(add[x], mu);
Mul(sum[x], mu);
} if(ad != 0){
Add(val[x], ad);
Add(add[x], ad);
Add(sum[x], (ll)size[x] * ad);
}
}
} namespace LCT{
void Access(int u){
int t = 0;
while(u){
splay::splay(u);
c[u][1] = t;
t = u;
splay::pushup(u);
u = fa[u];
}
} void Evert(int u){
Access(u);
splay::splay(u);
rev[u] ^= 1;
splay::pushdown(u);
} void link(int u, int v){
Evert(u);
fa[u] = v;
//Access(u); splay::splay(u);
} void cut(int u, int v){
Evert(u);
Access(v);
splay::splay(v);
c[v][0] = fa[u] = 0;
splay::pushup(v);
} void update(int u, int v, int mu, int ad){
Evert(u), Access(v);
splay::splay(v);
splay::update(v, mu, ad);
} int ask(int u, int v){
Evert(u), Access(v);
splay::splay(v);
return sum[v];
}
} int main(){
scanf("%d%d", &n, &q);
splay::init();
int u, v, u1, v1, c;
for(int i = 1; i < n; i ++){
scanf("%d%d", &u, &v);
LCT::link(u, v);
}
for(int i = 1; i <= q; i ++){
char ch = getchar();
for(; ch < '!'; ch = getchar());
if(ch == '+'){
scanf("%d%d%d", &u, &v, &c);
LCT::update(u, v, 1, c);
}
if(ch == '-'){
scanf("%d%d%d%d", &u, &v, &u1, &v1);
LCT::cut(u, v);
LCT::link(u1, v1);
}
if(ch == '*'){
scanf("%d%d%d", &u, &v, &c);
LCT::update(u, v, c, 0);
}
if(ch == '/'){
scanf("%d%d", &u, &v);
printf("%d\n", LCT::ask(u, v));
}
}
return 0;
}

[BZOJ 2631]tree的更多相关文章

  1. bzoj 2631: tree 动态树+常数优化

    2631: tree Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 1716  Solved: 576[Submit][Status] Descrip ...

  2. [BZOJ - 2631] tree 【LCT】

    题目链接:BZOJ - 2631 题目分析 LCT,像线段树区间乘,区间加那样打标记. 这道题我调了一下午. 提交之后TLE了,我一直以为是写错了导致了死循环. 于是一直在排查错误.直到.. 直到我看 ...

  3. BZOJ 2631: tree( LCT )

    LCT...略麻烦... -------------------------------------------------------------------------------- #inclu ...

  4. BZOJ 2631: tree [LCT splay区间]

    2631: tree Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 3854  Solved: 1292[Submit][Status][Discus ...

  5. BZOJ 2631 tree(动态树)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2631 [题目大意] 要求支持链修改,链查询,边断开,连接操作 [题解] 链修改分乘和加 ...

  6. BZOJ 2631 tree 动态树(Link-Cut-Tree)

    题目大意:维护一种树形数据结构.支持下面操作: 1.树上两点之间的点权值+k. 2.删除一条边.添加一条边,保证加边之后还是一棵树. 3.树上两点之间点权值*k. 4.询问树上两点时间点的权值和. 思 ...

  7. BZOJ 2631 Tree ——Link-Cut Tree

    [题目分析] 又一道LCT的题目,LCT只能维护链上的信息. [代码] #include <cstdio> #include <cstring> #include <cs ...

  8. bzoj 2631: tree link-cut-tree

    题目: Description 一棵n个点的树,每个点的初始权值为1.对于这棵树有q个操作,每个操作为以下四种操作之一: + u v c:将u到v的路径上的点的权值都加上自然数c: - u1 v1 u ...

  9. 洛谷 1501 [国家集训队]Tree II BZOJ 2631 Tree

    [题解] 维护乘法标记和加法标记的LCT #include<cstdio> #include<algorithm> #define Mod (51061) #define N ...

随机推荐

  1. BZOJ 2818

    2818:GCD Description 给定整数$N$,求$1\le x,y\le N$且$\gcd{x,y}$为素数的数对$(x,y)$有多少对. Input $N$ Output RT Samp ...

  2. Android中自定义Activity和Dialog的位置大小背景和透明度等demo

    1.自定义Activity显示样式 先在res/values下建colors.xml文件,写入: <?xml version="1.0" encoding="utf ...

  3. Java基础算法集50题

    最近因为要准备实习,还有一个蓝桥杯的编程比赛,所以准备加强一下算法这块,然后百度了一下java基础算法,看到的都是那50套题,那就花了差不多三个晚自习的时间吧,大体看了一遍,做了其中的27道题,有一些 ...

  4. 尖刀出鞘的display常用属性及css盒模型深入研究

    一:diplay:inline-block 含义:指元素创建了一个行级的块级元素,该元素内部(内容)被格式化成一个块级元素,同时元素本身则被格式化成一个行内元素.更简单的说就是说inline-bloc ...

  5. Rehashing

    The size of the hash table is not determinate at the very beginning. If the total size of keys is to ...

  6. Django的Context和RequestContext

    参考:http://www.dannysite.com/blog/38/ Django的模板渲染中,Context可以用来传递数据,一个Context是一系列变量和值的集合,它和Python的字典有点 ...

  7. 【转】利用mybatis-generator自动生成代码

    本文转自:http://www.cnblogs.com/yjmyzz/p/4210554.html mybatis-generator有三种用法:命令行.eclipse插件.maven插件.个人觉得m ...

  8. DP:Skiing(POJ 1088)

     北大教你怎么滑雪 题目是中文的,要读懂题目其实不难 其实这道题挺经典的,我们这样想,他最终要找到一个最大值,这个时候我们就想到要用动态规划 那怎么用呢?我们同时这样想,由于滑雪的最高点我们不能立马找 ...

  9. HDU 1848 Fibonacci again and again (斐波那契博弈SG函数)

    Fibonacci again and again Time Limit: 1000MS   Memory Limit: 32768KB   64bit IO Format: %I64d & ...

  10. WINDOWS xp用户账户怎么没了怎么办?

    这是因为系统的一个默认设置!新建用户会把管理员用户隐藏!只是修改了里面的注册表!在运行(windows徽标+R)里输入“regedit”可以打开注册表编辑器,定位到“HKEY_LOCAL_MACHIN ...