[bzoj1103][树链剖分]大都市

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

题解

树链剖分裸题吧。。线段树维护区间内土路个数，然后跑树剖就好。。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const int MAXN=410000;
struct node
{
    int x,y,next;
}a[1110000];int len,last[MAXN];
void ins(int x,int y)
{
    len++;
    a[len].x=x;a[len].y=y;
    a[len].next=last[x];last[x]=len;
}
struct trnode
{
    int lc,rc,l,r;
    int c,maxnum;
}tr[1110000];int trlen;
void bt(int l,int r)
{
    int now=++trlen;
    tr[now].l=l;tr[now].r=r;tr[now].c=0;tr[now].maxnum=0;
    tr[now].lc=tr[now].rc=-1;
    if(l<r)
    {
        int mid=(l+r)/2;
        tr[now].lc=trlen+1;bt(l,mid);
        tr[now].rc=trlen+1;bt(mid+1,r);
    }
}
int fa[MAXN],dep[MAXN],tot[MAXN],son[MAXN];
int n;
void pre_tree_node(int x)
{
    tot[x]=1;son[x]=0;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(y!=fa[x])
        {
            fa[y]=x;
            dep[y]=dep[x]+1;
            pre_tree_node(y);

            if(tot[y]>tot[son[x]])son[x]=y;
            tot[x]+=tot[y];
        }
    }
}
int top[MAXN],ys[MAXN],z;
void pre_tree_edge(int x,int tp)
{
    ys[x]=++z;top[x]=tp;
    if(son[x]!=0)pre_tree_edge(son[x],tp);
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(y!=fa[x] && y!=son[x])pre_tree_edge(y,y);
    }
}
void change(int now,int pp,int k)
{
    if(tr[now].l==tr[now].r){tr[now].maxnum=k;return ;}
    int lc=tr[now].lc,rc=tr[now].rc;
    int mid=(tr[now].l+tr[now].r)/2;
    if(pp<=mid)change(lc,pp,k);
    else change(rc,pp,k);
    tr[now].maxnum=tr[lc].maxnum+tr[rc].maxnum;
}
int findsum(int now,int l,int r)
{
    if(tr[now].l==l && tr[now].r==r)return tr[now].maxnum;
    int lc=tr[now].lc,rc=tr[now].rc;
    int mid=(tr[now].l+tr[now].r)/2;
    if(r<=mid)return findsum(lc,l,r);
    else if(mid+1<=l)return findsum(rc,l,r);
    else return (findsum(lc,l,mid)+findsum(rc,mid+1,r));
}
int ans;int p[MAXN];
int solvesum(int x,int y)
{
    ans=0;
    int tx=top[x],ty=top[y];
    while(tx!=ty)
    {
        if(dep[tx]>dep[ty])
        {
            swap(tx,ty);
            swap(x,y);
        }
        ans+=findsum(1,ys[ty],ys[y]);
        y=fa[ty];ty=top[y];
    }
    if(x==y)return ans+p[x];
    else 
    {
        if(dep[x]>dep[y])swap(x,y);
        ans+=findsum(1,ys[x],ys[y]);
        return ans;
    }
}
struct bian
{
    int x,y,c;
}e[MAXN];
char ss[10];int u;int v;

int main()
{
//  freopen("Test.txt","w",stdout);
    scanf("%d",&n);
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<n;i++)
    {
        scanf("%d%d",&e[i].x,&e[i].y);
        ins(e[i].x,e[i].y);
        ins(e[i].y,e[i].x);
    }   
    int x;  
    scanf("%d",&x);
    dep[1]=0;fa[1]=0;
    pre_tree_node(1);

    z=0;pre_tree_edge(1,1);

    trlen=0;bt(1,z);

    for(int i=1;i<n;i++)
    {
        if(dep[e[i].x]>dep[e[i].y])swap(e[i].x,e[i].y);
        change(1,ys[e[i].y],1); 
    }

    for(int T=1;T<=n+x-1;T++)
    {
        scanf("%s",ss+1);
        scanf("%d",&u);
        if(ss[1]=='A')
        {
            scanf("%d",&v);
            change(1,ys[v],0);
        }
        else printf("%d\n",solvesum(1,u));
    }
    return 0;
}