1. 程式人生 > >[動態樹 LCT] BZOJ 2594 [Wc2006]水管局長資料加強版

[動態樹 LCT] BZOJ 2594 [Wc2006]水管局長資料加強版

動態樹模板題

最開始加邊用Kruskal比較快,不用LCT,然後再用LCT維護最小生成樹

話說這道題小資料範圍的正解是什麼....

PS:以前LCT的求根都是打錯的,沒有pushdown()

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<map>
#include<algorithm>
using namespace std;
typedef pair<int,int> abcd;

inline char nc()
{
	static char buf[100000],*p1=buf,*p2=buf;
	if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
	return *p1++;
}

inline void read(int &x)
{
	char c=nc(),b=1;
	for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
	for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

int n,m;

int u[1000005],v[1000005],w[1000005],ivst[1000005];

struct Splay{
	struct node{
		int size,idx,rev;
		int val,maximum;
		node *p,*ch[2],*fat,*mpos;
		inline void setc(node *c,int d) { ch[d]=c; c->p=this; }
		inline bool dir() { return p->ch[1]==this; }
		inline void update() {
			size=ch[0]->size+ch[1]->size+1;
			maximum=val; mpos=this;
			if (ch[0]->maximum>maximum)
				maximum=ch[0]->maximum,mpos=ch[0]->mpos;
			if (ch[1]->maximum>maximum)
				maximum=ch[1]->maximum,mpos=ch[1]->mpos;
		}
		inline void reverse() { rev^=1; swap(ch[0],ch[1]); }
		inline void pushdown(node *null){
			if (rev){
				if (ch[0]!=null) ch[0]->reverse();
				if (ch[1]!=null) ch[1]->reverse();
				rev=0;
			}
		}
	}*null,Mem[1200005];
	Splay() { null=Mem; null->p=null->ch[0]=null->ch[1]=null->fat=null; null->size=0; null->maximum=-(1<<30)-(1<<29); null->mpos=null; }
	inline void rot(node *x){
		if (x==null || x->p==null) return ;
		bool d=x->dir(); node *p=x->p;
		if (p->p!=null) p->p->setc(x,p->dir()); else x->p=null;
		p->setc(x->ch[d^1],d); x->setc(p,d^1); p->update(); x->update(); swap(x->fat,p->fat);
	}
	node *sta[1200005];
	inline void splay(node *x){
		int pnt=0; node *y=x;
		while (y!=null) sta[++pnt]=y,y=y->p;
		for (int i=pnt;i;i--) sta[i]->pushdown(null);
		while (x->p!=null)
			if (x->p->p==null)
				rot(x);
			else
				x->dir()==x->p->dir()?(rot(x->p),rot(x)):(rot(x),rot(x));
	}
	inline node *Access(node *x){
		node *y=null;
		while (x!=null)
		{
			splay(x);
			x->ch[1]->p=null; x->ch[1]->fat=x;
			x->setc(y,1); y->fat=null;
			x->update();
			y=x; x=x->fat;
		}
		return y;
	}
	inline void Link(node *x,node *y){
		if (Jud(x,y)) return;
		Access(x)->reverse();
		splay(x);
		x->fat=y;
		Access(x);
	}
	inline void Cut(node *x){
		Access(x);
		splay(x);
		x->ch[0]->p=null; x->ch[0]=null;
		x->fat=null; x->update();
	}
	inline void Cut(node *x,node *y){
		Access(x)->reverse();
		Cut(y);
	}
	inline node *Root(node *x){
		Access(x);
		splay(x);
		node *y=x;
		while (y->ch[0]!=null) y->pushdown(null),y=y->ch[0];
		return y;
	}
	inline bool Jud(node *x,node *y){
		return Root(x)==Root(y);
	}
	inline int Query(node *x,node *y){
		Access(x)->reverse();
		return Access(y)->maximum;
	}
	inline node *Road(node *x,node *y){
		Access(x)->reverse();
		return Access(y)->mpos;
	}
}LCT;

Splay::node *pos[1200005];

inline void Init(){
	for (int i=1;i<=n+m;i++){
		pos[i]=LCT.Mem+i; 
		pos[i]->p=pos[i]->ch[0]=pos[i]->ch[1]=pos[i]->fat=LCT.null; 
		pos[i]->val=-1<<30; pos[i]->maximum=-1<<30; pos[i]->mpos=pos[i];
		pos[i]->idx=i; pos[i]->size=1;
	}
}

struct event{
	int f,u,v;
	int idx;
	bool operator < (const event &B) const{
		return w[idx]<w[B.idx];
	}
}eve[1200005];
int tot,icnt,pnt;

int ans[100005];
map<abcd,int> Map;

int fat[100005];
inline int getfat(int u){
	return u==fat[u]?u:fat[u]=getfat(fat[u]);
}

inline void Kru()
{
	int fx,fy,itot=0;
	for (int i=1;i<=n;i++) fat[i]=i;
	for (int i=1;i<=m;i++)
		if (!ivst[i])
			eve[++tot].f=2,eve[tot].u=u[i],eve[tot].v=v[i],eve[tot].idx=i;
	sort(eve+pnt+1,eve+tot+1);
	for (int i=pnt+1;i<=tot;i++)
	{
		fx=getfat(eve[i].u); fy=getfat(eve[i].v);
		if (fx!=fy)
		{
			fat[fx]=fy;
			LCT.Link(pos[eve[i].u],pos[eve[i].idx+n]);
			LCT.Link(pos[eve[i].v],pos[eve[i].idx+n]);
			itot++;
			if (itot==n-1) break;
		}
	}
}

int main()
{
	int iu,iv,Q;
	freopen("t.in","r",stdin);
	freopen("t.out","w",stdout);
	read(n); read(m); read(Q);
	Init();
	for (int i=1;i<=m;i++)
	{
		read(u[i]),read(v[i]),read(w[i]);
		pos[i+n]->val=w[i]; pos[i+n]->update();
		if (u[i]>v[i]) swap(u[i],v[i]);
		Map.insert(make_pair(abcd(u[i],v[i]),i));
	}
	for (int i=1;i<=Q;i++)
	{
		read(eve[++tot].f); read(eve[tot].u); read(eve[tot].v);
		if (eve[tot].u>eve[tot].v) swap(eve[tot].u,eve[tot].v);
		if (eve[tot].f==2)
			ivst[eve[tot].idx=Map[abcd(eve[tot].u,eve[tot].v)]]=1;
		else if (eve[tot].f==1)
			eve[tot].idx=++icnt;
	}
	pnt=tot;
	Kru();
	for (int i=pnt;i;i--)
		if (eve[i].f==1)
		{
			ans[eve[i].idx]=LCT.Query(pos[eve[i].u],pos[eve[i].v]);
		}
		else
		{
			iu=eve[i].u; iv=eve[i].v;
			Splay::node *road=LCT.Road(pos[iu],pos[iv]);
			if (road->val>w[eve[i].idx])
			{
				LCT.Cut(road,pos[u[road->idx-n]]);
				LCT.Cut(road,pos[v[road->idx-n]]);
				LCT.Link(pos[n+eve[i].idx],pos[iu]);
				LCT.Link(pos[n+eve[i].idx],pos[iv]);
			}
		}
	for (int i=1;i<=icnt;i++)
		printf("%d\n",ans[i]);
	return 0;
}


相關推薦

[動態 LCT] BZOJ 2594 [Wc2006]水管局長資料加強

動態樹模板題 最開始加邊用Kruskal比較快,不用LCT,然後再用LCT維護最小生成樹 話說這道題小資料範圍的正解是什麼.... PS:以前LCT的求根都是打錯的,沒有pushdown() #include<cstdio> #include<cstd

BZOJ 2594: [Wc2006]水管局長資料加強(lct)

傳送門 解題思路   一道\(lct\)維護動態最小生成樹。剛開始寫了一遍瘋狂\(Re\),冷靜了一下重新寫了一遍終於過了。首先題目中要求兩點之間最大值的最小值,其實就是維護一個最小生成樹,每次詢問最大值。要將刪邊轉化成加邊操作,就是倒著處理,這裡用\(set\)和\(map\)就比較方便。然後還要把邊權

BZOJ 2594: [Wc2006]水管局長資料加強(LCT+最小生成樹+離線)

題目戳我 Solution 題目大意:就是讓你維護一棵動態的最小生成樹, 並詢問兩點路徑中邊權最大值。 很明顯,用LCT來做這個。有幾個關鍵點: ①刪邊維護mst不好搞,我們離線然後刪邊變加邊。 ②一開始用kruskal搞出底圖的mst,然後加

BZOJ 2594 [Wc2006]水管局長資料加強 LCT

題意:連結 方法: LCT 解析: 搞了一個上午加1個小時的題,TM最後大錯誤居然是排序元素太多排不回原來的樣子! 我要重新學排序! 這題是用LCT維護動態最小生成樹,但是最小生成樹上刪邊應該是做不到的,所以我們可以離線操作,之後先把所

BZOJ 2594 [Wc2006]水管局長資料加強

由於題目有刪除操作,所以我們不妨倒敘列舉,把刪邊變成加邊。某人還起了個好聽的名字叫“時光倒流”…… 之後用動態樹維護一下動態最小生成樹就行了,每加入一條邊就用動態樹樹找出兩點見最長邊,比較大小看是否能替換就行了。 這裡還有一個小技巧,就是二分邊標號來找出每個

bzoj2594 [Wc2006]水管局長資料加強lct+kruskal+離線)

