BZOJ 2594 【WC2006】 水管局长数据加强版

题目链接:水管局长数据加强版

  好久没写博客了……

  上次考试的时候写了一发LCT,但是写挂了……突然意识到我已经很久没有写过LCT了,于是今天找了道题来练练手。

  首先,LCT这里不讲。这道题要求支持动态的查询两点间路径的最大值,还有删边的操作。很显然,把操作离线,倒过来做,删边变成了加边,那么只需要用LCT来动态维护最小生成树即可。和魔法森林那道题一样,把边看做点,丢到LCT中去,每次加边的时候如果形成了环,那么判断一下这条边是否需要替换掉环上的最大边,LCT即可。

  求删去的边的编号时可以通过给输入的边排序,每次在其中二分查找。注意边的标号不要弄混了。

  下面贴代码(自带巨大常数):

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
#define maxn 1100010

using namespace std;
typedef long long llg;

struct data{
	int u,v,c,b,id;
}a[maxn],b[maxn];
int s[maxn][2],addv[maxn],maxv[maxn],ff[maxn];
int val[maxn],d[maxn],n,m,q,fa[maxn],ans[maxn];
bool rev[maxn];

int getint(){
	int w=0;bool q=0;
	char c=getchar();
	while((c>'9'||c<'0')&&c!='-') c=getchar();
	if(c=='-') c=getchar(),q=1;
	while(c>='0'&&c<='9') w=w*10+c-'0',c=getchar();
	return q?-w:w;
}

bool cmp1(data a,data b){return a.u==b.u?a.v<b.v:a.u<b.u;}
bool cmp2(data a,data b){return a.b==b.b?a.c<b.c:a.b<b.b;}

bool isroot(int u){return u!=s[fa[u]][0] && u!=s[fa[u]][1];}
void update(int u){
	int l=s[u][0],r=s[u][1]; maxv[u]=u;
	if(val[maxv[l]]>val[maxv[u]]) maxv[u]=maxv[l];
	if(val[maxv[r]]>val[maxv[u]]) maxv[u]=maxv[r];
}

void rotate(int u){
	int p=fa[u],g=fa[p];
	bool l=(u==s[p][1]),r=!l;
	if(!isroot(p)) s[g][p==s[g][1]]=u;
	fa[s[u][r]]=p; s[p][l]=s[u][r];
	s[u][r]=p; fa[p]=u; fa[u]=g;
	update(p); update(u);
}

void splay(int u){
	d[d[0]=1]=u;
	for(int i=u;!isroot(i);i=fa[i]) d[++d[0]]=fa[i];
	for(int i=d[0],x;x=d[i],i;i--)
		if(rev[x]){
			swap(s[x][0],s[x][1]); rev[x]=0;
			rev[s[x][0]]^=1; rev[s[x][1]]^=1; 
		}
	while(!isroot(u)){
		int p=fa[u],g=fa[u];
		if(!isroot(p)){
			if((u==s[p][0])^(p==s[g][1])) rotate(p);
			else rotate(u);
		}
		rotate(u);
	}
}

int getf(int x){return ff[ff[x]]==ff[x]?ff[x]:ff[x]=getf(ff[x]);}
void access(int u){for(int t=0;u;t=u,u=fa[u]) splay(u),s[u][1]=t,update(u);}
void makert(int u){access(u);splay(u);rev[u]^=1;}
void link(int u,int v){makert(u); fa[u]=v;}
void cut(int u,int v){
	makert(u); access(v); splay(v);
	s[v][0]=fa[u]=0; update(v);
}

int query(int u,int v){
	makert(u); access(v); splay(v);
	return maxv[v];
}

int find(int u,int v){
	int l=1,r=m,mid;
	while(l!=r){
		mid=(l+r)>>1;
		if(u<a[mid].u || (u==a[mid].u && v<=a[mid].v)) r=mid;
		else l=mid+1;
	}
	return l;
}

int main(){
	File("a");
	n=getint(); m=getint(); q=getint();
	for(int i=1;i<=n;i++) ff[i]=i;
	for(int i=1;i<=m;i++){
		a[i].u=getint(),a[i].v=getint(),a[i].c=getint();
		if(a[i].u>a[i].v) swap(a[i].u,a[i].v);
	}
	sort(a+1,a+m+1,cmp1);
	for(int i=1;i<=m;i++) val[i+n]=a[i].c,a[i].id=i;
	for(int i=1;i<=q;i++){
		b[i].c=getint(),b[i].u=getint(),b[i].v=getint();
		if(b[i].u>b[i].v) swap(b[i].u,b[i].v);
		if(b[i].c==2) b[i].b=find(b[i].u,b[i].v),a[b[i].b].b=1;
	}
	sort(a+1,a+m+1,cmp2);
	for(int i=1,u,v;i<=m;i++){
		if(a[i].b) break;
		u=a[i].u; v=a[i].v;
		if(getf(u)!=getf(v)){
			ff[getf(u)]=getf(v);
			link(u,a[i].id+n);
			link(v,a[i].id+n);
		}
	}
	sort(a+1,a+m+1,cmp1);
	for(int i=q,u,v,t;i;i--){
		u=b[i].u,v=b[i].v;
		if(b[i].c==1) ans[i]=val[query(u,v)];
		else{
			if(getf(u)==getf(v)){
				t=query(u,v);
				if(val[t]>a[b[i].b].c) cut(a[t-n].u,t),cut(a[t-n].v,t);
				else continue;
			}
			link(u,b[i].b+n); link(v,b[i].b+n);
		}
	}
	for(int i=1;i<=q;i++) if(b[i].c==1) printf("%d\n",ans[i]);
	return 0;
}
posted @ 2017-02-02 20:47  lcf2000  阅读(171)  评论(0编辑  收藏  举报