|
|
@@ -0,0 +1,102 @@
|
|
|
++++
|
|
|
+date = "2021-08-03T19:45:00+08:00"
|
|
|
+draft = false
|
|
|
+tags = []
|
|
|
+title = "2021 杭电多校 (“中超联赛”) 第四场"
|
|
|
+authors = ["xry111"]
|
|
|
++++
|
|
|
+
|
|
|
+大概过程：
|
|
|
+
|
|
|
+* 开场签 1006，同时王老师签了 1003
|
|
|
+* 然后写 1007 错了 $114514$ 发，最后王老师改对了
|
|
|
+* 然后 1002 一开始觉得是压缩存储矩阵然后快速幂，发现事
|
|
|
+ $\mathcal{O}(n^4 \log L)$ 的。突然发现他有个傻了吧唧的条件是
|
|
|
+ $n | (M - 1)$，那就是裸的 Fourier 变换题了……
|
|
|
+ 结果被卡常数，加了个预处理单位根的幂才过。
|
|
|
+* 最后他俩在那自闭 1005 和 1009 不知道在干啥，我自闭 1004，
|
|
|
+ 最后没思路了开始用内嵌汇编乱卡，赛后卡过了
|
|
|
+
|
|
|
+# [String Mod](https://acm.hdu.edu.cn/showproblem.php?pid=7013)
|
|
|
+
|
|
|
+## 题意
|
|
|
+
|
|
|
+给定 $n$, $L$, $k$，设 $a^n = b^n = 1$，求
|
|
|
+
|
|
|
+$$(a + b + (k - 2))^n$$
|
|
|
+
|
|
|
+中 $a^i b^j$ 项的系数模 $M$。保证 $1 \leq n \leq 500$，$n | (M - 1)$。
|
|
|
+
|
|
|
+## 做法
|
|
|
+
|
|
|
+一开始想压缩存储一个矩阵然后快速幂，发现乘法的复杂度事 $\mathcal{O}(n^4)$，
|
|
|
+肯定过不了。然后突然发现一个怪异的条件 $n | (M - 1)$，这就意味着在模
|
|
|
+$M$ 意义下一定有 $n$ 阶单位根，
|
|
|
+再回头一看之前那个“矩阵乘”其实就是二维循环卷积，
|
|
|
+那么傻子都知道肯定是 Fourier 变换以后直接快速幂一下再变回去了
|
|
|
+(类似 [XDOJ 1085](https://acm.xidian.edu.cn/problem.php?id=1085)，
|
|
|
+只不过变成了二维)。
|
|
|
+
|
|
|
+因为 $n$ 才 $500$，不需要写 FFT (这个模数写 FFT 得看 $114514$ 篇论文，
|
|
|
+然后写 $1919810$ 行)，直接用 DFT 的定义 $\mathcal{O}(n^2)$
|
|
|
+地硬算每行每列就行了，这样二维 DFT 的复杂度就是 $\mathcal{O}(n^3)$。
|
|
|
+因此总的复杂度是 $\mathcal{O}(n^3 + n^2 \log L)$。
|
|
|
+
|
|
|
+本题需要对于每组读入预处理单位根的幂，否则会卡常。
|
|
|
+
|
|
|
+# [Another String](https://acm.hdu.edu.cn/showproblem.php?pid=7015)
|
|
|
+
|
|
|
+## 题意
|
|
|
+
|
|
|
+给定一个字符串 `S` 和一个数 $k$，对于每个 $i$，
|
|
|
+令 $U$ 是 `S[:i]` 所有子串的集合，$V$ 是 `S[i:]` 所有子串的集合，求
|
|
|
+
|
|
|
+$$\sum\_{u \in U} \sum\_{v \in V} \[|u| == |v|\]
|
|
|
+\\left \[\sum\_i \[u(i) \neq v(i)\] \leq k \\right \]$$
|
|
|
+
|
|
|
+## 做法 (?)
|
|
|
+
|
|
|
+很容易想到一个 $\mathcal{O}(n^3)$ 的做法：
|
|
|
+枚举 $a$, $t$, $b$，如果 `s[a:a+t]` 和 `s[b:b+t]`
|
|
|
+中字符不同的位置不超过 $k$ 个，就要对 $[a + t, b]$ 中的每个 $i$
|
|
|
+答案都加 $1$。因为只有一次查询，不用写 xx 树，最后求个前缀和即可。
|
|
|
+
|
|
|
+本地测了一下发现跑 $3$ 秒多，看了一下上面那个 $[a + t, b]$，
|
|
|
+发现对于固定的 $a$ 和 $b$，合法的 $t$ 肯定是从 $1$ 开始的一段连续的整数
|
|
|
+(除非没有合法的 $t$)。这样，如果我们对于特定的 $a$ 和 $b$
|
|
|
+能够较快地求出最大的合法的 $t$，然后把前缀和做两次，就能解决这个问题。
|
|
|
+最后实在没办法了，用汇编写了个求 $16$ 个字符中有多少不同的函数：
|
|
|
+
|
|
|
+```c++
|
|
|
+static inline int n_mismatch(const char *a, const char *b)
|
|
|
+{
|
|
|
+ int r;
|
|
|
+ static const long x = 1;
|
|
|
+ __asm__(
|
|
|
+ "vmovdqu (%1), %%xmm0\n"
|
|
|
+ "vpcmpeqb (%2), %%xmm0, %%xmm0\n"
|
|
|
+ "vpcmpeqd %%xmm1, %%xmm1, %%xmm1\n"
|
|
|
+ "vpxor %%xmm1, %%xmm0, %%xmm0\n"
|
|
|
+ "vpmovzxbw %%xmm0, %%ymm0\n"
|
|
|
+ "vextracti128 $1, %%ymm0, %%xmm1\n"
|
|
|
+ "vpbroadcastd %3, %%ymm2\n"
|
|
|
+ "vpmovzxwd %%xmm1, %%ymm1\n"
|
|
|
+ "vpand %%ymm2, %%ymm1, %%ymm1\n"
|
|
|
+ "vpmovzxwd %%xmm0, %%ymm0\n"
|
|
|
+ "vpand %%ymm2, %%ymm0, %%ymm0\n"
|
|
|
+ "vpaddd %%ymm1, %%ymm0, %%ymm0\n"
|
|
|
+ "vextracti128 $1, %%ymm0, %%xmm1\n"
|
|
|
+ "vpaddd %%xmm1, %%xmm0, %%xmm0\n"
|
|
|
+ "vpshufd $238, %%xmm0, %%xmm1\n"
|
|
|
+ "vpaddd %%xmm1, %%xmm0, %%xmm0\n"
|
|
|
+ "vpshufd $85, %%xmm0, %%xmm1\n"
|
|
|
+ "vpaddd %%xmm1, %%xmm0, %%xmm0\n"
|
|
|
+ "vmovd %%xmm0, %0"
|
|
|
+ : "=r"(r)
|
|
|
+ : "r"(a), "r"(b), "m"(x)
|
|
|
+ );
|
|
|
+ return r;
|
|
|
+}
|
|
|
+```
|
|
|
+
|
|
|
+赛后才写完卡了过去…… 正解是啥我暂时蒙在鼓里。
|
Xi Ruoyao