哈哈哈哈哈哈哈哈哈哈哈哈我终于过了！！！！！！！！！！！！！！！
从昨天上午就开始写了，下午回家之后调了一会儿没什么感觉，就删了重打了一遍，然后调了一晚上+今天半个上午......我*******终于过了......
一道极其恶心的动态树分治......
首先点分治，限制条件就变成了di+dj<=ri+rj，移项得rj-dj>=di-ri，对重心和子树开平衡树维护di-ri即可。
查询的时候先跳点分治树更新答案，然后逐层把di-ri插入平衡树，上跳的过程中判断是否有子树失衡(某个子树大于子树总大小的alpha倍)，有则提出来暴力重构成完全平衡的点分治树。
各种坑爹细节和脑残错误调了一天，不过还是过了，激动......
1 #include<cstdio>
2 #include<cstring>
3 #include<algorithm>
4 #include<vector>
5 using namespace std;
6 const int maxn=100010;
7 const double alpha=0.7;
8 struct node{
9 static int randint(){
10 static int a=1213,b=97818217,p=998244353,x=751815431;
11 x=a*x+b;x%=p;
12 return x<0?(x+=p):x;
13 }
14 int data,size,p;
15 node *ch[2];
16 node(int d):data(d),size(1),p(randint()){}
17 inline void refresh(){size=ch[0]->size+ch[1]->size+1;}
18 }*null=new node(0),*root[maxn],*root1[maxn][50];
19 void addnode(int,int);
20 void rebuild(int,int,int,int);
21 void dfs_getcenter(int,int,int&);
22 void dfs_getdis(int,int,int,int);
23 void dfs_destroy(int,int);
24 void insert(int,node*&);
25 int order(int,node*);
26 void destroy(node*&);
27 void rot(node*&,int);
28 vector<int>G[maxn],W[maxn];
29 int size[maxn]={0},siz[maxn][50]={0},son[maxn];
30 bool vis[maxn];
31 int depth[maxn],p[maxn],d[maxn][50],id[maxn][50];
32 int n,m,w[maxn],tmp;
33 long long ans=0;
34 int main(){
35 null->size=0;
36 null->ch[0]=null->ch[1]=null;
37 scanf("%*d%d",&n);
38 fill(vis,vis+n+1,true);
39 fill(root,root+n+1,null);
40 for(int i=0;i<=n;i++)fill(root1[i],root1[i]+50,null);
41 scanf("%*d%*d%d",&w[1]);
42 insert(-w[1],root[1]);
43 size[1]=1;
44 printf("0\n");
45 for(int i=2;i<=n;i++){
46 scanf("%d%d%d",&p[i],&tmp,&w[i]);
47 p[i]^=(ans%(int)1e9);
48 G[i].push_back(p[i]);
49 W[i].push_back(tmp);
50 G[p[i]].push_back(i);
51 W[p[i]].push_back(tmp);
52 addnode(i,tmp);
53 printf("%lld\n",ans);
54 }
55 return 0;
56 }
57 void addnode(int x,int z){//wj-dj>=di-wi
58 depth[x]=depth[p[x]]+1;
59 size[x]=1;
60 insert(-w[x],root[x]);
61 int rt=0;
62 for(int u=p[x],k=depth[p[x]];u;u=p[u],k--){
63 if(u==p[x]){
64 id[x][k]=x;
65 d[x][k]=z;
66 }
67 else{
68 id[x][k]=id[p[x]][k];
69 d[x][k]=d[p[x]][k]+z;
70 }
71 ans+=order(w[x]-d[x][k],root[u])-order(w[x]-d[x][k],root1[id[x][k]][k]);
72 insert(d[x][k]-w[x],root[u]);
73 insert(d[x][k]-w[x],root1[id[x][k]][k]);
74 size[u]++;
75 siz[id[x][k]][k]++;
76 if(siz[id[x][k]][k]>size[u]*alpha+5)rt=u;
77 }
78 id[x][depth[x]]=0;
79 d[x][depth[x]]=0;
80 if(rt){
81 dfs_destroy(rt,depth[rt]);
82 rebuild(rt,depth[rt],size[rt],p[rt]);
83 }
84 }
85 void rebuild(int x,int k,int s,int pr){
86 int u=0;
87 dfs_getcenter(x,s,u);
88 vis[x=u]=true;
89 p[x]=pr;
90 depth[x]=k;
91 size[x]=s;
92 d[x][k]=id[x][k]=0;
93 destroy(root[x]);
94 insert(-w[x],root[x]);
95 if(s<=1)return;
96 for(int i=0;i<(int)G[x].size();i++)if(!vis[G[x][i]]){
97 p[G[x][i]]=0;
98 d[G[x][i]][k]=W[x][i];
99 siz[G[x][i]][k]=p[G[x][i]]=0;
100 destroy(root1[G[x][i]][k]);
101 dfs_getdis(G[x][i],x,G[x][i],k);
102 }
103 for(int i=0;i<(int)G[x].size();i++)if(!vis[G[x][i]])rebuild(G[x][i],k+1,size[G[x][i]],x);
104 }
105 void dfs_getcenter(int x,int s,int &u){
106 size[x]=1;
107 son[x]=0;
108 for(int i=0;i<(int)G[x].size();i++)if(!vis[G[x][i]]&&G[x][i]!=p[x]){
109 p[G[x][i]]=x;
110 dfs_getcenter(G[x][i],s,u);
111 size[x]+=size[G[x][i]];
112 if(size[G[x][i]]>size[son[x]])son[x]=G[x][i];
113 }
114 if(!u||max(s-size[x],size[son[x]])<max(s-size[u],size[son[u]]))u=x;
115 }
116 void dfs_getdis(int x,int u,int rt,int k){
117 insert(d[x][k]-w[x],root[u]);
118 insert(d[x][k]-w[x],root1[rt][k]);
119 id[x][k]=rt;
120 siz[rt][k]++;
121 size[x]=1;
122 for(int i=0;i<(int)G[x].size();i++)if(!vis[G[x][i]]&&G[x][i]!=p[x]){
123 p[G[x][i]]=x;
124 d[G[x][i]][k]=d[x][k]+W[x][i];
125 dfs_getdis(G[x][i],u,rt,k);
126 size[x]+=size[G[x][i]];
127 }
128 }
129 void dfs_destroy(int x,int k){
130 vis[x]=false;
131 for(int i=0;i<(int)G[x].size();i++)if(depth[G[x][i]]>=k&&G[x][i]!=p[x]){
132 p[G[x][i]]=x;
133 dfs_destroy(G[x][i],k);
134 }
135 }
136 void insert(int x,node *&rt){
137 if(rt==null){
138 rt=new node(x);
139 rt->ch[0]=rt->ch[1]=null;
140 return;
141 }
142 int d=x>=rt->data;
143 insert(x,rt->ch[d]);
144 rt->refresh();
145 if(rt->ch[d]->p<rt->p)rot(rt,d^1);
146 }
147 int order(int x,node *rt){
148 int ans=0,d;
149 x++;
150 while(rt!=null){
151 if((d=x>rt->data))ans+=rt->ch[0]->size+1;
152 rt=rt->ch[d];
153 }
154 return ans;
155 }
156 void destroy(node *&x){
157 if(x==null)return;
158 destroy(x->ch[0]);
159 destroy(x->ch[1]);
160 delete x;
161 x=null;
162 }
163 void rot(node *&x,int d){
164 node *y=x->ch[d^1];
165 x->ch[d^1]=y->ch[d];
166 y->ch[d]=x;
167 x->refresh();
168 (x=y)->refresh();
169 }
170 /*
171 Wc2014 紫荆花之恋
172 di+dj<=wi+wj，移项得wj-dj>=di-wi，动态点分治维护即可。
173 每个重心和子树存一棵平衡树维护di-wi，加点的时候跳点分治树更新答案后加入新点，
174 加入之后判断子树是否失衡，找出最高的失衡节点，把整个子树暴力重构成完全平衡的点分治树。
175 */
时间： 2024-10-07 16:49:06