ACM-ICPC 2019南昌网络赛I题 Yukino With Subinterval

ACM-ICPC 2019南昌网络赛I题 Yukino With Subinterval

题目大意：给一个长度为n，值域为[1, n]的序列{a}，要求支持m次操作：

1. 单点修改
   1 pos val
2. 询问子区间中某个值域的数的个数，连续的相同数字只记为一个。（即统计数字段的个数）
   2 L R x y

数据范围：

1 ≤ n,m ≤ 2×10^5

1 ≤ a[i] ≤ n

---

解题思路：

• 连续重复的数字只记一次。所以考虑将每个数字段除第一个出现外的数字都删去（记为0）。在读入操作的时候暴力模拟，同时维护一个新序列，把原序列的连续重复数字在新序列都记为0。
• 询问在位置上和值域上都是可作差的。所以将原询问拆分为4个关于（t, pos, val）的子询问
• 将初值也转换为修改操作
• 转化成了一个带修改的二维偏序问题；算上时间序就是一个三维偏序

细节：

• 在修改a[i]时要分别考虑删除原数与添加新数对a[i - 1]和a[i + 1]的影响。删除和添加分开讨论会比较方便

if (pos != n && a[pos] == a[pos + 1])
	A[cnt++] = (ACT) {i, pos + 1, a[pos], 1, 1, 0};
if (pos == 1 || a[pos] != a[pos - 1])
	A[cnt++] = (ACT) {i, pos, a[pos], 1, -1, 0};
a[pos] = val;
if (pos != n && val == a[pos + 1])
	A[cnt++] = (ACT) {i, pos + 1, val, 1, -1, 0};
if (pos == 1 || val != a[pos - 1])
	A[cnt++] = (ACT) {i, pos, val, 1, 1, 0};

• 要特判每个询问的最左端点，如果是0，且会对答案有影响，（在值域范围内）要给答案+1

ans[Q] = (l == 1 || a[l] != a[l - 1]) ? 0 : 1;
if (a[l] < x || a[l] > y)	ans[Q] = 0; //无需修正

• 分治时，右半边的修改操作不需要再应用（其对右半边的查询已经被统计，对左半边的查询无影响，）
• 写分治的题时有几个不错的调试输出点，可以很快地找到错误：
  1. 树状数组的修改、查询函数开头处
  2. 查询操作的统计答案处
  3. 递归头部输出L、R，确定函数调用的状态

---

CDQ真的快。考场上队友用树套树，手写Treap + 读优还是被卡常，事后CDQ分治连cin都不用改就过了·_·！

---

AC代码：

#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std;

const int maxn = 2e5 + 1;
int cnt, ans[maxn], Q, a[maxn], n, sym[maxn << 3];
int _v[maxn];
struct ACT {
	int t, pos, val, opt, pls, id;
	bool operator<(ACT &rhs) const {
		return t < rhs.t || (t == rhs.t && \
			(pos < rhs.pos || (pos == rhs.pos && \
			(val < rhs.val || (val == rhs.val && \
			opt < rhs.opt)))));
	}
}A[maxn << 3], T[maxn << 3];

int lowbit(int x) {
	return x & (-x);
}
void add(int x, int v) {
	if (x == 0)
		return;
	//printf("add: %d, %d\n", x, v);
	while (x <= n) {
		_v[x] += v;
		x += lowbit(x);
	}
}
int query(int x) {
	int ret = 0;
	while (x) {
		ret += _v[x];
		x -= lowbit(x);
	}
	return ret;
}
void read() { //(A) {t, pos, val, opt, pls, id}
	int m, i, l, r, x, y, opt, pos, val;
	cin>>n>>m;
	for (i = 0; i < n; ++i) {
		cin>>a[i + 1];
		if (a[i + 1] != a[i])
			A[cnt++] = (ACT) {-1, i + 1, a[i + 1], 1, 1, 0};
	}
	for (i = 0; i < m; ++i) {
		cin>>opt;
		if (opt == 1) {
			cin>>pos>>val;
			if (pos != n && a[pos] == a[pos + 1])
				A[cnt++] = (ACT) {i, pos + 1, a[pos], 1, 1, 0};
			if (pos == 1 || a[pos] != a[pos - 1])
				A[cnt++] = (ACT) {i, pos, a[pos], 1, -1, 0};
			a[pos] = val;
			if (pos != n && val == a[pos + 1])
				A[cnt++] = (ACT) {i, pos + 1, val, 1, -1, 0};
			if (pos == 1 || val != a[pos - 1])
				A[cnt++] = (ACT) {i, pos, val, 1, 1, 0};
		} else {
			cin>>l>>r>>x>>y;
			ans[Q] = (l == 1 || a[l] != a[l - 1]) ? 0 : 1;
			if (a[l] < x || a[l] > y)
				ans[Q] = 0;
			A[cnt++] = (ACT) {i, r, y, 2, 1, Q};
			A[cnt++] = (ACT) {i, r, x - 1, 2, -1, Q};
			A[cnt++] = (ACT) {i, l - 1, y, 2, -1, Q};
			A[cnt++] = (ACT) {i, l - 1, x - 1, 2, 1, Q++};
		}
	}
	sort(A, A + cnt);
}
void mrdge(int l, int r, ACT* T) {
	if (r - l <= 1)	return;
	int mid = l + (r - l) / 2, p = l, q = mid, tot = 0;
	//printf("###l = %d, r = %d\n", l, r);
	mrdge(l, mid, T);	mrdge(mid, r , T);
	//printf("***l = %d, r = %d\n", l, r);
	if (query(n) != 0)
		cout<<"! BIT = "<<query(n)<<endl;
	while (p < mid || q < r) {
		if (q >= r || (p < mid && (A[p].pos < A[q].pos || (A[p].pos == A[q].pos && \
			(A[p].val < A[q].val || (A[p].val == A[q].val && \
			A[p].opt < A[q].opt)))))) {
			if (A[p].opt == 1) {
				add(A[p].val, A[p].pls);
				sym[tot] = 1;
			}
			T[tot++] = A[p++];
		}
		else {
			if (A[q].opt == 2) {
				ans[A[q].id] += query(A[q].val) * A[q].pls;
				//printf("querying[%d]...ans[%d] = %d\n", A[q].val, A[q].id, ans[A[q].id]);
			}
			T[tot++] = A[q++];
		}
	}
	for (int i = 0; i < tot; ++i) {
		if (T[i].opt == 1 && sym[i]) {
			add(T[i].val, -T[i].pls);
			sym[i] = 0;
		}
		A[l + i] = T[i];
	}
	return;
}
void print() {
	for (int i = 0; i < Q; ++i)
		cout<<ans[i]<<endl;
	return;
}
signed main(){
	read();
	mrdge(0, cnt, T);
	print();
	return 0;
}