BZOJ 1103: [POI2007]大都市meg
1103: [POI2007]大都市meg
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 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
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
1
0
1
HINT
Source
分析
求一个有根树上从根节点到某节点的链上的权值和,支持边权的修改。
求出dfs序后,可以将问题转化为区间加减,单点查询的问题,线段树或树状数组均可。
代码
1 #include <cmath> 2 #include <cstdio> 3 #include <cstring> 4 #include <cstdlib> 5 #include <iostream> 6 #include <algorithm> 7 8 using namespace std; 9 10 const int N = 250005; 11 const int M = 500005; 12 13 int n; 14 int m; 15 16 int tot; 17 int hd[N]; 18 int nt[M]; 19 int to[M]; 20 21 int tim; 22 int arv[N]; 23 int lev[N]; 24 25 void prework_dfs(int u) 26 { 27 arv[u] = ++tim; 28 29 for (int i = hd[u]; ~i; i = nt[i]) 30 if (to[i] > u)prework_dfs(to[i]); 31 32 lev[u] = tim; 33 } 34 35 int tree[N]; 36 37 int lowbit(int x) 38 { 39 return x & -x; 40 } 41 42 int query(int p) 43 { 44 int res = 0; 45 46 while (p) 47 { 48 res += tree[p]; 49 p -= lowbit(p); 50 } 51 52 return res; 53 } 54 55 void insert(int p, int val) 56 { 57 while (p <= n) 58 { 59 tree[p] += val; 60 p += lowbit(p); 61 } 62 } 63 64 char s[100]; 65 66 signed main(void) 67 { 68 scanf("%d", &n); 69 70 memset(hd, -1, sizeof(hd)); 71 72 for (int i = 1; i < n; ++i) 73 { 74 int x, y; scanf("%d%d", &x, &y); 75 76 nt[tot] = hd[x]; to[tot] = y; hd[x] = tot++; 77 nt[tot] = hd[y]; to[tot] = x; hd[y] = tot++; 78 } 79 80 prework_dfs(1); 81 82 memset(tree, 0, sizeof(tree)); 83 84 for (int i = 2; i <= n; ++i) 85 insert(arv[i], 1), insert(lev[i] + 1, -1); 86 87 scanf("%d", &m); m = n + m - 1; 88 89 for (int i = 1; i <= m; ++i) 90 { 91 scanf("%s", s); 92 93 if (s[0] == 'A') 94 { 95 int x, y; 96 scanf("%d%d", &x, &y); 97 98 int p = x < y ? y : x; 99 100 insert(arv[p], -1); 101 insert(lev[p] + 1, 1); 102 } 103 else 104 { 105 int p; scanf("%d", &p); 106 printf("%d\n", query(arv[p])); 107 } 108 } 109 }
@Author: YouSiki