dsu on tree总结

dsu on tree

树上启发式合并。我并不知道为什么要叫做这个名字。。。

干什么的

可以在\(O(n\log n)\)的时间内完成对子树信息的询问，可横向对比把树按\(dfs\)序转成序列问题的\(O(n\sqrt n)\)莫队算法。

怎么实现

当\(dfs\)到一个点\(u\)，执行以下操作：

看上去像是暴力？

复杂度分析

只有当存在一条轻边时，这条轻边所连接的子树才需要被重新计算一次贡献。那么一个点被重复计算的次数就等于这个点到根路径上经过的轻边条数。而根据轻重链剖分那套理论，一个点到根路径上的轻边不会超过\(\log n\)条，所以这个算法的复杂度为\(O(n\log n*t)\)，其中\(t\)为一次计算贡献的复杂度。

题目

懒得分开写博客了qaq

CF600E Lomsat gelral

luogu

一棵树有\(n\)个结点，每个结点都有一种颜色，求每个子树中出现次数最多的颜色(们)的编号之和。

sol：维护\(num_i\)表示有多少种颜色出现了\(i\)次，\(sum_i\)表示这些颜色的编号和。显然计算贡献是\(O(1)\)的。

code

#include<cstdio>

#include<algorithm>

#include<cstring>

using namespace std;

int gi(){

	int x=0,w=1;char ch=getchar();

	while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();

	if (ch=='-') w=0,ch=getchar();

	while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();

	return w?x:-x;

}

#define ll long long

const int N = 1e5+5;

int n,a[N],to[N<<1],nxt[N<<1],head[N],cnt;

int sz[N],son[N],vis[N],tot[N],num[N],Max;ll sum[N],ans[N];

void link(int u,int v){

	to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;

}

void dfs1(int u,int f){

	sz[u]=1;

	for (int e=head[u];e;e=nxt[e])

		if (to[e]!=f){

			dfs1(to[e],u);sz[u]+=sz[to[e]];

			if (sz[to[e]]>sz[son[u]]) son[u]=to[e];

		}

}

void update(int u,int f,int val){

	int &t=tot[a[u]];

	--num[t];sum[t]-=a[u];

	t+=val;

	++num[t];sum[t]+=a[u];

	if (val>0) Max=max(Max,t);

	else if (!num[Max]) --Max;

	for (int e=head[u];e;e=nxt[e])

		if (to[e]!=f&&!vis[to[e]])

			update(to[e],u,val);

}

void dfs(int u,int f,int keep){

	for (int e=head[u];e;e=nxt[e])

		if (to[e]!=f&&to[e]!=son[u])

			dfs(to[e],u,0);

	if (son[u]) dfs(son[u],u,1),vis[son[u]]=1;

	update(u,f,1);

	ans[u]=sum[Max];

	vis[son[u]]=0;

	if (!keep) update(u,f,-1);

}

int main(){

	n=gi();

	for (int i=1;i<=n;++i) a[i]=gi();

	for (int i=1;i<n;++i){

		int u=gi(),v=gi();

		link(u,v);link(v,u);

	}

	dfs1(1,0);dfs(1,0,1);

	for (int i=1;i<=n;++i) printf("%I64d ",ans[i]);

	puts("");return 0;

}

CF570D Tree Requests

luogu

给你一棵\(n\)个节点的树，每个节点上有一个小写英文字母。每次询问\(v_i,h_i\)，表示在\(v_i\)子树中所有深度为\(h_i\)的点上的字母拿出来重新组合能否形成一个回文串。

sol：一些字母重新组合后能够形成回文串当且仅当存在至多一个字母的出现次数为奇数。所以可以把每个小写英文字母当作是一个二进制位然后维护子树中每个深度的异或和即可。

code

#include<cstdio>

#include<algorithm>

#include<cstring>

using namespace std;

int gi(){

	int x=0,w=1;char ch=getchar();

	while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();

	if (ch=='-') w=0,ch=getchar();

	while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();

	return w?x:-x;

}

const int N = 5e5+5;

struct node{int h,nxt;}a[N];

int n,m,nxt[N],head[N],val[N],ft[N];char s[N];

