CF gym 100962D Deep Purple [后缀树，树链剖分，线段树]

Codeforces

思路

感觉这个离线的思路好神仙啊qwq

对于每个询问\([l,r]\)其实就是要求\(p_{max}\)，使得\(lcs(s[1,p],s[1,r])>p-l\)，也就是\(lcs(s[1,p],s[1,r])+l>p\)。

首先把询问离线按\(r\)排序，然后从右往左扫，每次

1. 处理之前已经被加进去的询问，看当前位置是否能被作为\(p\)，然后把已经处理完毕的询问给删掉。
2. 把当前询问塞进去。

建出反串的后缀树，那么不等式左边就是\(dep_{lca(p,r)}+l\)。

暴力的做法就是枚举\(lca\)，然后处理询问的时候查一下这个位置的\(l\)，塞询问的时候把\(l\)塞进去。

用树链剖分加速：对于每一条链把\(tag\)打在深度最深的地方。处理询问的时候就查一下\(p\)到根上的\(\max(dep_x+l)\)以及这个\(tag\)来自哪个询问，然后将它删掉。塞询问的时候就一路往上跳，打\(tag\)。

怎么查询？分最大的\(tag\)是在上面还是下面来讨论，如果是上面就直接查，否则查一下\(\max(l)\)再加上当前的\(dep\)。

由于每个点只能有一个最大的\(tag\)放在线段树上，所以还要用一个\(set\)记录这个点有哪些标记。

（可能讲得不是很清楚qwq）

代码

#include<bits/stdc++.h>
clock_t t=clock();
namespace my_std{
    using namespace std;
    #define pii pair<int,int>
    #define fir first
    #define sec second
    #define MP make_pair
    #define rep(i,x,y) for (int i=(x);i<=(y);i++)
    #define drep(i,x,y) for (int i=(x);i>=(y);i--)
    #define go(x) for (int i=head[x];i;i=edge[i].nxt)
    #define templ template<typename T>
    #define sz 402020
    typedef long long ll;
    typedef double db;
    mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
    templ inline T rnd(T l,T r) {return uniform_int_distribution<T>(l,r)(rng);}
    templ inline bool chkmax(T &x,T y){return x<y?x=y,1:0;}
    templ inline bool chkmin(T &x,T y){return x>y?x=y,1:0;}
    templ inline void read(T& t)
    {
        t=0;char f=0,ch=getchar();double d=0.1;
        while(ch>'9'||ch<'0') f|=(ch=='-'),ch=getchar();
        while(ch<='9'&&ch>='0') t=t*10+ch-48,ch=getchar();
        if(ch=='.'){ch=getchar();while(ch<='9'&&ch>='0') t+=d*(ch^48),d*=0.1,ch=getchar();}
        t=(f?-t:t);
    }
    template<typename T,typename... Args>inline void read(T& t,Args&... args){read(t); read(args...);}
    char __sr[1<<21],__z[20];int __C=-1,__zz=0;
    inline void Ot(){fwrite(__sr,1,__C+1,stdout),__C=-1;}
    inline void print(register int x)
    {
        if(__C>1<<20)Ot();if(x<0)__sr[++__C]='-',x=-x;
        while(__z[++__zz]=x%10+48,x/=10);
        while(__sr[++__C]=__z[__zz],--__zz);__sr[++__C]='\n';
    }
    void file()
    {
        #ifdef NTFOrz
        freopen("a.in","r",stdin);
        #endif
    }
    inline void chktime()
    {
        #ifndef ONLINE_JUDGE
        cout<<(clock()-t)/1000.0<<'\n';
        #endif
    }
    #ifdef mod
    ll ksm(ll x,int y){ll ret=1;for (;y;y>>=1,x=x*x%mod) if (y&1) ret=ret*x%mod;return ret;}
    ll inv(ll x){return ksm(x,mod-2);}
    #else
    ll ksm(ll x,int y){ll ret=1;for (;y;y>>=1,x=x*x) if (y&1) ret=ret*x;return ret;}
    #endif
//	inline ll mul(ll a,ll b){ll d=(ll)(a*(double)b/mod+0.5);ll ret=a*b-d*mod;if (ret<0) ret+=mod;return ret;}
}
using namespace my_std;

int n,Q;
char s[sz];
int pos[sz];
struct hh{int l,r,id;const bool operator < (const hh &x) const {return r>x.r;}}q[sz];
int ans[sz];

