BZOJ 1103: [POI2007]大都市meg

1103: [POI2007]大都市meg

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

 

[Submit][Status][Discuss]

分析

求一个有根树上从根节点到某节点的链上的权值和，支持边权的修改。

求出dfs序后，可以将问题转化为区间加减，单点查询的问题，线段树或树状数组均可。

代码

 #include <cmath>
 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <iostream>
 #include <algorithm>

 using namespace std;

 const int N = ;
 const int M = ;

 int n;
 int m;

 int tot;
 int hd[N];
 int nt[M];
 int to[M];

 int tim;
 int arv[N];
 int lev[N];

 void prework_dfs(int u)
 {
     arv[u] = ++tim;

     for (int i = hd[u]; ~i; i = nt[i])
         if (to[i] > u)prework_dfs(to[i]);

     lev[u] = tim;
 }

 int tree[N];

 int lowbit(int x)
 {
     return x & -x;
 }

 int query(int p)
 {
     int res = ;

     while (p)
     {
         res += tree[p];
         p -= lowbit(p);
     }

     return res;
 }

 void insert(int p, int val)
 {
     while (p <= n)
     {
         tree[p] += val;
         p += lowbit(p);
     }
 }

 char s[];

 signed main(void)
 {
     scanf("%d", &n);

     memset(hd, -, sizeof(hd));

     for (int i = ; i < n; ++i)
     {
         int x, y; scanf("%d%d", &x, &y);

         nt[tot] = hd[x]; to[tot] = y; hd[x] = tot++;
         nt[tot] = hd[y]; to[tot] = x; hd[y] = tot++;
     }

     prework_dfs();

     memset(tree, , sizeof(tree));

     for (int i = ; i <= n; ++i)
         insert(arv[i], ), insert(lev[i] + , -);

     scanf("%d", &m); m = n + m - ;

     for (int i = ; i <= m; ++i)
     {
         scanf("%s", s);

         if (s[] == 'A')
         {
             int x, y;
             scanf("%d%d", &x, &y);

             int p = x < y ? y : x;

             insert(arv[p], -);
             insert(lev[p] + , );
         }
         else
         {
             int p; scanf("%d", &p);
             printf("%d\n", query(arv[p]));
         }
     }
 }

BZOJ_1103.cpp

@Author: YouSiki