先将询问离线,把询问按照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;
}