BZOJ1103 [POI2007]大都市meg
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
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
1
0
1
HINT
正解:树链剖分
解题报告:
链剖裸题。dfs序+树状数组也可以做。
1 //It is made by jump~ 2 #include <iostream> 3 #include <cstdlib> 4 #include <cstring> 5 #include <cstdio> 6 #include <cmath> 7 #include <algorithm> 8 using namespace std; 9 typedef long long LL; 10 const int MAXN = 250011; 11 const int inf = (1<<30); 12 int n,total,ecnt,ans; 13 int id[MAXN],ql,qr; 14 int pre[MAXN],TT; 15 int top[MAXN],siz[MAXN],zhongerzi[MAXN],father[MAXN],deep[MAXN]; 16 int next[MAXN*2],to[MAXN*2],first[MAXN]; 17 char ch[12]; 18 19 struct node{ 20 int sum; 21 }a[MAXN*4]; 22 inline void link(int x,int y){ next[++ecnt]=first[x]; first[x]=ecnt; to[ecnt]=y; } 23 inline int getint() 24 { 25 int w=0,q=0; 26 char c=getchar(); 27 while((c<'0' || c>'9') && c!='-') c=getchar(); 28 if (c=='-') q=1, c=getchar(); 29 while (c>='0' && c<='9') w=w*10+c-'0', c=getchar(); 30 return q ? -w : w; 31 } 32 33 inline void build(int root,int l,int r){ 34 if(l==r){ if(l!=1) a[root].sum=1; else a[root].sum=0; return ; } 35 int lc=root*2,rc=lc+1; 36 int mid=(l+r)/2; build(lc,l,mid); build(rc,mid+1,r); 37 a[root].sum=a[lc].sum+a[rc].sum; 38 } 39 40 inline void dfs1(int u,int fa){ 41 siz[u]=1; 42 for(int i=first[u];i;i=next[i]){ 43 int v=to[i]; if(v==fa) continue; 44 father[v]=u;deep[v]=deep[u]+1; 45 dfs1(v,u);siz[u]+=siz[v]; 46 if(siz[v]>siz[ zhongerzi[u] ]) zhongerzi[u]=v; 47 } 48 } 49 50 inline void dfs2(int u,int fa){ 51 id[u]=++total; pre[total]=u; 52 if(zhongerzi[u]) { top[zhongerzi[u]]=top[u];dfs2(zhongerzi[u],u); } 53 for(int i=first[u];i;i=next[i]){ 54 int v=to[i]; if(v==fa || v==zhongerzi[u]) continue; 55 top[v]=v; dfs2(v,u); 56 } 57 } 58 59 inline void querysum(int root,int l,int r){ 60 if(ql<=l && r<=qr) { ans+=a[root].sum; return ; } 61 int mid=(l+r)/2; int lc=root*2,rc=lc+1; 62 if(ql<=mid) querysum(lc,l,mid); 63 if(qr>mid) querysum(rc,mid+1,r); 64 } 65 66 inline int findsum(int x){ 67 int f1=top[x]; ans=0; 68 while(f1!=0){ 69 ql=id[f1]; qr=id[x]; querysum(1,1,n); 70 x=father[f1]; f1=top[x]; 71 } 72 return ans; 73 } 74 75 inline void update(int root,int l,int r){ 76 if(l==r) { a[root].sum=0; return ; } 77 int mid=(l+r)/2; int lc=root*2,rc=lc+1; 78 if(TT<=mid) update(lc,l,mid); else update(rc,mid+1,r); 79 a[root].sum=a[lc].sum+a[rc].sum; 80 } 81 82 inline void work(){ 83 n=getint(); int x,y; 84 for(int i=1;i<n;i++){ 85 x=getint();y=getint(); if(x>y) swap(x,y); 86 link(x,y); 87 } 88 deep[1]=1; dfs1(1,0); 89 top[1]=1; dfs2(1,0); 90 build(1,1,n); 91 int Q=getint(); Q+=n-1; 92 93 for(int i=1;i<=Q;i++){ 94 scanf("%s",ch); 95 if(ch[0]=='W'){ 96 x=getint(); 97 printf("%d\n",findsum(x)); 98 } 99 else{ 100 x=getint(); y=getint(); if(x<y) swap(x,y); TT=id[x]; 101 update(1,1,n); 102 } 103 } 104 } 105 106 int main() 107 { 108 work(); 109 return 0; 110 }
本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!