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

30 篇文章 2 订阅

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 &lt; b a&lt;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 ] &gt; x [ k ] x[j]&gt;x[k] x[j]>x[k],我们不允许 x i &lt; x i + 1 x_i&lt;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 ] &lt; x [ k ] x[j]&lt;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<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;
}

  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值