传送门
我们令一个人开始收集的时间为 ti, 当前时间为t
那么答案就是路径上所有 ti < t - c 的人的个数, 如果没有修改的话, 可以用主席树
但是要一边修改一边询问, 我最开始想到的是树链剖分+线段树套权值线段树
后来发现, 如果我们一开始讲所有修改读入然后建主席树呢, 因为后面的修改一定不能作为当前的答案(ti 明显大于t)
于是就是裸的模板了
#include<bits/stdc++.h>
#define N 200050
using namespace std;
struct Query{int op,x,y,c;}q[N];
struct Node{int ls,rs,val;}t[N*40];
int n,m,val[N],rt[N],cnt;
int first[N],next[N],to[N],tot;
void add(int x,int y){next[++tot]=first[x],first[x]=tot,to[tot]=y;}
int fa[N][20],dep[N];
int read(){
int cnt=0,f=1;char ch=0;
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch))cnt=cnt*10+(ch-'0'),ch=getchar();
return cnt*f;
}
void dfs(int u,int f){
for(int i=1;i<=18;i++)
fa[u][i] = fa[fa[u][i-1]][i-1];
for(int i=first[u];i;i=next[i]){
int t=to[i]; dep[t] = dep[u] + 1;
fa[t][0] = u; dfs(t,u);
}
}
void Build(int &x,int l,int r){
x = ++cnt;
if(l==r) return;
int mid = (l+r) >> 1;
Build(t[x].ls, l, mid);
Build(t[x].rs, mid+1, r);
}
void Insert(int &x,int last,int l,int r,int pos){
x = ++cnt;
t[x] = t[last]; t[x].val++;
if(l==r) return;
int mid = (l+r) >> 1;
if(pos<=mid) Insert(t[x].ls, t[last].ls, l, mid, pos);
else Insert(t[x].rs, t[last].rs, mid+1, r, pos);
}
void dfs2(int u,int f){
val[u] = val[u] == 0 ? m + 1 : val[u];
Insert(rt[u], rt[f], 1, m+1, val[u]);
for(int i=first[u];i;i=next[i]){
int t=to[i]; dfs2(t,u);
}
}
int lca(int x,int y){
if(dep[x]<dep[y]) swap(x,y);
for(int i=18;i>=0;i--)
if(dep[fa[x][i]] >= dep[y]) x = fa[x][i];
if(x==y) return x;
for(int i=18;i>=0;i--)
if(fa[x][i] != fa[y][i]) x = fa[x][i], y = fa[y][i];
return fa[x][0];
}
int Quary(int a,int b,int c,int d,int l,int r,int L,int R){
if(L<=l && r<=R) return t[a].val + t[b].val - t[c].val - t[d].val;
int mid = (l+r) >> 1, ans = 0;
if(L<=mid) ans += Quary(t[a].ls, t[b].ls, t[c].ls, t[d].ls, l, mid, L, R);
if(R>mid) ans += Quary(t[a].rs, t[b].rs, t[c].rs, t[d].rs, mid+1, r, L, R);
return ans;
}
int main(){
n = read();
for(int i=1;i<=n;i++){
int x = read(); if(i!=1) add(x,i);
} dep[1] = 1, dfs(1,0);
m = read();
Build(rt[0], 1, m+1);
for(int i=1;i<=m;i++){
int op = read();
if(op==1) q[i] = (Query){1,read(),read(),read()};
if(op==2) val[read()] = i;
} dfs2(1,0);
for(int i=1;i<=m;i++){
int op = q[i].op; if(!op) continue;
int x = q[i].x, y = q[i].y, c = q[i].c;
int val = i - c - 1; int l = lca(x,y);
printf("%d ",dep[x] + dep[y] - dep[l] * 2 + 1);
if(val <= 0){printf("0\n"); continue;}
int ans = Quary(rt[x], rt[y], rt[l], rt[fa[l][0]], 1, m+1, 1, val);
printf("%d\n",ans);
} return 0;
}
