码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ1095 [ZJOI2007]Hide 捉迷藏

时间:2016-08-01 01:40:23      阅读:255      评论:0      收藏:0      [点我收藏+]

标签:

  动态树分治,用三个set分别维护每个重心到每一个子树的距离种类、每个重心所有子树的最大值和次大值、全局答案的最大值。复杂度O(nlogn^2)

  代码

  1 #include<cstdio>
  2 #include<algorithm>
  3 #include<vector>
  4 #include<set>
  5 #define pb push_back
  6 using namespace std;
  7 const int N = 501010;
  8 int size[N],flag[N],n,a,b,i,f[N],cnt,typ[N],deep[N];
  9 int jump[N][18];
 10 int q;
 11 char ch[5];
 12 vector<int> e[N];
 13 multiset<int> dist[N],ans[N],Ans;
 14 multiset<int>::iterator it;
 15 int getroot(int x,int fa,int n)
 16 {
 17     int i,tmp=0,rt=0;
 18     size[x]=1;
 19     for (i=0;i<e[x].size();i++)
 20     if ((e[x][i]!=fa)&&(!flag[e[x][i]]))
 21     {
 22         rt|=getroot(e[x][i],x,n);
 23         size[x]+=size[e[x][i]];
 24         if (size[e[x][i]]>n/2) tmp=1;
 25     }
 26     if (n-size[x]>n/2) tmp=1;
 27     if (tmp==1) return rt;else return x;
 28 }
 29 void gao(int x,int fa,int y,int dis)
 30 {
 31     dis++;
 32     int i;
 33     dist[y].insert(dis);
 34     for (i=0;i<e[x].size();i++)
 35     if ((!flag[e[x][i]])&&(e[x][i]!=fa))
 36     gao(e[x][i],x,y,dis);
 37 }
 38 void DFS(int x,int fa)
 39 {
 40     int i;
 41     deep[x]=deep[fa]+1;
 42     jump[x][0]=fa;
 43     for (i=1;i<=17;i++)
 44     jump[x][i]=jump[jump[x][i-1]][i-1];
 45     for (i=0;i<e[x].size();i++)
 46     if (e[x][i]!=fa) DFS(e[x][i],x);
 47 }
 48 int lca(int a,int b)
 49 {
 50     int i;
 51     if (deep[a]<deep[b]) a^=b^=a^=b;
 52     for (i=17;i>=0;i--)
 53     if (deep[jump[a][i]]>=deep[b]) a=jump[a][i];
 54     if (a==b) return a;
 55     for (i=17;i>=0;i--)
 56     if (jump[a][i]!=jump[b][i]) a=jump[a][i],b=jump[b][i];
 57     return jump[a][0];
 58 }
 59 int getdis(int x,int y)
 60 {
 61     int z=lca(x,y);
 62     return deep[x]+deep[y]-2*deep[z];
 63 }
 64 void del(int x)
 65 {
 66     int tmp=0;
 67     if (ans[x].size()>=2)
 68     {
 69         it=ans[x].end();
 70         --it;tmp+=*it;
 71         --it;tmp+=*it;
 72         it=Ans.lower_bound(tmp);
 73         Ans.erase(it); 
 74     }
 75 }
 76 void ins(int x)
 77 {
 78     int tmp=0;
 79     if (ans[x].size()>=2)
 80     {
 81         it=ans[x].end();
 82         --it;tmp+=*it;
 83         --it;tmp+=*it;
 84         Ans.insert(tmp); 
 85     }
 86 }
 87 void DEL(int y)
 88 {
 89     if (dist[y].size())
 90     {
 91         it=dist[y].end();--it;
 92         int tmp=*it;
 93         it=ans[f[y]].lower_bound(tmp);
 94         ans[f[y]].erase(it);
 95     }
 96 } 
 97 void INS(int y)
 98 {
 99     if (dist[y].size())
100     {
101         it=dist[y].end();--it;
102         int tmp=*it;
103         ans[f[y]].insert(tmp);
104     }
105 }
106 int change(int x,int y)
107 {
108     del(x);
109     if (typ[x]==1) 
110     {
111         it=ans[x].lower_bound(0);
112         ans[x].erase(it);
113     } 
114     else
115         ans[x].insert(0);
116     while (x)
117     {
118         ins(x);
119         if (f[x]) 
120         { 
121             del(f[x]);
122             DEL(x);
123             int dis=getdis(f[x],y);
124             if (typ[y]==1)
125             {
126                 it=dist[x].lower_bound(dis);
127                 dist[x].erase(it);
128             }
129             else
130                 dist[x].insert(dis);
131             INS(x);
132         } 
133         x=f[x];
134     }
135 }
136 int dfs(int x,int n,int fa)
137 {
138     int i,y;
139     x=getroot(x,0,n);
140     f[x]=fa;
141     flag[x]=1;
142     ans[x].insert(0); 
143     for (i=0;i<e[x].size();i++)
144     if (!flag[e[x][i]])
145     {
146         if (size[e[x][i]]>size[x])
147             y=dfs(e[x][i],n-size[x],x);
148         else
149             y=dfs(e[x][i],size[e[x][i]],x);
150         gao(e[x][i],0,y,0);
151         it=dist[y].end();
152         ans[x].insert(*(--it));
153     }
154     ins(x);
155     flag[x]=0;
156     return x;
157 }
158 int main()
159 {
160     scanf("%d",&n);
161     for (i=1;i<n;i++)
162     {
163         scanf("%d%d",&a,&b);
164         e[a].pb(b);
165         e[b].pb(a);
166     }
167     dfs(1,n,0);
168     DFS(1,0);
169     scanf("%d",&q);
170     for (i=1;i<=q;i++)
171     {
172         scanf("%s",ch+1);
173         if (ch[1]==G)
174         {
175             if (cnt==n)
176             printf("-1\n");
177             else
178             if (cnt==n-1)
179             printf("0\n");
180             else
181             {
182                 it=Ans.end();--it;
183                 printf("%d\n",*it);
184             }
185         }
186         else
187         if (ch[1]==C)
188         {
189             scanf("%d",&a);
190             if (typ[a]==0) 
191             cnt++,typ[a]=1;
192             else 
193             cnt--,typ[a]=0;
194             change(a,a);
195         }
196     }
197 } 

 

BZOJ1095 [ZJOI2007]Hide 捉迷藏

标签:

原文地址:http://www.cnblogs.com/fzmh/p/5724374.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!