【2022省选模拟】虚数之树——点分治、线段树

此题无链接

题目描述

原题面这里还有一个恶心的背景和图

题解

再一次抽自己：在没有先想log数据结构前，
不要想分块！
不要想分块！
不要想分块！

这题分块应该只需要对每个块建一个虚树，修改直接暴力更新虚树上的DP值，查询就整块+散点查即可。这个做法复杂度是常数较大的 $O(n\sqrt{n})$，不知道是不是对的，反正我考场上3.5h没调出来。

用点分治就很简单，直接对每个点分中心存一棵动态开点的线段树，加删点就直接在每层的点分中心的线段树上单点修改，值为到点分中心的距离，然后维护区间最小值。查询访问每层的线段树即可。

总复杂度是 $O(n{\log }^{2}n)$，空间消耗较大。

代码

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define uns unsigned
#define END putchar('\n')
#define IF (it->first)
#define IS (it->second)
#define lowbit(x) ((x)&-(x))
using namespace std;
const int MAXN=100003;
const ll INF=1e18;
inline ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[50],lpt;
inline void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}

struct itn{
	int ls,rs,a;itn(){}
	itn(int A){a=A,ls=rs=0;}
}t[MAXN*400];
int IN,UL=t[0].a=1e8;
inline void update(int x){t[x].a=min(t[t[x].ls].a,t[t[x].rs].a);}
inline int add(int x,int l,int r,int z,int a){
	if(!x)x=++IN,t[IN]=itn(1e8);
	if(l==r)return t[x].a=a,x;
	int mid=(l+r)>>1;
	if(z<=mid)t[x].ls=add(t[x].ls,l,mid,z,a);
	else t[x].rs=add(t[x].rs,mid+1,r,z,a);
	return update(x),x;
}
inline int query(int x,int l,int r,int a,int b){
	if(!x)return 1e8;
	if(a<=l&&b>=r)return t[x].a;
	int mid=(l+r)>>1,res=1e8;
	if(a<=mid)res=min(res,query(t[x].ls,l,mid,a,b));
	if(b>mid)res=min(res,query(t[x].rs,mid+1,r,a,b));
	return res;
}
#define pii pair<int,int>
#define fi first
#define se second
struct edge{
	int v,to;edge(){}
	edge(int V,int T){v=V,to=T;}
}e[MAXN<<1];
int EN,G[MAXN];
inline void addedge(int u,int v){
	e[++EN]=edge(v,G[u]),G[u]=EN;
	e[++EN]=edge(u,G[v]),G[v]=EN;
}
int n,m,ty,ct,mn,rt[MAXN];
bool dl[MAXN];
vector<pii>f[MAXN];
inline int pdfs(int x,int fa,int S){
	int siz=1,cp=0;
	for(int i=G[x];i;i=e[i].to){
		int v=e[i].v,sv;
		if(v==fa||dl[v])continue;
		sv=pdfs(v,x,S),siz+=sv,cp=max(cp,sv);
	}cp=max(cp,S-siz);
	if(cp<mn)mn=cp,ct=x;
	return siz;
}
inline void dfs1(int x,int fa,int c,int d){
	f[x].emplace_back(pii(c,d));
	for(int i=G[x];i;i=e[i].to){
		int v=e[i].v;
		if(v==fa||dl[v])continue;
		dfs1(v,x,c,d+1);
	}
}
inline void build(int x,int S){
	if(S<0)S=pdfs(x,0,S);
	ct=0,mn=1e8,pdfs(x,0,S),x=ct;
	dfs1(x,0,x,0),dl[x]=1;
	for(int i=G[x];i;i=e[i].to)
		if(!dl[e[i].v])build(e[i].v,-1);
}
int main()
{
	freopen("regression.in","r",stdin);
	freopen("regression.out","w",stdout);
	n=read(),m=read(),read(),ty=read();
	for(int i=1;i<n;i++)addedge(read(),read());
	build(1,n);
	int las=0;
	while(m--){
		int op=read();
		if(op==1){
			int x=(read()^(las*ty));
			for(pii y:f[x])rt[y.fi]=add(rt[y.fi],1,n,x,y.se);
		}else if(op==2){
			int x=(read()^(las*ty));
			for(pii y:f[x])rt[y.fi]=add(rt[y.fi],1,n,x,1e8);
		}else{
			int l=read(),r=read(),x=(read()^(las*ty));
			int ans=1e8;
			for(pii y:f[x])ans=min(ans,query(rt[y.fi],1,n,l,r)+y.se);
			if(ans>n)print(-1);
			else print(las=ans);
		}
	}
	return 0;
}