SPOJ8222

来自蒟蒻XXJ的做题记录

”PIA叽“ 哎呀我脸好疼……

其实这个就是在right数组上搞事情

然后注意right数组之间全都是包含关系

而每一个right数组的终点就是我现在插入进去的前缀

于是一开始计算的时候吧主链上的right数组全都开始++ 这样才能保证递推的开始是能够涵盖所有字串的

然后按照lth大小排序 把儿子的right数组全都加到父亲上面去

然后用f[i]去更新f[i-1]就可以了

#include<bits/stdc++.h>
#define mem(i,j) memset(i,j,sizeof(i))
#define mcy(i,j) memcpy(i,j,sizeof(i))
using namespace std;
const int MAXN=500010;
struct SAM{
    int trans[MAXN][26],r[MAXN],fa[MAXN],lth[MAXN];
    int last,cnt;
    SAM(){
        cnt=0;last=++cnt;
        mem(trans,0);mem(r,0);mem(fa,0);mem(lth,0);
    }
    void extend(int c){
        int np=++cnt,p=last;last=cnt;
        lth[np]=lth[p]+1;
        for(;p&&!trans[p][c];p=fa[p]) trans[p][c]=np;
        if(!p) fa[np]=1;
        else{
            int q=trans[p][c];
            if(lth[q]==lth[p]+1) fa[np]=q;
            else{
                int nq=++cnt;lth[nq]=lth[p]+1;
                mcy(trans[nq],trans[q]);
                fa[nq]=fa[q];
                fa[q]=fa[np]=nq;
                for(;p&&trans[p][c]==q;p=fa[p]) trans[p][c]=nq;
            }
        }
    }
    void build(char str[]){
        int len=strlen(str);
        for(int i=0;i<len;i++) extend(str[i]-'a');
    }
}sam;
char c[MAXN];
int rank[MAXN],n,ans[MAXN],f[MAXN];
bool cmp(int i,int j){
    return sam.lth[i]>sam.lth[j];
}
void init(){
    for(int i=1;i<=sam.cnt;++i) rank[i]=i;
    sort(rank+1,rank+sam.cnt+1,cmp);
    int p=1;
    for(int i=0;i<n;i++){
        int x=c[i]-'a';
        p=sam.trans[p][x];sam.r[p]++;
    }
    for(int i=1;i<=sam.cnt;++i){
        int ind=rank[i];
        if(!sam.fa[ind]) continue;
        sam.r[sam.fa[ind]]+=sam.r[ind];
    }
}
void input(){
    scanf("%s",c);n=strlen(c);
    sam.build(c);
    init();
}
void xxj(){
    for(int i=1;i<=sam.cnt;++i)
        ans[sam.lth[i]]=max(ans[sam.lth[i]],sam.r[i]);
    f[n]=ans[n];
    for(int i=n-1;i;--i) f[i]=max(max(f[i],f[i+1]),ans[i]);
}
void output(){
    for(int i=1;i<=n;i++) cout<<f[i]<<"\n";
}

int main(){
    input();
    xxj();
    output();
    return 0;
}

转载于:https://www.cnblogs.com/Xiaojianxiang/p/6598724.html

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值