[WC2006]水管局长 lct维护最小生成树

题目背景

SC 省 MY 市有着庞大的地下水管网络,嘟嘟是 MY 市的水管局长(就是管水管的啦)。

题目描述

每天供水公司可能要将一定量的水从 uu 处送往 vv 处,嘟嘟需要为供水公司找到一条从 uu 至 vv 的水管的路径,接着通过信息化的控制中心通知路径上的水管进入准备送水状态,等到路径上每一条水管都准备好了,供水公司就可以开始送水了。嘟嘟一次只能处理一项送水任务,等到当前的送水任务完成了,才能处理下一项。

在处理每项送水任务之前,路径上的水管都要进行一系列的准备操作,如清洗、消毒等等。嘟嘟在控制中心一声令下,这些水管的准备操作同时开始,但由于各条管道的长度、内径不同,进行准备操作需要的时间可能不同。

供水公司总是希望嘟嘟能找到这样一条送水路径,路径上的所有管道全都准备就绪所需要的时间尽量短。嘟嘟希望你能帮助他完成这样的一个选择路径的系统,以满足供水公司的要求。另外,由于 MY 市的水管年代久远,一些水管会不时出现故障导致不能使用,你的程序必须考虑到这一点。

不妨将 MY 市的水管网络看作一幅简单无向图(即没有自环或重边):水管是图中的边,水管的连接处为图中的结点。整张图共有 nn 个节点和 mm 条边,节点从 11 至 nn 编号。

输入格式

第一行有三个整数,分别表示管道连接处(结点)的数目 nn,目前水管(无向边)的数目 mm,以及你的程序需要处理的任务数目(包括寻找一条满足要求的路径和接受某条水管坏掉的事实)qq。

以下 mm 行,每行三个整数 u, v, tu,v,t,表示存在一条连接 (u, v)(u,v) 的水管,准备时间为 tt。

以下 qq 行,每行三个整数 k, u, vk,u,v,描述一项任务。其中 kk 表示任务类型:

  • 若 k = 1k=1,则表示你需要为供水公司寻找一条满足要求的从 uu 到 vv 的水管路径,满足准备时间最短;
  • 若 k = 2k=2,则表示直接连接 uu 和 vv 的水管宣布报废。

输出格式

对于每个 k = 1k=1 的任务,输出一行一个整数表示答案。

输入输出样例

输入 #1复制

4 4 3
1 2 2
2 3 3
3 4 2
1 4 2
1 1 4
2 1 4
1 1 4

输出 #1复制

2
3

由题意可知,任何时刻任何两点之间都是有路径的,显然这些路径可以构成一棵树。然后我们要求的是最大值最小,显然是一棵最小生成树。删边处理起来麻烦,所有我们反过来,进行加边。我们先生成一棵全局的最小生成树,这棵生成树一定是最终的图的最小生成树,然后对于删边操作,因为我们反过来了,就是加边,对于加入的这条边(u,v),它使得u到v变成了一个环,如果(u,v)这条边的权值大于等于(u,v)链上的最大值,那么它的加入没有意义。否则我们把这条边加入,并且去掉(u,v)链上的最大边,维护最小生成树。

另外lct我们一般用来维护点权,但是这题要维护边权。我们就把一条边化成三个点连两条  例如(u,v)这条边 就成 u,v,x+n三个点 其中x是这条边的编号

