如何用 Lyndon 薄纱 SA

没卡常,目前用时是次优解的一半左右。

推销 Lyndon 理论学习笔记(没写完)

P5108 仰望半月的夜空

给一个字符串 s,对 i[1,|s|] 问长度为 i 的最小子串的第一次出现。

先给复杂度:问最后一次出现可以 O(|s|),问第一次出现可以结合哈希二分做到 O(|s|log|s|)

字典序问题,可以在 Lyndon 分解上考虑,设分解得到 s=w1++wk

考虑长度为 i 的最小子串在哪里起头。首先根据 Lyndon 串的定义,在某个 wp 中间起头肯定不优,

而且分解出的 Lyndon 串字典序是单调不增的,所以应该在尽量靠后的 wp 起头,

所以可以得出结论:应该在 |wp+wp+1++wk|i 的最后一个 wp 处起头,

但这个子串不一定只在这里出现,考虑它还在哪里出现过。

分解出的 Lyndon 串字典序是单调不增的,所以这个子串只在 wp 及其前的若干个 Lyndon 串中出现。

二分这个子串最早在哪里出现即可。

#include <cstdio>
#include <cstring>
char _[300050];
int n, o, S, l[300050], a[300050];
unsigned long long p[300050], h[300050];
int Q(int l, int r) { return h[r] - h[l - 1] * p[r - l + 1]; }
int main()
{
    scanf("%d%d", &S, &n);
    if (S == 26)
    {
        scanf("%s", _ + 1);
        for (int i = 1; i <= n; ++i)
            a[i] = _[i];
    }
    else
    {
        for (int i = 1; i <= n; ++i)
            scanf("%d", a + i);
    }
    for (int i = p[0] = 1; i <= n; ++i)
        p[i] = p[i - 1] * 10000019, h[i] = h[i - 1] * 10000019 + a[i];
    int i = 1, j, k;
    while (i <= n)
    {
        j = i + 1, k = i;
        while (j <= n)
        {
            if (a[j] == a[k])
                ++j, ++k;
            else if (a[j] > a[k])
                ++j, k = i;
            else
                break;
        }
        while (i <= k)
            l[++o] = i, i += j - k;
    }
    for (int i = 1, j = o; i <= n; ++i)
    {
        if (n - l[j] + 1 < i)
            --j;
        int L = 1, R = j;
        while (L <= R)
        {
            int M = L + R >> 1;
            if (Q(l[M], l[M] + i - 1) == Q(l[j], l[j] + i - 1))
                R = M - 1;
            else
                L = M + 1;
        }
        printf("%d ", l[L]);
    }
    return 0;
}
posted @   Jijidawang  阅读(56)  评论(2编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 开源Multi-agent AI智能体框架aevatar.ai,欢迎大家贡献代码
· Manus重磅发布:全球首款通用AI代理技术深度解析与实战指南
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· AI技术革命,工作效率10个最佳AI工具
点击右上角即可分享
微信分享提示