Lyndon分解(后缀中的最小串)

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} xixi+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[ik] 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 kj),所以一节一节来就行。


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
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
#define fi first
#define se second
#define debug(x) cerr<<#x<<" = "<
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;
}

你可能感兴趣的:(字符串)