题目
思路
最简单的做法是,直接套用强化版问题,做到 O ( n log n ) \mathcal O(n\log n) O(nlogn) 。
但在这个过程中我们想到:不妨设原串是所有后缀中最大的(类似 L y n d o n \rm Lyndon Lyndon 串)那么 j j j 会因为考虑前缀变长而成为答案,只可能是 S [ j : ] S[j:] S[j:] 为 b o r d e r \rm border border 。
看上去没用。仔细想想,对于这种类 L y n d o n \rm Lyndon Lyndon 串,它满足 b o r d e r \rm border border 之后的字符随长度减小而增大。直接证明就是,记 S S S 的 b o r d e r \rm border border 为 T T T,则以 T T T 的后缀 b o r d e r \rm border border 作为起点,需要大于等于以前缀 b o r d e r \rm border border,即原串 S S S 的开头,作为起点。因此就能推出该结论。
事实上通过代码,归纳法即得。因为只需要暴力跳过 b o r d e r \rm border border 之后的字符过小的情况。时间复杂度 O ( n ) \mathcal O(n) O(n) 。
代码
#include <cstdio> // megalomaniac JZM yydJUNK!!!
#include <iostream> // Almighty XJX yyds!!!
#include <algorithm> // oracle: ZXY yydBUS!!!
#include <cstring> // DDG yyd Black Bridge!!!
#include <cctype> // decent XYX yydLONELY!!!
typedef long long llong;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
# define rep0(i,a,b) for(int i=(a); i!=(b); ++i)
inline int readint(){
int a = 0, c = getchar(), f = 1;
for(; !isdigit(c); c=getchar()) if(c == '-') f = -f;
for(; isdigit(c); c=getchar()) a = a*10+(c^48);
return a*f;
}
const int MAXN = 1000006;
char str[MAXN];
int a[MAXN], fail[MAXN];
int main(){
scanf("%s",str+1);
int n = int(strlen(str+1)), len = 1;
a[1] = str[1], fail[1] = 0; puts("1 1");
rep(i,2,n){
while(len && a[fail[len]+1] < str[i]) len = fail[len];
a[++ len] = str[i]; printf("%d %d\n",i-len+1,i);
int &j = fail[len] = fail[len-1]; // kmp
while(j && a[j+1] != a[len]) j = fail[j];
if(len != 1 && a[j+1] == a[len]) ++ j; // as normal
}
return 0;
}
后记
主要是没能发见那个性质。说起来, L y n d o n \rm Lyndon Lyndon 串本是比较奇特的玩意儿;据称,本题也有 L y n d o n \rm Lyndon Lyndon 串分解(基于 D u v a l \rm Duval Duval 算法)的做法。
我感觉 D u v a l \rm Duval Duval 肯定借用了 KMP \textit{KMP} KMP 的某些做法,但我找不出。