BZOJ2333 [SCOI2011]棘手的操作 【离线 + 线段树】

题目

有N个节点,标号从1到N,这N个节点一开始相互不连通。第i个节点的初始权值为a[i],接下来有如下一些操作:
U x y: 加一条边,连接第x个节点和第y个节点
A1 x v: 将第x个节点的权值增加v
A2 x v: 将第x个节点所在的连通块的所有节点的权值都增加v
A3 v: 将所有节点的权值都增加v
F1 x: 输出第x个节点当前的权值
F2 x: 输出第x个节点所在的连通块中,权值最大的节点的权值
F3: 输出所有节点中,权值最大的节点的权值

输入格式

输入的第一行是一个整数N,代表节点个数。
接下来一行输入N个整数,a[1], a[2], …, a[N],代表N个节点的初始权值。
再下一行输入一个整数Q,代表接下来的操作数。
最后输入Q行,每行的格式如题目描述所示。

输出格式

对于操作F1, F2, F3,输出对应的结果,每个结果占一行。

输入样例

3

0 0 0

8

A1 3 -20

A1 2 20

U 1 3

A2 1 10

F1 3

F2 3

A3 -10

F3

输出样例

-10

10

10

提示

对于30%的数据,保证 N<=100,Q<=10000

对于80%的数据,保证 N<=100000,Q<=100000

对于100%的数据,保证 N<=300000,Q<=300000

对于所有的数据,保证输入合法,并且 -1000<=v, a[1], a[2], …, a[N]<=1000

题解

据说此题很多人堆套堆,怎么这么难写
我那么弱当然是用线段树啦

我觉得线段树的确好写到不知哪里去
对于所有操作,似乎在线段树上都很好实现,唯一的难点就在于点的编号

那么问题就转化成了,给定一种编号方法,使任意时刻同一个联通块内的所有点编号连续
只需要分两种情况想就很容易实现了:

我们想象,一开始所有点相互独立,没什么关系

①当两个独立的点相连时,它们的编号一定是连续的,否则此时就不满足所需性质
那我们就先用链表将它们连起来,表示编号连续

②当两个联通块相连时,由我们维护的性质得:两个联通块内部的点编号一定是连续的,现在我们需要两个联通块编号连续,我们只需要将它们的编号衔接起来就好了,那么我们把其中一个联通块所对应的链 接到另一个联通块对应的链末尾就好了

可以发现,这样子操作之后,我们就会得出若干个链,表示链上的点编号必须按链上的顺序
所以我们按链的顺序标号,就能保证所有时刻联通块内部点的编号连续

取链头链尾用并查集实现
我们在询问的时候,也要用上并查集,并且链接顺序与标号的时候相同,保证每个联通块目前的代表元一定是标号时编号最小的点,所以我们再维护并查集的大小就可以轻松求出每次操作的区间啦~

数据结构部分就只用实现一个简单的线段树
比堆套堆不知道要好写到哪里去

丑丑的代码

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long int
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<' '; puts("");
#define ls (u << 1)
#define rs (u << 1 | 1)
using namespace std;
const int maxn = 300005,maxm = 100005,INF = 1000000000;
inline int read(){
	int out = 0,flag = 1; char c = getchar();
	while (c < 48 || c > 57) {if (c == '-') flag = -1; c = getchar();}
	while (c >= 48 && c <= 57) {out = (out << 3) + (out << 1) + c - '0'; c = getchar();}
	return out * flag;
}
struct Query{int opt,a,b,c;}q[maxn];
int n,m,pre[maxn],post[maxn],id[maxn],Hash[maxn],val[maxn],cnt;
int nxt[maxn],siz[maxn];
int mx[4 * maxn],tag[4 * maxn];
char opt[10];
int findu(int u){return u == pre[u] ? u : pre[u] = findu(pre[u]);}
int findd(int u){return u == post[u] ? u : post[u] = findd(post[u]);}
void build(int u,int l,int r){
	if (l == r) {mx[u] = val[Hash[l]]; return;}
	int mid = l + r >> 1;
	build(ls,l,mid);
	build(rs,mid + 1,r);
	mx[u] = max(mx[ls],mx[rs]);
}
void pd(int u){
	if (tag[u]){
		mx[ls] += tag[u]; tag[ls] += tag[u];
		mx[rs] += tag[u]; tag[rs] += tag[u];
		tag[u] = 0;
	}
}
void modify(int u,int l,int r,int L,int R,int v){
	if (l >= L && r <= R){mx[u] += v; tag[u] += v; return;}
	pd(u);
	int mid = l + r >> 1;
	if (mid >= L) modify(ls,l,mid,L,R,v);
	if (mid < R) modify(rs,mid + 1,r,L,R,v);
	mx[u] = max(mx[ls],mx[rs]);
}
int query(int u,int l,int r,int L,int R){
	if (l >= L && r <= R) return mx[u];
	pd(u);
	int mid = l + r >> 1;
	if (mid >= R) return query(ls,l,mid,L,R);
	else if (mid < L) return query(rs,mid + 1,r,L,R);
	else return max(query(ls,l,mid,L,R),query(rs,mid + 1,r,L,R));
}
int main(){
	n = read();
	for (int i = 1; i <= n; i++) val[i] = read(),pre[i] = post[i] = i;
	m = read();
	int fa,fb,sa,sb;
	for (int i = 1; i <= m; i++){
		scanf("%s",opt);
		if (opt[0] == 'U'){
			q[i].opt = 0,q[i].a = read(),q[i].b = read();
			fa = findu(q[i].a); fb = findu(q[i].b);
			if (fa == fb) continue;
			sa = findd(q[i].a); sb = findd(q[i].b);
			nxt[sa] = fb;
			pre[fb] = fa;
			post[sa] = sb;
		}
		else if (opt[0] == 'A'){
			q[i].a = read();
			if (opt[1] == '1') q[i].opt = 1,q[i].b = read();
			else if (opt[1] == '2') q[i].opt = 2,q[i].b = read();
			else q[i].opt = 3;
		}else {
			if (opt[1] == '1') q[i].opt = 4,q[i].a = read();
			else if (opt[1] == '2') q[i].opt = 5,q[i].a = read();
			else q[i].opt = 6;
		}
	}
	for (int i = 1; i <= n; i++){
		if (id[i]) continue;
		int u = findu(i);
		while (u) id[u] = ++cnt,Hash[cnt] = u,u = nxt[u];
	}
	build(1,1,n);
	for (int i = 1; i <= n; i++) pre[i] = i,siz[i] = 1;
	for (int i = 1; i <= m; i++){
		switch(q[i].opt){
			case 0:
				fa = findu(q[i].a); fb = findu(q[i].b);
				if (fa != fb){
					siz[fa] += siz[fb];
					pre[fb] = fa;
				}
				break;
			case 1:
				modify(1,1,n,id[q[i].a],id[q[i].a],q[i].b);
				break;
			case 2:
				fa = findu(q[i].a);
				modify(1,1,n,id[fa],id[fa] + siz[fa] - 1,q[i].b);
				break;
			case 3:
				modify(1,1,n,1,n,q[i].a);
				break;
			case 4:
				printf("%d
",query(1,1,n,id[q[i].a],id[q[i].a]));
				break;
			case 5:
				fa = findu(q[i].a);
				printf("%d
",query(1,1,n,id[fa],id[fa] + siz[fa] - 1));
				break;
			case 6:
				printf("%d
",mx[1]);
				break;
		}
	}
	return 0;
}

原文地址:https://www.cnblogs.com/Mychael/p/8490297.html