CF700E Cool Slogans

CF700E Cool Slogans

题目描述

给出一个长度为n的字符串\(s[1]\),由小写字母组成。定义一个字符串序列\(s[1....k]\),满足性质:\(s[i]\)\(s[i-1] (i>=2)\)中出现至少两次(位置可重叠),问最大的\(k\)是多少,使得从\(s[1]\)开始到\(s[k]\)都满足这样一个性质。

很妙的题啊。

首先\(s[i]\)一定是\(s[i+1]\)的后缀。因为如果不是,我们可以吧多余的部分删除,这样不影响答案。

我们建出后缀自动机,然后在\(fail\)树上\(DP\)。我们设\(f_v\)表示后缀自动机上\(v\)节点代表的子串作为最后一个串的答案。

更新的时候就判断\(v\)\(fail\)代表的节点是否在\(v\)代表的节点中出现了两次,如果是,\(f_v=f_{fail_v}+1\)。否则,\(v\)就没用了,我们更新\(v\)的儿子时仍然用\(fail_v\)来更新,因为他们\(f\)值相同,但是\(fail_v\)更短,所以更优。

判断一个串是否再另一个串中出现了两次可以用线段树合并搞一搞。

代码:

#include<bits/stdc++.h>
#define ll long long
#define N 400005

using namespace std;
inline int Get() {int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}while('0'<=ch&&ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;}

int n;
char s[N];
int mxlen[N<<1],fail[N<<1];
int ch[N<<1][26];
int last=1,cnt=1;
int pos[N<<1];
void Insert(int f,int Pos) {
    int p=last,v=++cnt;
    pos[v]=Pos;
    last=v;
    mxlen[v]=mxlen[p]+1;
    while(p&&!ch[p][f]) ch[p][f]=v,p=fail[p];
    if(!p) return fail[v]=1,void();
    int sn=ch[p][f];
    if(mxlen[sn]==mxlen[p]+1) return fail[v]=sn,void();
    int New=++cnt;
    memcpy(ch[New],ch[sn],sizeof(ch[sn]));
    mxlen[New]=mxlen[p]+1;
    fail[New]=fail[sn];
    fail[sn]=fail[v]=New;
    while(p&&ch[p][f]==sn) ch[p][f]=New,p=fail[p];
}

int rt[N<<1];
int lx,rx;
int tot;
int tag[N*40];
int ls[N*40],rs[N*40];
int mn[N<<1];
void Insert(int &v,int old,int lx,int rx,int p) {
    v=++tot;
    tag[v]=tag[old]+1;
    if(lx==rx) return ;
    int mid=lx+rx>>1;
    if(p<=mid) Insert(ls[v],ls[old],lx,mid,p);
    else Insert(rs[v],rs[old],mid+1,rx,p);
}

int Merge(int a,int b,int lx,int rx) {
    if(!a||!b) return a+b;
    int v=++tot;
    tag[v]=tag[a]+tag[b];
    if(lx==rx) return v;
    int mid=lx+rx>>1;
    ls[v]=Merge(ls[a],ls[b],lx,mid);
    rs[v]=Merge(rs[a],rs[b],mid+1,rx);
    return v;
}

int query(int v,int lx,int rx,int l,int r) {
    if(!v||lx>r||rx<l) return 0;
    if(l<=lx&&rx<=r) return tag[v];
    int mid=lx+rx>>1;
    return query(ls[v],lx,mid,l,r)+query(rs[v],mid+1,rx,l,r);
}

vector<int>e[N<<1];
void dfs(int v) {
    if(pos[v]) Insert(rt[v],rt[v],lx,rx,pos[v]);
    for(int i=0;i<e[v].size();i++) {
        int to=e[v][i];
        dfs(to);
        pos[v]=pos[to];
        rt[v]=Merge(rt[v],rt[to],lx,rx);
    }
}

int f[N],top[N];
int ans;
bool chk(int v,int f) {
    return query(rt[f],lx,rx,pos[v]-mxlen[v]+mxlen[f],pos[v])>=2;
}

void solve(int v) {
    ans=max(ans,f[v]);
    for(int i=0;i<e[v].size();i++) {
        int to=e[v][i];
        if(v==1) f[to]=1,top[to]=to;
        else {
            if(chk(to,top[v])) f[to]=f[v]+1,top[to]=to;
            else f[to]=f[v],top[to]=top[v];
        }
        solve(to);
    }
}

int main() {
    n=Get();
    scanf("%s",s+1);
    lx=1,rx=n;
    for(int i=1;i<=n;i++) Insert(s[i]-'a',i);
    for(int i=2;i<=cnt;i++) {
        e[fail[i]].push_back(i);
    }
    dfs(1);
    solve(1);
    cout<<ans;
    return 0;
}

转载于:https://www.cnblogs.com/hchhch233/p/10589048.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值