考慮離線操作,倒著做,就是維護不斷往裡加邊的最小生成樹。具體做法與魔法森林差不多。 至於如何找到一條邊的標號,可以二分查詢來做。 #include <cstdio> #include <cstring> #include <al

[BZOJ2594][Wc2006]水管局長資料加強(kruskal+lct

題目描述 傳送門 題解 寫lct就應該有那種誓死不看板子的氣魄。 這道題思路還是很清晰的,維護一棵最小生成樹,每一次找樹鏈上權值最大的邊 刪邊變成倒序加邊 最開始的時候用沒有刪的邊kruskal直接最小生成樹 動態的話就是維護一棵lct,每一次加邊

BZOJ2594 [Wc2006]水管局長資料加強LCT

Description SC省MY市有著龐大的地下水管網路,嘟嘟是MY市的水管局長(就是管水管的啦),嘟嘟作為水管局長的工作就是:每天供水公司可能要將一定量的水從x處送往y處,嘟嘟需要為供水公司找到一條從A至B的水管的路徑,接著通過資訊化的控制中心通知路徑上的

bzoj2594: [Wc2006]水管局長資料加強

題目大意:給定一個簡單圖,支援刪邊,每次詢問兩點間 最大邊權值最小的路徑。 思路:題目中權值在邊上,而LCT維護的全值在點上,所以要先轉化,對於連線x和y的第i條路徑,把邊變成第i+n個點,然後把權值放到邊點上,所有真正的點權值賦為0,這也是維護邊權的常用方法。接著就是一

【bzoj2594】[Wc2006]水管局長資料加強 link cut tree

好久的坑了,今天終於填上了。 一直想總結一下LCT,結果發現難題都不會做。 離線處理,因為刪邊比較難做,所以我們倒著做變成加邊操作。 問題轉化成加邊維護最小生成樹,用lct維護一下最大的邊是哪一條,一旦新加入的邊形成環了,那麼看一看環上最大的邊和新加入的邊哪個大,如果新加入

BZOJ 2594: [Wc2006]水管局長數據加強(lct)

-- for 維護 int void rotate roo ESS sin 傳送門 解題思路   一道\(lct\)維護動態最小生成樹。剛開始寫了一遍瘋狂\(Re\),冷靜了一下重新寫了一遍終於過了。首先題目中要求兩點之間最大值的最小值,其實就是維護一個最小生成樹,每次詢問

[BZOJ2594][Wc2006][LCT]水管局長資料加強

題意 給一張圖,每次操作或詢問所有u到v路徑上邊權最大值的最小值,或刪除一條邊 離線,反過來操作,刪邊變成加邊,用LCT維護圖的聯通就行。 #include <cstdio> #include <iostream> #inc

【BZOJ2594】【WC2006水管局長資料加強

【題目連結】【思路要點】若題目僅包含詢問操作,顯然我們只需保留一棵圖的最小生成樹,在樹上詢問兩點間邊權的最大值即可。而本題中多了一種刪邊操作,但沒有強制線上,因此,我們可以將操作離線,轉化為加邊操作。用LinkCutTree維護動態最小生成樹即可。時間複雜度\(O((M+Q)

洛谷.4172.[WC2006]水管局長(LCT Kruskal)

處理 pushd 重新 需要 樹邊 IT || 路徑 max 題目鏈接 洛谷(COGS上也有) 不想去做加強版了。。(其實處理一下矩陣就好了) 題意: 有一張圖,求一條x->y的路徑,使得路徑上最長邊盡量短並輸出它的長度。會有<=5000次刪邊。 這實際上就是動

BZOJ2594 [Wc2006]水管局長數據加強LCT

最短 一次 read amp std 之前 應輸入 NPU 進入 Description SC省MY市有著龐大的地下水管網絡,嘟嘟是MY市的水管局長(就是管水管的啦),嘟嘟作為水管局長的工作就是:每天供水公司可能要將一定量的水從x處送往y處,嘟嘟需要為供水公司找到一條從A至

Luogu4712 WC2006 水管局長 LCT

傳送門   老套路,刪邊換成加邊 然後就變成了$LCT$維護最小生成樹的裸題 化邊為點,對於每一個點記錄鏈上最大值和對應的邊的編號,每一次加入一條邊時考慮是否能通過割掉當前樹上路徑上的最大權值的邊獲得一個更小的生成樹。 1 #include<bits/stdc++.h&

洛谷 P4172 [WC2006]水管局長 LCT維護最小生成樹

題目描述 SC 省 MY 市有著龐大的地下水管網路,嘟嘟是 MY 市的水管局長(就是管水管的啦),嘟嘟作為水管局長的工作就是:每天供水公司可能要將一定量的水從 x 處送往 y 處,嘟嘟需要為供水公司找到一條從 A 至 B 的水管的路徑,接著通過資訊化的控制中心

bzoj2594: [Wc2006]水管局長數據加強

date get push roo 表示 sort 選擇 mit 管道 地址:http://www.lydsy.com/JudgeOnline/problem.php?id=2594 題目: 2594: [Wc2006]水管局長數據加強版 Time Limit:

BZOJ_2594_[Wc2006]水管局長數據加強_LCT

接受 無向圖 回車 des 操作 sample led rev 輸出 BZOJ_2594_[Wc2006]水管局長數據加強版_LCT Description SC省MY市有著龐大的地下水管網絡,嘟嘟是MY市的水管局長(就是管水管的啦),嘟嘟作為水管局長的工作就是:每

luogu4172 [WC2006]水管局長

als CI too cmp ID oid for air continue 就是用 lct 維護最小生成樹 ref #include <algorithm> #include <iostream> #include <cstdio> #