【LOJ6041】「雅礼集训 2017 Day7」事情的相似度（SAM）（LCT）（树状数组）

传送门

---

题面太暴力了

题解：

两个前缀的LCS就是SAM的fail树上的LCA。

将询问按照$r$排序。暴力处理询问的方式就是跳fail树看有没有之前打上的标记。然后求一个后缀最大值。

利用LCT的access来进行打标记，在切换儿子的时候将最大值扔到树状数组里面即可。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

namespace IO{
	inline char gc(){
		static cs int Rlen=1<<22|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	
	template<typename T>
	inline T get(){
		char c;T num;
		while(!isdigit(c=gc()));num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	}
	inline int gi(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

cs int N=2e5+15;

int n,m;
char s[N];int ans[N];
struct Query{int l,id;};
std::vector<Query> qry[N];

namespace SAM{
	int son[N][2],fa[N],len[N],now;
	inline void push_back(int i,int c){assert(0<=c&&c<=1);
		int cur=i,p=(i-1)?(i-1):(n+1);
		len[cur]=len[p]+1;
		for(;p&&!son[p][c];p=fa[p])son[p][c]=cur;
		if(!p)fa[cur]=n+1;
		else if(len[son[p][c]]==len[p]+1)fa[cur]=son[p][c];
		else {
			int nq=++now,q=son[p][c];
			memcpy(son[nq],son[q],sizeof son[q]);
			len[nq]=len[p]+1;fa[nq]=fa[q];
			fa[q]=fa[cur]=nq;
			for(;p&&son[p][c]==q;p=fa[p])son[p][c]=nq;
		}
	}
}

int tr[N];
inline void add(int p,int v){for(;p>=1;p^=p&-p)tr[p]=std::max(tr[p],v);}
inline int query(int p){int res=0;for(;p<=n;p+=p&-p)res=std::max(res,tr[p]);return res;}

namespace LCT{
	int fa[N],son[N][2],ps[N],val[N],tag[N];
	inline void push(int u,int v){ps[u]=v,tag[u]=v;}
	inline bool isrt(int u){return son[fa[u]][0]!=u&&son[fa[u]][1]!=u;}
	inline bool which(int u){return son[fa[u]][1]==u;}
	inline void Rotate(int u){
		int p=fa[u],pp=fa[p],d=which(u);
		if(!isrt(p))son[pp][which(p)]=u;
		son[p][d]=son[u][!d];if(son[u][!d])fa[son[p][d]]=p;
		son[u][!d]=p;fa[u]=pp,fa[p]=u;
	}
	inline void pushdown(int u){
		if(tag[u]){
			if(son[u][0])push(son[u][0],tag[u]);
			if(son[u][1])push(son[u][1],tag[u]);
			tag[u]=0;
		}
	}
	inline void Splay(int u){
		static int q[N],qn;q[qn=1]=u;
		for(int re p=u;!isrt(p);p=fa[p])q[++qn]=fa[p];
		while(qn)pushdown(q[qn--]);
		for(int re p=fa[u];!isrt(u);Rotate(u),p=fa[u])
		if(!isrt(p))Rotate(which(u)==which(p)?p:u);
	}
	inline void access(int u,int id){
		int ch=0;for(;u;u=fa[ch=u])
		Splay(u),son[u][1]=ch,add(ps[u],val[u]);
		push(ch,id);
	}
}

signed main(){
#ifdef zxyoi
	freopen("similar.in","r",stdin);
#endif
	scanf("%d%d%s",&n,&m,s+1);SAM::now=n+1;
	for(int re i=1;i<=n;++i)SAM::push_back(i,s[i]^48);
	for(int re i=1;i<=SAM::now;++i)LCT::val[i]=SAM::len[i],LCT::fa[i]=SAM::fa[i];
	for(int re i=1;i<=m;++i){
		int l=gi(),r=gi();
		qry[r].push_back((Query){l,i});
	}
	for(int re i=1;i<=n;++i){
		LCT::access(i,i);
		for(auto t:qry[i])ans[t.id]=query(t.l);
	}
	for(int re i=1;i<=m;++i)cout<<ans[i]<<"\n";
	return 0;
}