int TT;
//struct SegmentTree
//{
	pii tr[sz<<2][2];
	#define ls k<<1
	#define rs k<<1|1
	#define lson ls,l,mid
	#define rson rs,mid+1,r
	void pushup(int k){tr[k][TT]=max(tr[ls][TT],tr[rs][TT]);}
	void modify(int k,int l,int r,int x,pii w)
	{
		if (l==r) return (void)(tr[k][TT]=w);
		int mid=(l+r)>>1;
		if (x<=mid) modify(lson,x,w);
		else modify(rson,x,w);
		pushup(k);
	}
	pii query(int k,int l,int r,int x,int y)
	{
		if (x<=l&&r<=y) return tr[k][TT];
		int mid=(l+r)>>1;pii ret=MP(-1e9,0);
		if (x<=mid) chkmax(ret,query(lson,x,y));
		if (y>mid) chkmax(ret,query(rson,x,y));
		return ret;
	}
//}T1,T2;
// 1：维护一段里面最大的值是多少，以及它的位置。
// 2：维护一段里面加的l的最大值是多少，以及它的位置。 

namespace Tree
{
	struct hh{int t,nxt;}edge[sz];
	int head[sz],ecnt;
	void make_edge(int f,int t){edge[++ecnt]=(hh){t,head[f]};head[f]=ecnt;}
	int dfn[sz],top[sz],bot[sz],size[sz],son[sz],fa[sz],len[sz],T;
	#define v edge[i].t
	void dfs1(int x)
	{
		size[x]=1;
		go(x)
		{
			fa[v]=x;dfs1(v);
			size[x]+=size[v];
			if (size[v]>size[son[x]]) son[x]=v;
		}
	}
	void dfs2(int x,int tp)
	{
		dfn[x]=++T;
		top[x]=tp,bot[tp]=x;
		if (son[x]) dfs2(son[x],tp);
		go(x) if (v!=son[x]) dfs2(v,v);
	}
	#undef v
	set<pii>s[sz];
	// 维护这个地方打的标记有哪些
	pii query(int x)
	{
		pii ret=MP(0,0);
		while (x)
		{
			TT=0;pii w1=::query(1,1,T,dfn[top[x]],dfn[x]);
			TT=1;pii w2=::query(1,1,T,dfn[x],dfn[bot[top[x]]]);
			w2.fir+=len[x];
			chkmax(ret,w1),chkmax(ret,w2);
			x=fa[top[x]];
		}
		return ret;
	}
	void modify(int x,int w,int id,bool add)
	{
		while (x)
		{
			if (add) s[x].insert(MP(w,id));
			else s[x].erase(MP(w,id));
			if (s[x].size())
			{
				auto cur=*s[x].rbegin();
				TT=1;::modify(1,1,T,dfn[x],cur);
				cur.fir+=len[x];
				TT=0;::modify(1,1,T,dfn[x],cur);
			}
			else TT=0,::modify(1,1,T,dfn[x],MP(0,0)),TT=1,::modify(1,1,T,dfn[x],MP(0,0));
			x=fa[top[x]];
		}
	}
}

namespace SAM
{
	struct hh{int link,len,ch[28];}a[sz];
	int lst=1,cnt=1,root=1;
	int add(int c)
	{
		int cur=++cnt,p=lst;lst=cur;a[cur].len=a[p].len+1;
		while (p&&!a[p].ch[c]) a[p].ch[c]=cur,p=a[p].link;
		if (!p) return a[cur].link=1,cur;
		int q=a[p].ch[c];
		if (a[q].len==a[p].len+1) return a[cur].link=q,cur;
		int t=++cnt;a[t]=a[q];a[q].link=a[cur].link=t;a[t].len=a[p].len+1;
		while (p&&a[p].ch[c]==q) a[p].ch[c]=t,p=a[p].link;
		return cur;
	}
	void build(){rep(i,2,cnt) Tree::make_edge(a[i].link,i),Tree::len[i]=a[i].len;}
}

int main()
{
    file();
	read(n,Q);
	cin>>(s+1);
	rep(i,1,n) pos[i]=SAM::add(s[i]-'a'+1);
	SAM::build();
	Tree::dfs1(1),Tree::dfs2(1,1);
	rep(i,1,Q) read(q[i].l,q[i].r),q[i].id=i;
	sort(q+1,q+Q+1);
	int p=1;
	drep(i,n,1)
	{
		while (233)
		{
			pii w=Tree::query(pos[i]);int id=w.sec;
			if (w.fir<=i) break;
			ans[q[id].id]=i-q[id].l+1;
			Tree::modify(pos[q[id].r],q[id].l,id,0);
		}
		while (p<=Q&&q[p].r==i) Tree::modify(pos[q[p].r],q[p].l,p,1),++p;
	}
	rep(i,1,Q) printf("%d\n",ans[i]);
	return 0;
}