[WC2014] 紫荆花之恋 题解

啊啊啊啊啊啊啊啊啊啊啊我终于改完啦啊啊啊啊啊啊啊。

因为没有在最开始的时候将所有点设置为已经重构的，所以直接 \(R15-R70\) 间卡了两三天。

似乎也是我第一次大规模使用指针了。

---

这道题假如只有一次询问，就是一道简单淀粉质，直接在根节点建立平衡树，记录 \(r_x-dis(x,rt)\)，然后找到能和它配对的点对，再直接容斥即可。但是它不是。

考虑维护淀粉树（动态淀粉质），然后仍然使用上述操作，应该可以在离线情况下完成本题，但是它强制在线。

考虑维护动态淀粉树（动态动态淀粉质？？？），但似乎我们并没有学过动态淀粉树。我们借鉴平衡树的思路，在不平衡时使用一些方式，让它平衡。考虑替罪羊的 \(\alpha\) 判断是否平衡，再用普通的方式建立淀粉树即可。

平衡树我用的是替罪羊，简单好写常数小。

时间复杂度我猜可能是 \(O(n\log^2n)\) 吧，或者说是 \(O(n\sqrt n\log n)\)，本人时间复杂度分析小菜鸡，根本不会。

#include<bits/stdc++.h>
#define ls(x) sg[x].ls
#define rs(x) sg[x].rs
#define sz(x) sg[x].sz
#define vl(x) sg[x].vl
using namespace std;
const int N=1e5+5;
const int M=2e7+5,p=1e9;
const double alp=0.9;
int n,ndp[N],rb[N],cnt;
long long la;int* mdp;
namespace sad_goat_tree{
	struct sad_goat{
		int ls,rs,vl,sz;
	}sg[M];int tot,tp;
	int st[M],rts[N];
	int mk(){
		return tp?st[tp--]:++tot;
	}int fz(int k){
		int id=mk();
		sg[id]={0,0,k,1};
		return id;
	}void push_up(int x){
		sz(x)=sz(ls(x))+sz(rs(x))+1;
	}void build(int &x,int l,int r,int *a){
		if(l>r) return;
		int mid=(l+r)/2;x=fz(a[mid]);
		build(ls(x),l,mid-1,a);
		build(rs(x),mid+1,r,a);
		push_up(x);
	}void dfs(int x){
		if(ls(x)) dfs(ls(x));
		rb[++cnt]=vl(x);
		if(rs(x)) dfs(rs(x));
	}void clear(int &x){
		if(!x) return;
		clear(ls(x)),clear(rs(x));
		sg[x]={0,0,0,0},st[++tp]=x,x=0;
	}void rebuild(int &x){
		cnt=0,dfs(x),clear(x),build(x,1,cnt,rb);
	}int check(int x){
		return max(sz(ls(x)),sz(rs(x)))>alp*sz(x);
	}void add(int &x,int val){
		if(!x) return x=fz(val),void();
		sz(x)++,add((vl(x)<val)?rs(x):ls(x),val);
		if(check(x)) rebuild(x);
	}int que(int rtx,int va){
		int re=0;
		for(int x=rtx;x;){
			if(vl(x)>va) x=ls(x);
			else re+=sz(ls(x))+1,x=rs(x);
		}return re;
	}
}using namespace sad_goat_tree;
namespace starch_tree{
	struct edge{int v,w;};
	struct arc{int f,d,tr;};
	int t,siz[N],w[N],vis[N],num[N];
	vector<arc>ve[N];vector<edge>g[N];
	vector<int>so[N],nw[N];int top;
	int sum(int x,int fa){
		int re=1;
		for(auto y:g[x])
			if(y.v!=fa&&!vis[y.v])
				re+=sum(y.v,x);
		return re;
	}void get_rt(int x,int fa,int sm,int &rt){
		num[x]=0,siz[x]=1;
		for(auto y:g[x]){
			if(y.v==fa||vis[y.v]) continue;
			get_rt(y.v,x,sm,rt),siz[x]+=siz[y.v];
			num[x]=max(num[x],siz[y.v]);
		}num[x]=max(sm-siz[x],num[x]);
		if(!rt||num[x]<num[rt]) rt=x;
	}void dfs_tr(int x,int fa,int d,int a,int b){
		ve[x].push_back({a,d,b});
		so[a].push_back(x),mdp[++top]=d-w[x];
		for(auto y:g[x])
			if(y.v!=fa&&!vis[y.v])
				dfs_tr(y.v,x,d+y.w,a,b);
	}void solve(int x){
		int rt=0;get_rt(x,0,sum(x,0),rt);
		vis[rt]=1,so[rt].clear(),so[rt].push_back(rt);
		if(siz[x]==1) return clear(rts[rt]),add(rts[rt],-w[x]);
		mdp=ndp,top=-1,t=0;
		for(auto y:g[rt]){
			if(vis[y.v]) continue;
			mdp=mdp+top+1,t+=top+1,top=-1;int tr;
			dfs_tr(y.v,rt,y.w,rt,tp?st[tp]:tot+1);
			sort(mdp,mdp+top+1),build(tr,0,top,mdp);
			nw[rt].push_back(tr);
		}t+=top+1,ndp[t]=-w[rt],clear(rts[rt]);
		sort(ndp,ndp+t+1),build(rts[rt],0,t,ndp);
		for(auto y:g[rt]) if(!vis[y.v]) solve(y.v);
	}void rbuild(int x){
		for(auto y:so[x]){
			if(y!=x) while(ve[y].size()){
				if(ve[y].rbegin()->f==x){
					ve[y].pop_back();break;
				}else ve[y].pop_back();
			}for(auto it:nw[y]) clear(it);
			nw[y].clear(),vis[y]=0,clear(rts[y]);
		}solve(x);
	}int addn(int x,int fa,int va,int we){
		g[x].push_back({fa,va});
		g[fa].push_back({x,va});
		w[x]=we,ve[x]=ve[fa];
		for(auto it=ve[x].begin();it!=ve[x].end();it++)
			it->d+=va,add(it->tr,it->d-we),add(rts[it->f],it->d-we);
		int tr=0;add(tr,va-we);
		ve[x].push_back({fa,va,tr}),nw[fa].push_back(tr);
		add(rts[fa],va-we),ndp[0]=-we;int re=0;
		clear(rts[x]),add(rts[x],-we),so[x].push_back(x);
		for(auto y:ve[x]) so[y.f].push_back(x);
		re=que(rts[x],we)-1;
		for(auto y:ve[x])
			re+=que(rts[y.f],we-y.d)-que(y.tr,we-y.d);
		auto itl=ve[x].begin(),itr=itl;itr++;
		for(;itr!=ve[x].end();++itl,++itr)
			if(alp*so[itl->f].size()<=so[itr->f].size())
				{rbuild(itl->f);break;}
		return re;
	}
}using namespace starch_tree;
signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0),cout.tie(0);
	cin>>n>>n;int fa,va,we;
	for(int i=1;i<=n;i++) vis[i]=1;
	cin>>fa>>va>>we,w[1]=we,add(rts[1],-we);
	so[1].push_back(1),cout<<"0\n";
	for(int i=2;i<=n;i++){
		cin>>fa>>va>>we,fa^=la%p;
		la+=addn(i,fa,va,we),cout<<la<<"\n";
	}return 0;
}