Lyndon分解和最小循环表示学习

做CF594E涉及的两个知识点。以下字符串采用Python记法。

Lyndon分解

定义 $S$ 是Lyndon串，当且仅当对于任意有意义的正整数 $i$ 有 $S

定义 $S$ 的Lyndon分解是一个Lyndon串的序列 $s_1, s_2, \ldots, s_n$, 使得 $S=s_1s_2 \cdots s_n$ 并且 $s_1 \ge s_2 \ge \cdots \ge s_n$.

Lyndon分解存在且唯一。

不难发现，Lyndon分解可以这么得到：对于 $S$, 取最小的后缀 $S[i:]$, $S$ 的Lyndon分解，是在 $S[:i]$ 的Lyndon分解的最后加上一项 $S[i:]$ 得到的序列。

于是，Lyndon分解可以用后缀数组求，但是这样太复杂了。

Lyndon分解的Duval算法：

设原字符串为 $S$. 逐个加入字符，设已能确定 $S[:i]$ 的Lyndon分解 $s_1, s_2, \ldots, s_n$ 为 $S$ 的Lyndon分解的前缀，再尽可能地扩充字符串 $S[i:k]$, 使得 $S[i:k]$ 具有Lyndon周期 $t$, 也就是说，$S[i:k]=(m \times t)u$, 其中 $m$ 是正整数，$t$ 是Lyndon串，$u$ 是 $t$ 的（可以为空的）真前缀。考虑加入字符 $S[k]$, 若 $k=|S|$ 则令 $S[k]=-\infty$.

• 若 $S[k]
• 若 $S[k]=S[k-|t|]$, 则 $t$ 也是 $S[i:k+1]$ 的Lyndon周期。
• 若 $S[k]>S[k-|t|]$, 则 $S[i:k+1]$ 的Lyndon周期是其本身。

实现上我们只需要维护 $i, k$ 以及 $j=k-|t|$. 代码链接

扩充 $S[i:k]$ 的过程中 $i+k\le2|S|$ 且随扩充总轮数递增，因此该算法时间 $O(|S|)$, 除去输入输出只需要 $O(1)$ 额外空间，是一个非常简短而高效的算法。

最小循环表示

最小循环表示，也就是对于字符串 $S$, 求出最小的 $T_i=S[i:]S[:i]$.

我们维护两个决策 $i, j$, 满足 $i

初始时，设 $i=0, j=1$.

我们通过枚举比较求出 $T_i$ 和 $T_j$ 的最长公共前缀 $k$.

• 若 $T_i=T_j$, $T_i$ 即为最优解。这是因为 $S$ 以 $j-i$ 为循环节循环，而 $(i, j)$ 中任何整数都不是最优决策。
• 若 $T_i
• 若 $T_i>T_j$, 同理令 $i \gets i+k+1$, 若此后 $i \ge j$ 则再令 $j \gets i+1$, 上述性质仍未改变。

当 $j \ge |S|$ 时，我们发现 $[0, i)$ 和 $(i, |S|)$ 中的整数都不可能为最小决策了，所以 $T_i$ 是唯一的最优解；否则重复此过程。

$i+j+k \le 3|S|$ 随枚举总次数递增，因此该算法时间 $O(|S|)$, 除去输入输出只需要 $O(1)$ 额外空间，是一个非常简短而高效的算法。

代码：$n$（$n \le 3\times10^5$）项整数序列的最小循环表示。

不过这份代码由于在OJ上刷常数榜，不仅加入了输入输出优化，还破环为链以节省取模的时间，空间略大。

 1 #include 
 2 char r[1<<25], *rc=r, w[1<<25], *wc=w;
 3 int rd() {
 4     int a;
 5     bool b;
 6     unsigned char c;
 7     while(c=*rc++-48, c>9&c!=253);
 8     b=(c==253);
 9     for(a=b?0:c; (c=*rc++-48)<10; a=a*10+c);
10     return b?-a:a;
11 }
12 void wr(int v) {
13     char s[10], *p=s;
14     while(*p++=v%10+48, v/=10);
15     while(*wc++=*--p, p!=s);
16 }
17 const int N=3e5;
18 int n, x[2*N];
19 int main() {
20     int i, j, k;
21     fread(r, 1, 1<<25, stdin);
22     n=rd();
23     for(i=0; ird();
24     memcpy(x+n, x, n*sizeof(int));
25     for(i=0, j=1, k=0; jn; ) {
26         int xi=x[i+k], xj=x[j+k];
27         if(xi==xj) ++k;
28         else if(xi1, k=0;
29         else {
30             i+=k+1, k=0;
31             if(i>=j) j=i+1;
32         }
33     }
34     for(k=0; k' ';
35     wc[-1]='\n';
36     fwrite(w, 1, wc-w, stdout);
37     return 0;
38 }

最小循环表示