[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;
}

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章