[BZOJ2882][后缀自动机]工艺

LowestJN

于 2017-01-30 22:57:11 发布

阅读量550

点赞数

CC 4.0 BY-SA版权

分类专栏：后缀自动机文章标签：后缀自动机

本文链接：https://blog.youkuaiyun.com/Coldef/article/details/54790695

后缀自动机专栏收录该内容

14 篇文章

订阅专栏

本文介绍了一种利用后缀自动机解决特定字符串处理问题的方法。通过构建双倍序列的后缀自动机，并采用贪心策略选取路径，最终求解出最短匹配路径。文章包含完整的C++代码实现，适用于需要高效字符串搜索的应用场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

后缀自动机裸题
将序列复制两遍建后缀自动机
每次贪心地走最小的边
走N次得到答案

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <map>
#define N 1000010

using namespace std;

struct SAM_{
    int fail[N],stp[N],p,cnt;
    map<int,int> next[N];
    SAM_(){p=cnt=1;}
    void Extend(int x){
        int np=++cnt;stp[np]=stp[p]+1;
        while(p&&!next[p].count(x)) next[p][x]=np,p=fail[p];
        if(!p) fail[np]=1;
        else{
            int q=next[p][x];
            if(stp[q]==stp[p]+1) fail[np]=q;
            else{
                int nq=++cnt;stp[nq]=stp[q]+1;
                next[nq]=next[q];
                fail[nq]=fail[q];
                fail[q]=fail[np]=nq;
                while(p&&next[p][x]==q) next[p][x]=nq,p=fail[p];
            }
        }
        p=np;
    }
    void Explorer(int n){
        p=1;
        for(int i=1;i<=n;i++){
            printf("%d",next[p].begin()->first);
            p=next[p].begin()->second;
            if(i!=n)putchar(' ');else putchar('\n');
        }
    }
}SAM;

int n,A[N];

inline void reaD(int &x){
    char Ch=getchar();x=0;
    for(;Ch>'9'||Ch<'0';Ch=getchar());
    for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=getchar());
}

int main(){
    reaD(n);
    for(int i=1;i<=n;i++) reaD(A[i]);
    for(int i=1;i<=n;i++) SAM.Extend(A[i]);
    for(int i=1;i<=n;i++) SAM.Extend(A[i]);
    return SAM.Explorer(n),0;
}