题意
给一棵树和m个操作(持续m秒),操作如下
1.给一个点打上标记,被标记的点每秒加1危险度,(被标记时仍为0,后一秒为1)。
2.查询x,y的路径上有多少危险度大于k的点。
解析
这种每秒增加的标记显然是不好维护的。
考虑转化。
设某点在t时刻被标记,i时刻提出询问,那么此时它的危险度为i-t。
则有i-t>k,t
#include <cstdio>
#include <algorithm>
#define Rep( i , _begin , _end ) for(int i=(_begin),i##_END=(_end);i<=(i##_END);i++)
#define For( i , _begin , _end ) for(int i=(_begin),i##_END=(_end);i!=(i##_END);i++)
#define Lop( i , _begin , _end ) for(int i=(_begin),i##_END=(_end);i>=(i##_END);i--)
#define Dnt( i , _begin , _end ) for(int i=(_begin),i##_END=(_end);i!=(i##_END);i--)
using std :: max;
using std :: min;
const int maxx = 200000 + 25;
typedef int Array[maxx];
int n,m,x,y,z,num,cnt,tot,f,rt,t;
int head[maxx],to[maxx<<1],nxt[maxx<<1];
int T[maxx];
Array son,ftr,dpt,top,rnk,size,ans,ans1;
class Que{
public:
int l,r,id;
int v,pos,f;
}A[maxx];
bool cmp(Que a,Que b){
return a.v == b.v? a.f < b.f : a.v < b.v;
}
bool cmp1(Que a,Que b){
return a.id < b.id;
}
namespace Bit{
void Add(int x,int k){
for(int i=x;i<=n;i+=i&(-i))
T[i] += k;
}
int Query(int l,int r){
int ans = 0;
for(int i=r;i;i-=i&(-i))
ans += T[i];
for(int i=l-1;i;i-=i&(-i))
ans -= T[i];
return ans;
}
}
namespace Cute{
using namespace Bit;
void Ins(int x,int y){
if((!x) || (!y)) return;
to[++num] = y;nxt[num] = head[x];head[x] = num;
}
void Dfs1(int x){
size[x] = 1;
for(int i=head[x];i;i=nxt[i]){
if(to[i] == ftr[x]) continue;
dpt[to[i]] = dpt[x] + 1;ftr[to[i]] = x;
Dfs1(to[i]);size[x] += size[to[i]];
if(size[to[i]] > size[son[x]]) son[x] = to[i];
}
}
void Dfs2(int x,int brn){
rnk[x] = ++cnt;top[x] = brn;
if(son[x]) Dfs2(son[x],brn);
for(int i=head[x];i;i=nxt[i])
if(to[i] != ftr[x] && to[i] != son[x])
Dfs2(to[i],to[i]);
}
int Get(int x,int y,int &t){
int ans = 0;t = dpt[x] + dpt[y];
while(top[x] != top[y]){
if(dpt[top[x]] > dpt[top[y]]) std :: swap(x,y);
ans += Query(rnk[top[y]],rnk[y]);
y = ftr[top[y]];
}
if(rnk[x] > rnk[y]) std :: swap(x,y);
ans += Query(rnk[x],rnk[y]);
t -= 2*dpt[x];t ++;
return ans;
}
}
using namespace Cute;
int main(){
scanf("%d",&n);
Rep( i , 1 , n ){
scanf("%d",&x);
if(x == 0) rt = i;
Ins(i,x),Ins(x,i);
}
Dfs1(rt),Dfs2(rt,rt);
scanf("%d",&m);
Rep( i , 1 , m ){
scanf("%d",&f);
if(f == 2) scanf("%d",&x),A[i].pos = x,A[i].v = i;
if(f == 1){
scanf("%d%d%d",&x,&y,&z);
A[i].l = x;A[i].r = y;
A[i].v = i-z-1;
A[i].f = 1;A[i].id = i;
}
}
std :: sort(A+1,A+m+1,cmp);
Rep( i , 1 , m ){
if(A[i].f == 0) Add(rnk[A[i].pos],1);
if(A[i].f == 1){
Get(A[i].l,A[i].r,t);
ans[A[i].id] = t;
ans1[A[i].id] = Get(A[i].l,A[i].r,t);
//printf("%d %d\n",t,Get(A[i].l,A[i].r,t));
}
}
std :: sort(A+1,A+m+1,cmp1);
Rep( i , 1 , m ) if(A[i].f) printf("%d %d\n",ans[A[i].id],ans1[A[i].id]);
return 0;
}