CF1416D Graph and Queries

一、题目

点此看题

二、解法

删边有一个经典的套路就是变成加边离线,但是操作 1 1 1必须要按照时间顺序来,蛋疼。

这时候就要介绍一下 k r u s k a l \tt kruskal kruskal重构树(以前我还不知道有这个东西),在合并来个联通块的时候我们新建一个点连接这两个合并点的根,就建成了重构树。

离线的时候可以建出我们的重构树,因为重构树一个重要的性质就是新建点的子树代表了某个时刻这些点是连在一起的。询问的时候我们找到代表当时联通性的这个祖先,怎么找呢?我们把联通的时间放在这个祖先上,满足联通时间 ≥ \geq 当前时间就证明他已经联通了,倍增找到之后对子树进行询问和修改操作, d f n dfn dfn序加线段树就行了。

#include <cstdio>
#include <iostream>
using namespace std;
const int M = 500005;
#define pii pair<int,int>
int read()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,q,cnt,tot,f[M],fa[M],a[M],op[M],b[M],use[M],tim[M];
int Index,dfn[M],dl[M],dr[M],up[M][20];pii val[2*M];
struct node
{
	int u,v;
}s[M];
struct edge
{
	int v,next;
}e[2*M];
int find(int x)
{
	if(x!=fa[x]) fa[x]=find(fa[x]);
	return fa[x];
}
void add(int u,int v)
{
	//printf("nmsl %d %d\n",u,v);
	e[++tot]=edge{v,f[u]},f[u]=tot;
	e[++tot]=edge{u,f[v]},f[v]=tot;
}
void uni(int u,int v,int t)
{
	int x=find(u),y=find(v);
	if(x!=y)
	{
		cnt++;tim[cnt]=t;
		fa[x]=cnt;fa[y]=cnt;
		add(x,cnt);add(y,cnt);
	}
}
void dfs(int u,int p)
{
	dl[u]=1e9;
	if(a[u]) dfn[u]=dl[u]=dr[u]=++Index,tim[u]=tim[p];
	up[u][0]=p;
	for(int i=1;i<20;i++)
		up[u][i]=up[up[u][i-1]][i-1];
	for(int i=f[u];i;i=e[i].next)
	{
		int v=e[i].v;
		if(v==p) continue;
		dfs(v,u);
		dl[u]=min(dl[u],dl[v]);
		dr[u]=max(dr[u],dr[v]);
	}
}
void insert(int i,int l,int r,int id,int c)
{
	if(l==r)
	{
		val[i]=make_pair(c,id);
		return ;
	}
	int mid=(l+r)>>1;
	if(mid>=id) insert(i<<1,l,mid,id,c);
	else insert(i<<1|1,mid+1,r,id,c);
	val[i]=max(val[i<<1],val[i<<1|1]);
}
pii ask(int i,int l,int r,int L,int R)
{
	if(l>R || L>r) return make_pair(0,0);
	if(L<=l && r<=R) return val[i];
	int mid=(l+r)>>1;
	return max(ask(i<<1,l,mid,L,R),ask(i<<1|1,mid+1,r,L,R));
}
int jump(int u,int t)
{
	for(int i=19;i>=0;i--)
		if(up[u][i] && tim[up[u][i]]>=t)
			u=up[u][i];
	return u;
}
signed main()
{
	for(int i=1;i<M;i++) fa[i]=i;
	n=cnt=read();m=read();q=read();
	for(int i=1;i<=n;i++)
		a[i]=read();
	for(int i=1;i<=m;i++)
	{
		int u=read(),v=read();
		s[i]=node{u,v};
	}
	for(int i=1;i<=q;i++)
	{
		op[i]=read(),b[i]=read();
		if(op[i]==2) use[b[i]]=1;
	}
	for(int i=1;i<=m;i++)
		if(!use[i])
			uni(s[i].u,s[i].v,q);
	for(int i=q;i>=1;i--)
	{
		if(op[i]==2)
			uni(s[b[i]].u,s[b[i]].v,i);
	}
	for(int i=1;i<=cnt;i++)
		if(find(i)==i)
			dfs(i,0);
	for(int i=1;i<=n;i++)
	{
		//printf("%d\n",dfn[i]);
		insert(1,1,n,dfn[i],a[i]);
	}
	for(int i=1;i<=q;i++)
	{
		if(op[i]==2) continue;
		int x=jump(b[i],i);
		//printf("-> %d %d %d\n",b[i],dl[x],dr[x]);
		pii y=ask(1,1,n,dl[x],dr[x]);
		printf("%d\n",y.first);
		insert(1,1,n,y.second,0);
	}
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值