luogu P4172 [WC2006]水管局长

背景:

人生考试最低排名: R a n k = 91 Rank=91 Rank=91(年级 800 + 800+ 800+人)。
数学人生最低分数: 88 88 88(满分 120 120 120)。
靠语物化撑场。
纪念一下。

题目传送门:

https://www.luogu.org/problemnew/show/P4172

题意:

维护一棵支持删边的最小生成树。

思路:

显然最小生成树的做法无法支持删边。
不妨倒序来做,就变成了 L C T LCT LCT维护最小生成树的模板了。

代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#include<cstdlib>
using namespace std;
map<int,map<int,int> > MAP1,MAP2;
	struct node1{int d,fa,ma,pos,lazy,son[2];} tr[500010];
	struct node2{int x,y,z;} a[500010];
	struct node3{int t,x,y;} b[500010];
	int ans[500010];
	int n,m,q;
bool isroot(int x)
{
	return tr[tr[x].fa].son[0]!=x&&tr[tr[x].fa].son[1]!=x;
}
bool isson(int x)
{
	return x==tr[tr[x].fa].son[1];
}
void change(int x)
{
	if(!x)return;
	swap(tr[x].son[0],tr[x].son[1]);
	tr[x].lazy^=1;
}
void pushup(int x)
{
	int lc=tr[x].son[0],rc=tr[x].son[1];
	tr[x].ma=tr[x].d,tr[x].pos=x;
	if(tr[lc].ma>tr[x].ma) tr[x].ma=tr[lc].ma,tr[x].pos=tr[lc].pos;
	if(tr[rc].ma>tr[x].ma) tr[x].ma=tr[rc].ma,tr[x].pos=tr[rc].pos;
}
void pushdown(int x)
{
	if(!tr[x].lazy) return;
	change(tr[x].son[0]),change(tr[x].son[1]);
	tr[x].lazy=0;
}
void rot(int x)
{
	int w=isson(x),y=tr[x].fa,yy=tr[y].fa;
	tr[y].son[w]=tr[x].son[w^1];
	if(tr[y].son[w]) tr[tr[y].son[w]].fa=y;
	tr[x].fa=yy;
	if(!isroot(y)) tr[yy].son[isson(y)]=x;
	tr[x].son[w^1]=y;tr[y].fa=x;
	pushup(y);
}
int sta[300010];
int top;
void splay(int x)
{
	sta[top=1]=x;
	for(int i=x;!isroot(i);i=tr[i].fa)
		sta[++top]=tr[i].fa;
	while(top) pushdown(sta[top--]);
	for(int y=tr[x].fa;!isroot(x);rot(x),y=tr[x].fa)
		if(!isroot(y)) isson(x)^isson(y)?rot(x):rot(y);
	pushup(x);
}
void access(int x)
{
	for(int y=0;x;y=x,x=tr[x].fa)
		splay(x),tr[x].son[1]=y,pushup(x);
}
int findroot(int x)
{
	access(x);splay(x);
	while(tr[x].son[0]) x=tr[x].son[0];
	//splay(x);
	return x; 
}
void makeroot(int x)
{
	access(x);splay(x);change(x);
}
void link(int x,int y)
{
	makeroot(x);
	if(findroot(y)!=x) tr[x].fa=y;
}
void split(int x,int y)
{
	makeroot(x);access(y);splay(y);
}
void del(int x,int y)
{
	split(x,y);
	tr[x].fa=tr[y].son[0]=0,pushup(x),pushup(y);
}
node1 query(int x,int y)
{
	split(x,y);
	return tr[y];
}
void work(int id)
{
	if(findroot(a[id].x)!=findroot(a[id].y))
	{
		tr[id].d=a[id].z;
		link(a[id].x,id);link(a[id].y,id);
	}
	else
	{
		node1 MAX=query(a[id].x,a[id].y);
		if(MAX.pos>n&&a[id].z<MAX.ma)
		{
			del(a[MAX.pos].x,MAX.pos);del(a[MAX.pos].y,MAX.pos);
			tr[id].d=a[id].z;
			link(a[id].x,id);link(a[id].y,id);
		}
	}
}
int main()
{
	scanf("%d %d %d",&n,&m,&q);
	for(int i=0;i<=(n<<1);i++)
		tr[i]=(node1){0,0,0,i,0,0,0};
	for(int i=n+1;i<=n+m;i++)
	{
		scanf("%d %d %d",&a[i].x,&a[i].y,&a[i].z);
		MAP1[a[i].x][a[i].y]=MAP1[a[i].y][a[i].x]=i;
	}
	for(int i=1;i<=q;i++)
	{
		scanf("%d %d %d",&b[i].t,&b[i].x,&b[i].y);
		if(b[i].t==2) MAP2[b[i].x][b[i].y]=MAP2[b[i].y][b[i].x]=1;
	}
	for(int i=n+1;i<=n+m;i++)
		if(!MAP2[a[i].x][a[i].y]) work(i);
	int t=0;
	for(int i=q;i>=1;i--)
		if(b[i].t==1) ans[++t]=query(b[i].x,b[i].y).ma; else work(MAP1[b[i].x][b[i].y]);
	for(int i=t;i>=1;i--)
		printf("%d\n",ans[i]);
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值