【LOJ#6041】事情的相似度（后缀自动机）
题面
题解
(mbox{YCB})搬了这道题目。。。(mbox{QwQ})
还是用到(lcp)就是(parent)树上的(LCA)的(len)。
每次询问显然就是区间内点的贡献。
那么考虑所有可能出现的点对。
显然对于两个子串而言，只会匹配最靠近的两个。
那么用(set)维护(endpos)集合，每次合并的时候将两个最靠近的位置合并成为一个点对，其贡献就是当前点的(len)。
那么最终询问扫描线解决即可。
(zsy)还有一种用(LCT)的做法，大致口胡一下就是类似上面的操作，要求的就是一个链并。每次将当前点到根节点的路径染色，那么每次的交点就是一个(LCA)。发现染色操作类似(LCT)的(access)，因此直接这么模拟即可。
代码是第一种做法的。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
using namespace std;
#define ll long long
#define MAX 100100
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
int n,m;
char ch[MAX];
struct Node{int son[26],len,ff;}t[MAX<<1];
int last,tot,sum;
set<int> S[MAX<<1];
void extend(int id,int c)
{
int p=last,np=++tot;last=np;
t[np].len=t[p].len+1;
while(p&&!t[p].son[c])t[p].son[c]=np,p=t[p].ff;
if(!p)t[np].ff=1;
else
{
int q=t[p].son[c];
if(t[q].len==t[p].len+1)t[np].ff=q;
else
{
int nq=++tot;
t[nq]=t[q];t[nq].len=t[p].len+1;
t[np].ff=t[q].ff=nq;
while(p&&t[p].son[c]==q)t[p].son[c]=nq,p=t[p].ff;
}
}
S[np].insert(id);
}
struct Point{int u,v,len;}p[MAX<<5];
bool operator<(Point a,Point b){return (a.v!=b.v)?a.v<b.v:((a.u!=b.u)?a.u<b.u:a.len<b.len);}
struct Line{int v,next;}e[MAX<<2];
int h[MAX<<1],cnt=1;
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
void dfs(int u)
{
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;dfs(v);if(!t[u].len)continue;
if(S[u].size()<S[v].size())swap(S[u],S[v]);
set<int>::iterator it,nw,pre,nxt;
for(it=S[v].begin();it!=S[v].end();++it)
{
S[u].insert(*it);nw=pre=nxt=S[u].find(*it);++nxt;
if(pre!=S[u].begin())--pre,p[++sum]=(Point){*pre,*nw,t[u].len};
if(nxt!=S[u].end())p[++sum]=(Point){*nw,*nxt,t[u].len};
S[u].erase(*it);
}
for(it=S[v].begin();it!=S[v].end();++it)S[u].insert(*it);
}
}
struct Qry{int l,r,id;}q[MAX];
bool operator<(Qry a,Qry b){return (a.r!=b.r)?a.r<b.r:a.l<b.l;}
int ans[MAX];
int c[MAX];
int lb(int x){return x&(-x);}
void add(int x,int w){while(x)c[x]=max(c[x],w),x-=lb(x);}
int Query(int x){int ret=0;while(x<=n)ret=max(ret,c[x]),x+=lb(x);return ret;}
int main()
{
n=read();m=read();last=tot=1;scanf("%s",ch+1);
for(int i=1;i<=n;++i)extend(i,ch[i]-48);
for(int i=2;i<=tot;++i)Add(t[i].ff,i);
dfs(1);sort(&p[1],&p[sum+1]);
for(int i=1;i<=m;++i)q[i].l=read(),q[i].r=read(),q[i].id=i;
sort(&q[1],&q[m+1]);
for(int i=1,j=1;i<=m;++i)
{
while(j<=sum&&p[j].v<=q[i].r)add(p[j].u,p[j].len),++j;
ans[q[i].id]=Query(q[i].l);
}
for(int i=1;i<=m;++i)printf("%d
",ans[i]);
return 0;
}