[HNOI2016]网络 [树链剖分,可删除堆]
考虑在 |不在| 这条链上的所有点上放上一个 \(x\),删除也是,然后用可删除堆就随便草掉了。
// powered by c++11
// by Isaunoya
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
#include <bits/stdc++.h>
#define rep(i, x, y) for (register int i = (x); i <= (y); ++i)
#define Rep(i, x, y) for (register int i = (x); i >= (y); --i)
using namespace std;
using db = double;
using ll = long long;
using uint = unsigned int;
using pii = pair<int, int>;
#define ve vector
#define Tp template
#define all(v) v.begin(), v.end()
#define sz(v) ((int)v.size())
#define pb emplace_back
#define fir first
#define sec second
// the cmin && cmax
Tp<class T> void cmax(T& x, const T& y) {
if (x < y) x = y;
}
Tp<class T> void cmin(T& x, const T& y) {
if (x > y) x = y;
}
// sort , unique , reverse
Tp<class T> void sort(ve<T>& v) { sort(all(v)); }
Tp<class T> void unique(ve<T>& v) {
sort(all(v));
v.erase(unique(all(v)), v.end());
}
Tp<class T> void reverse(ve<T>& v) { reverse(all(v)); }
const int SZ = 0x191981;
struct FILEIN {
char qwq[SZ], *S = qwq, *T = qwq, ch;
char GETC() { return (S == T) && (T = (S = qwq) + fread(qwq, 1, SZ, stdin), S == T) ? EOF : *S++; }
FILEIN& operator>>(char& c) {
while (isspace(c = GETC()))
;
return *this;
}
FILEIN& operator>>(string& s) {
while (isspace(ch = GETC()))
;
s = ch;
while (!isspace(ch = GETC())) s += ch;
return *this;
}
Tp<class T> void read(T& x) {
bool sign = 1;
while ((ch = GETC()) < 0x30)
if (ch == 0x2d) sign = 0;
x = (ch ^ 0x30);
while ((ch = GETC()) > 0x2f) x = x * 0xa + (ch ^ 0x30);
x = sign ? x : -x;
}
FILEIN& operator>>(int& x) { return read(x), *this; }
FILEIN& operator>>(unsigned& x) { return read(x), *this; }
} in;
struct FILEOUT {
const static int LIMIT = 0x114514;
char quq[SZ], ST[0x114];
signed sz, O;
~FILEOUT() { flush(); }
void flush() {
fwrite(quq, 1, O, stdout);
fflush(stdout);
O = 0;
}
FILEOUT& operator<<(char c) { return quq[O++] = c, *this; }
FILEOUT& operator<<(string str) {
if (O > LIMIT) flush();
for (char c : str) quq[O++] = c;
return *this;
}
Tp<class T> void write(T x) {
if (O > LIMIT) flush();
if (x < 0) {
quq[O++] = 0x2d;
x = -x;
}
do {
ST[++sz] = x % 0xa ^ 0x30;
x /= 0xa;
} while (x);
while (sz) quq[O++] = ST[sz--];
return;
}
FILEOUT& operator<<(int x) { return write(x), *this; }
FILEOUT& operator<<(unsigned x) { return write(x), *this; }
} out;
int n, m;
struct P {
priority_queue<int> a, b;
inline void push(int x) { a.push(x); }
inline void erase(int x) { b.push(x); }
inline int top() {
while (!a.empty() && !b.empty() && a.top() == b.top()) a.pop(), b.pop();
if (a.empty()) return -1;
return a.top();
}
};
const int maxn = 1e5 + 51;
P s[maxn << 2];
void upd(int a, int b, int l, int r, int rt, int x) {
if (a <= l && r <= b) {
s[rt].push(x);
return;
}
int mid = l + r >> 1;
if (a <= mid) upd(a, b, l, mid, rt << 1, x);
if (b > mid) upd(a, b, mid + 1, r, rt << 1 | 1, x);
}
void erase(int a, int b, int l, int r, int rt, int x) {
if (a <= l && r <= b) {
s[rt].erase(x);
return;
}
int mid = l + r >> 1;
if (a <= mid) erase(a, b, l, mid, rt << 1, x);
if (b > mid) erase(a, b, mid + 1, r, rt << 1 | 1, x);
}
int ans = -1;
void qry(int l, int r, int x, int rt) {
cmax(ans, s[rt].top());
if (l == r) return;
int mid = l + r >> 1;
(x <= mid) ? qry(l, mid, x, rt << 1) : qry(mid + 1, r, x, rt << 1 | 1);
}
vector<int> g[maxn];
int sz[maxn], son[maxn], d[maxn], fa[maxn];
void dfs(int u) {
sz[u] = 1;
for (int v : g[u])
if (v ^ fa[u]) {
fa[v] = u, d[v] = d[u] + 1, dfs(v), sz[u] += sz[v];
if (sz[v] > sz[son[u]]) son[u] = v;
}
}
int top[maxn], dfn[maxn], idx = 0;
void dfs(int u, int t) {
dfn[u] = ++idx, top[u] = t;
if (son[u]) dfs(son[u], t);
for (int v : g[u])
if (v ^ fa[u] && v ^ son[u]) dfs(v, v);
}
vector<pii> v, vec;
void swap(int& x, int& y) { x ^= y ^= x ^= y; }
void upd(int x, int y, int c) {
v.clear(), vec.clear();
while (top[x] ^ top[y]) {
if (d[top[x]] < d[top[y]]) swap(x, y);
v.push_back({ dfn[top[x]], dfn[x] });
x = fa[top[x]];
}
if (d[x] > d[y]) swap(x, y);
v.push_back({ dfn[x], dfn[y] });
sort(v);
int las = 1;
for (pii x : v) {
if (x.first > las) vec.push_back({ las, x.first - 1 });
las = x.second + 1;
}
if (las <= n) vec.push_back({ las, n });
for (pii x : vec) upd(x.first, x.second, 1, n, 1, c);
}
void erase(int x, int y, int c) {
v.clear(), vec.clear();
while (top[x] ^ top[y]) {
if (d[top[x]] < d[top[y]]) swap(x, y);
v.push_back({ dfn[top[x]], dfn[x] });
x = fa[top[x]];
}
if (d[x] > d[y]) swap(x, y);
v.push_back({ dfn[x], dfn[y] });
sort(v);
int las = 1;
for (pii x : v) {
if (x.first > las) vec.push_back({ las, x.first - 1 });
las = x.second + 1;
}
if (las <= n) vec.push_back({ las, n });
for (pii x : vec) erase(x.first, x.second, 1, n, 1, c);
}
const int maxm = 2e5 + 52;
int A[maxm], B[maxm], V[maxm];
signed main() {
#ifdef _WIN64
freopen("testdata.in", "r", stdin);
#else
ios_base ::sync_with_stdio(false);
cin.tie(nullptr), cout.tie(nullptr);
#endif
// code begin.
in >> n >> m;
rep(i, 2, n) {
int u, v;
in >> u >> v, g[u].pb(v), g[v].pb(u);
}
dfs(1), dfs(1, 1);
rep(i, 1, m) {
int op;
in >> op;
if (op == 0) {
int x, y, c;
in >> x >> y >> c;
upd(x, y, c);
A[i] = x, B[i] = y, V[i] = c;
}
if (op == 1) {
int t;
in >> t, erase(A[t], B[t], V[t]);
}
if (op == 2) {
int u;
in >> u, ans = -1, qry(1, n, dfn[u], 1);
out << ans << '\n';
}
}
return 0;
// code end.
}
[HNOI2016]网络 [树链剖分,可删除堆]的更多相关文章
- BZOJ4538:[HNOI2016]网络(树链剖分,堆)
Description 一个简单的网络系统可以被描述成一棵无根树.每个节点为一个服务器.连接服务器与服务器的数据线则看做 一条树边.两个服务器进行数据的交互时,数据会经过连接这两个服务器的路径上的所有 ...
- [HNOI2016]网络 树链剖分,堆
[HNOI2016]网络 LG传送门 表示乱搞比正解难想. 整体二分很好想吧. 但是为了好写快乐,我们选择三个\(\log\)的乱搞. 先树剖,线段树套堆维护区间最大值.对于一次修改,如果是插入,就把 ...
- uoj30【CF Round #278】Tourists(圆方树+树链剖分+可删除堆)
- 学习了一波圆方树 学习了一波点分治 学习了一波可删除堆(巧用 ? STL) 传送门: Icefox_zhx 注意看代码看怎么构建圆方树的. tips:tips:tips:圆方树内存记得开两倍 CO ...
- luoguP3250 [HNOI2016]网络 树链剖分 + 堆
机房某大佬告诉我,一条链在全局线段树中的区间最多有$log$段 因此同样的,代表不在这条链上的区间同样只有$log$段 对这$log$段区间进行维护即可 为了能够删除,在线段树的每个节点暴力维护一个堆 ...
- BZOJ5210 最大连通子块和 【树链剖分】【堆】【动态DP】
题目分析: 解决了上次提到的<切树游戏>后,这道题就是一道模板题. 注意我们需要用堆维护子重链的最大值.这样不会使得复杂度变坏,因为每个重链我们只考虑一个点. 时间复杂度$O(nlog^2 ...
- 【BZOJ4538】【HNOI2016】网络(树链剖分,线段树,堆)
题目链接,我是真的懒得调题目的格式... 题解 树链剖分搞一下LCA 把线段树弄出来 这只是形式上的线段树 本质上是维护一段区间的一个堆 每次把堆插入节点, 询问的时候查询线段树上的堆的最大值就行了 ...
- 2019.01.13 bzoj4538: [Hnoi2016]网络(树链剖分)
传送门 树链剖分一眼题. 题意简述: 给定一棵树,有三种操作: 加入一条路径 删除一条已加入的路径 询问不过一个点x的路径的最大值. 思路: 直接树链剖分维护答案. 因为询问的事不过点xxx的最大值, ...
- 【bzoj5210】最大连通子块和 树链剖分+线段树+可删除堆维护树形动态dp
题目描述 给出一棵n个点.以1为根的有根树,点有点权.要求支持如下两种操作: M x y:将点x的点权改为y: Q x:求以x为根的子树的最大连通子块和. 其中,一棵子树的最大连通子块和指的是:该子树 ...
- 焦作网络赛E-JiuYuanWantstoEat【树链剖分】【线段树】
You ye Jiu yuan is the daughter of the Great GOD Emancipator. And when she becomes an adult, she wil ...
随机推荐
- t-io Java构建p2p网络
Java 构建p2p网络 这篇文章是一篇关于pbft算法实现的一篇补充文章,但是在这里不会涉及pbft的算法方面,所以可以当做一篇单独的文章食用.如果想查看关于区块链或者PBFT算法的文章,可以参考一 ...
- windows 安裝 gcc 編譯器 - MinGW
MinGW 介紹 MinGW 全稱 Minimalist GNU For Windows,是個精簡的Windows平台C/C++.ADA及Fortran編譯器,相比Cygwin而言,體積要小很多,使用 ...
- 深入理解JVM-类加载及类加载器
深入理解JVM 2020年02月06日22:43:09 - 记录学习过程 终于开始了.在学习这个之前,看了zhanglong老师的 java 8 和springboot 迫不及待了.先开始吧. 写在前 ...
- Sqli-Labs 闯关 less 42-49
Less 42 这一关一进去看着像前面的二次注入.发现也注入不了.. 我们观察代码发现这一关用的是堆叠注入. 登陆的这里可以看到login_password登陆的时候并没有使用mysqli_real_ ...
- XDOJ
1000.a+b. #include<bits/stdc++.h> using namespace std; int a,b; int main() { ios::sync_with_st ...
- WeChall_Encodings: URL (Training, Encoding)
Your task is to decode the following: %59%69%70%70%65%68%21%20%59%6F%75%72%20%55%52%4C%20%69%73%20%6 ...
- 给Hangfire的webjob增加callback和动态判断返回结果功能设计
背景介绍 通常业务中需要用到定时执行功能,我用hangfire搭建了一个调度服务,这个调度服务是独立于业务逻辑的,具体可以参考文章:https://github.com/yuzd/Hangfire.H ...
- LNK2019
原因:inline函数被外部文件的函数调用时,必须将inline函数定义在头文件中,不能定义在cpp文件中.
- VMware复制Linux虚拟机后网络配置
1.启动虚拟机,点击我已复制 点击已复制后,VMware将会为重置虚拟机网卡mac地址. 2.修改网卡mac地址 3.ifconfig显示网卡名称与配置不一致处理 Centos6中ifconfig显示 ...
- 20200104--python学习第六天
今日学习 集合 内存相关知识 深浅拷贝 内容回顾与补充 (1)列表: (a)reverse 反转 v1=[1,2,3111,32,13] print(v1) v1.reverse() print(v1 ...