CF1540D Inverse Inversions 题解

Description

传送门

Solution

算法一

考虑一个 $O(nq)$ 的算法。

我们动态地维护序列 $a$，其中 $a_i$ 表示 $p_i$ 是 $[1,i]$ 中第几小的。显然，有 $a_i=i-b_i$。

对于一次修改，直接 $O(1)$ 地更新 $a_i$。
对于一次查询，令 $ans=a_i$，并从 $i+1$ 一直扫描到 $n$。若当前的 $ans$ 满足 $ans\ge a_j$，那么就将 $ans$ 加上 $1$，否则不做操作。这里加 $1$ 的实质是，由于 $j>i$ 且 $p_j$ 在 $[1,j]$ 中从小到大的排名比 $[1,j-1]$ 中 $p_i$ 的排名要高，所以 $p_j<p_i$；此时，我们在 $[1,j-1]$ 的排名列表中在 $i$ 之前插入了一个位置，需将 $ans$ 加上 $1$ 表示 $i$ 在 $[1,j]$ 中的排名。

ans=a[i];
for (int j=i+1;j<=n;j++)
  if (ans>=a[j])  ans++;

最后输出 $ans$ 即可。

算法二

看到时限为 $\text{5.00s}$ 且 $n,q\le 10^5$，根据 $\text{Codeforces}$ 的做题经验，这显然是一道分块题（或者是一道树套树）。

考虑分块，令块长为 $B$。

对于每一个块，我们维护出一个列表 $v$。其实际意义是，若进入这个块时 $ans=x$，那么最终从这个块出来时 $ans=v_x+x$。

每次修改都是单点的，我们更新 $a_i$ 并重构该块的列表 $v$。
对于第一次查询，我们从左往右扫描，对于零散块直接在 $a$ 中扫描并更新 $ans$，对于整块直接调用列表更新 $ans$。

显然，每次查询都是 $O(B+\frac{n}{B})$ 的。关键在于，如何快速重构某一块的列表。

为方便叙述，令该块为 $[L,R]$。考虑从前往后扫描，实时维护一个数据结构，里面存储了若干个三元组 $(l,r,c)$，表示将 $v\in [l,r]$ 从块头推到当前位置时的值 $v+c$。每次更新的时候，我们通过二分找到 $\text{set}$ 里面第一个满足 $L\le a_i$ 的位置，并将该处的 $(l,r,c)$ 段成 $(l,L-1,c)$ 与 $(L,R,c+1)$，将后者向后面一个 $(l,r,c)$ 合并。接着，对于该位置往后的那一段后缀，将 $c$ 全部加一。显然，这个数据结构需要维护单点查询，单点插入与区间加，采用平衡树即可。由于最多有 $\log n$ 个区间，所以单次重构是 $O(B\log B)$ 的。

总复杂度 $O(q(B+\frac{n}{B}+B\log B))$。当 $n=10^5,q=10^5$ 时，取 $B=120$ 理论最优，实际应将 $B$ 调小（平衡树的常数很大）。

值得一提的是，若不维护若干个三元组，而用线段树维护长度为 $n$ 的序列，那么复杂度就是 $O(q(B+\frac{n}{B}+B\log n))$，理论劣于上面的解法。但是，线段树常数较小，所以其实二者差不太多。