【JZOJ 杂题选讲】【CF1270H】Number of Components

题目

对于一个序列a[i]，如果i<j且a[i]<a[j]，那么i和j有一条连边。

·有q次询问，每一次询问会修改一个a[i]的值，问每次询问之后连通块的个数。

·n,q<=5e5,ai<=1e6,

思路

基 础 线 段 树 题
首先我们可以发现，联通块一定是一段区间
我们只需找出有多少分界点。
如果一个点i是分界点，那么这个点i满足所有大于ai的在左边，所有小于ai的在右边
也就说只要呈现11111110000000的形势就有一个分界点
接下来只需要对所有的a[i]，维护它的f，求所有a[i]，f为1的个数。
问题在于修改，由于修改a[j]，只与a[j-1],a[j+1]以及对应的a[i]有关，对于a[i]的影响是区间性的，用权值线段树维护，区间修改，分类与a[j-1],a[j+1]的贡献就好了。
区间加减只需要记录最小值以及最小值的个数就好了。
注意到自己改自己不好搞，我们可以离线把询问中的q个a的值也当做线段树的一部分，同时维护这些还未生效的a的f

代码

#include<bits/stdc++.h>
#define ls(p) p << 1
#define rs(p) p << 1 | 1
const int N = 1e6 + 77;
int n,q,A[N],pos[N],val[N],lim;
int min[N * 4],cnt[N * 4],tag[N * 4];
void pushup(int p) 
{
	if (min[ls(p)] == min[rs(p)]) 
	{
		min[p] = min[ls(p)];
		cnt[p] = cnt[ls(p)] + cnt[rs(p)];
	} else if (min[ls(p)] < min[rs(p)]) {
		min[p] = min[ls(p)];
		cnt[p] = cnt[ls(p)];
	} else {
		min[p] = min[rs(p)];
		cnt[p] = cnt[rs(p)];
	}
}

void push(int p,int v) {
	min[p] += v;
	tag[p] += v;
}

void pushdown(int p) {
	if (tag[p]) {
		push(ls(p),tag[p]);
		push(rs(p),tag[p]);
		tag[p] = 0;
	}
}

void add(int p,int x,int v,int L,int R) {
	if (L == R) {
		cnt[p] += v;
		return;
	}
	pushdown(p);
	int mid = (L + R) >> 1;
	if (x <= mid)
		add(ls(p),x,v,L,mid);
	else
		add(rs(p),x,v,mid + 1,R);
	pushup(p);
}

void add(int p,int l,int r,int v,int L,int R) {
	if (l <= L && r >= R)
		return push(p,v);
	pushdown(p);
	int mid = (L + R) >> 1;
	if (l <= mid)
		add(ls(p),l,r,v,L,mid);
	if (r > mid)
		add(rs(p),l,r,v,mid + 1,R);
	pushup(p);
}

int query(int p,int l,int r,int L,int R) {
	if (l <= L && r >= R) {
		if (min[p] == 1)
			return cnt[p];
		return 0;
	}
	pushdown(p);
	int mid = (L + R) >> 1,ret = 0;
	if (l <= mid)
		ret = query(ls(p),l,r,L,mid);
	if (r > mid)
		ret += query(rs(p),l,r,mid + 1,R);
	return ret;
}

void update(int x,int y,int v) {
	if (x == y) return;
	if (x > y) std::swap(x,y);
	add(1,x,y - 1,v,0,lim);
}

int main() {
	scanf("%d%d",&n,&q);
	for(int i = 1; i <= n; ++i) {
		scanf("%d",&A[i]);
		lim = std::max(lim,A[i]);
	}
	for(int i = 1; i <= q; ++i) {
		scanf("%d%d",&pos[i],&val[i]);
		lim = std::max(lim,val[i]);
	}
	A[0] = ++lim;
	for(int i = 0; i <= n; ++i) {
		update(A[i],A[i + 1],1);
		add(1,A[i],1,0,lim);
	}
	for(int i = 1; i <= q; ++i) {
		int p = pos[i],v = val[i];
		add(1,A[p],-1,0,lim);
		update(A[p - 1],A[p],-1);
		update(A[p],A[p + 1],-1);
		add(1,v,1,0,lim);
		update(A[p - 1],v,1);
		update(v,A[p + 1],1);
		A[p] = v;
		printf("%d\n",query(1,1,lim - 1,0,lim));
	}
	return 0;
}