[HAOI2017]供给侧改革

最新推荐文章于 2019-07-31 21:59:00 发布

原创最新推荐文章于 2019-07-31 21:59:00 发布 · 1.5k 阅读

0 ·

CC 4.0 BY-SA版权

本文介绍了一种高效的字符串匹配算法，利用字典树结构和LCP（最长公共前缀）概念来解决批量字符串匹配的问题。通过预处理和巧妙的数据结构设计，算法能够在O(40N)的时间复杂度内完成所有询问，适用于需要频繁进行字符串比较的应用场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

先将询问离线，把询问按照r分类，从小到大回答每一个询问。
由随机的性质得知两个串的lcp不会特别长，我们令其最长为40
对于区间[1,R],[2,R],[3,R]…[R-1,R]，求其中最长的LCP一定是单调不增的
从左到右将后缀的前40位插入字典树，字典树路径上存下最后一次访问该点的时间，第二次访问这个深度为dep节点的时候就能得到一对时间x,y表示以x开头的后缀和以y开头的后缀的LCP长度大于等于dep
维护这40个值，对于每一个询问扫一遍这40个值即可
时空复杂度 $O(40N)$

#include <bits/stdc++.h>
#define N 100050
#define tp 50
using namespace std;
typedef long long LL;

int o[tp+5],F[tp*N],ch[tp*N][2],tot,n,q;
LL ans[N];
char s[N];
struct HbFS{int l,r,_;}Q[N];
inline int rd() {int r;scanf("%d",&r);return r;}

bool cmp(HbFS p1, HbFS p2) {return p1.r < p2.r;}

void add(int x) {
    int p = 0;
    for (int i=0;i+x<=n && i+1<=tp;i++) {
        int c = s[i+x]-'0';
        if (!ch[p][c]) {
            ch[p][c] = ++tot;
            F[tot] = x;
        } else {
            o[i+1] = max(F[ ch[p][c] ], o[i+1]);
            F[ ch[p][c] ] = x;
        }
        p = ch[p][c];
    }
}

int main() {
    freopen("b.in","r",stdin);
    freopen("b.out","w",stdout);
    n = rd(), q = rd();
    scanf("%s",s+1);
    for (int _=1;_<=q;_++) Q[_].l = rd(), Q[_].r = rd(), Q[_]._ = _;
    sort(Q+1,Q+q+1,cmp);

    int T = 1;
    for (int _=1;_<=n;_++) {
        add(_);
        while (T<=q && Q[T].r == _) {
            for (int i=1;i<=tp;i++)
                if (o[i] >= Q[T].l)
                    ans[ Q[T]._ ] += 1LL * i * (o[i] - max(Q[T].l-1, o[i+1]));
                else
                    break;
            T++;
        }
    }

    for (int _=1;_<=q;_++) printf("%lld\n",ans[_]);

    return 0;
}