神奇的树上启发式合并 (dsu on tree)
参考资料
https://www.cnblogs.com/zhoushuyu/p/9069164.html
https://www.cnblogs.com/candy99/p/dsuontree.html
https://www.cnblogs.com/zcysky/p/6822395.html
简介
树上启发式合并
用到了heavy−light decomposition树链剖分
把轻边子树的信息合并到重链上的点里
因为每次都是先dfs轻儿子再dfs重儿子,只有重儿子子树的贡献保留,所以可以保证dfs到每颗子树时当前全局维护的信息不会有别的子树里的,和莫队很像
算法实现
1.遍历轻儿子
2.遍历重儿子(保留数据)
3.计算所有轻儿子为根的子树
4.如果当前点是轻儿子,清空影响
复杂度分析
树链剖分后每个点到根的路径上有\(logn\)条轻边和\(logn\)条重链
每个点遇见轻边时合并一次,所以至多\(logn\)次
总复杂度\(O(nlogn)\)
例题
CF600E. Lomsat gelral
http://codeforces.com/contest/600/problem/E
题意:询问每颗子树中出现次数最多的颜色们编号和
板子题
Code
#include<bits/stdc++.h>
#define LL long long
#define RG register
using namespace std;
template<class T> inline void read(T &x) {
x = 0; RG char c = getchar(); bool f = 0;
while (c != '-' && (c < '0' || c > '9')) c = getchar(); if (c == '-') c = getchar(), f = 1;
while (c >= '0' && c <= '9') x = x*10+c-48, c = getchar();
x = f ? -x : x;
return ;
}
template<class T> inline void write(T x) {
if (!x) {putchar(48);return ;}
if (x < 0) x = -x, putchar('-');
int len = -1, z[20]; while (x > 0) z[++len] = x%10, x /= 10;
for (RG int i = len; i >= 0; i--) putchar(z[i]+48);return ;
}
const int N = 1e5+10;
int c[N], n;
struct node {
int to, nxt;
}g[N<<1];
int last[N], gl;
void add(int x, int y) {
g[++gl] = (node) {y, last[x]};
last[x] = gl;
}
int son[N], siz[N];
void dfs1(int u, int f) {
siz[u] = 1;
for (int i = last[u]; i; i = g[i].nxt) {
int v = g[i].to;
if (v == f) continue;
dfs1(v, u);
siz[u] += siz[v];
if (siz[son[u]] < siz[v]) son[u] = v;
}
return ;
}
int num[N], top;
LL sum[N], ans[N];
bool vis[N];
void calc(int u, int fa, int k) {
sum[num[c[u]]] -= c[u];
num[c[u]] += k;
sum[num[c[u]]] += c[u];
if (sum[top + 1]) top++;
else if (!sum[top]) top--;
for (int i = last[u]; i; i = g[i].nxt) {
int v = g[i].to;
if (v == fa || vis[v]) continue;
calc(v, u, k);
}
return ;
}
void dfs(int u, int fa, int op) {
for (int i = last[u]; i; i = g[i].nxt)
if (g[i].to != fa && g[i].to != son[u])
dfs(g[i].to, u, 0);
if (son[u])
dfs(son[u], u, 1), vis[son[u]] = 1;
calc(u, fa, 1); vis[son[u]] = 0;
ans[u] = sum[top];
if (!op) calc(u, fa, -1);
return ;
}
int main() {
read(n);
for (int i = 1; i <= n; i++) read(c[i]);
for (int i = 1; i < n; i++) {
int x, y; read(x), read(y);
add(x, y), add(y, x);
}
dfs1(1, 0);
dfs(1, 0, 1);
for (int i = 1; i <= n; i++)
printf("%I64d ", ans[i]);
return 0;
}
CF570D Tree Requests
http://codeforces.com/problemset/problem/570/D
https://www.luogu.org/problemnew/show/CF570D
构成回文串,奇数个的字母至多一个
用二进制状压判断即可
\(sum[x]\)表示深度为\(x\)构成的状态
#include<bits/stdc++.h>
#define LL long long
#define RG register
using namespace std;
template<class T> inline void read(T &x) {
x = 0; RG char c = getchar(); bool f = 0;
while (c != '-' && (c < '0' || c > '9')) c = getchar(); if (c == '-') c = getchar(), f = 1;
while (c >= '0' && c <= '9') x = x*10+c-48, c = getchar();
x = f ? -x : x;
return ;
}
template<class T> inline void write(T x) {
if (!x) {putchar(48);return ;}
if (x < 0) x = -x, putchar('-');
int len = -1, z[20]; while (x > 0) z[++len] = x%10, x /= 10;
for (RG int i = len; i >= 0; i--) putchar(z[i]+48);return ;
}
const int N = 500010;
struct node {
int to, nxt;
}g[N<<1], q[N];
int last[N], gl;
int n, m;
void add(int x, int y) {
g[++gl] = (node) {y, last[x]};
last[x] = gl;
g[++gl] = (node) {x, last[y]};
last[y] = gl;
}
char s[N];
int siz[N], son[N], val[N], dep[N], sum[N];
bool vis[N];
void dfs1(int u, int fa) {
siz[u] = 1;
for (int i = last[u]; i; i = g[i].nxt) {
int v = g[i].to;
if (v == fa) continue;
dep[v] = dep[u] + 1;
dfs1(v, u);
siz[u] += siz[v];
if (siz[son[u]] < siz[v]) son[u] = v;
}
}
void calc(int u, int fa) {
sum[dep[u]] ^= val[u];
for (int i = last[u]; i; i = g[i].nxt) {
int v = g[i].to;
if (v == fa || vis[v]) continue;
calc(v, u);
}
return ;
}
struct Node {
int h, nxt;
}a[N];
bool ans[N];
int S[N];
bool cal(int x) {
int tmp = 0;
while (x) {
tmp++;
x -= (x & (-x));
}
return tmp <= 1;
}
void dfs(int u, int fa, int op) {
for (int i = last[u]; i; i = g[i].nxt) {
int v = g[i].to;
if (v == fa || son[u] == v) continue;
dfs(v, u, 0);
}
if (son[u]) dfs(son[u], u, 1), vis[son[u]] = 1;
calc(u, fa); vis[son[u]] = 0;
for (int i = S[u]; i; i = a[i].nxt)
ans[i] = cal(sum[a[i].h]);
if (!op) calc(u, fa);
return ;
}
int main() {
read(n), read(m);
for (int i = 2; i <= n; i++) {
int x; read(x);
add(x, i);
}
scanf("%s", s+1);
for (int i = 1; i <= n; i++) val[i] = 1<<(s[i]-'a');
dep[1] = 1;
dfs1(1, 0);
for (int i = 1; i <= m; i++) {
int h, v;
read(v), read(h);
a[i].nxt = S[v];
S[v] = i; a[i].h = h;
}
dfs(1, 0, 1);
for (int i = 1; i <= m; i++) puts(ans[i] ? "Yes" : "No");
return 0;
}
神奇的树上启发式合并 (dsu on tree)的更多相关文章
- 树上启发式合并 (dsu on tree)
这个故事告诉我们,在做一个辣鸡出题人的比赛之前,最好先看看他发明了什么新姿势= =居然直接出了道裸题 参考链接: http://codeforces.com/blog/entry/44351(原文) ...
- 【CF600E】Lomset gelral 题解(树上启发式合并)
题目链接 题目大意:给出一颗含有$n$个结点的树,每个节点有一个颜色.求树中每个子树最多的颜色的编号和. ------------------------- 树上启发式合并(dsu on tree). ...
- dsu on tree 树上启发式合并 学习笔记
近几天跟着dreagonm大佬学习了\(dsu\ on\ tree\),来总结一下: \(dsu\ on\ tree\),也就是树上启发式合并,是用来处理一类离线的树上询问问题(比如子树内的颜色种数) ...
- 树上启发式合并(dsu on tree)学习笔记
有丶难,学到自闭 参考的文章: zcysky:[学习笔记]dsu on tree Arpa:[Tutorial] Sack (dsu on tree) 先康一康模板题吧:CF 600E($Lomsat ...
- dsu on tree (树上启发式合并) 详解
一直都没出过算法详解,昨天心血来潮想写一篇,于是 dsu on tree 它来了 1.前置技能 1.链式前向星(vector 建图) 2.dfs 建树 3.剖分轻重链,轻重儿子 重儿子 一个结点的所有 ...
- 【Luogu U41492】树上数颜色——树上启发式合并(dsu on tree)
(这题在洛谷主站居然搜不到--还是在百度上偶然看到的) 题目描述 给一棵根为1的树,每次询问子树颜色种类数 输入输出格式 输入格式: 第一行一个整数n,表示树的结点数 接下来n-1行,每行一条边 接下 ...
- CF741D Arpa’s letter-marked tree and Mehrdad’s Dokhtar-kosh paths 树上启发式合并(DSU ON TREE)
题目描述 一棵根为\(1\) 的树,每条边上有一个字符(\(a-v\)共\(22\)种). 一条简单路径被称为\(Dokhtar-kosh\)当且仅当路径上的字符经过重新排序后可以变成一个回文串. 求 ...
- 树上启发式合并(dsu on tree)
树上启发式合并属于暴力的优化,复杂度O(nlogn) 主要解决的问题特点在于: 1.对于树上的某些信息进行查询 2.一般问题的解决不包含对树的修改,所有答案可以离线解决 算法思路:这类问题的特点在于父 ...
- hdu6191(树上启发式合并)
hdu6191 题意 给你一棵带点权的树,每次查询 \(u\) 和 \(x\) ,求以 \(u\) 为根结点的子树上的结点与 \(x\) 异或后最大的结果. 分析 看到子树,直接上树上启发式合并,看到 ...
随机推荐
- 回顾2017系列篇(五):人工智能给UI/UX设计带来的影响
如今,我们正处于设计新纪年的转折点上,用机器人和人工智能方面的专家说法表达即“The end is near(终点近了)”.但这并不意味着世界末日,但未来机器人将毫无疑问地接管一部分目前被人类占领的工 ...
- /etc/profile、~/.bash_profile等几个文件的执行过程
/etc/profile.~/.bash_profile等几个文件的执行过程 摘自:http://blog.chinaunix.net/uid-14735472-id-3190130.html 分类 ...
- PBOC电子现金的交易流程
做一个电子现金的交易: 第一步当然是选中当前的应用, 方法是调用select命令, 传入当前的应用AID号, 如果卡片的状态码返回9000,则表示选中成功. 下面举一个例子: 发送: 00 a4 04 ...
- HUST软工1506班第2周作业成绩公布
说明 本次公布的成绩对应的作业为: 第2周个人作业:WordCount编码和测试 如果同学对作业成绩存在异议,在成绩公布的72小时内(截止日期4月26日0点)可以进行申诉,方式如下: 毕博平台的第二周 ...
- srand()、rand()、time()函数的用法
srand()就是给rand()提供种子seed. 如果srand每次输入的数值是一样的,那么每次运行产生的随机数也是一样的. 以一个固定的数值作为种子是一个缺点.通常的做法是 :以这样一句srand ...
- 使用 ServiceStack 构建跨平台 Web 服务(转)
出处:http://www.cnblogs.com/shanyou/p/3348347.html 本文主要来自MSDN杂志<Building Cross-Platform Web Service ...
- Memcached在Windows下的配置和使用(转)
出处:http://www.cnblogs.com/sunniest/p/4154209.html Memcached学习笔记---- 安装和配置 首先,下载Memcached相关文件. 打开控制台, ...
- 深入理解Cookie和Session机制
转载理解Cookie和Session机制 目录 Cookie机制什么是CookieCookie的不可跨域名性Unicode编码:保存中文BASE64编码:保存二进制图片设置Cookie的所有属性Coo ...
- 安卓利用Handlers,AsyncTask和Loaders运行后台程序
安卓的用户界面线程(user interface thread) 1.1 主线程 安卓修改用户界面并从一个单一用户界面线程中处理输入事件,这个线程也被称作主线程(main thread) Androi ...
- jmeter监控内存,CPU等方法
方法1: 使用插件来监控CPU,内存等的使用情况 1.需要的插件准备 JMeterPlugins-Standard-1.4.0.zip , JMeterPlugins-Extras-1.4.0.zip ...