After All
L y n d o n Lyndon Lyndon串: 对于字符串 x x x,如果 x x x的字典序严格小于 x x x的所有后缀的字典序,我们称 x x x是简单串,或者 L y n d o n Lyndon Lyndon串。
近似 L y n d o n Lyndon Lyndon串: 若 x x x为 L y n d o n Lyndon Lyndon串,则 x x x x ′ xxxx' xxxx′为近似 L y n d o n Lyndon Lyndon串, x ′ x' x′为 x x x的前缀。
L y n d o n Lyndon Lyndon分解: 将一个串 x x x分作 x 1 x 2 . . x k x_1x_2..x_k x1x2..xk,每一个部分都是 L y n d o n Lyndon Lyndon串,且 x i ≥ x i + 1 x_i\geq x_{i+1} xi≥xi+1。
定理:
-
L
y
n
d
o
n
Lyndon
Lyndon串是循环同构中字典序最小的一个。
-
L
y
n
d
o
n
Lyndon
Lyndon分解唯一。
- 两个
L
y
n
d
o
n
Lyndon
Lyndon串
a
,
b
a,b
a,b,若
a
<
b
a<b
a<b,有
a
b
ab
ab为
L
y
n
d
o
n
Lyndon
Lyndon串。
考虑 i i i之前的位置都已经被分解为 L y n d o n Lyndon Lyndon串,用指针 j j j指向 i i i, k k k指向 i + 1 i+1 i+1。
- 若
x
[
j
]
=
x
[
k
]
x[j]=x[k]
x[j]=x[k],加上
x
[
k
]
x[k]
x[k]不会影响
i
i
i开头的串为相似
L
y
n
d
o
n
Lyndon
Lyndon串。所以比较后一位,两个指针都往后移;
- 若
x
[
j
]
>
x
[
k
]
x[j]>x[k]
x[j]>x[k],我们不允许
x
i
<
x
i
+
1
x_i<x_{i+1}
xi<xi+1,而允许
L
y
n
d
o
n
Lyndon
Lyndon串后面跟一个小的串,所以包含这一位。这个串是一个整体,所以使
j
=
i
j=i
j=i,重新开始比较(
x
[
i
→
k
]
x[i\to k]
x[i→k]是
L
y
n
d
o
n
Lyndon
Lyndon串);
- 若
x
[
j
]
<
x
[
k
]
x[j]<x[k]
x[j]<x[k],就不行了,所以当前串就结束了;但是当前串可能是近似
L
y
n
d
o
n
Lyndon
Lyndon串(
L
y
n
d
o
n
Lyndon
Lyndon串的长度为
k
−
j
k-j
k−j),所以一节一节来就行。
Problem Link - https://loj.ac/problem/129
题意:
求 L y n d o n Lyndon Lyndon分解后所有串的右端点。
代码:
/*
* Author : Jk_Chen
* Date : 2019-09-05-10.09.48
*/
#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define rep(i,a,b) for(int i=(int)(a);i<=(int)(b);i++)
#define per(i,a,b) for(int i=(int)(a);i>=(int)(b);i--)
#define mmm(a,b) memset(a,b,sizeof(a))
#define pb push_back
#define pill pair<int, int>
#define fi first
#define se second
#define debug(x) cerr<<#x<<" = "<<x<<'\n';
const LL mod=1e9+7;
const int maxn=(1<<21);
LL rd(){ LL ans=0; char last=' ',ch=getchar();
while(!(ch>='0' && ch<='9'))last=ch,ch=getchar();
while(ch>='0' && ch<='9')ans=ans*10+ch-'0',ch=getchar();
if(last=='-')ans=-ans; return ans;
}
/*_________________________________________________________head*/
vector<int>Pos;
void deal(char *x,int n){
int i=0;
while(i<n){
int j=i,k=i+1;
while(k<n&&x[j]<=x[k]){
if(x[j]==x[k])
j++;
else
j=i;
k++;
}
while(i<=j){
i+=k-j;
Pos.push_back(i);
}
}
}
char x[maxn];
// every substring's endpos , index begin at 1
int main(){
gets(x);
deal(x,strlen(x));
for(auto P:Pos){
printf("%d ",P);
}
puts("");
return 0;
}