CodeForces - 1037H: Security（SAM+线段树合并）

题意：给定字符串S；  Q次询问，每次询问给出(L,R,T)，让你在S[L,R]里面找一个字典序最小的子串，其字典序比T大。 没有则输出-1；

思路：比T字典序大，而且要求字典最小，显然就是在T的尾巴加一个很小的字符，如果不存在，则依次删去尾巴，直到“存在”。

而“存在”是指，前缀lim+一个字符‘x’后存在于[L+lim,R]中， 所以我们需要预处理出所有点的endpos，这个用线段树合并就可以了。

注意这里是每一个节点都要记录，所以和普通的启发式合并不太一样，这里需要每次新开一个节点。

#include<bits/stdc++.h>
#define rep2(i,a,b) for(int i=a;i>=b;i--)
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
char c[maxn]; int rt[maxn],N,tot;
struct in{
    int L,R,sum;
}s[maxn*];
void ins(int &now,int L,int R,int p)
{
    if(!now) now=++tot; s[now].sum++;
    if(L==R) return ; int Mid=(L+R)>>;
    if(p<=Mid) ins(s[now].L,L,Mid,p);
    else ins(s[now].R,Mid+,R,p);
}
int merge(int x,int y)
{
    if(!x||!y) return x|y;
    int now=++tot;
    s[now].L=merge(s[x].L,s[y].L);
    s[now].R=merge(s[x].R,s[y].R);
    s[now].sum=s[s[now].L].sum+s[s[now].R].sum;
    return now;
}
bool query(int Now,int L,int R,int l,int r)
{
    if(!Now||l>r) return false;
    if(l<=L&&r>=R) return s[Now].sum;
    int Mid=(L+R)>>;
    if(l<=Mid&&query(s[Now].L,L,Mid,l,r)) return true;
    if(r>Mid) return query(s[Now].R,Mid+,R,l,r);
    return 0;
}
struct SAM{
    int ch[maxn][],fa[maxn],maxlen[maxn],cnt,last;
    int a[maxn],b[maxn];
    void init()
    {
        cnt=last=;
        memset(ch[],,sizeof(ch[]));
    }
    int add(int x,int id)
    {
       int np=++cnt,p=last; last=np;
       ins(rt[np],,N,id);
       maxlen[np]=maxlen[p]+; memset(ch[np],,sizeof(ch[np]));
       while(p&&!ch[p][x]) ch[p][x]=np,p=fa[p];
       if(!p) fa[np]=;
       else {
         int q=ch[p][x];
         if(maxlen[q]==maxlen[p]+) fa[np]=q;
         else {
            int nq=++cnt; maxlen[nq]=maxlen[p]+;
            fa[nq]=fa[q]; fa[q]=fa[np]=nq;
            memcpy(ch[nq],ch[q],sizeof(ch[q]));
            while(p&&ch[p][x]==q) ch[p][x]=nq,p=fa[p];
         }
       }
    }
    void sort()
    {
        rep(i,,cnt) a[i]=;
        rep(i,,cnt) a[maxlen[i]]++;
        rep(i,,cnt) a[i]+=a[i-];
        rep(i,,cnt) b[a[maxlen[i]]--]=i;
        for(int i=cnt;i>=;i--) {
            int x=b[i];
            if(fa[x]) rt[fa[x]]=merge(rt[fa[x]],rt[x]);
        }
    }
}T;
int pos[maxn],Q;
void solve()
{
    int L,R,len,fcy;
    scanf("%d%d%s",&L,&R,c+);
    len=strlen(c+);
    int now=,lim=len; c[len+]='a'-;
    rep(i,,len) {
        if(T.ch[now][c[i]-'a']) now=T.ch[now][c[i]-'a'];
        else { lim=i-; break;}
        pos[i]=now;
    }
    fcy=-; pos[]=;
    for(int i=lim;i>=&&fcy==-;i--){
        for(int j=c[i+]-'a'+;j<;j++){
            int v=T.ch[pos[i]][j];
            if(!v) continue;
            if(query(rt[v],,N,L+i,R)) {
                fcy=i+; c[fcy]='a'+j;
                break;
            }
        }
    }
    if(fcy==-) puts("-1");
    else { rep(i,,fcy) putchar(c[i]); puts(""); }
}
int main()
{
    T.init();
    scanf("%s",c+); N=strlen(c+);
    rep(i,,N) T.add(c[i]-'a',i);
    T.sort();
    scanf("%d",&Q);
    while(Q--) solve();
    return ;
}