【BZOJ2594】【WC2006】水管局长

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<algorithm>
  4 #include<cstring>
  5 #include<cmath>
  6 using namespace std;
  7 int rd(){int z=0;  char ch=getchar();
  8     while(ch<‘0‘||ch>‘9‘)  ch=getchar();
  9     while(ch>=‘0‘&&ch<=‘9‘){z=(z<<3)+(z<<1)+ch-‘0‘;  ch=getchar();}
 10     return z;
 11 }
 12 struct nds{int x,y,z,id;}e[1100000],q[110000];
 13 int n,m,o;
 14 int fth[1110000],chd[1110000][2],mxv[1110000],v[1110000],mxid[1110000],rvs[1110000];
 15 int stck[1110000],tp=0;
 16 bool flg[1100000];  int tl[1100000],tr[1100000];
 17 int ans[110000],ast=0;
 18 inline bool isrt(int x){  return (chd[fth[x]][0]!=x)&(chd[fth[x]][1]!=x);}
 19 inline void pshu(int x){
 20     if(!x)  return ;
 21     mxv[x]=v[x],mxid[x]=x;
 22     if(mxv[chd[x][0]]>mxv[x])  mxv[x]=mxv[chd[x][0]],mxid[x]=mxid[chd[x][0]];
 23     if(mxv[chd[x][1]]>mxv[x])  mxv[x]=mxv[chd[x][1]],mxid[x]=mxid[chd[x][1]];
 24 }
 25 inline void pshd(int x){
 26     if(!rvs[x])  return ;
 27     rvs[chd[x][0]]^=1,rvs[chd[x][1]]^=1,rvs[x]=0;
 28     swap(chd[x][0],chd[x][1]);
 29 }
 30 void rtt(int x){
 31     int y=fth[x],z=fth[fth[x]],l,r;
 32     r=(chd[y][0]==x);  l=r^1;
 33     if(!isrt(y))  chd[z][chd[z][1]==y]=x;
 34     fth[x]=z,fth[y]=x,fth[chd[x][r]]=y;
 35     chd[y][l]=chd[x][r],chd[x][r]=y;
 36     pshu(y),pshu(x);
 37 }
 38 void sply(int x){
 39     stck[tp=1]=x;
 40     for(int i=x;!isrt(i);i=fth[i])  stck[++tp]=fth[i];
 41     while(tp)  pshd(stck[tp--]);
 42     while(!isrt(x)){
 43         if(!isrt(fth[x]))  rtt((chd[fth[x]][0]==x)^(chd[fth[fth[x]]][0]==fth[x])?x:fth[x]);
 44         rtt(x);
 45     }
 46 }
 47 inline void accs(int x){  for(int i=0;x;sply(x),chd[x][1]=i,pshu(x),x=fth[i=x]);}
 48 inline void qdrt(int x){  accs(x),sply(x),rvs[x]^=1;}
 49 inline void lk(int x,int y){  qdrt(x),fth[x]=y,sply(x);}
 50 inline void ct(int x,int y){  qdrt(x),accs(y),sply(y),fth[x]=chd[y][0]=0,pshu(y);}
 51 inline int gtrt(int x){  while(fth[x])  x=fth[x];  return x;}
 52 inline int qry(int x,int y){  qdrt(x),accs(y),sply(y);
 53     return mxid[y]-n;}
 54 bool cmp1(nds x,nds y){  return x.x==y.x ? x.y<y.y : x.x<y.x;}
 55 bool cmp2(nds x,nds y){  return x.z<y.z;}
 56 bool cmp3(nds x,nds y){  return x.id<y.id;}
 57 int bnrsch(int x,int y){
 58     int l=tl[x],r=tr[x],md;
 59     while(l+1<r)  md=(l+r)>>1,(e[md].y<=y?l:r)=md;
 60     //while(l+1<r)  md=(l+r)>>1,(e[md].y>=y?l:r)=md;
 61     return e[l].y==y?l:r;
 62 }
 63 void krsk(){
 64     int bwl=0;
 65     sort(e+1,e+m+1,cmp2);
 66     for(int i=1;i<=m;++i)if((!flg[e[i].id])&(gtrt(e[i].x)!=gtrt(e[i].y))){
 67         lk(e[i].x,e[i].id+n),lk(e[i].y,e[i].id+n);
 68         if(++bwl==n-1)  break;
 69     }
 70 }
 71 void ist(int x){
 72     if(gtrt(e[x].x)!=gtrt(e[x].y)){
 73         lk(e[x].x,x+n),lk(e[x].y,x+n);
 74         return ;
 75     }
 76     int tmp=qry(e[x].x,e[x].y);
 77     if(e[tmp].z<e[x].z)  return ;
 78     ct(e[tmp].x,tmp+n),ct(e[tmp].y,tmp+n);
 79     mxv[tmp+n]=v[tmp+n]=0;
 80     lk(e[x].x,x+n),lk(e[x].y,x+n);
 81 }
 82 int main(){//freopen("ddd.in","r",stdin);
 83     //freopen("ddd.out","w",stdout);
 84     cin>>n>>m>>o;
 85     for(int i=1;i<=m;++i){
 86         e[i].x=rd(),e[i].y=rd(),e[i].z=rd();
 87         e[i].id=i,mxv[i+n]=v[i+n]=e[i].z,mxid[i+n]=i+n;
 88     }
 89     sort(e+1,e+m+1,cmp1);
 90     for(int i=1;i<=m+1;++i)if(e[i].x!=e[i-1].x)  tr[e[i-1].x]=i-1,tl[e[i].x]=i;
 91     for(int i=1;i<=o;++i){
 92         q[i].z=rd(),q[i].x=rd(),q[i].y=rd();
 93         if(q[i].z==2)  flg[q[i].id=e[bnrsch(q[i].x,q[i].y)].id]=1;
 94     }
 95     krsk();
 96     sort(e+1,e+m+1,cmp3);
 97     for(int i=o;i>=1;--i){
 98         if(q[i].z==2)  ist(q[i].id);
 99         else  ans[++ast]=e[qry(q[i].x,q[i].y)].z;
100     }
101     while(ast)  printf("%d\n",ans[ast--]);
102     return 0;
103 }