首先有个比较简单的 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 ) 。

题解做法简直抽象。
将两个操作变得更相似,可以发现
r i = 0 : c i − x + y 2 = c i − 1 + x − y 2 r i = 1 : c i − x + y 2 = − c i − 1 − x − y 2 ◊ r_i=0:c_i-\frac{x+y}2=c_{i-1}+\frac{x-y}2 \\
r_i=1:c_i-\frac{x+y}2=-c_{i-1}-\frac{x-y}2 \Large\color{red}\Diamond
r i = 0 : c i − 2 x + y = c i − 1 + 2 x − y r i = 1 : c i − 2 x + y = − c i − 1 − 2 x − y ◊
所以任意操作后都有 ( c i − x + y 2 ) 2 = ( c i − 1 + x − y 2 ) 2 (c_i-\frac{x+y}2)^2=(c_{i-1}+\frac {x-y}2)^2 ( c i − 2 x + y ) 2 = ( c i − 1 + 2 x − y ) 2 ,拆开得到 c i 2 − c i − 1 2 − ( x + y ) c i − ( x − y ) c i − 1 + x y = 0 c_i^2-c_{i-1}^2-(x+y)c_i-(x-y)c_{i-1}+xy=0 c i 2 − c i − 1 2 − ( x + y ) c i − ( x − y ) c i − 1 + x y = 0 。对于 i ∈ [ 1 , n ] i \in[1,n] i ∈ [ 1 , n ] 求和后就得到了 c n 2 − y c n + n x y = x ∑ ( c i + c i + 1 ) c_n^2-yc_n+nxy=x\sum(c_i+c_{i+1}) c n 2 − y c n + n x y = x ∑ ( c i + c i + 1 ) ,由于 f ( r ) = ∑ c − i f(r)=\sum c-i f ( r ) = ∑ c − i ,所以 f ( r ) = c n 2 + ( x − y ) c n + n x y 2 x f(r)=\frac{c_n^2+(x-y)c_n+nxy}{2x} f ( r ) = 2 x c n 2 + ( x − y ) c n + n x y ,所以知道 c n c_n c n 就能直接求出 f ( r ) f(r) f ( r ) 了。由于每次操作二都会将之前的 y y y 的系数 × ( − 1 ) + 1 \times(-1)+1 × ( − 1 ) + 1 ,由归纳法可得 y y y 的系数只会是 0 / 1 0/1 0 / 1 ,所以只会有至多 2 n 2n 2 n 种 c n c_n c n ,求解直接设 f i , j ∈ { 0 , 1 } , k f_{i,j\in\{0,1\},k} f i , j ∈ { 0 , 1 } , k 表示前 i i i 个数 y y y 的系数为 j j j ,x x x 的系数为 k k k 是否可行,转移用 bitset 优化,复杂度 O ( n 2 ω ) \mathcal O(\frac{n^2}\omega) O ( ω n 2 ) ,要滚动数组。
*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 ) 。
由于每次询问只会删除一个点,所以可以考虑缺一分治。◊ \Large\color{red}\Diamond ◊ 同时由于 n ≤ 300 n\le 300 n ≤ 3 0 0 ,所以考虑用 floyd 做最短路。具体的令 s o l v e ( l , r ) solve(l,r) s o l v e ( l , r ) 表示进入时 [ 1 , l ) [1,l) [ 1 , l ) 和 ( r , n ] (r,n] ( r , n ] 都已经加入图中并处理完最短路,则答案直接离线后在 s o l v e ( p i , p i ) solve(p_i,p_i) s o l v e ( p i , p i ) 中求即可。每次加入一段区间 [ l , r ] [l,r] [ l , r ] 只要将其设为 floyd 中的转移点即可,因为不是它作为转移点的在之前已经处理过了,大概过程就是 i ∈ [ l , r ] , j ∈ [ 1 , n ] , k ∈ [ 1 , n ] : f j , k ← f j , i + f j , k i\in[l,r],j\in [1,n],k\in [1,n] :f_{j,k}\gets f_{j,i}+f_{j,k} i ∈ [ l , r ] , j ∈ [ 1 , n ] , k ∈ [ 1 , n ] : f j , k ← f j , i + f j , k ,这个过程单次复杂度为 O ( n 2 s i z ) \mathcal O(n^2siz) O ( n 2 s i z ) ,由于在分治过程中 ∑ s i z = O ( n log n ) \sum siz=\mathcal O(n\log n) ∑ s i z = O ( n log n ) ,所以总复杂度为 O ( n 3 log n + q ) \mathcal O(n^3\log n+q) O ( n 3 log n + q ) 。由于 floyd 无法删除,撤销只能记录当时的最短路数组后拷贝,所以空间复杂度是 O ( n 3 ) \mathcal O(n^3) O ( n 3 ) 。
每次新到达的点 x x x 都要满足 t x ≥ m a x c t_x\ge maxc t x ≥ m a x c ,所以考虑从小到大枚举 m a x c maxc m a x c 并转移能到达的最大 s s s 。若 t i < m a x c t_i<maxc t i < m a x c 则现在无法被加入和到达,直接全部删除不考虑。剩余的若 c i ≤ m a x c c_i\le maxc c i ≤ m a x c 则当前点可以转移到其他点也可以被转移到,将其称为 1 类点,否则只能被其他点到达转移称为 2 类点。若一条边连接了两个 1 类点,则可以互相转移,令 b e l ( x ) bel(x) b e l ( x ) 为 x x x 所在的极大 1 类点连通块。而 1 类点连向 2 类点的边只能单向转移,并且 2 类点当前无法再向外转移,所以最后的图就是若干个 1 类连通块像菊花一样向外连若干条单向边。
对于一条边 ( u , v ) (u,v) ( u , v ) 会在 m a x c ∈ [ max ( c u , c v ) , min ( t u , t v ) ] maxc\in [\max(c_u,c_v),\min(t_u,t_v)] m a x c ∈ [ max ( c u , c v ) , min ( t u , t v ) ] 时称为 1 类边,在 m a x c ∈ [ c u , min ( t u , c v − 1 ) ] maxc\in [c_u,\min(t_u,c_v-1)] m a x c ∈ [ c u , min ( t u , c v − 1 ) ] 时成为 2 类边,注意反向边也要建。考虑如何顺着这个图进行转移,由于连通块之间的信息是相同的,而且不好撤销,而一条边会在一个区间出现,所以套路地使用线段树分治,1 类连通块用可撤销并查集维护。复杂度为 O ( n log 2 n ) \mathcal O(n\log^2n) O ( n log 2 n ) 。
其实到这一步时可以求出每个时刻每个连通块的答案,但是还无法将连通块的贡献传到每个点。优化可以发现对于一条无向边 ( u , v ) (u,v) ( u , v ) ,其产生的 1 类边时间一定完全晚于 2 类边时间,所以需要在线段树上倒着扫 m a x c maxc m a x c ,否则就会导致先扫到的 2 类点在被删掉后更新 1 类连通块。
容斥一下,方案数为
a n s = ∑ i + j > 0 ( − 1 ) i + j − 1 ( n i ) ( n j ) f ( i , j ) ans = \sum_{i+j>0}(-1)^{i+j-1}\binom{n}{i}\binom{n}{j}f(i,j)
a n s = i + j > 0 ∑ ( − 1 ) i + j − 1 ( i n ) ( j n ) f ( i , j )
其中 f ( i , j ) f(i,j) f ( i , j ) 表示至少有 i i i 列同色和 j j j 行同色的方案数,则 f ( 0 , i ) = f ( i , 0 ) = 3 i + n ( n − i ) f(0,i)=f(i,0)=3^{i+n(n-i)} f ( 0 , i ) = f ( i , 0 ) = 3 i + n ( n − i ) ,f ( i , j ) = 3 n 2 − ( i + j ) n + i j + 1 f(i,j)=3^{n^2-(i+j)n+ij+1} f ( i , j ) = 3 n 2 − ( i + j ) n + i j + 1 。i = 0 i=0 i = 0 和 j = 0 j=0 j = 0 分开讨论,剩下是 i , j > 0 i,j>0 i , j > 0 。
∑ i , j > 0 ( − 1 ) i + j − 1 ( n i ) ( n j ) 3 n 2 − ( i + j ) n + i j + 1 = 3 n 2 + 1 ∑ i = 1 n ( − 1 ) i − 1 3 − i n ( n i ) ∑ j = 1 n ( − 1 ) j ( n j ) 3 − j n + i j = − 3 n 2 + 1 ∑ i = 1 n ( − 1 ) i 3 − i n ( n i ) ∑ j = 1 n ( n j ) ( − ( 3 − n + i ) ) j \sum_{i,j>0}(-1)^{i+j-1}\binom ni\binom nj3^{n^2-(i+j)n+ij+1} \\
=3^{n^2+1}\sum_{i=1}^n (-1)^{i-1}3^{-in}\binom ni\sum_{j=1}^n(-1)^j\binom nj 3^{-jn+ij}\\
=-3^{n^2+1}\sum_{i=1}^n (-1)^{i}3^{-in}\binom ni\sum_{j=1}^n\binom nj (-(3^{-n+i}))^j
i , j > 0 ∑ ( − 1 ) i + j − 1 ( i n ) ( j n ) 3 n 2 − ( i + j ) n + i j + 1 = 3 n 2 + 1 i = 1 ∑ n ( − 1 ) i − 1 3 − i n ( i n ) j = 1 ∑ n ( − 1 ) j ( j n ) 3 − j n + i j = − 3 n 2 + 1 i = 1 ∑ n ( − 1 ) i 3 − i n ( i n ) j = 1 ∑ n ( j n ) ( − ( 3 − n + i ) ) j
看到 j j j 式子的很像二项式定义,考虑
( − ( 3 − n + i ) + 1 ) n = ∑ j = 0 n ( n j ) ( − 3 − n + i ) j ∑ j = 1 n ( n j ) ( − ( 3 − n + i ) ) j = ( − ( 3 − n + i ) + 1 ) n − 1 (-(3^{-n+i})+1)^n=\sum_{j=0}^n\binom nj(-3^{-n+i})^j \\
\sum_{j=1}^n\binom nj (-(3^{-n+i}))^j=(-(3^{-n+i})+1)^n-1
( − ( 3 − n + i ) + 1 ) n = j = 0 ∑ n ( j n ) ( − 3 − n + i ) j j = 1 ∑ n ( j n ) ( − ( 3 − n + i ) ) j = ( − ( 3 − n + i ) + 1 ) n − 1
然后把式子带回去就可以 O ( n log n ) \mathcal O(n\log n) O ( n log n ) 做完了。
\begin{align}
\sum F^k &= \sum_{i=0}^k{k\brace i}i!\sum \binom{F}{i} \\
&= \sum_{i=0}^k{k\brace i}i!\binom{n}{i}\lceil\frac{m}2\rceil^i m^{n-i}
\end{align}
第二个式子是因为考虑 ( F i ) \binom Fi ( i F ) 的组合意义就是在所有的奇数编号球中选择 i i i 个的方案数,所以可以先钦定被选择的 i i i 个位置选奇数,然后剩下的随便选。代码实现时将 i ! ( n i ) i!\binom ni i ! ( i n ) 转成了 n i ‾ n^{\underline i} n i 求解,复杂度 O ( T k 2 ) \mathcal O(Tk^2) O ( T k 2 ) 。枚举 i i i 时动态维护 n i ‾ n^{\underline i} n i 做到 O ( T k ) \mathcal O(Tk) O ( T k ) 。
神秘困难图染色构造题,这个 ∣ S ∣ = 13 |S|=13 ∣ S ∣ = 1 3 甚至有用。
先把限制 a n s u ≠ a n s v ans_u\neq ans_v a n s u = a n s v 的边 ( u , v ) (u,v) ( u , v ) 先建出来,无向边数只有 6 n 6n 6 n 条,若有自环则显然无解。考虑稀疏图的性质,由于有向边数只有 12 n 12n 1 2 n ,而一共只有 n n n 个点,由鸽巢原理得至少有一个点的度数 ≤ 12 ◊ \le 12\Large\color{red}\Diamond ≤ 1 2 ◊ ,否则至少会有 13 n 13n 1 3 n 条有向边,矛盾。然后删除一个点后一定还有至少一个点度数 ≤ 12 \le 12 ≤ 1 2 ,所以可以考虑做一个拓扑排序类似的事情。取到拓扑序后倒着染色,每个点此时至多会有 12 12 1 2 条出度,但是有 13 13 1 3 种颜色,所以可以染上。
考虑将一个在二维上的点 ( x , y ) (x,y) ( x , y ) 转化成一维上的一段区间 [ x , n − y ] [x,n-y] [ x , n − y ] ◊ \Large\color{red}\Diamond ◊ ,由于所有的点在直角三角形中所以一定有 x + y ≤ n x+y\le n x + y ≤ n 。然后发现 H 操作就是将 r ≥ n − L , l ← max ( l , n − L ) r\ge n-L,l\gets \max(l,n-L) r ≥ n − L , l ← max ( l , n − L ) ,V 操作就是将 l ≤ L , r ← min ( r , L ) l\le L,r\gets \min(r,L) l ≤ L , r ← min ( r , L ) 。考虑建立猫树(线段树)后将所有区间挂上去,下面只考虑操作 H 对于 r ≥ x , l ← max ( l , x ) r\ge x,l\gets \max(l,x) r ≥ x , l ← max ( l , x ) 。
若 x ≤ m i d x\le mid x ≤ m i d ,则需要将所有区间按照 l l l 排序后取一段前缀推平成 x x x ,直接维护 l l l 相同的连续段即可。
若 x > m i d x>mid x > m i d ,则所有当前线段树节点上的区间都会被移动到右儿子,这个过程对于每个区间显然只会出现 O ( log n ) \mathcal O(\log n) O ( log n ) 次,所以直接暴力将当前节点的区间取出后全部重新插入即可。
维护这个过程应该直接对每个线段树节点都使用堆和 odt 维护就可以了。总复杂度为 O ( ( n + q ) log 2 n ) \mathcal O((n+q)\log^2n) O ( ( n + q ) log 2 n ) ,难写+难调+速度慢+空间大。
题面太臭了,其实就是要求:
n ← ∣ A ∣ , m ← ∣ B ∣ max i = 1 P LCP ( A [ i , P ] , B [ m − S + 1 , m ] ) = max i = 1 P min ( P − i + 1 , LCP ( A [ i , n ] , B [ m − S + 1 , m ] ) ) n\gets |A|,m\gets |B|\\
\max_{i=1}^P \operatorname{LCP}(A[i,P],B[m-S+1,m]) \\
=\max_{i=1}^P \min(P-i+1,\operatorname{LCP}(A[i,n],B[m-S+1,m])) \\
n ← ∣ A ∣ , m ← ∣ B ∣ i = 1 max P L C P ( A [ i , P ] , B [ m − S + 1 , m ] ) = i = 1 max P min ( P − i + 1 , L C P ( A [ i , n ] , B [ m − S + 1 , m ] ) )
求两个串后缀的 LCP,套路地对串 A+#+B 做 SA◊ \Large\color{red}\Diamond ◊ ,则答案就变成了
max i = 1 P min ( P − i + 1 , min j = r k i + 1 r k l e n − S + 1 h t j ) \max_{i=1}^P\min(P-i+1,\min_{j=rk_i+1}^{rk_{len-S+1}} ht_j)
i = 1 max P min ( P − i + 1 , j = r k i + 1 min r k l e n − S + 1 h t j )
后面的 min \min min 实际上是上下界可能会相反,大概意思知道就行了。
这个区间 min \min min 还是不太好求,考虑将 r k i rk_i r k i 排序后建立对应的 h t ht h t 的小根笛卡尔树,然后就将区间 min \min min 转化成求笛卡尔树上的 LCA◊ \Large\color{red}\Diamond ◊ 。考虑枚举 r k l e n − Q + 1 rk_{len-Q+1} r k l e n − Q + 1 的 LCA 为 u u u ,则对答案贡献 min ( P − m n u + 1 , h t u ) \min(P-mn_u+1,ht_u) min ( P − m n u + 1 , h t u ) ,其中 m n u mn_u m n u 表示 u u u 子树内的编号最小值。发现 u u u 越往上 m n u mn_u m n u 越小,P − m n u + 1 P-mn_u+1 P − m n u + 1 越大,h t u ht_u h t u 越小,所以可以在链上二分求解两个函数交点就是答案。复杂度 O ( ∑ ( n + m + q ) log ( n + m ) ) \mathcal O(\sum (n+m+q)\log (n+m)) O ( ∑ ( n + m + q ) log ( n + m ) ) 。
令 f i , j f_{i,j} f i , j 表示长度为 i i i 的排列,当前的权值和为 j j j 的方案数。
f i , j = ∑ k = 0 i − 1 ∑ l = 0 j − i ( i − 1 k ) f k , l f i − 1 − k , j − i − l f_{i,j}=\sum_{k=0}^{i-1}\sum_{l=0}^{j-i} \binom{i-1}k f_{k,l}f_{i-1-k,j-i-l}
f i , j = k = 0 ∑ i − 1 l = 0 ∑ j − i ( k i − 1 ) f k , l f i − 1 − k , j − i − l
由于第二维的大小为 O ( n 2 ) \mathcal O(n^2) O ( n 2 ) 的,所有直接做复杂度是 O ( n 6 ) \mathcal O(n^6) O ( n 6 ) ,即使用 FFT 也只能做到 O ( n 4 log n ) \mathcal O(n^4\log n) O ( n 4 log n ) 。
看到第二维是卷积形式,考虑对序列 f i f_i f i 做生成函数◊ \Large\color{red}\Diamond ◊ 。令 F i ( x ) = ∑ j = i n 2 f i , j x j F_i(x)=\sum_{j=i}^{n^2}f_{i,j} x^j F i ( x ) = ∑ j = i n 2 f i , j x j ,则转移式可以写成
F i = ∑ k = 0 i − 1 ( i − 1 k ) x i F k F i − 1 − k F_i=\sum_{k=0}^{i-1}\binom{i-1}k x^i F_kF_{i-1-k}
F i = k = 0 ∑ i − 1 ( k i − 1 ) x i F k F i − 1 − k
由于 F i F_i F i 是一个至多 n 2 n^2 n 2 次的多项式,所以只要求出 O ( n 2 ) \mathcal O(n^2) O ( n 2 ) 个点值后拉差即可复原系数。◊ \Large\color{red}\Diamond ◊ 对于每个 i i i 维护 g i , j = F i ( j ) g_{i,j}=F_i(j) g i , j = F i ( j ) ,卷积时直接对位乘即可,最后要给 n 2 n^2 n 2 次多项式做拉差复杂度为 O ( ( n 2 ) 2 ) = O ( n 4 ) \mathcal O((n^2)^2)=\mathcal O(n^4) O ( ( n 2 ) 2 ) = O ( n 4 ) 。
这题能对 FFT 加速的主要原因是需要做 O ( n 2 ) \mathcal O(n^2) O ( n 2 ) 次卷积,若只维护点值单次做卷积的复杂度就是 O ( n 2 ) \mathcal O(n^2) O ( n 2 ) ,最后在拉差还原时也是 O ( n 4 ) \mathcal O(n^4) O ( n 4 ) 的。若只做一次,拉差还原部分就已经是 O ( n 2 ) \mathcal O(n^2) O ( n 2 ) 的,不如暴力。
已知 k k k 次多项式的 k + 1 k+1 k + 1 个点值,O ( k 2 ) \mathcal O(k^2) O ( k 2 ) 还原多项式。
将拉格朗日插值式子展开:
f ( x ) = ∑ i f ( x i ) ∏ j ≠ i x − x j x i − x j = ∑ i f ( x i ) ∏ j ≠ i ( x i − x j ) ∏ j ≠ i ( x − x j ) f(x)=\sum_if(x_i)\prod_{j\neq i}\frac{x-x_j}{x_i-x_j} \\
=\sum_{i}\frac{f(x_i)}{\prod_{j\neq i}(x_i-x_j)} \prod_{j\neq i} (x-x_j) \\
f ( x ) = i ∑ f ( x i ) j = i ∏ x i − x j x − x j = i ∑ ∏ j = i ( x i − x j ) f ( x i ) j = i ∏ ( x − x j )
前半部分的值是已知的,直接 O ( k 2 ) \mathcal O(k^2) O ( k 2 ) 枚举 i , j i,j i , j 求即可,后半部分就是多个一次函数相乘,直接每次暴力展开是 O ( k 3 ) \mathcal O(k^3) O ( k 3 ) 的,但是可以写成 ∏ ( x − x j ) x − x i \frac{\prod(x-x_j)}{x-x_i} x − x i ∏ ( x − x j ) ,然后预处理 ∏ ( x − x i ) \prod(x-x_i) ∏ ( x − x i ) 的系数后每次做大除法即可,复杂度为 O ( k 2 ) \mathcal O(k^2) O ( k 2 ) 。
假的追忆,但更困难了。
首先这个题的空间是开不下 O ( n 2 ) \mathcal O(n^2) O ( n 2 ) 的 bitset 的,但是由于只需要求可达的集合大小,所以可以类似将 [ 1 , n ] [1,n] [ 1 , n ] 按照 B B B 分块◊ \Large\color{red}\Diamond ◊ ,每次只求到达一个块内元素的可达性,空间为 O ( n B ) \mathcal O(nB) O ( n B ) ,时间为 O ( n 2 ω + n 2 B ) \mathcal O(\frac{n^2}\omega+\frac{n^2}B) O ( ω n 2 + B n 2 ) ,只要将 B B B 开到刚好卡到空间限制就基本不会对时间复杂度有影响。
其实这个矩阵内的点是假的,就算开到空间内的点/k 维偏序 也只是常数上的区别,因为都已经有 O ( n 2 ω ) \mathcal O(\frac{n^2}\omega) O ( ω n 2 ) 的复杂度了就不用考虑怎么拆多维偏序了,直接对于每一维求出在询问区间内的点的集合,然后全部 & 起来就是满足的集合。需要求满足的点的可达集合,就可以先对每个点 u u u 记录 b s u , i bs_{u,i} b s u , i 表示 u u u 是否在第 i i i 个询问中,然后转移时跑拓扑排序,将 u u u 可达的点 t t t 都做 b s t ← b s t ∨ b s u bs_t\gets bs_t \vee bs_u b s t ← b s t ∨ b s u 表示 t t t 会贡献到询问 i i i ◊ \Large\color{red}\Diamond ◊ 。
现在有一个 n × q n\times q n × q 的 bitset 记录了 b s u , i bs_{u,i} b s u , i 表示 u u u 是否会贡献到询问 i i i ,若把其当成一个 01 矩阵,则第 i i i 个询问的答案就是第 i i i 列的和。考虑 诡异操作 的 trick,记录 c n t i , j cnt_{i,j} c n t i , j 表示第 i i i 列的和的第 j j j 个二进制位,然后将 c n t cnt c n t 转置后就变成了一个 log 2 n × q \log_2n \times q log 2 n × q 的矩阵,每次做对两个矩阵的每一行依次做二进制加法,单次复杂度为 O ( q log n ω ) \mathcal O(\frac{q\log n}\omega) O ( ω q log n ) ,做 n n n 次加法复杂度就是 O ( n q log n ω ) \mathcal O(\frac{nq\log n}{\omega}) O ( ω n q log n ) 。但是可以使用分治将 log n \log n log n 进一步减小,具体的,每次递归求解 [ l , m i d ] [l,mid] [ l , m i d ] 和 [ m i d + 1 , r ] [mid+1,r] [ m i d + 1 , r ] 的和,再做一次加法,复杂度 T ( n ) = 2 T ( n 2 ) + O ( q log n ω ) T(n)=2T(\frac n2)+\mathcal O(\frac{q\log n}\omega) T ( n ) = 2 T ( 2 n ) + O ( ω q log n ) ,由主定理可以得到 T ( n ) = O ( n q ω ) T(n)=\mathcal O(\frac{nq}\omega) T ( n ) = O ( ω n q ) ,虽然都是做 n − 1 n-1 n − 1 次加法,但分治可以将 log n \log n log n 不卡满,所以复杂度少一个 log \log log 。总复杂度为 O ( n q ω ) \mathcal O(\frac{nq}\omega) O ( ω n q ) ,需要使用第一段的卡空间方法,将 q q q 次询问分块处理,平衡后时间复杂度做到 O ( n q ω + q B ) \mathcal O(\frac{nq}\omega+qB) O ( ω n q + q B ) ,空间 O ( n B ω ) \mathcal O(\frac{nB}\omega) O ( ω n B ) 。