12.2
P7706 「Wdsr-2.7」文文的摄影布置 FST=1 \color{Blue}\text{P7706 「Wdsr-2.7」文文的摄影布置\ \ \ \ \ }\color{Red}\text{FST=1} P7706 「 Wdsr-2.7 」文文的摄影布置 FST=1
题目中的信息可以直接在区间上合并,于是用线段树维护即可。
12.3
P5278 算术天才9与等差数列 FST=19 \color{Purple}\text{P5278 算术天才9与等差数列\ \ \ \ \ }\color{Red}\text{FST=19} P5278 算术天才 9 与等差数列 FST=19
每次询问的等差数列是可以求出的,于是可以数列 Hash \text{Hash} Hash ,因为数据水,所以可以判最小最大+线性和+平方和水过。
12.5
被锁在机房外\kk
12.6
P9073 [WC/CTS2023] 楼梯 FST=1 \color{Black}\text{P9073 [WC/CTS2023] 楼梯\ \ \ \ \ }\color{Red}\text{FST=1} P9073 [WC/CTS2023] 楼梯 FST=1
神仙题,没看懂第一篇题解的特殊性质。
所形成的阶梯可能很大,所以考虑维护其边界,令竖直的边界为 1 1 1 ,水平的为 0 0 0 ,那么楼梯就变成了 01 01 01 序列,还有一个神奇的性质:任何一段连续的 1...0 1...0 1...0 都是一个子楼梯!接着发现每次询问的 q q q 很特殊,都是边界数 s u m − 1 sum-1 s u m − 1 的因子,显然不是用来根号分治的, 神仙的点在于,设边界为 B 1 ∼ s u m B_{1\sim sum} B 1 ∼ s u m ,单独拎出 B 1 , B 1 + q , B 1 + 2 q , . . . , B s u m B_1,B_{1+q},B_{1+2q},...,B_{sum} B 1 , B 1 + q , B 1 + 2 q , ... , B s u m ,若有相邻的 B 1 + i × q , B 1 + i × q + q B_{1+i\times q},B_{1+i\times q+q} B 1 + i × q , B 1 + i × q + q 为 1 , 0 1,0 1 , 0 ,那么就能找出一组解,而 B 1 = 1 , B s u m = 0 B_1=1,B_{sum}=0 B 1 = 1 , B s u m = 0 ,必然可以用二分找出一组解!
题目中有回退操作,可以用可持久化线段树维护,复杂度 O ( n log 2 n ) O(n\log^2n) O ( n log 2 n ) 。
P8557 炼金术(Alchemy) No FST \color{Yellow}\text{P8557 炼金术(Alchemy)\ \ \ \ \ }\color{Green}\text{No FST} P8557 炼金术( Alchemy ) No FST
开始水题。
分别考虑每种金属,方案数显然为 2 n − 1 2^n-1 2 n − 1 ,然后乘法原理即可。
P5303 [GXOI/GZOI2019] 逼死强迫症 FST=1 \color{Purple}\text{P5303 [GXOI/GZOI2019] 逼死强迫症\ \ \ \ \ }\color{Red}\text{FST=1} P5303 [GXOI/GZOI2019] 逼死强迫症 FST=1
记 f i f_i f i 为斐波那契第 i i i 项(f 0 = f 1 = 1 f_0=f_1=1 f 0 = f 1 = 1 ),s f i sf_i s f i 为 ∑ j = 0 i f j \sum_{j=0}^if_j ∑ j = 0 i f j ,那么答案为 ∑ i = 0 n − 3 f i × s f n − i − 3 \sum_{i=0}^{n-3}f_i\times sf_{n-i-3} ∑ i = 0 n − 3 f i × s f n − i − 3 ,经推式子可知能用距乘求解,用到的一个可能只有我不知道的结论就是 s f i = f i + 2 − 1 sf_i=f_{i+2}-1 s f i = f i + 2 − 1 。
P5239 回忆京都 No FST \color{Green}\text{P5239 回忆京都\ \ \ \ \ }\color{Green}\text{No FST} P5239 回忆京都 No FST
前缀和即可。
12.7
CF13C Sequence P4597 序列 sequence CF713C Sonya and Problem Wihtout a Legend No FST \color{Blue}\text{CF13C Sequence\ }\color{Purple}\text{P4597 序列 sequence\ CF713C Sonya and Problem Wihtout a Legend\ \ \ \ \ }\color{Green}\text{No FST} CF13C Sequence P4597 序列 sequence CF713C Sonya and Problem Wihtout a Legend No FST
感谢 @ xhgua @\color{Red}\text{xhgua} @ xhgua 送的经验。
先考虑单调不降,递增就只需再减去下标。
如果一段区间要修改为同一个数,显然要取中位数,如果有两个相邻区间中位数逆序就直接合并即可,要用可并堆。
12.9
P8576 「DTOI-2」星之界 FST=50+ \color{Purple}\text{P8576 「DTOI-2」星之界\ \ \ \ \ }\color{Red}\text{FST=50+} P8576 「 DTOI-2 」星之界 FST=50+
一眼切,但是调两天。
修改操作很套路,直接分块+并查集即可,对于询问,可以发现这就是多重集的排列,预处理阶乘即可。
但是还没有做完,这样只有 70 p t s 70pts 70 pt s ,需要时空都进一步优化:
并查集尽量少新建节点,能替代就替代。
可以预处理快速幂,不然会被卡。
P8360 [SNOI2022] 军队 By Mouth \color{Purple}\text{P8360 [SNOI2022] 军队\ \ \ \ \ }\color{Blue}\text{By Mouth} P8360 [SNOI2022] 军队 By Mouth
做法同上题,只是维护的标记变了。
12.10
P9361 [ICPC2022 Xi’an R] Contests FST=1 \color{Purple}\text{P9361 [ICPC2022 Xi'an R] Contests\ \ \ \ \ }\color{Red}\text{FST=1} P9361 [ICPC2022 Xi’an R] Contests FST=1
如果某次在第 i i i 行跳并且下次在第 j j j 行跳的话,显然要跳可以跳到的在第 j j j 行排名最小的才最优,并且对于一个选手,他能调到的第 j j j 行的最优点是确定的(总共就 m m m 种情况),于是设 f t , k , j f_{t,k,j} f t , k , j 表示选手 j j j 跳 2 t 2^t 2 t 步到第 k k k 行的最优点,询问只要倍增查询即可。
CF1511G Chips on a Board No FST \color{Purple}\text{CF1511G Chips on a Board\ \ \ \ \ }\color{Green}\text{No FST} CF1511G Chips on a Board No FST
显然询问是个 Nim \text{Nim} Nim 博弈,我们只需求出每个棋子能走的步数的异或和即可。
根据标签可知这个题关于 2 2 2 的次幂的合并有优美的性质,例如合并两个长度为 2 x , y 2^x,y 2 x , y 的区间(y ≤ 2 x y\le 2^x y ≤ 2 x ),右区间增加的 2 x 2^x 2 x 贡献只与右区间的点数有关,很方便维护,于是倍增即可。
P3643 [APIO2016] 划艇 No FST \color{Purple}\text{P3643 [APIO2016] 划艇\ \ \ \ \ }\color{Green}\text{No FST} P3643 [APIO2016] 划艇 No FST
想写拉插,结果写成了这个题。
直接 d p dp d p ,设 f i , j f_{i,j} f i , j 表示第 i i i 个学校出 j j j 个划艇的方案数,但这样的话时空均爆炸,考虑将所有的 [ a i , b i + 1 ) [a_i,b_i+1) [ a i , b i + 1 ) 离散化为 O ( n ) O(n) O ( n ) 个区间,将状态改为第 i i i 个学校出的游艇数在第 j j j 个区间内的方案数,那么 f i , j f_{i,j} f i , j 就会等于 f p , j − 1 f_{p,j-1} f p , j − 1 乘上一些学校艇数在第 j j j 个区间里的方案数,显然可以算出 p + 1 ∼ i p+1\sim i p + 1 ∼ i 中艇数能在第 j j j 个区间的学校数 m m m ,问题就转化为了在 l e n len l e n 个数中选不超过 m m m 个数的方案数。
有一个巧妙的构造,即该问题等价于在 m m m 个 0 0 0 和 1 ∼ l e n 1\sim len 1 ∼ l e n 组成的序列(即 0 , 0 , 0 , . . . , 0 , 0 , 1 , 2 , 3 , . . . , l e n 0,0,0,...,0,0,1,2,3,...,len 0 , 0 , 0 , ... , 0 , 0 , 1 , 2 , 3 , ... , l e n )中选 m m m 个数的方案数,选择第 i i i 个 0 0 0 等价于第 i i i 个学校不参与,于是答案为 C m + l e n m C_{m+len}^m C m + l e n m 。但是本题强制第 i i i 个学校选,所以应为 C m + l e n − 1 m C_{m+len-1}^m C m + l e n − 1 m 。
如此 d p dp d p 即可做到 O ( n 3 ) O(n^3) O ( n 3 ) 。
12.12
P8290 [省选联考 2022] 填树 WA&TLE \color{Purple}\text{P8290 [省选联考 2022] 填树\ \ \ \ \ }\color{Black}\text{WA\&TLE} P8290 [ 省选联考 2022] 填树 WA&TLE
调两天终究没调出来的题。
对于第一问,显然可以枚举取值区间 [ x , x + k ] [x,x+k] [ x , x + k ] ,但要求必须有最小值,于是可以容斥,转化为 w ( x , x + k ) − w ( x + 1 , x + k ) w(x,x+k)-w(x+1,x+k) w ( x , x + k ) − w ( x + 1 , x + k ) 的形式,于是可以树上换根 d p dp d p 做到 O ( n w ) O(nw) O ( n w ) ,但值域过大,可以发现 x x x 在一定范围内时,每个点的种类数可以表示为一次函数,而这样的范围离散化后只有 O ( n ) O(n) O ( n ) 种,于是扫描线,可以求出对应的多项式再用拉格朗日插值优化求出。
对于第二问也同理。
12.13
P3306 [SDOI2013] 随机数生成器 FST=6 \color{Blue}\text{P3306 [SDOI2013] 随机数生成器\ \ \ \ \ }\color{Red}\text{FST=6} P3306 [SDOI2013] 随机数生成器 FST=6
关于 b b b 的项是一个等比数列,于是化简可得:a x − 1 ≡ a t − a + b x 1 t − x 1 + b ( m o d p ) a^{x-1}\equiv \frac{at-a+b}{x_1t-x_1+b}\pmod p a x − 1 ≡ x 1 t − x 1 + b a t − a + b ( mod p ) ,直接 BSGS \text{BSGS} BSGS 即可。
CF1106F Lunar New Year and a Recursive Sequence No FST \color{Purple}\text{CF1106F Lunar New Year and a Recursive Sequence\ \ \ \ \ }\color{Green}\text{No FST} CF1106F Lunar New Year and a Recursive Sequence No FST
乘积又乘方不好处理,这种时候不妨取离散对数,以 998244353 998244353 998244353 的原根 3 3 3 为底数,就能出现一个齐次线性递推的柿子,另 f k = g 1 f_k=g^1 f k = g 1 ,那么可以距乘求出 f n = g a f_n=g^a f n = g a (注意此时模数为 998244352 998244352 998244352 ),故 f k a ≡ m ( m o d 998244353 ) f_k^a\equiv m\pmod {998244353} f k a ≡ m ( mod 998244353 ) ,仍使用原根,设 f k = g x , m = g b f_k=g^x,m=g^b f k = g x , m = g b ,那么 b b b 可通过 BSGS \text{BSGS} BSGS 求出,最终转化为 a x ≡ b ( m o d 998244352 ) ax\equiv b\pmod{998244352} a x ≡ b ( mod 998244352 ) ,利用 exgcd \text{exgcd} exgcd 求解即可。
可以发现,原根的好处就是可以把一个未知数表示为 g x g^x g x ,实现降幂,方便求解。
12.15
P3723 [AH2017/HNOI2017] 礼物 FST=1 \color{Purple}\text{P3723 [AH2017/HNOI2017] 礼物\ \ \ \ \ }\color{Red}\text{FST=1} P3723 [AH2017/HNOI2017] 礼物 FST=1
c c c 取非负与取任意整数是等价的,于是转化为求 ∑ i = 1 n ( a i + c − b i ) 2 = ∑ i = 1 n a i 2 + ∑ i = 1 n b i 2 + n c 2 + 2 c ∑ i = 1 n ( a i − b i ) − 2 ∑ i = 1 n a i b i \sum_{i=1}^n(a_i+c-b_i)^2=\sum_{i=1}^na_i^2+\sum_{i=1}^nb_i^2+nc^2+2c\sum_{i=1}^n(a_i-b_i)-2\sum_{i=1}^na_ib_i ∑ i = 1 n ( a i + c − b i ) 2 = ∑ i = 1 n a i 2 + ∑ i = 1 n b i 2 + n c 2 + 2 c ∑ i = 1 n ( a i − b i ) − 2 ∑ i = 1 n a i b i 。
前面的部分可以取对称轴,问题转化为求 ∑ i = 1 n a i b i \sum_{i=1}^na_ib_i ∑ i = 1 n a i b i 的最大值,将 b b b 倒序后就变成 ∑ i = 1 n a i b n + 1 − i \sum_{i=1}^na_ib_{n+1-i} ∑ i = 1 n a i b n + 1 − i ,是一个卷积的形式,于是可以将 b b b 复制两倍,将多项式 a , b a,b a , b 相乘,取 n + 1 ∼ 2 n n+1\sim2n n + 1 ∼ 2 n 项系数的最大值即可。
12.19
P6406 [COCI2014-2015#2] Norma FST=2 \color{Purple}\text{P6406 [COCI2014-2015\#2] Norma\ \ \ \ \ }\color{Red}\text{FST=2} P6406 [COCI2014-2015#2] Norma FST=2
数据结构不好维护这个柿子,发现没有修改,于是考虑 cdq \text{cdq} cdq 分治,对于跨过 m i d mid mi d 的部分,枚举左端点,右端点可按照贡献种类分为 3 3 3 段,分别讨论即可。
12.20
P3350 [ZJOI2016] 旅行者 FST=8 \color{Purple}\text{P3350 [ZJOI2016] 旅行者\ \ \ \ \ }\color{Red}\text{FST=8} P3350 [ZJOI2016] 旅行者 FST=8
对于每个询问,我们可以选择一个行 s s s ,求出 s s s 上的每个点到询问点的最短路,再合并就能求出一个可能答案,对于两点在 s s s 异侧的点,显然这就是答案,否则在两侧递归处理即可。
P8231 [AGM 2022 资格赛] 农场 No FST \color{Purple}\text{P8231 [AGM 2022 资格赛] 农场\ \ \ \ \ }\color{Green}\text{No FST} P8231 [AGM 2022 资格赛 ] 农场 No FST
整体分治板题。
CF1217F Forced Online Queries Problem WA&Fake \color{Purple}\text{CF1217F Forced Online Queries Problem\ \ \ \ \ }\color{Black}\text{WA\&Fake} CF1217F Forced Online Queries Problem WA&Fake
如果没有强制在线,直接线段树分治+可撤销并查集即可。
如果有强制在线,第一篇题解告诉我不影响答案,于是写了并且假了。
注意起影响作用的 l s t ∈ { 0 , 1 } lst\in\{0,1\} l s t ∈ { 0 , 1 } ,所以只有两类情况,一条边的影响可以分解为若干区间来做。明天月考,就不补了。
P1641 [SCOI2010] 生成字符串 No FST \color{Blue}\text{P1641 [SCOI2010] 生成字符串\ \ \ \ \ }\color{Green}\text{No FST} P1641 [SCOI2010] 生成字符串 No FST
一个经典的模型转化:从 ( 0 , 0 ) (0,0) ( 0 , 0 ) 出发,每次向上或向右走,不穿过 y = x y=x y = x 的直线到 ( n , m ) (n,m) ( n , m ) 的方案数。
不考虑限制,方案数为 C ( n + m , m ) C(n+m,m) C ( n + m , m ) 。
对于经过直线 y = x y=x y = x 的方案,我们考虑第一步直接走到 ( 0 , 1 ) (0,1) ( 0 , 1 ) ,那么到 ( n , m ) (n,m) ( n , m ) 就必定经过 y = x y=x y = x ,如果不是一开始就走 ( 0 , 1 ) (0,1) ( 0 , 1 ) 的方案,可以翻转图案得到这个结果,这样从 ( 0 , 1 ) (0,1) ( 0 , 1 ) 出发的方案就对应了所有不合法的情况,此时情况数为 C ( n + m , m − 1 ) C(n+m,m-1) C ( n + m , m − 1 ) 。
二者相减即为答案。
12.24
P3206 [HNOI2010] 城市建设 FST=6 \color{Purple}\text{P3206 [HNOI2010] 城市建设\ \ \ \ \ }\color{Red}\text{FST=6} P3206 [HNOI2010] 城市建设 FST=6
神仙抽象题。
考虑分治,也就是说我们要使处理一个区间的复杂度在 O ( l e n ) O(len) O ( l e n ) 级别,当分治到 l = r l=r l = r 的区间时,我们希望除了第 l l l 次修改的边其他边都能被更新好。
对于区间 [ l , r ] [l,r] [ l , r ] ,不妨设 [ l , r ] [l,r] [ l , r ] 的修改牵扯到的边为动态边,其余的为静态边。
先强制选择所有动态边,再跑 Kruskal \text{Kruskal} Kruskal ,则静态边可分为两类:必选边与可选可不选边。
把必选边都计入答案,并将此时的图(不包含动态边)缩点,这样就保证了点数在 l e n len l e n 级别。
强制不选择所有动态边,再跑 Kruskal \text{Kruskal} Kruskal ,则可选可不选边中会出现一类无用边,直接删去即可。这样也保证的边数在 O ( l e n ) O(len) O ( l e n ) 级别。
这样就满足了上面的要求。只需要在点数 = 1 =1 = 1 或 l = r l=r l = r 时特判即可。
P8321 『JROI-4』沈阳大街 2 No FST \color{Purple}\text{P8321 『JROI-4』沈阳大街 2\ \ \ \ \ }\color{Green}\text{No FST} P8321 『 JROI-4 』沈阳大街 2 No FST
把 a , b a,b a , b 两数列合并,从大到小排序,属于 a a a 的染白, b b b 的染黑,问题转化为求将一个数组中白黑元素两两配对,一对的权值为靠后的元素的权值,一种方案的权值为 n n n 对权值之积,求所有方案权值之和。
直接 d p dp d p 即可。
P8688 [蓝桥杯 2019 省 A] 组合数问题 P6669 [清华集训2016] 组合数问题 FST=3 \color{Purple}\text{P8688 [蓝桥杯 2019 省 A] 组合数问题\ P6669 [清华集训2016] 组合数问题\ \ \ \ \ }\color{Red}\text{FST=3} P8688 [ 蓝桥杯 2019 省 A] 组合数问题 P6669 [ 清华集训 2016] 组合数问题 FST=3
用卢卡斯定理,可将 ( n m ) \dbinom{n}{m} ( m n ) 拆解为若干 ( x y ) ( x < k , y < k ) \dbinom{x}{y}(x<k,y<k) ( y x ) ( x < k , y < k ) 的乘积,( n m ) \dbinom{n}{m} ( m n ) 为 0 0 0 当且仅当存在 ( x y ) \dbinom{x}{y} ( y x ) 且 x < y x<y x < y ,考虑容斥,用总方案数减去全部满足 x ≥ y x\ge y x ≥ y 的方案数,直接数位 d p dp d p 即可。
P8594 「KDOI-02」一个仇的复 FST=1 \color{Purple}\text{P8594 「KDOI-02」一个仇的复\ \ \ \ \ }\color{Red}\text{FST=1} P8594 「 KDOI-02 」一个仇的复 FST=1
先从简单的情况考虑,如果要用 k k k 个矩阵填满一个 1 × n 1\times n 1 × n 的网络,那么可以令第 x x x 个元素选的意义为方案中存在以 x x x 为左端点的矩阵,第一个元素是必选的,故方案数为 ( n − 1 k − 1 ) \dbinom{n-1}{k-1} ( k − 1 n − 1 ) ,两行的方案数为 ∑ i ( n − 1 i − 1 ) ( n − 1 k − i − 1 ) \sum_{i}\dbinom{n-1}{i-1}\dbinom{n-1}{k-i-1} ∑ i ( i − 1 n − 1 ) ( k − i − 1 n − 1 ) ,利用范德蒙德卷积可得这可以化简为 ( 2 n − 2 k − 2 ) \dbinom{2n-2}{k-2} ( k − 2 2 n − 2 ) 。
对于题目,我们假设两个并排的 1 × x 1\times x 1 × x 矩阵有 i i i 段,2 × 1 2\times 1 2 × 1 的有 j j j 个,那么:
1 × x 1\times x 1 × x 矩阵的方案数可以用上述方法合并为 ( 2 n − 2 i − 2 j k − i − 2 j ) \dbinom{2n-2i-2j}{k-i-2j} ( k − i − 2 j 2 n − 2 i − 2 j ) ,但由于要分成 i i i 段,还要乘上 ( n − j − 1 i − 1 ) \dbinom{n-j-1}{i-1} ( i − 1 n − j − 1 ) 。
再考虑把 j j j 个 2 × 1 2\times1 2 × 1 矩阵插入 i + 1 i+1 i + 1 个空位(首尾还有 2 2 2 个)中,其中首尾的可以不放,中间的至少要放一个,方案数为 ( j + 1 i ) \dbinom{j+1}{i} ( i j + 1 ) 。
最终答案为 ∑ i = 1 k ∑ j = 0 k ( 2 n − 2 i − 2 j k − i − 2 j ) ( n − j − 1 i − 1 ) ( j + 1 i ) + [ n = k ] \sum_{i=1}^k\sum_{j=0}^k\dbinom{2n-2i-2j}{k-i-2j}\dbinom{n-j-1}{i-1}\dbinom{j+1}{i}+[n=k] ∑ i = 1 k ∑ j = 0 k ( k − i − 2 j 2 n − 2 i − 2 j ) ( i − 1 n − j − 1 ) ( i j + 1 ) + [ n = k ] 。
12.25
P8367 [LNOI2022] 盒 FST=50+ \color{Black}\text{P8367 [LNOI2022] 盒\ \ \ \ \ }\color{Red}\text{FST=50+} P8367 [LNOI2022] 盒 FST=50+
神仙题,写了个假做法调了一节课。
先把几个 trick \text{trick} trick 写在前面:
j ( i + j − 1 j ) = i ( i + j − 1 i ) j\dbinom{i+j-1}{j}=i\dbinom{i+j-1}{i} j ( j i + j − 1 ) = i ( i i + j − 1 ) 。展开即可证得。
( n + m n ) = ∑ j = 0 m ( i + j j ) ( n − i − 1 + m − j m − j ) \dbinom{n+m}{n}=\sum_{j=0}^m\dbinom{i+j}{j}\dbinom{n-i-1+m-j}{m-j} ( n n + m ) = ∑ j = 0 m ( j i + j ) ( m − j n − i − 1 + m − j ) 。考虑组合意义:从 ( 0 , 0 ) (0,0) ( 0 , 0 ) 走到 ( n , m ) (n,m) ( n , m ) 的方案数相当于走到第 i i i 行并从 ( i , j ) (i,j) ( i , j ) 走到 ( i + 1 , j ) (i+1,j) ( i + 1 , j ) 的方案数。
∑ i = 0 q ( p + i i ) ( n − p − 1 + m − i m − i ) = ∑ i = p + 1 n ( i + q i ) ( n − i + m − q − 1 n − i ) \sum_{i=0}^q\dbinom{p+i}{i}\dbinom{n-p-1+m-i}{m-i}=\sum_{i=p+1}^n\dbinom{i+q}{i}\dbinom{n-i+m-q-1}{n-i} ∑ i = 0 q ( i p + i ) ( m − i n − p − 1 + m − i ) = ∑ i = p + 1 n ( i i + q ) ( n − i n − i + m − q − 1 ) 。考虑组合意义:前者相当于从 ( 0 , 0 ) (0,0) ( 0 , 0 ) 走到 ( n , m ) (n,m) ( n , m ) 时,从 ( p , i ) (p,i) ( p , i ) 到 ( p + 1 , i ) (p+1,i) ( p + 1 , i ) 的 i i i 不能超过 q q q ,即第 p p p 行到 p + 1 p+1 p + 1 行时列不超过 q q q 的方案数;后者相当于从第 q q q 列到第 q + 1 q+1 q + 1 列时行已超过 p p p 的方案数,二者显然是等价的。
i i i 个非负元素相加等于 j j j 的方案数:( i + j − 1 i − 1 ) \dbinom{i+j-1}{i-1} ( i − 1 i + j − 1 ) ,由隔板法可得。
令 c i c_i c i 为 a i a_i a i 的前缀和,那么有:
a n s = ∑ i = 1 n − 1 w i ∑ j = 0 S ∣ c i − j ∣ ( i + j − 1 i − 1 ) ( n − i − 1 + S − j n − i − 1 ) ans=\sum_{i=1}^{n-1}w_i\sum_{j=0}^S|c_i-j|\dbinom{i+j-1}{i-1}\dbinom{n-i-1+S-j}{n-i-1} an s = i = 1 ∑ n − 1 w i j = 0 ∑ S ∣ c i − j ∣ ( i − 1 i + j − 1 ) ( n − i − 1 n − i − 1 + S − j )
令 r e s i = ∑ j = 0 S ∣ c i − j ∣ ( i + j − 1 i − 1 ) ( n − i − 1 + S − j n − i − 1 ) {res}_i=\sum_{j=0}^S|c_i-j|\dbinom{i+j-1}{i-1}\dbinom{n-i-1+S-j}{n-i-1} res i = ∑ j = 0 S ∣ c i − j ∣ ( i − 1 i + j − 1 ) ( n − i − 1 n − i − 1 + S − j ) ,拆去不等号并化简,得到从 0 0 0 开始的形式:
∑ j = 0 S ( j − c i ) ( i + j − 1 j ) ( n − i − 1 + S − j S − j ) + 2 ∑ j = 0 c i ( c i − j ) ( i + j − 1 j ) ( n − i − 1 + S − j S − j ) \sum_{j=0}^S(j-c_i)\dbinom{i+j-1}{j}\dbinom{n-i-1+S-j}{S-j}+2\sum_{j=0}^{c_i}(c_i-j)\dbinom{i+j-1}{j}\dbinom{n-i-1+S-j}{S-j} j = 0 ∑ S ( j − c i ) ( j i + j − 1 ) ( S − j n − i − 1 + S − j ) + 2 j = 0 ∑ c i ( c i − j ) ( j i + j − 1 ) ( S − j n − i − 1 + S − j )
先拆开左式,右式同理解决:
∑ j = 0 S j ( i + j − 1 j ) ( n − i − 1 + S − j S − j ) − c i ∑ j = 0 S ( i + j − 1 j ) ( n − i − 1 + S − j S − j ) \sum_{j=0}^Sj\dbinom{i+j-1}{j}\dbinom{n-i-1+S-j}{S-j}-c_i\sum_{j=0}^S\dbinom{i+j-1}{j}\dbinom{n-i-1+S-j}{S-j} j = 0 ∑ S j ( j i + j − 1 ) ( S − j n − i − 1 + S − j ) − c i j = 0 ∑ S ( j i + j − 1 ) ( S − j n − i − 1 + S − j )
利用 trick \text{trick} trick 可得:
i ∑ j = 0 S ( i + j − 1 j − 1 ) ( n − i − 1 + S − j S − j ) − c i ∑ j = 0 S ( i + j − 1 j ) ( n − i − 1 + S − j S − j ) i\sum_{j=0}^S\dbinom{i+j-1}{j-1}\dbinom{n-i-1+S-j}{S-j}-c_i\sum_{j=0}^S\dbinom{i+j-1}{j}\dbinom{n-i-1+S-j}{S-j} i j = 0 ∑ S ( j − 1 i + j − 1 ) ( S − j n − i − 1 + S − j ) − c i j = 0 ∑ S ( j i + j − 1 ) ( S − j n − i − 1 + S − j )
为了使左右两边的形式相同,我们把左式的上界向下挪一位:
i ∑ j = 0 S − 1 ( i + j j ) ( n − i − 1 + S − j + 1 S − j + 1 ) − c i ∑ j = 0 S ( i + j − 1 j ) ( n − i − 1 + S − j S − j ) i\sum_{j=0}^{S-1}\dbinom{i+j}{j}\dbinom{n-i-1+S-j+1}{S-j+1}-c_i\sum_{j=0}^S\dbinom{i+j-1}{j}\dbinom{n-i-1+S-j}{S-j} i j = 0 ∑ S − 1 ( j i + j ) ( S − j + 1 n − i − 1 + S − j + 1 ) − c i j = 0 ∑ S ( j i + j − 1 ) ( S − j n − i − 1 + S − j )
i ∑ j = 0 S − 1 ( ( i + 1 ) + j − 1 ( i + 1 ) − 1 ) ( ( n + 1 ) − ( i + 1 ) − 1 + S + 1 − j ( n + 1 ) − ( i + 1 ) − 1 ) − c i ∑ j = 0 S ( i + j − 1 i − 1 ) ( n − i − 1 + S − j n − i − 1 ) i\sum_{j=0}^{S-1}\dbinom{(i+1)+j-1}{(i+1)-1}\dbinom{(n+1)-(i+1)-1+S+1-j}{(n+1)-(i+1)-1}-c_i\sum_{j=0}^S\dbinom{i+j-1}{i-1}\dbinom{n-i-1+S-j}{n-i-1} i j = 0 ∑ S − 1 ( ( i + 1 ) − 1 ( i + 1 ) + j − 1 ) ( ( n + 1 ) − ( i + 1 ) − 1 ( n + 1 ) − ( i + 1 ) − 1 + S + 1 − j ) − c i j = 0 ∑ S ( i − 1 i + j − 1 ) ( n − i − 1 n − i − 1 + S − j )
令 f ( n , s , i , k ) = ∑ j = 0 k ( i + j − 1 i − 1 ) ( n − i − 1 + s − j n − i − 1 ) f(n,s,i,k)=\sum_{j=0}^k\dbinom{i+j-1}{i-1}\dbinom{n-i-1+s-j}{n-i-1} f ( n , s , i , k ) = ∑ j = 0 k ( i − 1 i + j − 1 ) ( n − i − 1 n − i − 1 + s − j ) ,那么上式等于 i × f ( n + 1 , S − 1 , i + 1 , S − 1 ) − c i × f ( n , S , i , S ) i\times f(n+1,S-1,i+1,S-1)-c_i\times f(n,S,i,S) i × f ( n + 1 , S − 1 , i + 1 , S − 1 ) − c i × f ( n , S , i , S ) ,另一式也同理即可得出:
r e s i = i × f ( n + 1 , S − 1 , i + 1 , S − 1 ) − s i × f ( n , S , i , S ) + 2 × f ( n , S , i , s i ) − 2 i × f ( n + 1 , S − 1 , i + 1 , s i − 1 ) {res}_i=i\times f(n+1,S-1,i+1,S-1)-s_i\times f(n,S,i,S)+2\times f(n,S,i,s_i)-2i\times f(n+1,S-1,i+1,s_i-1) res i = i × f ( n + 1 , S − 1 , i + 1 , S − 1 ) − s i × f ( n , S , i , S ) + 2 × f ( n , S , i , s i ) − 2 i × f ( n + 1 , S − 1 , i + 1 , s i − 1 )
虽然单次计算是 O ( k ) O(k) O ( k ) 的,但是根据公式,可以由 f ( n , s , i , k ) O ( 1 ) f(n,s,i,k)\ \text{O}(1) f ( n , s , i , k ) O ( 1 ) 转移到 f ( n , s , i , k + 1 ) f(n,s,i,k+1) f ( n , s , i , k + 1 ) ,可 i i i 的转移就不好办了,这个时候就可以用上 trick \text{trick} trick :
∑ j = 0 k ( i + j − 1 i − 1 ) ( n − i − 1 + s − j n − i − 1 ) = ∑ j = i + 1 n ( j + k − 1 k ) ( n − j + s − k − 1 n − j ) \sum_{j=0}^k\dbinom{i+j-1}{i-1}\dbinom{n-i-1+s-j}{n-i-1}=\sum_{j=i+1}^n\dbinom{j+k-1}{k}\dbinom{n-j+s-k-1}{n-j} j = 0 ∑ k ( i − 1 i + j − 1 ) ( n − i − 1 n − i − 1 + s − j ) = j = i + 1 ∑ n ( k j + k − 1 ) ( n − j n − j + s − k − 1 )
就也可以 O ( 1 ) \text{O}(1) O ( 1 ) 实现 i i i 到 i + 1 i+1 i + 1 的转移了。
12.26
P4562 [JXOI2018] 游戏 FST=4 \color{Blue}\text{P4562 [JXOI2018] 游戏\ \ \ \ \ }\color{Red}\text{FST=4} P4562 [JXOI2018] 游戏 FST=4
显然如果选择了 x x x ,其因数不会被覆盖,那么选 x x x 就是无用的,不妨找出最大的因数(不包括自己)不在 [ l , r ] [l,r] [ l , r ] 的数,特别地, 1 1 1 也算,记总个数为 c n t cnt c n t ,记这些数为关键数,那么 t ( p ) t(p) t ( p ) 为最后一个关键数的位置,直接计数即可。
P2606 [ZJOI2010] 排列计数 FST=3 \color{Blue}\text{P2606 [ZJOI2010] 排列计数\ \ \ \ \ }\color{Red}\text{FST=3} P2606 [ZJOI2010] 排列计数 FST=3
将 ⌊ p 2 ⌋ \left\lfloor\frac{p}{2}\right\rfloor ⌊ 2 p ⌋ 向 p p p 连边,就构成了一颗以 1 1 1 为根的二叉树,在树上 d p dp d p 即可。
然后 WA 100 pts \text{WA 100 pts} WA 100 pts ,被 hack \text{hack} hack 了。
n n n 可能超过模数,逆元就不存在,这时套路的另开一个数组存模数的次幂即可。
12.27
P8863 「KDOI-03」构造数组 No FST \color{Purple}\text{P8863 「KDOI-03」构造数组\ \ \ \ \ }\color{Green}\text{No FST} P8863 「 KDOI-03 」构造数组 No FST
记 ∑ i = 1 n b i = s u m , ∑ j = 1 i b j = a i \sum_{i=1}^nb_i=sum,\sum_{j=1}^ib_j=a_i ∑ i = 1 n b i = s u m , ∑ j = 1 i b j = a i ,若 s u m ≡ 1 ( m o d 2 ) sum\equiv1\pmod 2 s u m ≡ 1 ( mod 2 ) 那么无解。
否则,问题等价于有 n n n 种元素,每种有 b i b_i b i 个,组成 s u m 2 \frac{sum}{2} 2 s u m 个无序数对且同种元素不能在一对的方案数的方案数。
按编号依次考虑每种元素,这样就可以不用考虑无序关系。
对于第 i i i 种元素,前面的元素占用配对的情况一定只有三种:2 2 2 个,1 1 1 个,0 0 0 个。记其个数分别为 k k , k , j kk,k,j kk , k , j ,那么有:
k k + k + j = s u m 2 kk+k+j=\frac{sum}{2} kk + k + j = 2 s u m
2 k k + k = a i 2kk+k=a_i 2 kk + k = a i
所以 k , j k,j k , j 一定都可以用 k k kk kk 表示,故设 f i , k k f_{i,kk} f i , kk 为第 i i i 种元素放完后有 k k kk kk 个满配对的方案数,那么:
f i , k k = ∑ p = 0 max { k k , b i } f i − 1 , k k − p ( j + b i − p b i − p ) ( k + p − ( b i − p ) p ) f_{i,kk}=\sum_{p=0}^{\max\{kk,b_i\}}f_{i-1,kk-p}\dbinom{j+b_i-p}{b_i-p}\dbinom{k+p-(b_i-p)}{p} f i , kk = p = 0 ∑ m a x { kk , b i } f i − 1 , kk − p ( b i − p j + b i − p ) ( p k + p − ( b i − p ) )
时间复杂度 O ( s u m 2 ) O(sum^2) O ( s u m 2 ) ,但是能过(?
P4456 [CQOI2018] 交错序列 FST=8 \color{Purple}\text{P4456 [CQOI2018] 交错序列\ \ \ \ \ }\color{Red}\text{FST=8} P4456 [CQOI2018] 交错序列 FST=8
显然 x + y = n x+y=n x + y = n ,所以一种情况的贡献值为 ( n − y ) a y b (n-y)^ay^b ( n − y ) a y b ,展开后可得:
x a y b = ∑ i = 0 a ( a i ) n i ( − 1 ) a − i y a + b − i x^ay^b=\sum_{i=0}^a\dbinom{a}{i}n^i(-1)^{a-i}y^{a+b-i} x a y b = i = 0 ∑ a ( i a ) n i ( − 1 ) a − i y a + b − i
( a i ) n i ( − 1 ) a − i \dbinom{a}{i}n^i(-1)^{a-i} ( i a ) n i ( − 1 ) a − i 这部分显然可以提出来,故只需算出 y a + b − i y^{a+b-i} y a + b − i 的系数即可。
设 f i , j , 0 / 1 f_{i,j,0/1} f i , j , 0/1 表示第 i i i 位为 0 / 1 0/1 0/1 的 y a + b − i y^{a+b-i} y a + b − i 项的系数,那么:
f i , j , 0 = f i − 1 , j , 0 + f i − 1 , j , 1 f_{i,j,0}=f_{i-1,j,0}+f_{i-1,j,1} f i , j , 0 = f i − 1 , j , 0 + f i − 1 , j , 1
f i , j , 1 = ∑ k = 0 j ( j k ) f i − 1 , k , 0 f_{i,j,1}=\sum_{k=0}^{j}\dbinom{j}{k}f_{i-1,k,0} f i , j , 1 = k = 0 ∑ j ( k j ) f i − 1 , k , 0
矩阵快速幂即可。
然后获得了比暴力多 10 10 10 分的好成绩
卡常技巧:
矩乘时先乘后统一取模。
在交换循环顺序的情况下,进入第三层循环前判断 a [ i ] [ k ] a[i][k] a [ i ] [ k ] 是否非 0 0 0 。
P5163 WD与地图 By mouth \color{Black}\text{P5163 WD与地图\ \ \ \ \ }\color{Blue}\text{By mouth} P5163 WD 与地图 By mouth
显然要先倒序操作,把删边变为加边。
如果边是无向的,只要动态合并连通分量即可,可以用合并权值线段树解决。
但这题边是有向的,可以发现发生合并的时间就是某条边第一次在一个连通分量内的时间,并且这个时间是可二分的,于是整体二分即可。
具体地,用可撤销并查集维护联通情况,跑只与新加的边有关的点的 tarjan \text{tarjan} tarjan 。
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步