ACGO 第7次排位赛题解
写在前面
本次排位赛D题 火星背包\tt{火星背包}火星背包 测试数据已加强。
非常感谢大家参加本场排位赛！
为了大家更好的参赛体验，本场比赛：
1. 所有题目采用全新的排版风格；
新的排版风格本着更方便参赛选手审题的原则，对“时空限制”，“题目内容”，“数据范围”，“输入输出格式”，等做了重新排版。
* 每道题目，均将时间限制与内存限制放置于题目开头。
* 题目内容中不再在变量后标注数据范围。
* 统一将题目中出现的所有变量的数据范围，归置于 数据范围\tt{数据范围}数据范围 板块，在题目内容下。
* 对于 多组测试用例\tt{多组测试用例}多组测试用例 的题目在 数据范围\tt{数据范围}数据范围 板块上方作强调说明。
* 输入格式部分，不再使用生硬，不易理解的文字来说明，而是使用更为直接的 数据样式示例\tt{数据样式示例}数据样式示例 说明。
2. 所有题目判题采用 CHECKER(SPECIAL JUDGE)\TT{CHECKER(SPECIAL\ JUDGE)}CHECKER(SPECIAL JUDGE)
* 对于存在多组合法方案的题目，不使用 checker 的题目，往往要求按照一定的顺序（例如字典序）输出，这往往会给解题带来一些不必要的麻烦，亦或是限制了参赛者的解题思路。而 checker 判题，则给了解题代码更大的包容性，输出任意顺序或者满足要求的答案即可。
* 对于答案为浮点数的题目，要求保留若干位小数输出，在不使用 checker 判题的情况下，往往可能会导致一些不必要的精度问题，
* 对于输出 Yes\tt{Yes}Yes 或者 No\tt{No}No 此类型的题目，使用 cheker 可以忽略解题代码输出答案的大小写，即：你可以选择你喜欢的形式，输出 Yes\tt{Yes}Yes 或者 No\tt{No}No。
3. 所有题目给出样例说明，并尽可能给出较好的样例
* 我们希望题目的 样例\tt{样例}样例 和 样例说明\tt{样例说明}样例说明 拥有让参赛者能够更好理解题目意思的作用，而不是仅仅因为需要有一个样例而存在。（比如本场比赛的 AAA 题。）
* 我们希望题目给出的样例能够在解题程序出现很大错误的时候是无法得到正确答案的，而不是在提交后发现 答案错误\tt{答案错误}答案错误。（比如本场比赛的 DDD 题的样例2，如果你采用一般的贪心的错误解法，是无法通过这个样例的。）
* 在不透露题目解题思路的情况下，尽可能给出一个数据略大的样例。（比如本场比赛的 B,D,FB, D, FB,D,F 题。）
4. 简洁清晰的题面
* 选手应该具备一定的阅读理解能力，但是我们认为竞赛重心仍然是放在解题上，让每位选手能够读清楚，读懂题目，然后思考问题。而不是故意编造一些和题目关联性很低的故事，或者是在题目中玩 文字游戏\tt{文字游戏}文字游戏 埋坑，或者使用一些可能会存在歧义的描述。
* 对于题目中提到的特殊概念如 好数组\tt{好数组}好数组，好数字\tt{好数字}好数字 等，做特殊标注并做相应解释；对于题目中的重点信息，使用加粗标注。
总而言之，我们希望大家更好地享受算法竞赛，减少一些其他原因影响解题的情况。
大家如果有其他的建议，欢迎在评论区提出哦～
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
简易题解
详细题解的AC代码包括复杂度分析，请点击题目标题查看。
A - 统计线段上格点的数量
令 X=∣xA−xB∣,Y=∣yA−yB∣X = \vert x_A - x_B \vert, Y = \vert y_A - y_B \vertX=∣xA −xB ∣,Y=∣yA −yB ∣，那么我们可以将题目转化为：求最多可以将 XXX 和 YYY 均分为多少段，使得每段的长度均为整数。
令均分的段数为 kkk，那么我们需要求的结果就是找到一个可以整除 XXX 和 YYY 的最大的 kkk。
也就是 XXX 和 YYY 的最大公约数，即 gcd(X,Y)\gcd(X, Y)gcd(X,Y)。
答案即为 gcd(X,Y)+1\gcd(X, Y) + 1gcd(X,Y)+1。
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
B - 大数求和
对于题目中的式子，我们有 ∑k=1Nkmod M=N×(N+1)/2mod M\sum_{k=1}^{N}k \mod M = N \times (N + 1) / 2 \mod M∑k=1N kmodM=N×(N+1)/2modM。
但是由于本题的 N(1≤N≤104×106)N (1 \le N \le 10^{4 \times 10^6})N(1≤N≤104×106) 非常大，直接计算上面的式子时间复杂度为 O(log2N)\mathrm{\mathrm{O}}(\log^2{N})O(log2N)，计算量大概为 1.6×10131.6 \times 10^{13}1.6×1013，会超出时间限制。
我们考虑 1+2+3+⋯+2×M−1+2×M=M×(2×M+1)1 + 2 + 3 + \cdots + 2 \times M - 1 + 2 \times M = M \times (2 \times M + 1)1+2+3+⋯+2×M−1+2×M=M×(2×M+1)。
而 M×(2×M+1)∣MM \times (2 \times M + 1) \vert MM×(2×M+1)∣M，即每 2×M2 \times M2×M 个数字的和对 MMM 取余的值为 000。
令 K=Nmod (M×2)K = N \mod (M \times 2)K=Nmod(M×2)，我们只需要计算最后的 KKK 个数字的和即可。
而最后的 KKK 个数字的形式为 i×M×2+1,i×M×2+2,i×M×2+3,⋯+i×M×2+Ki \times M \times 2 + 1, i \times M \times 2 + 2, i \times M \times 2 + 3, \cdots + i \times M \times 2 + Ki×M×2+1,i×M×2+2,i×M×2+3,⋯+i×M×2+K。
每个数字对 MMM 取余得到 1+2+3+⋯+K=K×(K+1)/21 + 2 + 3 + \cdots + K = K \times (K + 1) / 21+2+3+⋯+K=K×(K+1)/2。
对于本题 M(1≤M≤109),K(0≤K<M×2)M(1 \le M \le 10^9), K(0 \le K < M \times 2)M(1≤M≤109),K(0≤K<M×2)，可以直接做整形运算。
对于 NNN，令其长度为 nnn， 我们可以将其写成 ∑i=0na[i]×10i\sum_{i=0}^{n}a[i] \times 10^i∑i=0n a[i]×10i，其中 a[i]a[i]a[i] 为其从低位到高位上的数字。那么我们便可以从逐位处理求得 Nmod (M×2)N \mod (M \times 2)Nmod(M×2)。
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
C - 黑白方格
使用 BFS/DFS\tt{BFS/DFS}BFS/DFS 或者 并查集\tt{并查集}并查集 所有不同的连通块进行标记。
统计网格中连通块的数量记为 kkk，白色方格的数量为 www。
遍历所有 白色\tt{白色}白色 方格，统计每个白色方格周围不同的连通块数量，记为 cic_ici ，那么如果将当前方格染为 黑色\tt{黑色}黑色，网格上的连通块的数量变为为 k−ci+1k - c_i + 1k−ci +1。求出所有的 k−ci+1k - c_i + 1k−ci +1，最后的答案为 ∑i=1wk−ci+1w\frac{{\sum_{i=1}^{w}{k - c_i + 1}}}{w}w∑i=1w k−ci +1 。
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
D - 火星背包
本题的背包的容量，以及物品的大小都非常大，无法使用常规的01背包解法。
另外一方面 N(1≤N≤40)N(1 \le N \le 40)N(1≤N≤40)的范围不支持我们使用”搜索+剪枝“来在时限范围内求解。
但是我们可以把所有矿石分为两个部分，分别求解，最后将两个部分的答案合并。
这种方法叫做 折半枚举（meet in middle）\tt{折半枚举（meet \ in \ middle）}折半枚举（meet in middle）。
拆分成两个部分后，每个部分的矿石数量为不超过 202020，可以在 O(2N/2N)\mathrm{O}(2^{N/2}N)O(2N/2N) 的复杂度下求得两部分的答案。
我们可以先预处理前半或者后半部分中的选取方法对应的重量和价值总和记为 w1,v1w_1, v_1w1 ,v1 。这样在另一部分寻找总重 w2≤M−w1w_2 \le M - w_1w2 ≤M−w1 时使 v2v_2v2 最大的选取方法即可。
接下来思考如何从枚举得到的 (w2,v2)(w_2, v_2)(w2 ,v2 ) 集合中高效寻找 max{v2∣w2≤M′}\max\{v_2|w2\le M'\}max{v2 ∣w2≤M′} 的方法。
首先，可以排除所有 w2[i]≤w[j]w_2[i] \le w[j]w2 [i]≤w[j] 并且 v2[i]≥v2[j]v_2[i] \ge v_2[j]v2 [i]≥v2 [j] 的 jjj。这一点可以通过按照 w2,v2w_2, v_2w2 ,v2 的字典序排序后，处理得到。然后剩余的元素都满足 w2[i]<w2[j]⇔v2[i]<v2[j]w_2[i] < w_2[j] \Leftrightarrow v_2[i] < v_2[j]w2 [i]<w2 [j]⇔v2 [i]<v2 [j]，要计算 max{v2∣w2≤M′}\max\{v_2|w2\le M'\}max{v2
∣w2≤M′}，只需要找到满足 w2[i]≤M′w_2[i] \le M'w2 [i]≤M′ 的最大的 iii 即可。这一点可以通过二分查找高效实现。
令一半矿石的答案个数为 M(M≤2N/2)M(M \le 2^{N/2})M(M≤2N/2)，一次搜索需要
O(logM)\mathrm{O}(\log{M})O(logM) 的时间。这个算法总的时间复杂度为 O(2N/2N)\mathrm{O}(2^{N/2}N)O(2N/2N).
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
E - 连通三元组
我们使用 AiA_iAi 的二进制位来表示 iii 和其他点之间是否有边。
为了方便运算，我们可以只存储编号较小的点到编号较大点的边。
当 iii 和 jjj 固定不变时，满足条件的 kkk 的个数即 Ai,k=Aj,k=1A_{i,k}=A_{j,k}=1Ai,k =Aj,k =1 的个数，即 Ai AND AjA_i\ \mathrm{AND}\ A_jAi AND Aj 中 111 的个数。
因此，我们可以使用 std::bitset 来实现
A1,A2,…,ANA_1,A_2,\dots,A_NA1 ,A2 ,…,AN 并枚举 (i,j)(i,j)(i,j) 计算 Ai AND AjA_i\ \mathrm{AND}\ A_jAi AND Aj 中的 111 个数，总共需要 O(Nw)\mathrm{O}(\frac{N}{w})O(wN ) 时间其中 w=64w = 64w=64，因此总的时间复杂度为 O(N3w)\mathrm{\mathrm{O}}(\frac{N^3}{w})O(wN3 )。
由于 3000364=421875000≃4×108\frac{3000^3}{64}=421875000 \simeq 4 \times 10^86430003 =421875000≃4×108 ，这个解法已经足够快了。我们只需枚举 i,j(i<j)i,j(i < j)i,j(i<j) 就能进一步将执行时间减半。
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
F - 统计含K个非零数字的整数
考虑 数位DP\bf{数位DP}数位DP。
将数字 nnn 转化为字符串 sss 定义 f(i,cnt,isLim,isNum)f(i, cnt, isLim, isNum)f(i,cnt,isLim,isNum) 表示构造第 iii 位及其之后数位的合法方案数，其余参数含义为：
* cntcntcnt 表示前面选择的非零数字的数量。
* isLimisLimisLim 表示当前是否受到了 nnn 的约束（注意要构造的数字不能超过 nnn）。若为真，则第 iii 位填入的数字至多为 s[i]s[i]s[i]，否则可以是 999。如果在受到约束的情况下填了 s[i]s[i]s[i]，那么后续填入的数字仍会受到 nnn 的约束。例如 n=1234n=1234n=1234，那么 i=0i=0i=0 填的是 111 的话，i=2i=2i=2 的这一位至多填 222。
* isNumisNumisNum 表示 iii 前面的数位是否填了数字。若为假，则当前位可以跳过（不填数字），或者要填入的数字至少为 111；若为真，则要填入的数字可以从 000 开始。例如
n=123n=123n=123，在 i=0i=0i=0 时跳过的话，相当于后面要构造的是一个 999999 以内的数字了，如果 i=1i=1i=1 不跳过，那么相当于构造一个 101010 到 999999 的两位数，如果 i=1i=1i=1 跳过，相当于构造的是一个 999 以内的数字。
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
