【bzoj 1103】[POI2007]大都市meg

1103: [POI2007]大都市meg

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 2324   Solved: 1231
[ 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

本来以为树链剖分就完了，做完翻题解dfs序？？？抽空补一下

#include<bits/stdc++.h>
using namespace std;
const int maxn = 250050;
int h[maxn],nx[maxn<<1],to[maxn<<1],len;
int fa[maxn],top[maxn],sz[maxn],ch[maxn],rk[maxn],dep[maxn],cnt;
int val[maxn],seg[maxn<<2];
int n,m;
void insert(int a,int b)
{
    len++;nx[len]=h[a];to[len]=b;h[a]=len;
    len++;nx[len]=h[b];to[len]=a;h[b]=len;
}
void pushup(int rt)
{
    seg[rt]=seg[rt<<1]+seg[rt<<1|1];
}
void build(int rt,int l,int r)
{
    if(l==r){seg[rt]=1;if(l==1){seg[rt]=0;}return;}
    int mid=(l+r)>>1;
    build(rt<<1,l,mid);
    build(rt<<1|1,mid+1,r);
    pushup(rt);
}
int query(int rt,int L,int R,int l,int r)
{
    if(l<=L && r>=R)return seg[rt];
    int mid=(L+R)>>1;
    if(r<=mid)return query(rt<<1,L,mid,l,r);
    else if(l>mid)return query(rt<<1|1,mid+1,R,l,r);
    return query(rt<<1,L,mid,l,mid) + query(rt<<1|1,mid+1,R,mid+1,r);
}
void modify(int rt,int l,int r,int loc,int val)
{
    int mid=(l+r)>>1;
    if(l==r){seg[rt]=0; return;}
    if(loc<=mid)modify(rt<<1,l,mid,loc,val);
    else modify(rt<<1|1,mid+1,r,loc,val);
    pushup(rt);
}
int lca(int a,int b)
{
    while(top[a] != top[b])
    {
    //  cout<<a<<" "<<b<<endl;
        dep[top[a]] > dep[top[b]] ? a=fa[top[a]] : b=fa[top[b]];  
    }
    return dep[a] < dep[b] ? a : b;
}
void dfs1(int nd,int pr)
{
    fa[nd]=pr;
    dep[nd]=dep[pr]+1;
    sz[nd]=1;
    for(int i=h[nd];i;i=nx[i])
    {
        if(to[i]!=fa[nd])
        {
            dfs1(to[i],nd);
            sz[nd]+=sz[to[i]];
            if(sz[to[i]] > sz[ch[nd]])
            {
                ch[nd]=to[i];
            }
        }
    }
}
void dfs2(int nd,int tp)
{
    top[nd]=tp;
    rk[nd]=++cnt;
    if(ch[nd]!=0)dfs2(ch[nd],tp);
    for(int i=h[nd];i;i=nx[i])
    {
        if(to[i]!=fa[nd] && to[i]!=ch[nd])
        {
            dfs2(to[i],to[i]);
        }
    }
}
int getans(int nd)
{
    int ret=0;
    while(top[nd]!=1)
    {
        ret+=query(1,1,cnt,rk[top[nd]],rk[nd]);
         
        nd=fa[top[nd]];
    }//cout<<rk[nd]<<" "<<rk[top[nd]]<<" "<<ret<<endl;
    ret+=query(1,1,cnt,rk[top[nd]],rk[nd]);
    return ret;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n-1;i++)
    {
        int _a,_b;
        scanf("%d%d",&_a,&_b);
        insert(_a,_b);
    }
     
    dfs1(1,0);
    dfs2(1,1);
    build(1,1,cnt);
    scanf("%d",&m);
    for(int i=1;i<=n+m-1;i++)
    {
        char _c[60];int _a,_b;
        scanf("%s",_c);
        if(_c[0]=='W')
        {
            scanf("%d",&_a);
            printf("%d\n",getans(_a));
        }
        if(_c[0]=='A')
        {
            scanf("%d%d",&_a,&_b);
            modify(1,1,cnt,rk[_b],0);
        }
    }
    return 0;
}