POJ 3321 Apple Tree 树状数组+建树

古人学问无遗力,少壮工夫老始成。这篇文章主要讲述POJ 3321 Apple Tree 树状数组+建树相关的知识,希望能为你提供帮助。
题目链接:http://poj.org/problem?id=3321
Apple Tree

Time Limit:  2000MS   Memory Limit:  65536K
Total Submissions:  34812   Accepted:  10469
Description
There is an apple tree outside of kaka\'s house. Every autumn, a lot of apples will grow in the tree. Kaka likes apple very much, so he has been carefully nurturing the big apple tree.
The tree has  N  forks which are connected by branches. Kaka numbers the forks by 1 to  N  and the root is always numbered by 1. Apples will grow on the forks and two apple won\'t grow on the same fork. kaka wants to know how many apples are there in a sub-tree, for his study of the produce ability of the apple tree.
The trouble is that a new apple may grow on an empty fork some time and kaka may pick an apple from the tree for his dessert. Can you help kaka?
 
POJ 3321 Apple Tree 树状数组+建树

文章图片
 
Input
The first line contains an integer  N  (N  ≤ 100,000) , which is the number of the forks in the tree.
The following  N  - 1 lines each contain two integers  u  and  v, which means fork  u  and fork  v  are connected by a branch.
The next line contains an integer  M  (M  ≤ 100,000).
The following  M  lines each contain a message which is either
"C  x" which means the existence of the apple on fork  x  has been changed. i.e. if there is an apple on the fork, then Kaka pick it; otherwise a new apple has grown on the empty fork.
or
"Q  x" which means an inquiry for the number of apples in the sub-tree above the fork  x, including the apple (if exists) on the fork x
Note the tree is full of apples at the beginning
Output
For every inquiry, output the correspond answer per line.Sample Input
3 1 2 1 3 3 Q 1 C 2 Q 1

Sample Output
3 2

Source
POJ Monthly--2007.08.05, Huang, Jinsong 
 
题意概括:
给一个 N 阶无向图,要求把这个无向图转换为一颗以 1 为根的树后进行两种操作:
C x :x结点上如果有苹果就摘掉,如果没有苹果就长出一个新的来
【POJ 3321 Apple Tree 树状数组+建树】Q x:    查询 x 结点与它的所有后代分支一共有几个苹果
解题思路:
因为原本苹果之间的关系我们不好处理,所以就给他们这些结点重新编号。
并且为了方便子树的苹果数求和(即区间查询)我们在给结点增加一个左值一个右值表示当前结点管辖的范围。
如何建树?
根据题目一开始给的边与边的关系我们可以通过 dfs 进行建树。
L[ x ] (X的左值): X结点的新编号就是该结点的左值,表示它管辖范围的下限
R[ x ] (X的右值):而子树的编号最大值为该结点的右值,表示它管辖范围的上限。
至于dfs顺序,如下
举个栗子:
N = 6;
1 - 2
1 - 6
2 - 5
2 - 3
6 - 4
 
POJ 3321 Apple Tree 树状数组+建树

文章图片

POJ 3321 Apple Tree 树状数组+建树

文章图片

POJ 3321 Apple Tree 树状数组+建树

文章图片

