BZOJ1103: [POI2007]大都市meg
1103: [POI2007]大都市meg
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1145 Solved: 569
[Submit][Status]
Description
在经济全球化浪潮的影响下,习惯于漫步在清晨的乡间小路的邮递员Blue Mary也开始骑着摩托车传递邮件了。不过,她经常回忆起以前在乡间漫步的情景。昔日,乡下有依次编号为1..n的n个小村庄,某些村庄之间有一些双向的土路。从每个村庄都恰好有一条路径到达村庄1(即比特堡)。并且,对于每个村庄,它到比特堡的路径恰好只经过编号比它的编号小的村庄。另外,对于所有道路而言,它们都不在除村庄以外的其他地点相遇。在这个未开化的地方,从来没有过高架桥和地下铁道。随着时间的推移,越来越多的土路被改造成了公路。至今,Blue Mary还清晰地记得最后一条土路被改造为公路的情景。现在,这里已经没有土路了——所有的路都成为了公路,而昔日的村庄已经变成了一个大都市。 Blue Mary想起了在改造期间她送信的经历。她从比特堡出发,需要去某个村庄,并且在两次送信经历的间隔期间,有某些土路被改造成了公路.现在Blue Mary需要你的帮助:计算出每次送信她需要走过的土路数目。(对于公路,她可以骑摩托车;而对于土路,她就只好推车了。)
Input
第一行是一个数n(1 < = n < = 2 50000).
以下n-1行,每行两个整数a,b(1 < = a以下一行包含一个整数m(1 < = m < = 2 50000),表示Blue Mary曾经在改造期间送过m次信。
以下n+m-1行,每行有两种格式的若干信息,表示按时间先后发生过的n+m-1次事件:
若这行为 A a b(a若这行为 W a, 则表示Blue Mary曾经从比特堡送信到村庄a。
Output
有m行,每行包含一个整数,表示对应的某次送信时经过的土路数目。
Sample Input
1 2
1 3
1 4
4 5
4
W 5
A 1 4
W 5
A 4 5
W 5
W 2
A 1 2
A 1 3
Sample Output
1
0
1
HINT
Source
题解:
神奇的DFS序!
先DFS,记录前后时间戳l[i],r[i]则区间(l[i],r[i])就表示它的子树
我们add(l[i],1) add (r[i],-1) 则对 l[i]求前缀和就是从这个点到根节点的土路数目
因为如果一进一出 +1 -1正好抵消了!而如果在 i 的子树中就只有+1,就算上这条边了
可以理解为如果从i的父节点进入i,则土路+1,如果从i 回到i 的父节点,土路-1
下面考虑修改,直接add(l[i],-1) add(r[i],1)即可
这正好与之前的操作相反,这样无论在不在 i 的子树内都是 i 节点的贡献都是0,也就代表了这场为了一条公路
需要注意的几个细节:
1.根节点1不用add(l[i],1) add(r[i],-1),因为1没有父节点,而且从1出发,初始土路数肯定为0
2.修改(x,y)这条路要看谁是谁的儿子,我们可以记录dep来判断,也可以记录fa,总之只需要改掉儿子的信息即可
ps:这题用树链剖分可以骗到100分(14个测试点),其他略有超时。。。
代码:
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 {$inline on} 2 {$M 10000000,0,maxlongint} 3 const maxn=300000+100; 4 type node=record 5 go,next:longint; 6 end; 7 8 var e:array[0..2*maxn] of node; 9 head,dep,l,r:array[0..maxn] of longint; 10 s:array[0..2*maxn] of longint; 11 i,n,x,y,tot,clock:longint; 12 ch:char; 13 procedure swap(var x,y:longint);inline; 14 var t:longint; 15 begin 16 t:=x;x:=y;y:=t; 17 end; 18 procedure insert(x,y:longint);inline; 19 begin 20 inc(tot); 21 e[tot].go:=y;e[tot].next:=head[x];head[x]:=tot; 22 end; 23 procedure dfs(x:longint);inline; 24 var i,y:longint; 25 begin 26 inc(clock);l[x]:=clock; 27 i:=head[x]; 28 while i<>0 do 29 begin 30 y:=e[i].go; 31 if dep[y]=0 then 32 begin 33 dep[y]:=dep[x]+1; 34 dfs(y); 35 end; 36 i:=e[i].next; 37 end; 38 inc(clock);r[x]:=clock; 39 end; 40 procedure add(x,y:longint);inline; 41 begin 42 while x<=2*n do 43 begin 44 inc(s[x],y); 45 inc(x,x and (-x)); 46 end; 47 end; 48 function sum(x:longint):longint;inline; 49 begin 50 sum:=0; 51 while x>0 do 52 begin 53 inc(sum,s[x]); 54 dec(x,x and (-x)); 55 end; 56 end; 57 procedure init; 58 begin 59 readln(n); 60 for i:=1 to n-1 do begin readln(x,y);insert(x,y);insert(y,x);end; 61 dep[1]:=1; 62 dfs(1); 63 //for i:=1 to n do writeln(i,' ',l[i],' ',r[i],' ',dep[i]); 64 for i:=2 to n do begin add(l[i],1);add(r[i],-1);end; 65 end; 66 procedure solveask;inline; 67 begin 68 readln(x);writeln(sum(l[x])); 69 end; 70 procedure solvechange;inline; 71 begin 72 readln(x,y); 73 if dep[x]<dep[y] then swap(x,y); 74 add(l[x],-1);add(r[x],1); 75 end; 76 procedure main; 77 begin 78 readln; 79 while not eof do 80 begin 81 read(ch); 82 case ch of 83 'W':solveask; 84 'A':solvechange; 85 end; 86 end; 87 end; 88 begin 89 assign(input,'input.txt');assign(output,'output.txt'); 90 reset(input);rewrite(output); 91 init; 92 main; 93 close(input);close(output); 94 end.