dsu on tree题表
dsu on tree,又名树上启发式合并、重链剖分,是一类十分实用的trick,它常常可以作为一些正解的替代算法:
1.DFS序+线段树/主席树/线段树合并
2.对DFS序分块的树上莫队
3.长链剖分(但复杂度会多一个log)
4.点分治(通常可以做有根树的点分治)
重链剖分的概念,用一个DFS找到每个点最大的一个儿子,作为它的重儿子,并将它标记。则从上到下一段连续的标记点就成为一条重链。
重链剖分有一个常用的性质:每个点到根的路径上,至多经过$O(\log n)$条重链。点分治、树链剖分都用到了这个性质。
dsu on tree 是一个优化后的暴力,主要优化的地方在于它先递归轻子树并消除影响,后递归重子树并保留影响。之后再计算该节点需要的信息。
它可以解决大部分无修改子树查询问题,需要问题满足以下几个条件:
1.有一个树上的$O(n^2)$暴力算法。
2.从轻子树合并上来的复杂度是线性的。
3.从重子树合并上来的复杂度是$O(1)$的。
4.可以在O(子树大小)时间内清空递归后的数组(也就是线性撤销所有影响)
dsu on tree的流程:
1.递归到所有轻儿子并消除影响。
2.递归到重儿子并保留影响。
3.递归所有轻儿子计算子树内除重子树之外的点对当前点答案的影响。
4.若此点不是父亲的重儿子则消除子树内所有影响(即将数组清空)。
另外,dsu on tree的题,只要使用的暴力数组是以深度为下标的,几乎都可以被长链剖分替代且复杂度少一个log。如下面的例2,3。
例一:[CF600E]询问每棵子树中出现次数最多的颜色(可能不只一个)的编号和。
首先考虑暴力算法,对每个点x DFS下去,得到一个计数器数组co[i]表示x的子树内颜色i的点的个数,同时维护x的答案。
然后重链剖分,按流程做即可,具体模板见代码。
#include<cstdio>
#include<algorithm>
#include<iostream>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
#define For(i,x) for (int i=h[x],k; i; i=nxt[i])
typedef long long ll;
using namespace std; const int N=;
int n,u,v,cnt,tot[N],mx,col[N],sz[N],son[N],h[N],to[N],nxt[N];
ll ans[N],sm;
bool skip[N]; void add(int u,int v){ to[++cnt]=v; nxt[cnt]=h[u]; h[u]=cnt; } void get(int x,int fa){
sz[x]=;
For(i,x) if ((k=to[i])!=fa){
get(k,x); sz[x]+=sz[k];
if (sz[k]>sz[son[x]]) son[x]=k;
}
} void dfs(int x,int fa,int op){
tot[col[x]]+=op;
if (op> && tot[col[x]]>=mx){
if (tot[col[x]]>mx) sm=,mx=tot[col[x]];
sm+=col[x];
}
For(i,x) if ((k=to[i])!=fa && !skip[k]) dfs(k,x,op);
} void work(int x,int fa,bool cl){
For(i,x) if ((k=to[i])!=fa && k!=son[x]) work(k,x,);
if (son[x]) work(son[x],x,),skip[son[x]]=;
dfs(x,fa,); ans[x]=sm; skip[son[x]]=;
if (cl) dfs(x,fa,-),mx=sm=;
} int main(){
scanf("%d",&n);
rep(i,,n) scanf("%d",&col[i]);
rep(i,,n) scanf("%d%d",&u,&v),add(u,v),add(v,u);
get(,); work(,,);
rep(i,,n) cout<<ans[i]<<' ';
return ;
}
CF600E
例二:[CF570D]一棵树,询问某棵子树指定深度的点能否构成回文。
用二进制记录每个深度的每个字母出现次数的奇偶性即可。
#include<cstdio>
#include<vector>
#include<algorithm>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
#define For(i,x) for (int i=h[x],k; i; i=nxt[i])
using namespace std; const int N=;
char s[N];
int n,m,x,y,cnt,skip,a[N],d[N],ans[N],fa[N],sz[N],son[N],to[N],nxt[N],h[N];
struct P{ int x,y; };
vector<P>ve[N]; void add(int u,int v){ to[++cnt]=v; nxt[cnt]=h[u]; h[u]=cnt; } void dfs(int x){
d[x]=d[fa[x]]+; sz[x]=;
For(i,x){
dfs(k=to[i]); sz[x]+=sz[k];
if (sz[k]>sz[son[x]]) son[x]=k;
}
} void get(int x){
a[d[x]]^=<<(s[x]-'a');
For(i,x) if ((k=to[i])!=skip) get(k);
} void work(int x,int op){
For(i,x) if ((k=to[i])!=son[x]) work(k,);
if (son[x]) work(son[x],),skip=son[x];
get(x); skip=; int ed=ve[x].size()-;
rep(i,,ed) ans[ve[x][i].x]=__builtin_popcount(a[ve[x][i].y])<=;
if (op) get(x);
} int main(){
freopen("570D.in","r",stdin);
freopen("570D.out","w",stdout);
scanf("%d%d",&n,&m);
rep(i,,n) scanf("%d",&fa[i]),add(fa[i],i);
rep(i,,n) scanf(" %c",&s[i]);
rep(i,,m) scanf("%d%d",&x,&y),ve[x].push_back((P){i,y});
dfs(); work(,);
rep(i,,m) if (ans[i]) puts("Yes"); else puts("No");
return ;
}
CF570D
例三:
[CF246E]一个森林,求k级后代中多少种不同的权值。
[CF208E]给出一个森林,求和一个点有相同k级祖先的点有多少。
第二个问题倍增找到k级祖先后就是弱化版的第一个问题了。
与例二类似,对每个深度维护一个set记录出现过的权值数目。
注意由于流程的最后一步是清空整个数组,所以删除的时候不需要考虑记录次数的问题。
#include<set>
#include<map>
#include<string>
#include<cstdio>
#include<vector>
#include<algorithm>
#include<iostream>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
#define For(i,x) for (int i=h[x],k; i; i=nxt[i])
using namespace std; const int N=;
string s;
map<string,int>mp;
int n,Q,x,y,tot,skip,cnt,rts,id[N],fa[N],rt[N],ans[N],sz[N];
int son[N],d[N],h[N],to[N<<],nxt[N<<];
set<int>S[N];
struct P{ int x,y; };
vector<P>ve[N];
void add(int u,int v){ to[++cnt]=v; nxt[cnt]=h[u]; h[u]=cnt; } void dfs(int x){
d[x]=d[fa[x]]+; sz[x]=;
For(i,x){
dfs(k=to[i]); sz[x]+=sz[k];
if (sz[k]>sz[son[x]]) son[x]=k;
}
} void get(int x,int op){
if (op) S[d[x]].insert(id[x]); else S[d[x]].erase(id[x]);
For(i,x) if ((k=to[i])!=skip) get(k,op);
} void work(int x,int op){
For(i,x) if ((k=to[i])!=son[x]) work(k,);
if (son[x]) work(son[x],),skip=son[x];
get(x,); skip=; int ed=ve[x].size()-;
rep(i,,ed) ans[ve[x][i].x]=S[d[x]+ve[x][i].y].size();
if (op) get(x,);
} int main(){
scanf("%d",&n);
rep(i,,n){
cin>>s>>fa[i];
if (mp.find(s)!=mp.end()) id[i]=mp[s]; else mp[s]=id[i]=++tot;
if (!fa[i]) rt[++rts]=i; else add(fa[i],i);
}
scanf("%d",&Q);
rep(i,,Q) scanf("%d%d",&x,&y),ve[x].push_back((P){i,y});
rep(i,,rts) dfs(rt[i]),work(rt[i],);
rep(i,,Q) printf("%d\n",ans[i]);
return ;
}
CF246E
例四:[CF741D]一棵有根树,边上有字母a~v,求每个子树中最长的边,满足这个边上的所有字母重排后可以构成回文。
很容易想到点分治,但由于是有根树,所以点分治是做不了的。
先对没给点记录w[x]表示x到根的路径的各个字母出现奇偶性,然后p[S]记录满足w[x]=S的所有点x的最大深度。
类似点分治,一棵子树一棵子树地处理,每次先对子树内每个点查找之前子树中是否有能与它拼成回文串的点与它的最大深度,再更新p[]。
注意查找和更新不能同时进行,要整棵子树都查询完毕之后再进行更新操作。注意查找与更新当前子树的根。
#include<cstdio>
#include<algorithm>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
#define For(i,x) for (int i=h[x],k; i; i=nxt[i])
using namespace std; const int N=,inf=1e9;
char s[N];
int n,cnt,d[N],fa[N],ans[N],w[N],sz[N],son[N],p[<<],h[N],to[N<<],nxt[N<<];
void add(int u,int v){ to[++cnt]=v; nxt[cnt]=h[u]; h[u]=cnt; } void dfs(int x){
d[x]=d[fa[x]]+; sz[x]=;
For(i,x){
k=to[i]; w[k]=w[x]^(<<(s[k]-'a')); dfs(k); sz[x]+=sz[k];
if (sz[k]>sz[son[x]]) son[x]=k;
}
} void get(int rt,int x){
ans[rt]=max(ans[rt],p[w[x]]+d[x]-*d[rt]);
rep(i,,) ans[rt]=max(ans[rt],p[w[x]^(<<i)]+d[x]-*d[rt]);
For(i,x) get(rt,k=to[i]);
} void get2(int x){ p[w[x]]=max(p[w[x]],d[x]); For(i,x) get2(k=to[i]); } void del(int x){ p[w[x]]=-inf; For(i,x) del(k=to[i]); } void work(int x,int op){
For(i,x) if ((k=to[i])!=son[x]) work(k,),ans[x]=max(ans[x],ans[k]);
if (son[x]) work(son[x],),ans[x]=max(ans[x],ans[son[x]]);
For(i,x) if ((k=to[i])!=son[x]) get(x,k),get2(k);
ans[x]=max(ans[x],p[w[x]]-d[x]);
rep(i,,) ans[x]=max(ans[x],p[w[x]^(<<i)]-d[x]);
p[w[x]]=max(p[w[x]],d[x]);
if (op) del(x);
} int main(){
freopen("741D.in","r",stdin);
freopen("741D.out","w",stdout);
scanf("%d",&n);
rep(i,,n) scanf("%d %c",&fa[i],&s[i]),add(fa[i],i);
rep(i,,(<<)-) p[i]=-inf;
dfs(); work(,);
rep(i,,n) printf("%d ",ans[i]);
return ;
}
CF741D
dsu on tree题表的更多相关文章
- [CSP-S模拟测试]:射手座之日(dsu on tree)
题目传送门(内部题103) 输入格式 第一行一个数$n$,表示结点的个数. 第二行$n–1$个数,第$i$个数是$p[i+1]$.$p[i]$表示结点$i$的父亲是$p[i]$.数据保证$p[i]&l ...
- CF 741D. Arpa’s letter-marked tree and Mehrdad’s Dokhtar-kosh paths [dsu on tree 类似点分治]
D. Arpa's letter-marked tree and Mehrdad's Dokhtar-kosh paths CF741D 题意: 一棵有根树,边上有字母a~v,求每个子树中最长的边,满 ...
- [dsu on tree]【学习笔记】
十几天前看到zyf2000发过关于这个的题目的Blog, 今天终于去学习了一下 Codeforces原文链接 dsu on tree 简介 我也不清楚dsu是什么的英文缩写... 就像是树上的启发式合 ...
- dsu on tree入门
先瞎扯几句 说起来我跟这个算法好像还有很深的渊源呢qwq.当时在学业水平考试的考场上,题目都做完了不会做,于是开始xjb出题.突然我想到这么一个题 看起来好像很可做的样子,然而直到考试完我都只想出来一 ...
- 洛谷P4482 [BJWC2018]Border 的四种求法 字符串,SAM,线段树合并,线段树,树链剖分,DSU on Tree
原文链接https://www.cnblogs.com/zhouzhendong/p/LuoguP4482.html 题意 给定一个字符串 S,有 q 次询问,每次给定两个数 L,R ,求 S[L.. ...
- BZOJ.4182.Shopping(点分治/dsu on tree 树形依赖背包 多重背包 单调队列)
BZOJ 题目的限制即:给定一棵树,只能任选一个连通块然后做背包,且每个点上的物品至少取一个.求花费为\(m\)时最大价值. 令\(f[i][j]\)表示在点\(i\),已用体积为\(j\)的最大价值 ...
- Codeforces 741D Arpa’s letter-marked tree and Mehrdad’s Dokhtar-kosh paths(dsu on tree)
感觉dsu on tree一定程度上还是与点分类似的.考虑求出跨过每个点的最长满足要求的路径,再对子树内取max即可. 重排后可以变成回文串相当于出现奇数次的字母不超过1个.考虑dsu on tree ...
- [学习笔记]Dsu On Tree
[dsu on tree][学习笔记] - Candy? - 博客园 题单: 也称:树上启发式合并 可以解决绝大部分不带修改的离线询问的子树查询问题 流程: 1.重链剖分找重儿子 2.sol:全局用桶 ...
- Codeforces.600E.Lomsat gelral(dsu on tree)
题目链接 dsu on tree详见这. \(Description\) 给定一棵树.求以每个点为根的子树中,出现次数最多的颜色的和. \(Solution\) dsu on tree模板题. 用\( ...
随机推荐
- C语言复习---打印菱形
#define _CRT_SECURE_NO_WARNINGS #include <stdio.h> #include <stdlib.h> #include <math ...
- Spark记录-Scala异常处理与文件I/O
Scala的异常处理类似许多其他语言(如Java).它不是以正常方式返回值,方法可以通过抛出异常来终止. 但是,Scala实际上并没有检查异常. 当您想要处理异常时,要像Java一样使用try {.. ...
- hdu 3022 Sum of Digits
http://acm.hdu.edu.cn/showproblem.php?pid=3022 题意: 最多不超过10000组数据,每组数据给定两个数n,m,求一个最小的数,使得该数每一位之和等于n,每 ...
- Spring RedisTemplate操作-String操作(2)
@Autowired @Resource(name="redisTemplate") private RedisTemplate<String, String> rt; ...
- Jad查看源码
需要者两个文件:下载地址:http://pan.baidu.com/s/11qq4I 1,解压jar包 有两个文件 分别是net.sf.jadclipse_3.3.0 jad.exe 文件 2,找到m ...
- Tju_Oj_3988Password
这个题是给树的前序和中序,输出后序. 做法是根据前序找根,根据根在中序中找中序的左右子树,根据左右子树长度找前序的左右子树,依此递归. 做过之后感觉还是比较基础的,废话不多说,上题上代码. Bob w ...
- Chrome插件笔记之content_scripts
一.概论 说这个之前先看一个段子,讲的是甲方有一奇葩客户,这客户看一网站某些样式很别扭不得劲,非要让乙方修改,乍一听没毛病,但关键是这网站不是乙方家的,根本没有修改权限,怎么办,客户就是上帝,上帝的要 ...
- 分模块开发创建dao子模块——(七)
1.选中父工程右键新建maven module
- js实现避免浏览器拦截弹出新页面的方法
1 问题描述 点击button按钮,提交页面的form表单,后台执行完毕后返回参数,前台页面需要该参数实现跳转,如何实现保留该原来的页面,并在浏览器选项卡新建一个页面,且不被浏览器拦截? 2 方法及问 ...
- Linux本地解析文件/etc/hosts说明【原创】
windows的域名本地解析文件hosts是可以一个域名对多个IP,如果有一个IP有问题,可以去解析到其他IP Linux的本地解析文件/etc/hosts,是否也可以这样呢.下面做了个测试 先看一下 ...