bzoj 2631
lct 基础(' ' ) 就当个纪念吧(' ' ) 毕竟写了4h, cut 部分一直naive 总是想找谁是儿子,然后最后发现直接提根就好了啊(' ' )
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std; typedef long long ll;
const ll mod = 51061;
const ll maxn = 100010; ll pl(ll a, ll b) {
ll ret = a + b;
if(ret >= mod) ret %= mod;
return ret;
} ll mul(ll a, ll b) {
return a * b % mod;
} struct node {
ll ans, lm, lp, lr, size, data, p;
node *son[2], *fa;
}e[maxn]; ll ne = 0; void test(node* x) {
if(!x) return;
cout << x-> data <<" "<<x-> size <<" "<< x-> p << endl;
for(ll i = 0; i < 2; ++ i) test(x-> son[i]);
} void update(node* x) {
x-> ans = x-> data, x-> size = 1;
for(ll i = 0; i < 2; ++ i)
if(x-> son[i])
x-> ans = pl(x-> ans, x-> son[i]-> ans), x-> size = pl(x-> son[i]-> size, x-> size);
} void swap(node* &a, node* &b) {
node* mid = a; a = b, b = mid;
} void pushdown(node* x) {
if(!x || (!x-> lp && !x-> lr && x-> lm == 1)) return;
if(x-> lr) {
swap(x-> son[0], x-> son[1]);
for(ll i = 0; i < 2; ++ i) if(x-> son[i]) x-> son[i]-> lr ^= 1;
x-> lr = 0;
}
for(ll i = 0; i < 2; ++ i) {
if(x-> son[i]) {
x-> son[i]-> ans = pl(mul(x-> son[i]-> ans, x-> lm), mul(x-> son[i]-> size, x-> lp));
x-> son[i]-> data = pl(mul(x-> son[i]->data, x-> lm), x-> lp);
x-> son[i]-> lp = pl(mul(x-> son[i]-> lp, x-> lm), x-> lp);
x-> son[i]-> lm = mul(x-> son[i]-> lm, x-> lm);
}
}
x-> lm = 1, x-> lp = 0;
} void rotate(node* x, ll f) {
node* y = x-> fa;
if(y-> fa) {
if(y-> fa-> son[0] == y) y-> fa-> son[0] = x;
else y-> fa-> son[1] = x;
}
x-> fa = y-> fa; x-> size = y-> size; y-> fa = x; x-> p = y-> p;
x-> ans = y-> ans;
y-> son[f] = x-> son[!f];
if(x-> son[!f]) x-> son[!f]-> fa = y;
x-> son[!f] = y;
update(y);
} void splay(node* x, node* f) {
pushdown(x);
while(x-> fa != f) {
if(x-> fa-> fa == f) {
pushdown(x-> fa-> fa), pushdown(x-> fa), pushdown(x);
ll a = x-> fa-> son[0] == x ? 0 : 1;
rotate(x, a);
}
else {
node *y = x-> fa, *z = y-> fa;
pushdown(z), pushdown(y), pushdown(x);
ll a = z-> son[0] == y ? 0 : 1;
ll b = y-> son[0] == x ? 0 : 1;
if(a == b) rotate(y, a), rotate(x, b);
else rotate(x, b), rotate(x, a);
}
}
} void access(ll cur) {
node* x = e + cur; pushdown(x);
node* y; splay(x, NULL);
if(x-> son[1]) x-> son[1]-> p = cur, x-> son[1]-> fa = NULL, x-> son[1] = NULL;
update(x);
ll pp;
while((pp = x-> p)) {
y = e + pp;
splay(y, NULL);
if(y-> son[1]) y-> son[1]-> p = pp, y-> son[1]-> fa = NULL, y-> son[1] = NULL;
y-> son[1] = x;
x-> fa = y;
update(y);
splay(x, NULL);
}
} void reserve(ll x) {
access(x);
//test(x + e); cout << endl;
(e + x)-> lr ^= 1;
} ll n, m; void link(ll a, ll b) {
access(a);
reserve(b);
(e + b)-> p = a;
update(e + a), update(e + b);
} void cut(ll a, ll b) {
reserve(a), access(b);
//test(b + e); cout << endl;
(e + a)-> p = 0, (e + a)-> fa = NULL, (e + b)-> son[0] = NULL;
update(e + a), update(e + b);
} ll ll_get() {
ll x = 0; char c = (char)getchar();
while(!isdigit(c) && c != '-') c = (char)getchar();
bool f = 0; if(c == '-') f = 1, c = (char)getchar();
while(isdigit(c)) {
x = x * 10 + (ll)(c - '0');
c = (char)getchar();
}
if(f) x = -x;
return x;
} struct edge {
ll t; edge* next;
}se[maxn * 2], *head[maxn]; ll oe = 0; void addedge(ll f, ll t) {
se[oe].t = t, se[oe].next = head[f], head[f] = se + oe ++;
} bool vis[maxn]; void build(ll x, ll pre) {
vis[x] = 1;
(e + x)-> p = pre; (e + x)-> data = (e + x)-> ans = 1; (e + x)-> size = 1; (e + x)-> lm = 1;
for(edge* p = head[x]; p; p = p-> next) {
if(!vis[p-> t]) build(p-> t, x);
}
} void read() {
n = ll_get(), m = ll_get();
for(ll i = 1; i < n; ++ i) {
ll f = ll_get(), t = ll_get();
addedge(f, t), addedge(t, f);
}
memset(vis, 0, sizeof(vis));
build(1, 0);
} void sov() {
char s[10];
// for(int i = 1; i <= n; ++ i) test(e + i), cout << endl;
//cout << "****\n";
while(m --) {
scanf("%s", s + 1);
if(s[1] == '+') {
ll a = ll_get(), b = ll_get(), c = ll_get();
reserve(a); access(b);
node* x = (e + b);
x-> lp = pl(c, x-> lp), x-> ans = pl(x-> ans, mul(c, x-> size)), x-> data = pl(x-> data, c);
}
if(s[1] == '*') {
ll a = ll_get(), b = ll_get(), c = ll_get();
reserve(a);
access(b);
node* x = (e + b);
x-> lm = mul(x-> lm, c), x-> ans = mul(x-> ans, c), x-> data = mul(x-> data, c), x-> lp = mul(x-> lp, c);
}
if(s[1] == '-') {
ll a, b, c, d; a = ll_get(), b = ll_get(), c = ll_get(), d = ll_get();
cut(a, b);
//for(int i = 1; i <= n; ++ i) cout << (e + i)-> p << endl<< endl;
//for(int i = 1; i <= n; ++ i) test(e + i), cout << endl;
//cout << "****\n";
link(c, d); //for(int i = 1; i <= n; ++ i) test(e + i), cout << endl; }
if(s[1] == '/') {
ll a, b; a = ll_get(), b = ll_get();
reserve(a), access(b);
printf("%lld\n", (e + b)-> ans);
}
//for(int i = 1; i <= n; ++ i) test(e + i), cout << endl;
//cout << "****\n";
//cout << (e + 1)-> son[0] << endl;
}
} int main() {
//freopen("test.in", "r", stdin);
//freopen("test.out", "w", stdout);
read();
sov();
}
bzoj 2631的更多相关文章
- [BZOJ - 2631] tree 【LCT】
题目链接:BZOJ - 2631 题目分析 LCT,像线段树区间乘,区间加那样打标记. 这道题我调了一下午. 提交之后TLE了,我一直以为是写错了导致了死循环. 于是一直在排查错误.直到.. 直到我看 ...
- bzoj 2631: tree 动态树+常数优化
2631: tree Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 1716 Solved: 576[Submit][Status] Descrip ...
- BZOJ 2631: tree( LCT )
LCT...略麻烦... -------------------------------------------------------------------------------- #inclu ...
- BZOJ 2631: tree [LCT splay区间]
2631: tree Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 3854 Solved: 1292[Submit][Status][Discus ...
- BZOJ 2631 tree(动态树)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2631 [题目大意] 要求支持链修改,链查询,边断开,连接操作 [题解] 链修改分乘和加 ...
- BZOJ 2631 [国家集训队]Tree II (LCT)
题目大意:给你一棵树,让你维护一个数据结构,支持 边的断,连 树链上所有点点权加上某个值 树链上所有点点权乘上某个值 求树链所有点点权和 (辣鸡bzoj又是土豪题,洛谷P1501传送门) LCT裸题, ...
- BZOJ 2631 Tree ——Link-Cut Tree
[题目分析] 又一道LCT的题目,LCT只能维护链上的信息. [代码] #include <cstdio> #include <cstring> #include <cs ...
- [BZOJ 2631]tree
裸LCT..QAQ写了三遍没写对 真是老了..QAQ 主要错的地方是 init: size[i] = sum[i] = val[i] = mul[i] = 1; pushdown: 注意判断左右儿子是 ...
- bzoj 2631: tree link-cut-tree
题目: Description 一棵n个点的树,每个点的初始权值为1.对于这棵树有q个操作,每个操作为以下四种操作之一: + u v c:将u到v的路径上的点的权值都加上自然数c: - u1 v1 u ...
- 洛谷 1501 [国家集训队]Tree II BZOJ 2631 Tree
[题解] 维护乘法标记和加法标记的LCT #include<cstdio> #include<algorithm> #define Mod (51061) #define N ...
随机推荐
- Maven介绍及安装
1.maven是一个管理第三方库的jar package 2.从该页面下载相应的Maven jar包(http://maven.apache.org/download.cgi),linux OS下 ...
- 【python3】基于scrapyd + scrapydweb 的可视化部署
一.部署组件概览 该部署方式适用于 scrapy项目.scrapy-redis的分布式爬虫项目 需要安装的组件有: 1.scrapyd 服务端 [运行打包后的爬虫代码](所有的爬虫机器都要安 ...
- Java排序算法 [选择、冒泡、快排]
选择排序: 简述:从数组的第一个元素开始,依次与其他所有的元素对比,如果比自身大或小(取决于升序或降序)交换位置. package com.sort; import java.util.Arrays; ...
- spring的 onApplicationEvent方法被执行两次问题
原文地址:http://www.cnblogs.com/a757956132/p/5039438.html 在做web项目开发中,尤其是企业级应用开发的时候,往往会在工程启动的时候做许多的前置检查. ...
- 升级ceph
参考文档 https://blog.51cto.com/jerrymin/2140258 https://www.virtualtothecore.com/en/upgrade-ceph-cluste ...
- Mybatis基于接口注解配置SQL映射器(二)
Mybatis之增强型注解 MyBatis提供了简单的Java注解,使得我们可以不配置XML格式的Mapper文件,也能方便的编写简单的数据库操作代码.但是注解对动态SQL的支持一直差强人意,即使My ...
- Kubernetes v1.16 对API的更改
前段时间安装Kubernetes v1.16.2,然后从v1.14版本的拷贝yaml文件过来执行,很多都报没有相应的api,查看一下新版本的api admissionregistration.k8s. ...
- Nacos Config客户端与Spring Boot、Spring Cloud深度集成
目录 Nacos与Spring Boot集成 @NacosPropertySource和@NacosValue com.alibaba.nacos.spring.core.env.NacosPrope ...
- HTML5: HTML5 介绍
ylbtech-HTML5: HTML5 介绍 1. 什么是 HTML5?返回顶部 HTML5 是下一代 HTML 标准. HTML,HTML 4.01的上一个版本诞生于1999年.自从那儿以后, ...
- mybatis plus的条件构造器
我们在使用条件构造器的时候要使用QueryWrapper或者UpdateWrapper来充当条件语句来进行构造 QueryWrapper(LambdaQueryWrapper) 和 UpdateWra ...