◊ \Large\color{red}\Diamond ◊ 是关键转化。
首先有个比较简单的 dp,令 f u , i , j f_{u,i,j} f u , i , j 表示 u u u 子树内删了 i i i 条边,目前 u u u 连通块的大小为 j j j 是否可行。
转移就类似树上背包合并,考虑如何优化 j j j 这一维。
显然 j j j 的最大值是 T u − i L T_u-iL T u − i L ,最小值是 T u − i R T_u-iR T u − i R ,其中 T u T_u T u 表示 u u u 子树内的 A A A 和,所以 j j j 的极差 ≤ i ( R − L ) \le i(R-L) ≤ i ( R − L ) 。
接着证明对于 j = x , y , z ( x < y < z , z − x < R − L ) j=x,y,z(x<y<z,z-x<R-L) j = x , y , z ( x < y < z , z − x < R − L ) 时,状态 y y y 可以省略掉◊ \Large\color{red}\Diamond ◊ ,设选择了 y y y 后最后形成的连通块大小为 s z ∈ [ l , r ] sz\in[l,r] s z ∈ [ l , r ] ,则新增了 s z − y sz-y s z − y ,分讨后可以证明 s z − y + x sz-y+x s z − y + x 和 s z − y + z sz-y+z s z − y + z 至少一个会 ∈ [ l , r ] \in [l,r] ∈ [ l , r ] ,所以 y y y 状态可以不计。此时所有记录的递增状态 j 1 , j 2 … j s j_1,j_2\dots j_s j 1 , j 2 … j s 都满足 j i + 2 − j i ≥ R − L j_{i+2}-j_i\ge R-L j i + 2 − j i ≥ R − L ,而又O ( k ) \mathcal O(k) O ( k ) 个。所以状态只有 O ( n k 2 ) \mathcal O(nk^2) O ( n k 2 ) 了,总复杂度为 O ( n k 3 ) \mathcal O(nk^3) O ( n k 3 ) 。具体为什么不是 O ( n k 4 ) \mathcal O (nk^4) O ( n k 4 ) ,是因为 f u , i f_{u,i} f u , i 和 f v , j f_{v,j} f v , j 卷的时候由于树上背包的证明是 O ( n k ) \mathcal O(nk) O ( n k ) 的,剩余两个卷积是 O ( k 2 ) \mathcal O(k^2) O ( k 2 ) 的。
令 w ( l , r ) w(l,r) w ( l , r ) 表示区间 [ l , r ] [l,r] [ l , r ] 中所有子串中合法括号序列数量,可以证明 w w w 满足四边形不等式。令 f i , j f_{i,j} f i , j 表示前 i i i 个字符分成了 j j j 段的方案数,则 f i , j = min k < i f k − 1 , j − 1 + w ( k , i ) f_{i,j}=\min_{k<i}f_{k-1,j-1}+w(k,i) f i , j = min k < i f k − 1 , j − 1 + w ( k , i ) ,所以 F ( x ) = f n , x F(x)=f_{n,x} F ( x ) = f n , x 为下凸函数。
套路地,显示使用 wqs 二分,然后需要求 g i = min j < i g j − 1 + w ( j , i ) − m i d g_i=\min_{j<i}g_{j-1}+w(j,i)-mid g i = min j < i g j − 1 + w ( j , i ) − m i d 。因为 w w w 满足四边形不等式,所以 g g g 满足决策单调性,但 w w w 不好直接求,并且 g g g 求解过程需要在线。
发现 w ( l , r ) w(l,r) w ( l , r ) 在做 l → l ± 1 l\to l\pm 1 l → l ± 1 或 r → r ± 1 r\to r\pm 1 r → r ± 1 时可以直接 O ( 1 ) \mathcal O(1) O ( 1 ) 求解,所以可以用 cdq 分治套决策单调性分治解决,复杂度为 O ( n log 3 n ) \mathcal O(n\log^3n) O ( n log 3 n ) 。但是可以使用LARSCH/简化 LARSCH (感觉像科技)算法优化至 dp 部分只要 O ( n ) / O ( n log n ) \mathcal O(n)/\mathcal O(n\log n) O ( n ) / O ( n log n ) ,总复杂度为 O ( n log n / n log 2 n ) \mathcal O(n\log n/n\log^2n) O ( n log n / n log 2 n ) 。◊ \Large\color{red}\Diamond ◊
简化 LARSCH 算法其实类似于 cdq 分治套决策单调性分治,令 o p t t ( x ) opt_t(x) o p t t ( x ) 表示只考虑从 [ 1 , t ] [1,t] [ 1 , t ] 的转移到 x x x 的最优决策,o p t ( x ) = o p t x − 1 ( x ) opt(x)=opt_{x-1}(x) o p t ( x ) = o p t x − 1 ( x ) 。具体过程是定义 s o l v e ( l , r ) solve(l,r) s o l v e ( l , r ) 表示已知 [ 1 , l ) [1,l) [ 1 , l ) 的 g , o p t g,opt g , o p t 和 o p t l − 1 ( r ) opt_{l-1}(r) o p t l − 1 ( r ) ,求解 [ 1 , r ] [1,r] [ 1 , r ] 的 g , o p t g,opt g , o p t 。
用 i ∈ [ o p t l − 1 , o p t l − 1 , r ] i\in [opt_{l-1},opt_{l-1,r}] i ∈ [ o p t l − 1 , o p t l − 1 , r ] 求出 o p t l − 1 ( m i d ) opt_{l-1}(mid) o p t l − 1 ( m i d ) ,因为 o p t t ( x ) ≤ o p t t ( y ) ≤ o p t t ( z ) , x < y < z opt_{t}(x)\le opt_{t}(y)\le opt_{t}(z),x<y<z o p t t ( x ) ≤ o p t t ( y ) ≤ o p t t ( z ) , x < y < z 。
调用 s o l v e ( l , m i d ) solve(l,mid) s o l v e ( l , m i d ) 。
用 o p t ( i ) , i ∈ [ l , m i d ] opt(i),i\in [l,mid] o p t ( i ) , i ∈ [ l , m i d ] 和 o p t l − 1 ( r ) opt_{l-1}(r) o p t l − 1 ( r ) 求出 o p t m i d ( r ) opt_{mid}(r) o p t m i d ( r ) 。
调用 s o l v e ( m i d + 1 , r ) solve(mid+1,r) s o l v e ( m i d + 1 , r ) 。
注意 1 操作中要移动 w 1 w1 w 1 的指针,3 操作中要移动 w 2 w2 w 2 的指针,可以证明整个过程中 w 1 w1 w 1 和 w 2 w2 w 2 的指针一共只会移动 O ( n log n ) \mathcal O(n\log n) O ( n log n ) 次,若两种操作移动同一对指针会退化到 O ( n 2 ) \mathcal O(n^2) O ( n 2 ) 次。
考虑对于每个点的度数和儿子子树大小根号分治,对于度数 ≤ B \le B ≤ B 的点的询问考虑暴力枚举每个儿子,> B >B > B 的点只有 O ( n ) \mathcal O(\sqrt n) O ( n ) 个再考虑优化。◊ \Large\color{red}\Diamond ◊
记 f ( u , l , r ) f(u,l,r) f ( u , l , r ) 表示 u u u 子树中编号在 [ l , r ] [l,r] [ l , r ] 的点的数量,则对于询问 ( l , r , x ) (l,r,x) ( l , r , x ) 答案就是 ( f ( x , l , r ) 2 ) − ∑ v is son of x ( f ( v , l , r ) 2 ) \binom{f(x,l,r)}2-\sum_\text{v is son of x}\binom{f(v,l,r)}{2} ( 2 f ( x , l , r ) ) − ∑ v is son of x ( 2 f ( v , l , r ) ) 。对于 d e g x ≤ n deg_x\le \sqrt n d e g x ≤ n ,直接枚举 v v v 后就会形成 O ( n n ) \mathcal O(n\sqrt n) O ( n n ) 个求子树内编号在 [ l , r ] [l,r] [ l , r ] 的数量的询问,把子树用 dfs 序拍掉就是 O ( n ) \mathcal O(n) O ( n ) 个点 O ( n n ) \mathcal O(n\sqrt n) O ( n n ) 次询问二维数点,离线后就是 O ( n ) \mathcal O(n) O ( n ) 次单点加,O ( n n ) \mathcal O(n\sqrt n) O ( n n ) 次区间求和,可以使用分块内记录块内和块间前缀和做到 O ( n ) − O ( 1 ) \mathcal O(\sqrt n)-\mathcal O(1) O ( n ) − O ( 1 ) 修改/查询,复杂度为 O ( n n ) \mathcal O(n\sqrt n) O ( n n ) 。但是空间还是 O ( n n ) \mathcal O(n\sqrt n) O ( n n ) ,因为要存下 O ( n n ) \mathcal O(n\sqrt n) O ( n n ) 个询问,可以发现 x x x 的儿子的询问区间互不相交,所以可以扫描求完最小的区间后再把下一个加入,空间就是 O ( n ) \mathcal O(n) O ( n ) 的了。
剩下考虑 d e g x > n deg_x >\sqrt n d e g x > n 的怎么做。再次考虑分成两种情况,对于 x x x 的儿子 v v v 的子树大小 > n >\sqrt n > n 的,x x x 只会有 n \sqrt n n 个这种儿子,可以理解为度数只有 n \sqrt n n ,所以直接套用上面的做法即可。◊ \Large\color{red}\Diamond ◊ 对于 v v v 子树大小 ≤ n \le \sqrt n ≤ n ,直接暴力 O ( s i z v 2 ) \mathcal O(siz_v^2) O ( s i z v 2 ) 枚举出所有二元组 ( x , y ) , x < y (x,y),x<y ( x , y ) , x < y ,若满足 l ≤ x < y ≤ r l\le x< y\le r l ≤ x < y ≤ r 就会贡献答案,这也是一个二维偏序(数点),但是现在有 O ( n n ) \mathcal O(n\sqrt n) O ( n n ) 对修改,O ( m ) \mathcal O(m) O ( m ) 个询问,使用 O ( 1 ) − O ( n ) \mathcal O(1)-\mathcal O(\sqrt n) O ( 1 ) − O ( n ) 的分块处理即可。这部分的复杂度还是 O ( n n ) \mathcal O(n\sqrt n) O ( n n ) 的,修改简单处理就可以做到线性空间。
总结一下就是直接对于每个点的儿子子树中前 n \sqrt n n 大的做 O ( n ) − O ( 1 ) \mathcal O(\sqrt n) -O(1) O ( n ) − O ( 1 ) 的二维数点,对剩余的做 O ( 1 ) − O ( n ) \mathcal O(1)-\mathcal O(\sqrt n) O ( 1 ) − O ( n ) 的二维数点。
首先将二维数点离线扫描线改成对 a a a 做区间 ± 1 \pm 1 ± 1 ,由于 a i ≤ n a_i\le n a i ≤ n ,所以 m m m 的倍数只会有 O ( n m ) \mathcal O(\frac nm) O ( m n ) 种。考虑以 B B B 为块长做序列分块,区间加是好做的,由于初始时 a a a 的差分数组都是 c i = 0 c_i=0 c i = 0 ,每次区间加会让 c l ← c l ± 1 , c r + 1 ← c r + 1 ∓ 1 c_l\gets c_{l}\pm 1,c_{r+1}\gets c_{r+1}\mp 1 c l ← c l ± 1 , c r + 1 ← c r + 1 ∓ 1 ,只会有 O ( n ) \mathcal O(n) O ( n ) 次区间加,所以 ∑ ∣ c i ∣ ≤ 2 n \sum |c_i|\le 2n ∑ ∣ c i ∣ ≤ 2 n ,所以每个块内的 max − min \max-\min max − min 的和是 O ( n ) \mathcal O(n) O ( n ) 级别的。询问时对每个块内枚举 m m m 的倍数并查询出现次数即可,m m m 的倍数一共只会有 O ( max − min m ) \mathcal O(\frac {\max-\min} m) O ( m max − min ) 个,所以查询复杂度为 O ( n m ) \mathcal O(\frac n m) O ( m n ) ,区间加复杂度为 O ( n B + B ) \mathcal O(\frac nB+B) O ( B n + B ) ,平衡后整体复杂度为 O ( n n ) \mathcal O(n\sqrt n) O ( n n ) 。注意重构的时候应该对每个块内的数删除而不是 O ( max − min ) \mathcal O(\max-\min) O ( max − min ) 的清空,否则复杂度会退化到 O ( n 2 ) \mathcal O(n^2) O ( n 2 ) 。
这也太困难了。
p ← 317 p\gets 317 p ← 3 1 7 ,固定 c c c 数组后,分配宝石的方案数就是 ( n c 1 , c 2 , … c m ) \binom {n}{c_1,c_2,\dots c_m} ( c 1 , c 2 , … c m n ) ,由 Lucas 定理的过程可以得到 ( n c 1 , c 2 … c n ) ≡ 0 ( m o d p ) \binom{n}{c_1,c_2\dots c_n}\not\equiv 0(\bmod p) ( c 1 , c 2 … c n n ) ≡ 0 ( m o d p ) 的必要条件是 ∀ k , ∑ b i t k ( c i ) = b i t k ( n ) \forall k,\sum bit_k(c_i)=bit_k(n) ∀ k , ∑ b i t k ( c i ) = b i t k ( n ) ,其中 b i t k ( x ) bit_k(x) b i t k ( x ) 表示 x x x 在 p p p 进制下的第 k k k 位,即需要满足 ∑ c i \sum c_i ∑ c i 求和过程中不能在 p p p 进制下进位。◊ \Large\color{red}\Diamond ◊
先考虑 n < p n<p n < p ,此时就只有一层。考虑直接枚举排序后 c i c_i c i 之间的偏序关系,即相邻的两个位置之间的关系是 < 还是 =,一共只有 O ( 2 m − 1 ) \mathcal O(2^{m-1}) O ( 2 m − 1 ) 种。然后在确定相对偏序关系后,( n c 1 , c 2 … c m ) \binom n {c_1,c_2\dots c_m} ( c 1 , c 2 … c m n ) 的贡献和 A , B A,B A , B 造成的贡献就独立了,分别考虑。◊ \Large\color{red}\Diamond ◊
先考虑求出 A , B A,B A , B 的贡献。令 d p S 1 , S 2 dp_{S1,S2} d p S 1 , S 2 表示已经填入了集合 S 1 S1 S 1 的选手,前面的偏序关系是 S 2 S2 S 2 的贡献和,转移考虑枚举一个新的极长的相等段,选手编号集合为 S 3 S3 S 3 ,然后可以提前预处理出 w S 1 , S 3 w_{S1,S3} w S 1 , S 3 造成的贡献进行转移。复杂度为 O ( ∑ i = 1 m ( m i ) × 2 i − 1 × 2 m − i ) = O ( 4 m ) \mathcal O(\sum_{i=1}^m \binom{m}{i}\times 2^{i-1}\times 2^{m-i})=\mathcal O(4^m) O ( ∑ i = 1 m ( i m ) × 2 i − 1 × 2 m − i ) = O ( 4 m ) 。
然后考虑求 ( n c 1 , c 2 … c m ) \binom{n}{c_1,c_2\dots c_m} ( c 1 , c 2 … c m n ) 的贡献。由于 c i c_i c i 的顺序不影响,所以这里只计算 c c c 有序时的贡献和。令 f i , j , S f_{i,j,S} f i , j , S 表示前 i i i 个数和为 j j j ,相对偏序关系为 S S S 的总贡献,转移时讨论 c i + 1 c_{i+1} c i + 1 与 c i c_i c i 的偏序关系。复杂度为 O ( ∑ i = 1 m 2 i − 1 p 2 ) = O ( 2 m p 2 ) \mathcal O(\sum_{i=1}^m 2^{i-1}p^2)=\mathcal O(2^mp^2) O ( ∑ i = 1 m 2 i − 1 p 2 ) = O ( 2 m p 2 ) 。
拓展到 n ≥ p n\ge p n ≥ p ,现在一共有 log p n \log_p n log p n 层。类似上面的做法,将两种贡献分开计算,A , B A,B A , B 造成的贡献直接 O ( 4 m ) \mathcal O(4^m) O ( 4 m ) 预处理即可。考虑如何适应层数对偏序关系的影响。从高位到低位,原本相等的一段区间可能会分裂成几个有序的区间,考虑用 dp 记录这个过程。令 f i , S f_{i,S} f i , S 表示到达 i i i 层时的偏序关系为 S S S 的方案数(b i t 2 ( S ) = [ c i ≠ c i + 1 ] bit_2(S)=[c_i \ne c_{i+1}] b i t 2 ( S ) = [ c i = c i + 1 ] ),考虑从 f i , S f_{i,S} f i , S 转移到 f i − 1 , S ′ f_{i-1,S'} f i − 1 , S ′ (S ⊆ S ′ S\sube S' S ⊆ S ′ ) 的转移系数 g S , S ′ , b i t i − 1 ( n ) g_{S,S',bit_{i-1}(n)} g S , S ′ , b i t i − 1 ( n ) ,实际上 g S , S ′ , i g_{S,S',i} g S , S ′ , i 可以将这个过程分成 S S S 的每个 0 区间(等价类)变成 S ′ S' S ′ 对应区间后对 ∑ b i t k ( c ) \sum bit_k(c) ∑ b i t k ( c ) 维做卷积得来的。具体的,令 h S , i h_{S,i} h S , i 表示相对偏序关系为 S S S ,∑ c = i \sum c=i ∑ c = i 的方案数,则 ( S , S ′ ) (S,S') ( S , S ′ ) 的转移可以通过 h S , i h_{S,i} h S , i 的第二维 i i i 的和卷积得来。h h h 的转移和上一段落中的类似,g g g 的状态中有许多可以通过记忆化解决◊ \Large\color{red}\Diamond ◊ ,比如下图中三个绿色段的相对顺序不会影响 g g g 的值,所以可以将其排序后存入map进行记忆化。令map的大小为 F ( m ) F(m) F ( m ) ,则预处理 g g g 的复杂度为 O ( F ( m ) p 2 ) \mathcal O(F(m)p^2) O ( F ( m ) p 2 ) 。若没有记忆化,F ( m ) = 3 m F(m)=3^m F ( m ) = 3 m (即枚举 S ⊆ S ′ S\sube S' S ⊆ S ′ ),但是显然记忆化后会大幅减少状态数,实际上 F ( 12 ) = 17547 \color{red}F(12)=17547 F ( 1 2 ) = 1 7 5 4 7 。对于每次询问跑 f f f 转移后记得最后与前面算的 A , B A,B A , B 造成的贡献 d p dp d p 拼起来。总复杂度为 O ( 4 m + F ( m ) p 2 + q 3 m log p n ) \mathcal O(4^m+F(m)p^2+q3^m\log_p n) O ( 4 m + F ( m ) p 2 + q 3 m log p n ) 。
倒闭了,写不了一点,还要卡常,下播了。upd. 改完了。
由于 A , B A,B A , B 很小,所以考虑直接枚举其中一个。令 f u , i f_{u,i} f u , i 表示到 u u u 时经过的 ∑ a = i \sum a=i ∑ a = i 的最小的 ∑ b \sum b ∑ b ,转移不难发现由于 a i > 0 a_i>0 a i > 0 ,所以这个转移是分层图,可以从小到大枚举 i i i 后枚举每条边转移,复杂度为 O ( ( n + m ) m V ) \mathcal O((n+m)mV) O ( ( n + m ) m V ) ,若没有发现分层图直接跑 dijkstra 会多一只 log 无法通过。
由于 k k k 和 ∑ k \sum k ∑ k 很小,可以考虑模拟费用流。首先有个显然的费用流模型,对每个 i i i 连边 ( S , i , 1 , 0 ) , ( i , i + 1 , 1 , a i ) , ( i , T , 1 , 0 ) (S,i,1,0),(i,i+1,1,a_i),(i,T,1,0) ( S , i , 1 , 0 ) , ( i , i + 1 , 1 , a i ) , ( i , T , 1 , 0 ) ,还有 k k k 的限制所以答案就是在这个图中增广 k k k 次的最大花费。每次增广就是找到当前的最大区间和,然后将经过的链上的边反向并将费用变成相反数 ◊ \Large\color{red}\Diamond ◊ 。这个过程可以使用线段树模拟,就是做最大子段和还有区间取相反数,复杂度为 O ( q log n + ∑ k n log n ) \mathcal O(q\log n+\sum kn\log n) O ( q log n + ∑ k n log n ) 。
其实也可以从反悔贪心的方面思考刚刚的过程,就是先取当前和最大的区间 [ l , r ] [l,r] [ l , r ] ,然后反悔再选取 [ l ′ , r ′ ] ⊆ [ l , r ] [l',r'] \sube [l,r] [ l ′ , r ′ ] ⊆ [ l , r ] 并减去 s u m [ l ′ , r ′ ] sum[l',r'] s u m [ l ′ , r ′ ] ,这样区间就断成了 [ l , l ′ ) [l,l') [ l , l ′ ) 和 ( r ′ , r ] (r',r] ( r ′ , r ] ,完成了反悔操作,也就对应模拟费用流的反向弧的操作。
依旧是模拟费用流。首先这显然是一个带权二分图匹配问题,考虑建立费用流后如何优化。
如果手玩这个退流和走反向弧的过程不难发现每次实际上就是在右边选出一个序列 p 1 , p 2 , … p m p_1,p_2,\dots p_m p 1 , p 2 , … p m 满足 p [ 1 , m ) p[1,m) p [ 1 , m ) 都曾经被流满了(匹配完 c p i c_{p_i} c p i 个左部点了),然后不断进行退流操作(即走反向弧)直到最后走到了一个没有流满的点 p m p_m p m 后走到 T T T ,下图红色的边就体现了退流的过程。由于 m ≤ k ≤ 10 m\le k\le 10 m ≤ k ≤ 1 0 ,所以考虑对右部点做一些本质的操作。◊ \Large\color{red}\Diamond ◊ 实际上每次寻找最短路时可以看作从 p i p_i p i 走到 p i + 1 p_{i+1} p i + 1 ,距离就是 p i → x → p i + 1 p_i\to x \to p_{i+1} p i → x → p i + 1 两条边的距离之和,这个可以提前预处理出任意两个右部点的 x x x 的集合处理,所以就是已知任意两个右部点之间的距离求一条最短路 r s ⇝ r t rs\rightsquigarrow rt r s ⇝ r t 且保证 r t rt r t 没有流满。直接 O ( k 3 ) \mathcal O(k^3) O ( k 3 ) 跑 floyd 即可,由于要模拟退流,所以要记录转移点来复原路径。复杂度为 O ( n k 3 + n k 2 log n ) \mathcal O(nk^3+nk^2\log n) O ( n k 3 + n k 2 log n ) 。
*3500。这种题写还是算了,太恐怖了一点。Ж ( l , r ) → f ( l , r ) \operatorname{Ж}(l,r)\to f(l,r) Ж ( l , r ) → f ( l , r )
f ( l , r ) = ∑ i = l r ∣ lcp ( s [ l , r ] , s [ i , r ] ) ∣ = ∑ i = l r ∣ min ( lcp ( s [ l , n ] , s [ i , n ] ) , r − i + 1 ) ∣ f(l,r)=\sum_{i=l}^r|\operatorname{lcp}(s[l,r],s[i,r])|=\sum_{i=l}^r|\min(\operatorname{lcp}(s[l,n],s[i,n]),r-i+1)|
f ( l , r ) = i = l ∑ r ∣ l c p ( s [ l , r ] , s [ i , r ] ) ∣ = i = l ∑ r ∣ min ( l c p ( s [ l , n ] , s [ i , n ] ) , r − i + 1 ) ∣
由于 lcp ( s [ i , n ] , s [ j , n ] ) = min k = r k i + 1 r k j h t k \operatorname{lcp}(s[i,n],s[j,n])=\min_{k=rk_{i}+1}^{rk_{j}}ht_k l c p ( s [ i , n ] , s [ j , n ] ) = min k = r k i + 1 r k j h t k ,所以考虑将 i ∈ [ l , r ] i\in [l,r] i ∈ [ l , r ] 分成 r k i < r k l rk_i< rk_l r k i < r k l 和 r k i > r k l rk_i>rk_l r k i > r k l 和 i = l i=l i = l ,然后将询问挂在 r k l rk_l r k l 上。后面其实就是数据结构的事了。
r k i > r k l rk_i>rk_l r k i > r k l 和 r k i < r k l rk_i<rk_l r k i < r k l 是本质相同的,这里只讨论前者。区间最小值的拆贡献考虑离线后对 r k rk r k 分治,令 s u f i = min j = i m i d h t j suf_i=\min_{j=i}^{mid} ht_j s u f i = min j = i m i d h t j ,p r e i = min j = m i d i h t j pre_i=\min_{j=mid}^i ht_j p r e i = min j = m i d i h t j 。从 m i d → L mid\to L m i d → L 扫描线,当遇到询问 [ l , r ] [l,r] [ l , r ] 时,令 m n = s u f r k l mn=suf_{rk_l} m n = s u f r k l ,满足∀ m i d < i ≤ p , p r e i ≥ m n \forall mid<i\le p,pre_i\ge mn ∀ m i d < i ≤ p , p r e i ≥ m n 的 p p p 最大,m n mn m n 单调不降,p r e pre p r e 单调不降,因此 p p p 单调不降,考虑拆开贡献成 ∑ i = m i d + 1 p [ s a i ∈ [ l , r ] ] min ( m n , r − s a i + 1 ) + ∑ i = p + 1 R [ s a i ∈ [ l , r ] min ( p r e i , r − s a i + 1 ) ] \sum_{i=mid+1}^p [sa_i\in [l,r]]\min(mn,r-sa_i+1)+\sum_{i=p+1}^R[sa_i\in [l,r]\min(pre_i,r-sa_i+1)] ∑ i = m i d + 1 p [ s a i ∈ [ l , r ] ] min ( m n , r − s a i + 1 ) + ∑ i = p + 1 R [ s a i ∈ [ l , r ] min ( p r e i , r − s a i + 1 ) ] 。前半部分是讨论对 ( s a i , i ) (sa_i,i) ( s a i , i ) 的二维数点,扫 p p p 的时候加入树状数组即可。后半部分当 p r e i ≤ r − s a i + 1 pre_i\le r-sa_i+1 p r e i ≤ r − s a i + 1 时还需要满足 s a i ∈ [ l , r ] , i ∈ [ p + 1 , R ] sa_i\in [l,r],i\in [p+1,R] s a i ∈ [ l , r ] , i ∈ [ p + 1 , R ] ,所以是对 ( i , s a i , s a i + p r e i ) (i,sa_i,sa_i+pre_i) ( i , s a i , s a i + p r e i ) 的三维数点,考虑优化。题解做法是将限制拆成 m n ≤ r − s a i + 1 mn\le r-sa_i+1 m n ≤ r − s a i + 1 加上 p r e i ≤ r − s a i + 1 < m n pre_i\le r-sa_i+1< mn p r e i ≤ r − s a i + 1 < m n ,因为 p r e i < m n pre_i< mn p r e i < m n 。第一个限制是二维数点,乍一看第二个限制还是三维数点,但是由于对于 i ∈ ( m i d , p ] i\in (mid,p] i ∈ ( m i d , p ] 的都满足 p r e i ≥ m n pre_i\ge mn p r e i ≥ m n ,所以就忽略掉了对 i ∈ [ p + 1 , R ] i\in [p+1,R] i ∈ [ p + 1 , R ] 的限制也不会多算,所以就改成了二维数点。◊ \Large\color{red}\Diamond ◊ 其他的讨论都类似就算了。
复杂度为 O ( ( n + q ) log 2 n ) \mathcal O((n+q)\log^2n) O ( ( n + q ) log 2 n ) 。