在经济全球化浪潮的影响下,习惯于漫步在清晨的乡间小路的邮递员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
5 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
2 1 0 1
Hint
题解:这道题目,我看到就想到了树链剖分,但是根据标签是树状数组的题目( ⊙ o ⊙ )啊!
然后看了hzw的blog,发现,为什么求dfs序还不用递归的,这个服了,愣是看了半天。
然后还是自己理解打了一个,最后还是A了。
1 #include<cstdio> 2 #include<cstring> 3 #include<cmath> 4 #include<iostream> 5 #include<algorithm> 6 using namespace std; 7 const int N=250007; 8 9 int n,m,tim; 10 int cnt,head[N],next[N*2],rea[N*2]; 11 int l[N*2],r[N*2],tree[N*2]; 12 char ch[2]; 13 14 void add(int u,int v){next[++cnt]=head[u],head[u]=cnt,rea[cnt]=v;} 15 int lowbit(int x){return x&(-x);} 16 void dfs(int u,int fa) 17 { 18 l[u]=++tim; 19 for (int i=head[u];i!=-1;i=next[i]) 20 { 21 int v=rea[i]; 22 if (v==fa) continue; 23 dfs(v,u); 24 } 25 r[u]=++tim; 26 } 27 void ab(int x,int z) 28 { 29 for (int i=x;i<=2*n;i+=lowbit(i)) tree[i]+=z; 30 } 31 int query(int x) 32 { 33 int res=0; 34 for (int i=x;i>=1;i-=lowbit(i)) 35 res+=tree[i]; 36 return res; 37 } 38 int main() 39 { 40 memset(head,-1,sizeof(head)); 41 scanf("%d",&n); 42 int x,y; 43 for (int i=1;i<n;i++) 44 { 45 scanf("%d%d",&x,&y); 46 add(x,y),add(y,x); 47 } 48 dfs(1,-1); 49 for (int i=2;i<=n;i++) 50 ab(l[i],1),ab(r[i],-1); 51 scanf("%d",&m); 52 for (int i=1;i<=m+n-1;i++) 53 { 54 scanf("%s",ch); 55 if (ch[0]=='W') 56 { 57 scanf("%d",&x); 58 printf("%d ",query(l[x])); 59 } 60 else 61 { 62 scanf("%d%d",&x,&y); 63 x=max(x,y); 64 ab(l[x],-1),ab(r[x],1); 65 } 66 } 67 }
就是差分一下就可以了,=。