#include<bits/stdc++.h>
#define R register int 
#define I inline void
#define lc c[x][0]
#define rc c[x][1]
using namespace std;
const int N = 1e6+10;
int c[N][2],st[N],f[N],mx[N],fa[N],val[N],r[N],n,m,Q;
inline int in(){
	int w=0,x=0;char c=0;
	while(c<'0'||c>'9') w|=c=='-',c=getchar();
	while(c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return w?-x:x;
}
inline bool nroot(R x){
	return c[f[x]][0]==x||c[f[x]][1]==x;
} 
I pushup(R x){
	mx[x]=x;
	if(val[mx[lc]]>val[mx[x]]) mx[x]=mx[lc];
	if(val[mx[rc]]>val[mx[x]]) mx[x]=mx[rc];  
}
I pushr(R x){
	swap(lc,rc);r[x]^=1;
}
I pushdown(R x){
	if(r[x]){
		if(lc) pushr(lc);
		if(rc) pushr(rc);
		r[x]=0; 
	}
}
I rotate(R x){
	R y=f[x],z=f[y],k=c[y][1]==x,w=c[x][k^1];
	if(nroot(y)) c[z][c[z][1]==y]=x; c[x][k^1]=y;c[y][k]=w; 
	if(w) f[w]=y; f[y]=x;f[x]=z;
	pushup(y);
}
I splay(R x){
	R y=x,z=0;
	st[++z]=y;
	while(nroot(y)) st[++z]=y=f[y];
	while(z) pushdown(st[z--]);
	while(nroot(x)){
		y=f[x],z=f[y];
		if(nroot(y)) rotate((c[y][1]==x)^(c[z][1]==y)?x:y);
		rotate(x);
	}
	pushup(x);
}
I access(R x){
	for(R y=0;x;x=f[y=x])
	splay(x),rc=y,pushup(x);
}
I makeroot(R x){
	access(x);splay(x);
	pushr(x);
}
inline int findroot(R x){
	access(x);splay(x);
	while(lc) pushdown(x),x=lc;
	splay(x);
	return x; 
}
I split(R x,R y){
	makeroot(x);
	access(y);splay(y);
}
I link(R x,R y){
	makeroot(x);
	if(findroot(y)!=x) f[x]=y; 
}
I cut(R x,R y){
	makeroot(x);
	if(findroot(y)==x&&f[y]==x&&!c[y][0]){
		c[x][1]=f[y]=0;
		pushup(x);
	}	
}
int query(int u,int v){
	split(u,v);return mx[v];
}
struct edge{
	int u,v,w,id,d;
	bool operator < (const edge &a)const{
		if(w==a.w) return id<a.id;//此处不写会出大问题,会导致两次排序不一样(因为边权相同的边是随机排的)
		return w<a.w;
	}
}e[N];
bool cmp(edge a,edge b){
	return a.u<b.u||(a.u==b.u&&a.v<b.v);
}
struct que{
	int u,v,id,op,ans;
}q[N];
int get(R x){
	return x==fa[x]?x:fa[x]=get(fa[x]);
}
inline int find(R u,R v){
	int l=1,r=m;
	while(l<=r){
		int mid = l+r>>1;
		if(e[mid].u<u||(e[mid].u==u&&e[mid].v<v)) l=mid+1;
		else if(e[mid].u==u&&e[mid].v==v) return mid;
		else r=mid-1;
	}
}
int main(){
	n=in(),m=in(),Q=in();
	for(int i = 1; i <= n; i++) fa[i]=i;
	for(int i = 1; i <= m; i++){
		e[i].u=in(),e[i].v=in(),e[i].w=in();
		if(e[i].u>e[i].v) swap(e[i].u,e[i].v);
	}
//	printf("1\n");
	sort(e+1,e+1+m);
	for(int i = 1; i <= m; i++){
		val[i+n]=e[i].w;
		mx[i+n]=i+n;
		e[i].id=i;	
	}
	//printf("2\n");
	sort(e+1,e+1+m,cmp);
	for(int i = 1; i <= Q; i++){
		q[i].op=in(),q[i].u=in(),q[i].v=in();
		if(q[i].op==2){
			if(q[i].u>q[i].v)swap(q[i].u,q[i].v);
			int t=find(q[i].u,q[i].v);
			q[i].id=e[t].id;e[t].d=1;
		}
	}
	//printf("3\n");
	sort(e+1,e+1+m);
	int tot=0;
	for(int i = 1; i <= m; i++){
		if(!e[i].d){
			int x=get(e[i].u),y=get(e[i].v);
			if(x!=y){
				tot++;
				fa[x]=y;
				link(e[i].u,i+n);link(e[i].v,i+n);
				if(tot==n-1) break;
			}
		}
	}
	//printf("4\n");
	for(int i = Q; i >= 1; i--){
		if(q[i].op==1){
			q[i].ans=val[query(q[i].u,q[i].v)];
		}else{
			int t=query(q[i].u,q[i].v),k=q[i].id;
			if(e[k].w<val[t]){
				cut(e[t-n].u,t);cut(e[t-n].v,t);
				link(q[i].u,k+n);link(q[i].v,k+n);
			}
		}
	}
	for(int i = 1; i <= Q; i++) if(q[i].op==1) printf("%d\n",q[i].ans);
	return 0;
}

猜你喜欢

转载自blog.csdn.net/weixin_43824564/article/details/106504602