【BZOJ-1103】大都市meg 树状数组 + DFS序

1103: [POI2007]大都市meg

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 2009  Solved: 1056
[Submit][Status][Discuss]

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

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

Source

Solution

一开始想出来一种树链剖分的做法，不过$nlog^{2}n$的复杂度显然承受不来

后来一想，其实很简单

因为一个点所能影响的只是他的子树，所以我们用DFS序处理一下

查分一下扔到树状数组上就好了，时间复杂度$O(nlogn)$

Code

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdio>
using namespace std;
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
#define MAXN 250010
int N,M;
struct EdgeNode{int next,to,from;}edge[MAXN<<1];
int head[MAXN],cnt=1;
void AddEdge(int u,int v) {cnt++; edge[cnt].next=head[u]; head[u]=cnt; edge[cnt].to=v; edge[cnt].from=u;}
void InsertEdge(int u,int v) {AddEdge(u,v); AddEdge(v,u);}
namespace BIT
{
int tree[MAXN];
inline int lowbit(int x) {return x&-x;}
inline void Change(int pos,int del)
{
    for (int i=pos; i<=N; i+=lowbit(i))
        tree[i]+=del;
}
inline int Query(int pos)
{
    int re=0;
    for (int i=pos; i; i-=lowbit(i))
        re+=tree[i];
    return re;
}
}
int pl[MAXN],pre[MAXN],pr[MAXN],dfn;
void DFS(int now,int last)
{
    pl[now]=++dfn; pre[dfn]=now;
    for (int i=head[now]; i; i=edge[i].next)
        if (edge[i].to!=last)
            DFS(edge[i].to,now);
    pr[now]=dfn;
}
int main()
{
    N=read();
    for (int u,v,i=1; i<=N-1; i++) u=read(),v=read(),InsertEdge(u,v);
    DFS(1,0);
    for (int i=2; i<=N; i++) BIT::Change(pl[i],1),BIT::Change(pr[i]+1,-1);
    M=read()+N-1;
    while (M--)
        {
            char opt[5]; scanf("%s",opt); int a,b;
            if (opt[0]=='A') a=read(),b=read(),BIT::Change(pl[b],-1),BIT::Change(pr[b]+1,1);
            if (opt[0]=='W') a=read(),printf("%d
",BIT::Query(pl[a]));
        }
    return 0;
}