LOJ3298 封印

题目传送门

分析：
对于每一个左端点\(l\)，我们都能找到一个\(r\)，使得\(l \leq j \leq r\)的\(S_{l,j}\)为\(T\)的子串，\(j \geq r\)的\(S_{l,j}\)不为\(T\)的子串
当\(l\)增大时，对应的\(r\)单调不降，尝试求出每一个\(l\)对应的\(r\)
考虑建一个SAM，目前的\(S_{l,r}\)对应的子串所在节点为\(p\)，\(r\)加一位，相当于沿着字符转移边跳一步
如果目前节点子串集合最短长度（即父亲的\(len\)加一）大于了\(r-l+1\)，则\(p\)往后缀树上的父亲跳一步
复杂度是\(O(n)\)的
找到了每一个位置\(i\)对应的\(r_i\)，考虑询问\(L,R\)
需要分\(r_i \leq R\)和\(r_i>R\)分类讨论
我们对询问离线，两种情况分两棵线段树，单点修改区间取最大值可以简单维护

#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm>
#include<queue>
#include<bitset>
#include<map>
#include<set>

#define maxn 400005
#define INF 0x3f3f3f3f

using namespace std;

inline int getint()
{
	int num=0,flag=1;char c;
	while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
	while(c>='0'&&c<='9')num=num*10+c-48,c=getchar();
	return num*flag;
}

int n;
struct Sam{
	int fa,len,nxt[26];
}t[maxn];
char S[maxn],T[maxn];
int lst,tot,cnt;
int ans[maxn];
struct node{
	int l,r,op,id;
}q[maxn];
inline bool cmp(node x,node y)
{return x.r==y.r?x.op<y.op:x.r<y.r;}

inline void ins(int c)
{
	int p=lst,np=lst=++tot;
	t[np].len=t[p].len+1;
	while(p&&!t[p].nxt[c])t[p].nxt[c]=np,p=t[p].fa;
	if(!p)t[np].fa=1;
	else
	{
		int q=t[p].nxt[c];
		if(t[q].len==t[p].len+1)t[np].fa=q;
		else
		{
			int nq=++tot;
			memcpy(t[nq].nxt,t[q].nxt,sizeof t[q].nxt);
			t[nq].fa=t[q].fa,t[nq].len=t[p].len+1;
			t[q].fa=t[np].fa=nq;
			while(p&&t[p].nxt[c]==q)t[p].nxt[c]=nq,p=t[p].fa;
		}
	}
}

struct SegmentTree{
	int mx[maxn<<2];
	inline void update(int i,int l,int r,int p,int x)
	{
		if(l==r){mx[i]=x;return;}
		int mid=(l+r)>>1;
		if(p<=mid)update(i<<1,l,mid,p,x);
		else update(i<<1|1,mid+1,r,p,x);
		mx[i]=max(mx[i<<1],mx[i<<1|1]);
	}
	inline int getmx(int i,int l,int r,int ql,int qr)
	{
		if(qr<l||r<ql)return -INF;
		if(ql<=l&&r<=qr)return mx[i];
		int mid=(l+r)>>1;
		return max(getmx(i<<1,l,mid,ql,qr),getmx(i<<1|1,mid+1,r,ql,qr));
	}
}T1,T2;

int main()
{
	scanf("%s%s",S+1,T+1);lst=tot=1;
	n=strlen(T+1);
	for(int i=1;i<=n;i++)ins(T[i]-97);
	n=strlen(S+1);int r=0,p=1;
	for(int i=1;i<=n;i++)
	{
		r=max(i-1,r);
		if(r==i-1)p=1;
		while(p!=1&&t[t[p].fa].len+1>r-i+1)p=t[p].fa;
		while(r<n&&t[p].nxt[S[r+1]-97])p=t[p].nxt[S[++r]-97];
		q[++cnt]=(node){i,r,0,0};
	}
	for(int i=1;i<=n;i++)T1.update(1,1,n,i,-INF),T2.update(1,1,n,i,-i);
	int m=getint();
	for(int i=1;i<=m;i++)
	{
		int l=getint(),r=getint();
		q[++cnt]=(node){l,r,1,i};
	}
	sort(q+1,q+cnt+1,cmp);
	for(int i=1;i<=cnt;i++)
		if(!q[i].op)
		{
			T1.update(1,1,n,q[i].l,q[i].r-q[i].l+1);
			T2.update(1,1,n,q[i].l,-INF);
		}
		else ans[q[i].id]=max(T1.getmx(1,1,n,q[i].l,q[i].r),q[i].r+T2.getmx(1,1,n,q[i].l,q[i].r)+1);
	for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
}