ACM-ICPC 2019南昌网络赛I题 Yukino With Subinterval
题目大意：给一个长度为n，值域为[1, n]的序列{a}，要求支持m次操作：
- 单点修改
1 pos val - 询问子区间中某个值域的数的个数，连续的相同数字只记为一个。（即统计数字段的个数）
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; //无需修正- 分治时，右半边的修改操作不需要再应用（其对右半边的查询已经被统计，对左半边的查询无影响，）
- 写分治的题时有几个不错的调试输出点，可以很快地找到错误：
- 树状数组的修改、查询函数开头处
- 查询操作的统计答案处
- 递归头部输出L、R，确定函数调用的状态
CDQ真的快。考场上队友用树套树，手写Treap + 读优还是被卡常，事后CDQ分治连cin都不用改就过了·_·！
AC代码：
#include
#include
#include
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 = "<= 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<