bzoj 2594 [Wc2006]水管局长数据加强版（LCT+最小生成树）

【深坑勿入】

【给个链接】

　　http://blog.csdn.net/popoqqq/article/details/41348549

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #define FOR(a,b,c) for(int a=b;a<=c;a++)
 #define trav(u,i) for(int i=front[u];i;i=e[i].nxt)
 using namespace std;

 typedef long long ll;
 typedef unsigned int ul;
 const int N = 4e5+;

 ll read() {
     char c=getchar();
     ll f=,x=;
     while(!isdigit(c)) {
         if(c=='-') f=-; c=getchar();
     }
     while(isdigit(c))
         x=x*+c-'',c=getchar();
     return x*f;
 }

 struct Edge {
     int u,v,w;
     bool operator < (const Edge& rhs) const {
         return u<rhs.u||(u==rhs.u&&v<rhs.v);
     }
 }e[N<<];

 namespace LCT {

     struct Node {
         Node *ch[],*fa;
         int rev,id,maxe;
         Node() ;
         void reverse() {
             rev^=;
             swap(ch[],ch[]);
         }
         void up_push() {
             if(fa->ch[]==this||fa->ch[]==this)
                 fa->up_push();
             if(rev) {
                 ch[]->reverse();
                 ch[]->reverse();
                 rev=;
             }
         }
         void maintain() {
             int _max=-;
             if(e[ch[]->maxe].w>_max)
                 _max=e[ch[]->maxe].w,maxe=ch[]->maxe;
             if(e[ch[]->maxe].w>_max)
                 _max=e[ch[]->maxe].w,maxe=ch[]->maxe;
             if(e[id].w>_max) maxe=id;
         }
     } *null=new Node,T[N],E[N<<];
     Node::Node() {
         id=maxe=rev=;
         fa=ch[]=ch[]=null;
     }

     void rot(Node* o,int d) {
         Node *p=o->fa;
         p->ch[d]=o->ch[d^];
         o->ch[d^]->fa=p;
         o->ch[d^]=p;
         o->fa=p->fa;
         if(p==p->fa->ch[])
             p->fa->ch[]=o;
         else if(p==p->fa->ch[])
             p->fa->ch[]=o;
         p->fa=o;
         p->maintain();
     }
     void splay(Node *o) {
         o->up_push();
         Node *nf,*nff;
         while(o->fa->ch[]==o||o->fa->ch[]==o) {
             nf=o->fa,nff=nf->fa;
             if(o==nf->ch[]) {
                 if(nf==nff->ch[]) rot(nf,);
                 rot(o,);
             } else {
                 if(nf==nf->ch[]) rot(nf,);
                 rot(o,);
             }
         }
         o->maintain();
     }
     void Access(Node* o) {
         Node *son=null;
         while(o!=null) {
             splay(o);
             o->ch[]=son;
             o->maintain();
             son=o; o=o->fa;
         }
     }
     void evert(Node* o) {
         Access(o);
         splay(o);
         o->reverse();
     }
     void Link(Node* u,Node* v) {
         evert(u);
         u->fa=v;
     }
     void Cut(Node* u,Node* v) {
         evert(u);
         Access(v); splay(v);
         u->fa=v->ch[]=null;
         v->maintain();
     }
     Node* find(Node* o) {
         while(o->fa!=null) o=o->fa;
         return o;
     }

 }
 using namespace LCT;

 int n,m,q,flag[N],ans[N];

 struct Q { int op,x,y;
 }que[N];

 int query(Node* u,Node* v)
 {
     if(find(u)!=find(v)) return ;
     evert(u);
     Access(v),splay(v);
     return v->maxe;
 }
 void insert(Node* u,Node* v,int id)
 {
     if(find(u)==find(v)) {
         int maxe=query(u,v);
         if(e[maxe].w<=e[id].w) return ;
         Cut(&E[maxe],&T[e[maxe].u]);
         Cut(&E[maxe],&T[e[maxe].v]);
     }
     Link(&E[id],u);
     Link(&E[id],v);
 }
 int search(Edge x)
 {
     int l=,r=m,mid;
     while(l<=r) {
         mid=l+r>>;
         if((!(e[mid]<x))&&(!(x<e[mid]))) return mid;
         if(x<e[mid]) r=mid-;
         else l=mid+;
     }
     return ;
 }

 int main()
 {
     freopen("in.in","r",stdin);
     freopen("out.out","w",stdout);
     n=read(),m=read(),q=read();
     FOR(i,,m) E[i].maxe=E[i].id=i;
     FOR(i,,m) {
         e[i].u=read(),e[i].v=read(),e[i].w=read();
         if(e[i].u>e[i].v) swap(e[i].u,e[i].v);
     }
     e[].w=;
     sort(e+,e+m+);
     FOR(i,,q) {
         que[i].op=read(),que[i].x=read(),que[i].y=read();
         if(que[i].x>que[i].y) swap(que[i].x,que[i].y);
         if(que[i].op==) {
             int x=search((Edge){que[i].x,que[i].y,});
             if(x==) puts("error");
             flag[x]=;
         }
     }
     FOR(i,,m) if(!flag[i]) {
         insert(&T[e[i].u],&T[e[i].v],i);
     }
     for(int i=q;i;i--) {
         if(que[i].op==) {
             int x=search((Edge){que[i].x,que[i].y,});
             insert(&T[e[x].u],&T[e[x].v],x);
         } else {
             ans[i]=e[query(&T[que[i].x],&T[que[i].y])].w;
         }
     }
     FOR(i,,q) if(que[i].op==)
         printf("%d\n",ans[i]);
     return ;
 }