bzoj4154

一开始读错题，各种不会做，后来发现染色只是染孩子……

那不就简单了吗……注意这题是允许离线的

染色如果没有距离限制，它就是个dfs序

距离限制怎么做呢？我们考虑扩展一维变成二维的问题，将每个点变为二维平面上的点(x,y)，y=d[x]表示x的深度

染色a，距离限制l实际上就是对x∈[l,r],y<=d[a]+l的二维点染色([l,r]表示dfs序对应的区间）

这是一个非常经典的思路，我们可以扩展一维方便表示限制

注意染色的贡献是独立的，考虑离线的cdq分治，对于所有操作，我们把l变成d[a]+l;

然后以l为关键字降序排序，然后对时间线分治,设cdq(l,r)表示解决的第[l,r]的所有操作

考虑时间线[l,m]的染色操作对[m+1,r]的查询操作的影响，我们已经按l为关键字排序了，

对询问点产生影响的一定是区间包含这个点的最后一次操作，因此我们维护以时间为关键字的线段树

这样就是维护最大值，区间覆盖，单点查询，用线段树即可

然后递归下去处理即可

  1 const mo=1000000007;
  2 type node=record
  3        l,c,a:longint;
  4      end;
  5      way=record
  6        po,next:longint;
  7      end;
  8
  9 var  a:array[0..100010] of node;
 10      e:array[0..100010] of way;
 11      v:array[0..100010*4] of boolean;
 12      ans,fa,d,ld,rd,p,q,tq:array[0..100010] of longint;
 13      tag,st:array[0..100010*4] of longint;
 14      tot,k,i,n,m,tt,t,x,y,len:longint;
 15      s:int64;
 16
 17 procedure add(x,y:longint);
 18   begin
 19     inc(len);
 20     e[len].po:=y;
 21     e[len].next:=p[x];
 22     p[x]:=len;
 23   end;
 24
 25 procedure dfs(x:longint);
 26   var i,y:longint;
 27   begin
 28     inc(t);
 29     ld[x]:=t;
 30     i:=p[x];
 31     while i<>0 do
 32     begin
 33       y:=e[i].po;
 34       d[y]:=d[x]+1;
 35       dfs(y);
 36       i:=e[i].next;
 37     end;
 38     rd[x]:=t;
 39   end;
 40
 41 procedure swap(var a,b:longint);
 42   var c:longint;
 43   begin
 44     c:=a;
 45     a:=b;
 46     b:=c;
 47   end;
 48
 49 function cmp(i,j:longint):boolean;
 50   begin
 51     if a[i].l=a[j].l then exit(a[i].c>a[j].c);
 52     exit(a[i].l>a[j].l);
 53   end;
 54
 55 procedure sort(l,r:longint);
 56   var i,j,x:longint;
 57   begin
 58     i:=l;
 59     j:=r;
 60     x:=q[(l+r) shr 1];
 61     repeat
 62       while cmp(q[i],x) do inc(i);
 63       while cmp(x,q[j]) do dec(j);
 64       if not(i>j) then
 65       begin
 66         swap(q[i],q[j]);
 67         inc(i);
 68         dec(j);
 69       end;
 70     until i>j;
 71     if l<j then sort(l,j);
 72     if i<r then sort(i,r);
 73   end;
 74
 75 procedure deal(i,x:longint);
 76   begin
 77     if (tag[i]<x) then tag[i]:=x;
 78     if not v[i] then
 79     begin
 80       inc(tot);
 81       st[tot]:=i;
 82       v[i]:=true;
 83     end;
 84   end;
 85
 86 procedure push(i:longint);
 87   begin
 88     if tag[i]<>0 then
 89     begin
 90       deal(i*2,tag[i]);
 91       deal(i*2+1,tag[i]);
 92       tag[i]:=0;
 93     end;
 94   end;
 95
 96 procedure work(i,l,r,x,y,z:longint);
 97   var m:longint;
 98   begin
 99     if (x<=l) and (y>=r) then
100       deal(i,z)
101     else begin
102       m:=(l+r) shr 1;
103       push(i);
104       if x<=m then work(i*2,l,m,x,y,z);
105       if y>m then work(i*2+1,m+1,r,x,y,z);
106     end;
107   end;
108
109 function ask(i,l,r,x:longint):longint;
110   var m:longint;
111   begin
112     if (l=r) then exit(tag[i])
113     else begin
114       m:=(l+r) shr 1;
115       push(i);
116       if x<=m then exit(ask(i*2,l,m,x))
117       else exit(ask(i*2+1,m+1,r,x));
118     end;
119   end;
120
121 procedure cdq(l,r:longint);
122   var i,t,x,f1,f2,h:longint;
123   begin
124     if l=r then exit;
125     m:=(l+r) shr 1;
126     t:=0;
127     tot:=0;
128     h:=0;
129     for i:=l to r do
130     begin
131       if (q[i]<=m) and (a[q[i]].c<>0) or (q[i]>m) and (a[q[i]].c=0) then
132       begin
133         inc(t);
134         tq[t]:=q[i];
135       end;
136       if q[i]<=m then inc(h);
137     end;
138     for i:=1 to t do
139     begin
140       x:=tq[i];
141       if a[x].c=0 then
142       begin
143         f1:=ask(1,1,n,ld[a[x].a]);
144         if f1>0 then ans[x]:=a[f1].c;
145       end
146       else work(1,1,n,ld[a[x].a],rd[a[x].a],x);
147     end;
148     for i:=1 to tot do
149     begin
150       tag[st[i]]:=0;
151       v[st[i]]:=false;
152     end;
153     f1:=l;
154     f2:=l+h;
155     for i:=l to r do
156       if q[i]<=m then
157       begin
158         tq[f1]:=q[i];
159         inc(f1);
160       end
161       else begin
162         tq[f2]:=q[i];
163         inc(f2);
164       end;
165     for i:=l to r do
166       q[i]:=tq[i];
167     cdq(l,f1-1);
168     cdq(f1,r);
169   end;
170
171 begin
172   readln(tt);
173   while tt>0 do
174   begin
175     dec(tt);
176     fillchar(p,sizeof(p),0);
177     readln(n,m,k);
178     len:=0;
179     for i:=2 to n do
180     begin
181       read(fa[i]);
182       add(fa[i],i);
183     end;
184     t:=0;
185     dfs(1);
186     for i:=1 to k do
187     begin
188       readln(a[i].a,a[i].l,a[i].c);
189       if a[i].c<>0 then
190       begin
191         a[i].l:=a[i].l+d[a[i].a];
192         ans[i]:=0;
193       end
194       else begin
195         a[i].l:=d[a[i].a];
196         ans[i]:=1;
197       end;
198       q[i]:=i;
199     end;
200       sort(1,k);
201     cdq(1,k);
202     s:=0;
203     for i:=1 to k do
204       s:=(s+int64(ans[i])*int64(i)) mod mo;
205     writeln(s);
206   end;
207 end.