bzoj2865 字符串识别

Description

XX在进行字符串研究的时候，遇到了一个十分棘手的问题。

在这个问题中，给定一个字符串S，与一个整数K，定义S的子串T=S(i, j)是关于第K位的识别子串，满足以下两个条件：

1、i≤K≤j。

2、子串T只在S中出现过一次。

例如，S="banana"，K=5，则关于第K位的识别子串有"nana"，"anan"，"anana"，"nan"，"banan"和"banana"。

现在，给定S，XX希望知道对于S的每一位，最短的识别子串长度是多少，请你来帮助他。

Input

仅一行，输入长度为N的字符串S。

Output

输出N行，每行一个整数，第i行的整数表示对于第i位的最短识别子串长度。

建SAM，只有right集大小为1的节点对答案有贡献，

若其出现位置右端点为r，此节点可接受的最短串长为x，最长串长为y，

则对(r-x,r]用x更新最小值，对r-k (y<k≤x)则用k更新最小值

用两棵线段树维护答案，分别处理以上两种情况

此题卡内存...

#include<cstdio>
#include<cstring>
const int N=;
int nx[N][],fa[N],l[N],r[N],t[N],q[N],ql=,qr=,d[N],ptr=,pv=,len=;
char s[];
void ins(int w,int pos){
    int p=pv,np=++ptr;
    t[np]=;
    r[np]=pos;
    l[np]=l[p]+;
    while(p&&!nx[p][w])nx[p][w]=np,p=fa[p];
    if(!p)fa[np]=;
    else{
        int q=nx[p][w];
        if(l[q]==l[p]+)fa[np]=q;
        else{
            int nq=++ptr;
            memcpy(nx[nq],nx[q],sizeof(nx[]));
            l[nq]=l[p]+;
            fa[nq]=fa[q];
            fa[q]=fa[np]=nq;
            while(p&&nx[p][w]==q)nx[p][w]=nq,p=fa[p];
        }
    }
    pv=np;
}
void build(){
    for(int i=;i<=ptr;i++)++d[fa[i]];
    for(int i=;i<=ptr;i++)if(!d[i])q[qr++]=i;
    while(ql!=qr){
        int w=q[ql++];
        int u=fa[w];
        if(!u)continue;
        t[u]+=t[w];
        r[u]=r[w];
        if(!--d[u])q[qr++]=u;
    }
}
int mn[],mn2[];
inline void mins(int&a,int b){if(a>b)a=b;}
void setmin(int l,int r,int x,int*mn){
    for(l+=,r+=;l^r^;l>>=,r>>=){
        if(~l&)mins(mn[l^],x);
        if(r&)mins(mn[r^],x);
    }
}
int ans[];
int main(){
    memset(mn,,sizeof(mn));
    memset(mn2,,sizeof(mn2));
    scanf("%s",s);
    for(len=;s[len];++len)ins(s[len]-'a',len+);
    build();
    for(int i=;i<=ptr;i++)if(t[i]==){
        int ml=l[fa[i]];
        setmin(r[i]-ml+,r[i],ml+,mn);
        setmin(r[i]-l[i]+,r[i]-ml,r[i]+,mn2);
    }
    for(int i=;i<;i++){
        int lc=i<<,rc=lc^;
        mins(mn[lc],mn[i]),mins(mn[rc],mn[i]);
        mins(mn2[lc],mn2[i]),mins(mn2[rc],mn2[i]);
    }
    for(int i=;i<=len;i++){
        ans[i]=mn[i+];
        mins(ans[i],mn2[i+]-i);
    }
    for(int i=;i<=len;i++)printf("%d\n",ans[i]);
    return ;
}