int dep[N],sz[N],son[N],vis[N],sum[N],ans[N];

void add(int v,int h,int id){

	a[id]=(node){h,ft[v]};ft[v]=id;

}

void dfs1(int u,int f){

	dep[u]=dep[f]+1;sz[u]=1;

	for (int v=head[u];v;v=nxt[v]){

		dfs1(v,u),sz[u]+=sz[v];

		if (sz[v]>sz[son[u]]) son[u]=v;

	}

}

void update(int u){

	sum[dep[u]]^=val[u];

	for (int v=head[u];v;v=nxt[v])

		if (!vis[v]) update(v);

}

bool cal(int x){

	int res=0;

	while (x) ++res,x^=x&-x;

	return res<=1;

}

void dfs(int u,int f,int keep){

	for (int v=head[u];v;v=nxt[v])

		if (v!=son[u]) dfs(v,u,0);

	if (son[u]) dfs(son[u],u,1),vis[son[u]]=1;

	update(u);

	for (int i=ft[u];i;i=a[i].nxt) ans[i]=cal(sum[a[i].h]);

	vis[son[u]]=0;

	if (!keep) update(u);

}

int main(){

	n=gi();m=gi();

	for (int i=2,f;i<=n;++i)

		f=gi(),nxt[i]=head[f],head[f]=i;

	scanf("%s",s+1);

	for (int i=1;i<=n;++i) val[i]=1<<s[i]-'a';

	for (int i=1,v,h;i<=m;++i) v=gi(),h=gi(),add(v,h,i);

	dfs1(1,0);dfs(1,0,1);

	for (int i=1;i<=m;++i) puts(ans[i]?"Yes":"No");

	return 0;

}

CF208E Blood Cousins

luogu

给你一片森林，每次询问和一个点有相同的\(k\)次祖先的点有多少个。

sol：倍增跳到那个祖先上然后就是只需要知道该深度有多少个点就行了。

code

#include<cstdio>

#include<algorithm>

#include<cstring>

using namespace std;

int gi(){

	int x=0,w=1;char ch=getchar();

	while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();

	if (ch=='-') w=0,ch=getchar();

	while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();

	return w?x:-x;

}

const int N = 1e5+5;

int n,m,nxt[N],head[N],root[N],fa[18][N],ft[N];

int dep[N],sz[N],son[N],vis[N],tot[N],ans[N];

struct node{int d,nxt;}a[N];

void add(int x,int d,int id){

	a[id]=(node){d,ft[x]};ft[x]=id;

}

void dfs1(int u,int f){

	fa[0][u]=f;dep[u]=dep[f]+1;sz[u]=1;

	for (int i=1;i<=17;++i) fa[i][u]=fa[i-1][fa[i-1][u]];

	for (int v=head[u];v;v=nxt[v]){

		dfs1(v,u);sz[u]+=sz[v];

		if (sz[v]>sz[son[u]]) son[u]=v;

	}

}

void update(int u,int val){

	tot[dep[u]]+=val;

	for (int v=head[u];v;v=nxt[v])

		if (!vis[v]) update(v,val);

}

void dfs(int u,int keep){

	for (int v=head[u];v;v=nxt[v])

		if (v!=son[u]) dfs(v,0);

	if (son[u]) dfs(son[u],1),vis[son[u]]=1;

	update(u,1);

	for (int i=ft[u];i;i=a[i].nxt) ans[i]=tot[a[i].d]-1;

	vis[son[u]]=0;

	if (!keep) update(u,-1);

}

int main(){

	n=gi();

	for (int i=1;i<=n;++i){

		int f=gi();

		if (!f) root[++root[0]]=i;

		else nxt[i]=head[f],head[f]=i;

	}

	for (int i=1;i<=root[0];++i) dfs1(root[i],0);

	m=gi();

	for (int i=1;i<=m;++i){

		int x=gi(),k=gi(),d=dep[x];

		for (int j=0;j<=17;++j) if (k&(1<<j)) x=fa[j][x];

		add(x,d,i);

	}

	for (int i=1;i<=root[0];++i) dfs(root[i],0);

	for (int i=1;i<=m;++i) printf("%d ",ans[i]);

	puts("");return 0;

}