到这里感觉有一点点线段树的味道了。。。
但这里要用树状数组对这些区间进行维护和查询。
例如:我要 Q 2;
由上面建好的新树我们可以知道 结点2这颗子树包含了结点2、结点5、结点3,而结点2所管辖的区间【2,4】刚好把它们包括进去了。
如果我们要求 这可子树的苹果数量,其实就是用树状数组求区间【2,4】的苹果总数
也就是右值的前缀和减去左值的前缀和(但不包括左值,因为结点本身也要算进去)即:ans = sum( R[ 2 ] ) - sum( L[ 2 ] - 1);
而C操作其实就是树状数组的单点更新而已。
注意事项:
一开始用 stl 的 vector 存无向图,虽然过了但效率不咋地。想想好久没用静态邻接表了,换了一下,快了好多。
AC code:
1 ///建树+树状数组(stl)8080k 1047ms 2 /* 3 #include < cstdio> 4 #include < iostream> 5 #include < cstring> 6 #include < vector> 7 #include < algorithm> 8 #define INF 0x3f3f3f3f 9 using namespace std; 10 typedef vector< int> ve; 11 12 const int MAXN= 1e5+10; 13 int t[MAXN], l[MAXN], r[MAXN]; 14 bool vis[MAXN]; 15 int N, M, cnt; 16 17 vector< ve> node(MAXN); 18 19 int lowbit(int x) 20 { 21return x& (-x); 22 } 23 void add(int x, int value) 24 { 25for(int i = x; i < = N; i+=lowbit(i)) 26t[i]+=value; 27 } 28 int sum(int x) 29 { 30int res = 0; 31for(int i = x; i > 0; i-=lowbit(i)) 32res+=t[i]; 33return res; 34 } 35 void dfs(int k) 36 { 37l[k] = cnt; 38for(int i = 0; i < node[k].size(); i++) 39{ 40cnt++; 41dfs(node[k][i]); 42} 43r[k] = cnt; 44return; 45 } 46 int main() 47 { 48scanf("%d", & N); 49int a, b; 50for(int i = 1; i < = N; i++) 51{ 52vis[i] = 1; 53add(i, 1); 54} 55for(int i = 1; i < N; i++) 56{ 57scanf("%d%d", & a, & b); 58node[a].push_back(b); 59} 60cnt = 1; 61dfs(1); 62scanf("%d", & M); 63while(M--) 64{ 65char com[3]; 66int px; 67scanf("%s", & com); 68if(com[0] == \'C\') 69{ 70scanf("%d", & px); 71if(vis[px]) add(l[px], -1); 72else add(l[px], 1); 73vis[px] = !vis[px]; 74} 75else 76{ 77scanf("%d", & px); 78int ans = sum(r[px]) - sum(l[px]-1); 79printf("%d\\n", ans); 80} 81} 82return 0; 83 } 84 */ 85 86 ///建树+树状数组(静态连接表) 4940k 454ms 87 #include < cstdio> 88 #include < iostream> 89 #include < cstring> 90 #include < vector> 91 #include < algorithm> 92 #define INF 0x3f3f3f3f 93 using namespace std; 94 95 const int MAXN = 1e5+10; 96 const int MAXM = 1e5+10; 97 98 struct node 99 { 100int to, next; 101 }edge[MAXM< < 1]; 102 103 int head[MAXN]; 104 bool vis[MAXN]; 105 int t[MAXN], l[MAXN], r[MAXN]; 106 int N, M, cnt, key; 107 108 int lowbit(int x) 109 { 110return x& (-x); 111 } 112 void add(int x, int value) 113 { 114for(int i = x; i < = N; i+=lowbit(i)) 115t[i]+=value; 116 } 117 int sum(int x) 118 { 119int res = 0; 120for(int i = x; i > 0; i-=lowbit(i)) 121res+=t[i]; 122return res; 123 } 124 125 void dfs(int k) 126 { 127l[k] = key; 128for(int i = head[k]; i != -1; i = edge[i].next) 129{ 130if(l[edge[i].to]) continue; 131key++; 132dfs(edge[i].to); 133} 134r[k] = key; 135 } 136 137 void init() 138 { 139memset(head, -1, sizeof(head)); 140cnt = 1; 141 } 142 143 void add_e(int u, int v) 144 { 145edge[cnt].to = v; 146edge[cnt].next = head[u]; 147head[u] = cnt++; 148 } 149 150 int main() 151 { 152scanf("%d", & N); 153init(); 154int a, b; 155for(int i = 1; i < = N; i++) 156{ 157vis[i] = 1; //每个节点一开始都有苹果 158add(i, 1); //初始化树状数组 159} 160for(int i = 1; i < N; i++) 161{ 162scanf("%d%d", & a, & b); ///建无向图 163add_e(a, b); 164add_e(b, a); 165} 166key = 1; dfs(1); //建树 167scanf("%d", & M); 168int px; 169char com[3]; 170while(M--) 171{ 172scanf("%s", & com); 173if(com[0] == \'C\') 174{ 175scanf("%d", & px); 176if(vis[px]) add(l[px], -1); 177else add(l[px], 1); 178vis[px] = !vis[px]; 179} 180else 181{ 182scanf("%d", & px); 183printf("%d\\n", sum(r[px])-sum(l[px]-1)); 184} 185} 186return 0; 187 }

 

    推荐阅读