2021年07月训练

一、练习题分析

CF 1547F

其实我们已经想到了本质是求最长的一段不是全局公约数的子序列长度，而且$O(\sum n)$的提示也容易使人明白需要用类似尺取的方法。但差点意思的是怎么尺取呢？

之前失败的思路，问题在于要考虑的质因数太多，并不能动态维护当前段的所有质因数。所以这时候我们需要估计对于每一个数，其质因数个数的上限——其实就是8个，所以一个一个考虑是足够的。那么问题需要转化，也就是首先所有的数都除掉全局公约数，这样就是求最长的一段$gcd>1$的子序列。

那么我们按照线性筛得到的筛法树的边权**（注意严格保存最小质因子）**，容易得到每个数的因数分解，那么对于每一个数$a_i$,枚举它的所有质因数$p$，然后向左向右扩展，每扩展一个数，就把它的$p$从其质因数集合中消去，这样就可以不断更新最大长度。

算法标签 数论 线性筛 双指针/尺取法

CF 1547G

基本思路是正确的，就是要用SCC来统一处理环。然后1是唯一确定的源点，这个条件很好，使得DFS和Tarjan过程仅从1开始就行。然后从1开始到某点的访问次数，由于只考虑2次以上有限次，所以DFS就可以统计，并不需要DP。最后凡是从1到某个点路径经过SCC（或者自环）的，$an{s}_i=-1$

算法标签 图论 强连通分量

CF 1391D

这个题其实还是一个奇偶的构造。

首先行数超过4行和行数只有1行的情况很容易处理。那么先考虑两行的情况。这时的偶方阵就是$2\times 2$的方阵，所以如果有奇数个1，相邻$2\times 1$的列的1个数奇偶性应该不同。所以整个矩阵最终效果，就是“奇偶奇偶”的情形，或者"偶奇偶奇"的情形。按照这两种情况扫一遍即可。

那么行数为3的情况怎么样呢？其实这里就是化归，前两行和后两行分别实现$n=2$时奇偶分布，所以是四种情况，同样扫就行。

算法标签 奇偶分析

CF 1550D

这题核心在于一个偏移量的思想。设置$a_i=i+k_i$然后分析$k_i$的性质，来计数。具体思路见题解。

CF1550F

这个题目很综合，也体现建图和图的处理的一些技巧。

首先怎么考虑每一次$[d-k,d+k]$这个范围的跳跃的构建。其实最朴素的想法是枚举两个节点，连接一条边，边权就是坐标差的绝对值。那么适配这个题目，就再把$d$考虑进去，变成$w=|d-|a_u-a_v||$。所以如果从$s$到$i$经过的若干边，他的边权$w$都不超过$k$，那么就是满足要求的。

考虑到每一次的$k$是变动的，$k$越大，相当于跳跃的约束越小。因此我们就考虑尽可能用$w$的边来构成从$s$到$i$的路径，这样就可以更大化地满足要求苛刻的$k$。于是我们可以考虑求一棵最小生成树。

但是个$n^2$量级的完全图，因此普通的Kruskal和Prim是不行的。我们用一个新算法——$Boruvka$。

Boruvka 算法的一句话思想便是：

“从所有当前的连通块向其他连通块扩展出最小边，直到只剩一个连通块”，其中取最小边的贪心思想是 Kruskal 的主体，而向外扩展又是 Prim 的思想 —— 基于另外两种生成树算法，Boruvka 的正确性显然。

其代码框架如下：

while 连通块个数>1
    for 每个连通块 i
        mn[i] = 连接 i 与其他连通块的最小边
    for 每个连通块 i
        if mn[i] 连接两个不同的连通块
            ans += mn[i]
            Merge( mn[i] 连接的连通块 )
            连通块个数 --

而不同的题目难点就在于快速计算$m{n}_i$。对于这道题而言，其实$m{n}_i$就是要找到可以向其他连通块扩展的权值在$[a_i-d,a_i+d]$范围内的最小边。那么我们可以利用set+lower_bound来求这个最小边。

所以求得最小边后，就用vector维护的并查集来合并（小并入大）维护连通块，同时在新图上连接无向边。

运行完Boruvka算法后，再从$s$出发DFS遍历整个MST，求得它到其他点路径上的权值最大值。最后处理询问时只需要考虑这个权值最大值是不是比k来的小就行了。

复杂度$O(nlo{g}^{2}n)$

算法标签 最小生成树 启发式合并 图论建图

CF1550E

这道题核心在于模拟匹配。

从Frozen_Guardian的博客上学到了关于这种状压的一些真知灼见。

状压$dp$就是全排列问题的一种优化，所以状压$dp$可以解决的问题，全排列暴力枚举肯定也是可以解决的。所以先将问题简化，假设本题放开全排列该如何去做？

我们要全排列什么？其实相当于从左到右依次放置待判定长度的相同字母子串。所以问题的关键就在于判别不同排列时能不能恰好匹配？所以这里就借助了序列自动机的概念。我们记$nnext[i][j]$表示从字符串第$i$个字符开始，其最早完全匹配第$j$种字符的待判定长度的字串时，下一个位置是多少？那么计算就是从后往前计算，如果借助？和已有的字母能够新增出一个完全的匹配，那么就赋值，否则照搬$nnext[i+1][j]$。

再考虑$dp$的转移，很容易想到$dp[mask]$就是已经使用字母的状态，这个dp值的含义是完成状态内这些字母的匹配后，最小的下一个位置（相当于一个$end$指针）。所以说我们其实并不关心从左到右到底是如何安排的，只需要按照$nnext$来辅助转移就可以了。

所以最后的判定标准就是如果$dp[2^k-1]$在$n-1$范围内，就是可以的。主函数二分答案即可。

算法标签 二分答案 序列自动机 状压DP

Vector的简易写法和for的C++0x用法

vector<int> a;
a.resize(n);  //使得其大小为n，访问[0,n-1]
vector<int> a(n); //大小为n的vector，访问同上
vector<int> a(1,7);//一个元素，值为7
for(auto x : comp)   //x为comp容器的遍历器，下面直接访问x如何
{
    xxx
}
for(int i:{-d,d})//数值范围内遍历
//如果是vector<vector <xxx> > Q；这样的二维vector情形，则考虑
for(const vector<int> &comp : comps) {  //不过遇事不决用auto 就可以了

计算几何学习

见手写笔记