【学习笔记】CF700E Cool Slogans

仰望星空的蚂蚁

已于 2023-09-29 15:02:59 修改

阅读量361

点赞数

文章标签：学习笔记算法

于 2023-06-30 11:18:21 首次发布

本文链接：https://blog.csdn.net/cqbzlydd/article/details/131365755

版权

文章讨论了如何利用SAM（SuffixTree）结构解决字符串问题，强调了后缀的性质在最优解中的重要性。通过在后缀树上进行动态规划，可以实现线性时间复杂度的解决方案。文章提到了在SAM上维护信息的难点，并给出了一段相关的C++代码示例。

摘要由CSDN通过智能技术生成

服了，做这道题还要重新学一遍 $S A M$ ，泪目了。

先不考虑交叉的情况，简单画一个图可以发现最优方案中 $s_{i-1}$ 一定是 $s_i$ 的后缀。因此，考虑在后缀树上统计答案。

比较让人纠结的地方在于，同一等价类对应的答案相同吗？这就涉及到我们对 $S A M$ 基本性质的了解了。抛出一个结论：设 $s$ 是节点 $u$ 对应的最长串， $v$ 是 $u$ 的祖先，那么 $v$ 对应的所有字符串在 $s$ 中的匹配情况相同。证明的关键在于利用 $\text{endpos}$ 集合相同。感觉这个证明脑补一下应该不难就略过了。这也说明了一件事，对后缀性质的充分理解才是解决问题的关键。

最后是关于实现。直接暴力二分的复杂度是 $O(n\log^2 n)$ 。但是因为是在树上所以能做的操作还是比较多的，比如说这道题就可以根据父亲的信息来计算，注意到一条链上的 $D P$ 值是连续的，并且相邻两个位置的 $D P$ 值不超过 $1$ ，因此可以做到线性。这个地方的处理还是比较巧妙的。复杂度 $O(n\log n)$ 。

就是模板有点长。

打了一发代码，感觉最恶心的地方还是在于在 $S A M$ 上维护一些比较繁琐的信息。如果对 $S A M$ 不熟悉的话还是很容易搞错的。

#include<bits/stdc++.h>
#define ll long long
#define fi first
#define se second
#define pb push_back
#define db double
#define inf 0x3f3f3f3f3f3f3f3f
#define uint unsigned int
using namespace std;
const int N=4e5+5;
const int M=N*20;
int n,tot,last,rt[N],nums[N],cnt,degree[N],ls[M],rs[M],totnodes,dp[N],fr[N],res;
string s;
queue<int>Q;
struct node{
    int link,len,ch[26],pos;
}t[N];
void ins(int &rt,int l,int r,int x){
    if(!rt)rt=++totnodes;
    if(l==r)return;
    int mid=l+r>>1;
    x<=mid?ins(ls[rt],l,mid,x):ins(rs[rt],mid+1,r,x);
}
void extend(int c,int x){
    int cur=++tot;
    t[cur].len=t[last].len+1;
    t[cur].pos=x;
    int p=last;
    while(p!=-1&&!t[p].ch[c]){
        t[p].ch[c]=cur;
        p=t[p].link;
    }
    if(p==-1){
        t[cur].link=0;
    }
    else{
        //assert(t[p].ch[c]);
        int q=t[p].ch[c];
        assert(q);
        if(t[q].len==t[p].len+1){
            t[cur].link=q;
        }
        else{
            int clone=++tot;
            t[clone].link=t[q].link,t[clone].len=t[p].len+1;for(int i=0;i<26;i++)t[clone].ch[i]=t[q].ch[i];t[clone].pos=x;
            t[cur].link=clone;
            t[q].link=clone;
            while(p!=-1&&t[p].ch[c]==q){
                t[p].ch[c]=clone;
                p=t[p].link;
            } 
        }
    }
    last=cur;
    ins(rt[cur],0,n-1,x);
}
int merge(int p,int q,int l,int r){
    if(l>r)return 0;
    if(!p||!q)return p+q;
    int o=++totnodes,mid=l+r>>1;
    ls[o]=merge(ls[p],ls[q],l,mid);
    rs[o]=merge(rs[p],rs[q],mid+1,r);
    return o;
}
int query(int p,int l,int r,int ql,int qr){
    if(!p)return 0;
    if(ql<=l&&r<=qr)return 1;
    int mid=l+r>>1;
    if(ql<=mid&&query(ls[p],l,mid,ql,qr))return 1;
    if(mid<qr&&query(rs[p],mid+1,r,ql,qr))return 1;
    return 0;
}
signed main(){
    //freopen("data.in","r",stdin);
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
	cin>>n>>s;t[0].link=-1;
	for(int i=0;i<n;i++){
        extend(s[i]-'a',i);
    }
    for(int i=1;i<=tot;i++)degree[t[i].link]++;
    for(int i=1;i<=tot;i++){
        if(degree[i]==0){
            Q.push(i);
        }
    }
    while(Q.size()){
        int u=Q.front(),v=t[u].link;Q.pop();
        nums[++cnt]=u;
        if(v){
            rt[v]=merge(rt[v],rt[u],0,n-1);
            if(--degree[v]==0)Q.push(v);
        }
    }for(int i=1;i<=tot;i++)assert(degree[i]==0);
    assert(cnt==tot);
    for(int i=cnt;i>=1;i--){
        int x=nums[i],y=t[x].link;
        if(fr[y]==0||query(rt[fr[y]],0,n-1,t[x].pos-t[x].len+t[fr[y]].len,t[x].pos-1)){
            dp[x]=dp[fr[y]]+1;
            fr[x]=x;
        }
        else{
            dp[x]=dp[fr[y]];
            fr[x]=fr[y];
        }
        res=max(res,dp[x]);
    }
    cout<<res;
}