[BZOJ 1095] [ZJOI 2007]Hide 捉迷藏

在BZ上连续MLE n次后,终于A了.

自己YY的动态点分写法,思路还是很清楚的,但是比较卡内存.

用到了MAP导致复杂度比其他的代码多了一个log,看来需要去借鉴一下别人怎么写的.

updata in 2017-05-25:

发现了一些没必要储存的东西.

1. 存储当前重心某子树堆的位置的MAP可以利用在重心树上的儿子记录此信息.

2.求树上距离可以用RMQ,没必要直接拿MAP记下来.

改完这两个令我颇不舒服的地方,复杂度就正常多了.

 /*
 动态树分治
 BZOJ1095
 全局路径查询
 单点修改,全局查询
 */
 #include<cstdio>
 #include<iostream>
 #include<cstdlib>
 #include<cstring>
 #include<ctime>
 #include<cmath>
 #include<vector>
 #include<map>
 #include<set>
 #include<bitset>
 #include<algorithm>
 #include<iomanip>
 #include<queue>
 #include<stack>
 using namespace std;
 #define ll long long
 #define up(i,j,n) for(int i=j;i<=n;++i)
 #define db double
 #define piii pair< pair<int,int>,int>
 #define pb push_back
 #define FILE "dealing"
 #define eps 1e-8
 template<class T> inline bool cmin(T&a,T b){return a>b?a=b,true:false;}
 template<class T> inline bool cmax(T&a,T b){return a<b?a=b,true:false;}
 template<class T> inline T squ(T a){return a*a;}
 const int maxn=+,base=,limit=3e8,inf=;
 int read(){
     int x=,f=,ch=getchar();
     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>=''&&ch<='')x=(x<<)+(x<<)+ch-'',ch=getchar();
     return x*f;
 }
 int n,Q;
 vector<int>E[maxn];
 struct dui{
     priority_queue<int> q,q2;
     void pop(){
         while(!q.empty()&&!q2.empty()&&q.top()==q2.top())q.pop(),q2.pop();
         if(!q.empty())q.pop();
     }
     int top(){
         while(!q.empty()&&!q2.empty()&&q.top()==q2.top())q.pop(),q2.pop();
         if(!q.empty())return q.top();
         else return -;
     }
     void erase(int x){
         q2.push(x);
         while(!q.empty()&&!q2.empty()&&q.top()==q2.top())q.pop(),q2.pop();
     }
     void push(int x){
         q.push(x);
         while(!q.empty()&&!q2.empty()&&q.top()==q2.top())q.pop(),q2.pop();
     }
     int getans(){
         int s1=top();pop();
         int s2=top();push(s1);
         if(s1==-)return -;
         if(s2==-)return -;
         return s1+s2;
     }
 }All,DA[maxn],MA[maxn<<];int cnt=;
 int siz[maxn],vis[maxn],Mx[maxn],dep[maxn];
 int getsize(int x,int fa){
     siz[x]=;
     for(int i=;i<E[x].size();i++)if(E[x][i]!=fa&&!vis[E[x][i]]){
         getsize(E[x][i],x);
         siz[x]+=siz[E[x][i]];
     }
 }
 void getroot(int x,int fa,int& rt,int S){
     Mx[x]=;
     for(int i=;i<E[x].size();i++)if(E[x][i]!=fa&&!vis[E[x][i]]){
         getroot(E[x][i],x,rt,S);
         cmax(Mx[x],siz[E[x][i]]);
     }
     cmax(Mx[x],S-siz[x]);
     if(Mx[x]<Mx[rt])rt=x;
 }
 int fa[maxn];
 vector<int> F[maxn];
 map<int,int>A[maxn],B[maxn];
 void DFS(int x,int fa,int Cnt,int FA){
     A[FA][x]=Cnt;
     B[FA][x]=dep[x];
     MA[Cnt].push(dep[x]);
     for(int i=;i<E[x].size();i++)if(E[x][i]!=fa&&!vis[E[x][i]]){
         dep[E[x][i]]=dep[x]+;
         DFS(E[x][i],x,Cnt,FA);
     }
 }
 int work(int x,int S){
     int rt=;
     getroot(x,,rt,siz[x]);
     vis[rt]=;
     getsize(rt,);x=rt;
     cnt++;F[x].push_back(cnt);MA[cnt].push();
     A[x][x]=cnt;B[x][x]=;
     for(int i=;i<E[x].size();i++)if(!vis[E[x][i]]){
         cnt++;
         F[x].push_back(cnt);
         dep[E[x][i]]=;
         DFS(E[x][i],,cnt,x);
     }
     for(int i=;i<F[x].size();i++)
         DA[x].push(MA[F[x][i]].top());
     int val=DA[x].getans();
     if(val!=-&&val!=-)All.push(val);
     for(int i=;i<E[x].size();i++)if(!vis[E[x][i]]){
         int y=work(E[x][i],siz[E[x][i]]);
         fa[y]=x;
     }
     return x;
 }
 char s[];
 int f[maxn],sum;
 void motify(int x,int FA){
     while(x){
         int val=DA[x].getans();
         if(val!=-&&val!=-)All.erase(val);

         int id=A[x][FA];
         val=MA[id].top();if(val!=-)DA[x].erase(val);
         if(!f[FA])MA[id].push(B[x][FA]);
         else MA[id].erase(B[x][FA]);
         val=MA[id].top();if(val!=-)DA[x].push(val);

         val=DA[x].getans();
         if(val!=-&&val!=-)All.push(val);
         x=fa[x];
     }
 }
 int main(){
     freopen(FILE".in","r",stdin);
     freopen(FILE".out","w",stdout);
     n=read();Mx[]=inf;sum=n;
     up(i,,n){
         int x=read(),y=read();
         E[x].push_back(y);
         E[y].push_back(x);
     }
     getsize(,);
     work(,siz[]);
     Q=read();
     while(Q--){
         scanf("%s",s+);
         if(s[]=='G'){
             if(sum==)printf("-1\n");
             else if(sum==)printf("0\n");
             else printf("%d\n",All.top());
         }
         if(s[]=='C'){
             int x=read();
             f[x]^=;
             if(!f[x])sum++;
             else sum--;
             motify(x,x);
         }
     }
     return ;
 }