标签：last fa Reincarnation memset len 后缀 int HDU4622 sizeof

Description

给定一个长度不超过 \(2 \times 10^3\) 的字符串，有 \(q \le 10^4\) 个询问，每个询问给出一个区间，问区间中共有多少个不同的子串。

Solution

枚举一个后缀，增量构造 SAM，即可得到它所有子串中含有的不同子串的数目

在一次 sam.extend() 过程中，新增的本质不同子串数目为 \(len[last]-len[fa[last]]\)，其中 \(cnt\) 是 \(endpos\) 集合的大小

#include <bits/stdc++.h>
using namespace std;
const int N = 4005;

struct SAM
{
    int len[N], ch[N][26], fa[N], ind, last;
    int t[N], a[N], cnt[N], f[N];
    SAM()
    {
        ind = last = 1;
    }
    void clear()
    {
        memset(len,0,sizeof len);
        memset(ch,0,sizeof ch);
        memset(fa,0,sizeof fa);
        memset(t,0,sizeof t);
        memset(a,0,sizeof a);
        memset(cnt,0,sizeof cnt);
        memset(f,0,sizeof f);
        ind = last = 1;
    }
    inline void extend(int id)
    {
        int cur = (++ ind), p;
        len[cur] = len[last] + 1;
        cnt[cur] = 1;
        for (p = last; p && !ch[p][id]; p = fa[p]) ch[p][id] = cur;
        if (!p) fa[cur] = 1;
        else
        {
            int q = ch[p][id];
            if (len[q] == len[p] + 1) fa[cur] = q;
            else
            {
                int tmp = (++ ind);
                len[tmp] = len[p] + 1;
                for(int i=0; i<26; i++) ch[tmp][i] = ch[q][i];
                fa[tmp] = fa[q];
                for (; p && ch[p][id] == q; p = fa[p]) ch[p][id] = tmp;
                fa[cur] = fa[q] = tmp;
            }
        }
        last = cur;
    }
    int extend(char ch)
    {
        extend(ch-'a');
        return len[last]-len[fa[last]];
    }
} sam;

int ans[N/2][N/2];

void solve()
{
    sam.clear();

    string str;
    cin>>str;

    memset(ans,0,sizeof ans);

    int m=str.length();

    for(int i=1;i<=m;i++)
    {
        sam.clear();
        for(int j=i;j<=m;j++)
        {
            ans[i][j]=ans[i][j-1]+sam.extend(str[j-1]);
        }
    }

    int n;
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        int p,q;
        cin>>p>>q;
        cout<<ans[p][q]<<endl;
    }
}

signed main()
{
    ios::sync_with_stdio(false);

    int t;
    cin>>t;
    while(t--) solve();
}

标签：last,fa,Reincarnation,memset,len,后缀,int,HDU4622,sizeof
来源： https://www.cnblogs.com/mollnn/p/13583496.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。