「Luogu P3273 [SCOI2011]棘手的操作」

题目大意

给出一些点和一些操作,资瓷链接两点,某个点所在的联通块都加上一个数,某个点加上一个数,全部加上一个数,查询单点的值,查询某个点所在的连通块中的最大值,查询所有值中的最大值.

分析

可并堆裸题.因为连边这个东西很难处理,所以考虑离线,考虑对每个点记录一下 (tim_i) 表示当前节点连接其他节点的时间.再考虑建树,但是如果直接在需要连接的两点直接连边会出问题,没有办法处理某个时间的连通块,所以考虑用并查集维护一下,这样就可以保证对于某个节点的子节点连接父亲的时间都是在自己连接父亲之前,考虑当前修改的连通块在 DFS 序中必定是连续的一段,且左端点必定是当前节点的某个祖先,那么寻找祖先就可以用倍增来维护,在知道祖先后考虑修改的部分必定包含当前节点的若干子树,且在儿子中链接当前节点的时间是单调递增的,所以可以用二分来快速查询,然后就是一个线段树模板维护区间 (max) 就好了.

代码

#include<bits/stdc++.h>
#define REP(i,first,last) for(int i=first;i<=last;++i)
#define DOW(i,first,last) for(int i=first;i>=last;--i)
using namespace std;
const int MAXN=3e5+5;
const int INF=1e9;
int n,m;
vector<int>edge[MAXN];
char opt[MAXN][2];
int x[MAXN],y[MAXN],v[MAXN];
int father[MAXN][21];
int Find(int now)//并查集部分
{
	if(now==father[now][0])
	{
		return now;
	}
	return father[now][0]=Find(father[now][0]);
}
int tim[MAXN]={MAXN};//记录链接父亲的时间
int val[MAXN];
int id[MAXN];
int dfs[MAXN];
int dfs_cnt=0;
int siz[MAXN];//记录子树大小
void DFS(int now)//遍历所有树
{
	REP(i,1,19)//倍增预处理
	{
		father[now][i]=father[father[now][i-1]][i-1];
	}
	id[now]=++dfs_cnt;//当前节点在dfs序的位置
	dfs[dfs_cnt]=now;
	siz[now]=1;
	if(edge[now].size())
	{
		REP(i,0,edge[now].size()-1)//遍历子节点
		{
			father[edge[now][i]][0]=now;
			DFS(edge[now][i]);
			siz[now]+=siz[edge[now][i]];
		}
	}
}
int l,r;
void Getlr(int now,int now_time)//得到对于某个时间时某个节点的联通快
{
	if(tim[now]<now_time)//如果当前节点在当前时间已经连接了父亲,那么就倍增找祖先
	{
		DOW(i,19,0)
		{
			if(tim[father[now][i]]<now_time)
			{
				now=father[now][i];
			}
		}
		now=father[now][0];
	}
	int left=0,right=edge[now].size()-1,middle;
	l=id[now];
	r=id[now];
	while(left<=right)//在自己的儿子内二分
	{
		middle=(left+right)>>1;
		if(tim[edge[now][middle]]<now_time)
		{
			r=id[edge[now][middle]];
			left=middle+1;
		}
		else
		{
			right=middle-1;
		}
	}
	if(r^id[now])//如果自己的子树有修改了那么这整颗子树都会被修改
	{
		r+=siz[dfs[r]]-1;
	}
}
//以下为一颗区间加,区间查询max的线段树,过于模板就不多讲了
struct LazyTag
{
	int add;
	void Clean()
	{
		add=0;
	}
}for_make;
LazyTag MakeTag(int add)
{
	for_make.add=add;
	return for_make;
}
struct SegmentTree
{
	long long sum;
	long long max;
	LazyTag tag;
}sgt[MAXN*4];
#define LSON (now<<1)
#define RSON (now<<1|1)
#define MIDDLE ((left+right)>>1)
#define LEFT LSON,left,MIDDLE
#define RIGHT RSON,MIDDLE+1,right
#define NOW now_left,now_right
void PushUp(int now)
{
	sgt[now].sum=sgt[LSON].sum+sgt[RSON].sum;
	sgt[now].max=max(sgt[LSON].max,sgt[RSON].max);
}
void Build(int now=1,int left=1,int right=n)
{
	if(left==right)
	{
		sgt[now].max=sgt[now].sum=val[dfs[left]];
		return;
	}
	Build(LEFT);
	Build(RIGHT);
	PushUp(now);
}
void Down(LazyTag tag,int now,int left,int right)
{
	sgt[now].sum+=1ll*tag.add*(right-left+1);
	sgt[now].max+=tag.add;
	sgt[now].tag.add+=tag.add;
}
void PushDown(int now,int left,int right)
{
	Down(sgt[now].tag,LEFT);
	Down(sgt[now].tag,RIGHT);
	sgt[now].tag.Clean();
}
void Updata(int now_left,int now_right,int add,int now=1,int left=1,int right=n)
{
	if(now_right<left||right<now_left)
	{
		return;
	}
	if(now_left<=left&&right<=now_right)
	{
		Down(MakeTag(add),now,left,right);
		return;
	}
	PushDown(now,left,right);
	Updata(NOW,add,LEFT);
	Updata(NOW,add,RIGHT);
	PushUp(now);
}
int QueryMax(int now_left,int now_right,int now=1,int left=1,int right=n)
{
	if(now_right<left||right<now_left)
	{
		return -INF;
	}
	if(now_left<=left&&right<=now_right)
	{
		return sgt[now].max;
	}
	PushDown(now,left,right);
	return max(QueryMax(NOW,LEFT),QueryMax(NOW,RIGHT));
}
#undef LSON
#undef RSON
#undef MIDDLE
#undef LEFT
#undef RIGHT
#undef NOW
int main()
{
	scanf("%d",&n);
	REP(i,1,n)
	{
		scanf("%d",&val[i]);
		father[i][0]=i;
	}
	scanf("%d",&m);
	REP(i,1,m)//离线
	{
		cin>>opt[i];
		if(opt[i][0]=='U')//对于连边操作连边
		{
			scanf("%d%d",&x[i],&y[i]);
			if(Find(y[i])^Find(x[i]))
			{
				tim[Find(y[i])]=i;
				edge[Find(x[i])].push_back(Find(y[i]));
				father[Find(y[i])][0]=Find(x[i]);
			}
		}
		if(opt[i][0]=='A')
		{
			if(opt[i][1]=='1')
			{
				scanf("%d%d",&x[i],&v[i]);
			}
			if(opt[i][1]=='2')
			{
				scanf("%d%d",&x[i],&v[i]);
			}
			if(opt[i][1]=='3')
			{
				scanf("%d",&v[i]);
			}
		}
		if(opt[i][0]=='F')
		{
			if(opt[i][1]=='1')
			{
				scanf("%d",&x[i]);
			}
			if(opt[i][1]=='2')
			{
				scanf("%d",&x[i]);
			}
			if(opt[i][1]=='3')
			{

			}
		}
	}
	REP(i,1,n)//获得这个森林的dfs序
	{
		if(!tim[i])
		{
			tim[i]=INF;
			DFS(i);
		}
	}
	Build();//建树
	REP(i,1,m)
	{
		if(opt[i][0]=='A')//修改操作
		{
			if(opt[i][1]=='1')
			{
				Updata(id[x[i]],id[x[i]],v[i]);
			}
			if(opt[i][1]=='2')
			{
				Getlr(x[i],i);
				Updata(l,r,v[i]);
			}
			if(opt[i][1]=='3')
			{
				Updata(1,n,v[i]);
			}
		}
		if(opt[i][0]=='F')//查询操作
		{
			if(opt[i][1]=='1')
			{
				printf("%d
",QueryMax(id[x[i]],id[x[i]]));
			}
			if(opt[i][1]=='2')
			{
				Getlr(x[i],i);
				printf("%d
",QueryMax(l,r));
			}
			if(opt[i][1]=='3')
			{
				printf("%d
",QueryMax(1,n));
			}
		}
	}
	return 0;
}

原文地址:https://www.cnblogs.com/Sxy_Limit/p/12825806.html