hdu_5221_Occupation(树剖)

题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=5221

题意:给你一棵树,每个节点有一定的值,有三种操作: 1 x y 表示占领树上x-y的所有节点,2 x 表示去掉这点,3 x表示占领以这个节点为根的所有节点

题解:1 2 都是树剖操作,3则要记录树剖的时间戳,一个根节点下的点的新编号都是连续的,所以只要记录该根节点的最大时间戳就知道这个区间了

 1 #include<cstdio>
 2 #pragma comment(linker, "/STACK:102400000,102400000")
 3 #define F(i,a,b) for(int i=a;i<=b;i++)
 4 #define root 1,n,1
 5 #define ls l,m,rt<<1
 6 #define rs m+1,r,rt<<1|1
 7 
 8 const int N=1e5+7;
 9 int dep[N],sz[N],hs[N],top[N],tid[N],fid[N],fa[N],idx,mp[N];
10 int t,n,m,pp,x,y,a[N],sum[N*4],vis[4*N],ans[N*4],g[N],nxt[2*N],v[2*N],ed;
11 
12 inline int max(int a,int b){return a>b?a:b;}
13 inline void adg(int x,int y){v[++ed]=y,nxt[ed]=g[x],g[x]=ed;}
14 //线段树
15 void build(int l,int r,int rt){
16     vis[rt]=ans[rt]=0;
17     if(l==r){sum[rt]=a[fid[l]];return;}
18     int m=(l+r)>>1;
19     build(ls),build(rs);
20     sum[rt]=sum[rt<<1]+sum[rt<<1|1];
21 }
22 
23 inline void pd(int rt){
24     if(vis[rt])
25     vis[rt<<1]=vis[rt<<1|1]=1,ans[rt<<1]=sum[rt<<1],ans[rt<<1|1]=sum[rt<<1|1],vis[rt]=0;
26 }
27 void update(int L,int R,int op,int l,int r,int rt){
28     if(L<=l&&r<=R){
29         vis[rt]=op,ans[rt]=op?sum[rt]:0;
30         return;
31     }
32     pd(rt);
33     int m=(l+r)>>1;
34     if(L<=m)update(L,R,op,ls);
35     if(R>m)update(L,R,op,rs);
36     ans[rt]=ans[rt<<1]+ans[rt<<1|1];
37 }
38 //树剖
39 void dfs1(int u,int pre){
40     dep[u]=dep[pre]+1,hs[u]=0,fa[u]=pre,sz[u]=1;
41     for(int i=g[u];i;i=nxt[i])if(v[i]!=pre)
42     dfs1(v[i],u),sz[u]+=sz[v[i]],hs[u]=(sz[v[i]]>sz[hs[u]])?v[i]:hs[u];
43 }
44 void dfs2(int u,int tp){
45     top[u]=tp,tid[u]=++idx,fid[idx]=u,mp[u]=idx;
46     if(hs[u])dfs2(hs[u],tp),mp[u]=max(mp[u],mp[hs[u]]);
47     for(int i=g[u];i;i=nxt[i])
48     if(v[i]!=fa[u]&&v[i]!=hs[u])dfs2(v[i],v[i]),mp[u]=max(mp[u],mp[v[i]]);
49 }
50 
51 void up(int x,int y){
52     int fx=top[x],fy=top[y];
53     while(fx!=fy){
54         if(dep[fx]>dep[fy])update(tid[fx],tid[x],1,root),x=fa[fx],fx=top[x];
55         else update(tid[fy],tid[y],1,root),y=fa[fy],fy=top[y];
56     }
57     if(dep[x]>dep[y])x^=y,y^=x,x^=y;
58     update(tid[x],tid[y],1,root);
59 }
60 
61 int main(){
62     scanf("%d",&t);
63     while(t--){
64         scanf("%d",&n);
65         F(i,1,n)scanf("%d",a+i);
66         F(i,1,n)g[i]=0;ed=0;
67         F(i,1,n-1)scanf("%d%d",&x,&y),adg(x,y),adg(y,x);
68         dfs1(1,0),idx=0,dfs2(1,1),build(root);
69         scanf("%d",&m);
70         F(i,1,m){
71             scanf("%d",&pp);
72             if(pp==1)scanf("%d%d",&x,&y),up(x,y);
73             else if(pp==2)scanf("%d",&x),update(tid[x],tid[x],0,root);
74             else scanf("%d",&x),update(tid[x],mp[x],1,root);
75             printf("%d\n",ans[1]);
76         }
77     }
78     return 0;
79 }
View Code

 

posted @ 2016-06-28 18:32  bin_gege  阅读(153)  评论(0编辑  收藏  举报