BZOJ 3230 后缀数组+ST

思路：
首先我们已经会了后缀数组求本质不同的子串个数
这道题跟那个差不多
首先我们可以知道按字典序排好的每个后缀之前包含多少本质不同的字串
就是sigma(n-sa[i]+1-ht[i]+bi[i-1])
在这上面二分就可以求出来原串是什么了
就把两个后缀的LCP和原串的长度取个min即可
再把串倒过来也差不多这么搞一遍就好啦

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define int long long
const int N=100050;
int n,q,A[N],B[N],cntA[N],cntB[N],rk[N],ht[N],sa[N],tsa[N],f[N][20],bi[N],base[N];
char s[N];
struct Ask{int x,y,bx,by,backx,backy,lim,ans;}ask[N];
void SA(){
    memset(cntA,0,sizeof(cntA));
    for(int i=1;i<=n;i++)cntA[s[i]]++;
    for(int i='b';i<='z';i++)cntA[i]+=cntA[i-1];
    for(int i=n;i;i--)sa[cntA[s[i]]--]=i;
    rk[sa[1]]=1;
    for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);
    for(int l=1;rk[sa[n]]<n;l<<=1){
        memset(cntA,0,sizeof(cntA));
        memset(cntB,0,sizeof(cntB));
        for(int i=1;i<=n;i++)
            cntA[A[i]=rk[i]]++,
            cntB[B[i]=(i+l<=n?rk[i+l]:0)]++;
        for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1];
        for(int i=n;i;i--)tsa[cntB[B[i]]--]=i;
        for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i];
        rk[sa[1]]=1;
        for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]);
    }
    for(int i=1,j=0;i<=n;i++){
        j=j?j-1:0;
        while(s[i+j]==s[sa[rk[i]-1]+j])j++;
        ht[rk[i]]=j;
    }
}
void init_rmq(){
    for(int i=1;i<=n;i++)f[i][0]=ht[i];
    for(int j=1;j<=18;j++)
        for(int i=1;i+(1<<(j-1))<=n;i++)
            f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
int LCP(int x,int y){
    int t=base[y-x+1];
    return min(f[x][t],f[y-(1<<t)+1][t]);
}
signed main(){
    scanf("%lld%lld%s",&n,&q,s+1);
    base[0]=-1;
    for(int i=1;i<=n;i++)base[i]=base[i>>1]+1;
    SA(),init_rmq();
    for(int i=1;i<=n;i++)bi[i]=n-sa[i]+1-ht[i]+bi[i-1];
    for(int i=1;i<=q;i++){
        scanf("%lld%lld",&ask[i].x,&ask[i].y);
        if(ask[i].x>ask[i].y)swap(ask[i].x,ask[i].y);
        ask[i].bx=lower_bound(bi+1,bi+1+n,ask[i].x)-bi,
        ask[i].by=lower_bound(bi+1,bi+1+n,ask[i].y)-bi;
        int lena=n-(bi[ask[i].bx]-ask[i].x)-sa[ask[i].bx]+1,
            lenb=n-(bi[ask[i].by]-ask[i].y)-sa[ask[i].by]+1;
        ask[i].backx=n-(lena+sa[ask[i].bx]-1)+1;
        ask[i].backy=n-(lenb+sa[ask[i].by]-1)+1;
        if(ask[i].by>n)ask[i].lim=-1;
        else ask[i].lim=min(lena,lenb);
        if(~ask[i].lim){
            int t=min(ask[i].lim,ask[i].bx!=ask[i].by?LCP(ask[i].bx+1,ask[i].by):N);
            ask[i].ans+=t*t;
        }
    }
    for(int i=1;i<=n/2;i++)swap(s[i],s[n-i+1]);
    SA(),init_rmq();
    for(int i=1;i<=n;i++)bi[i]=n-sa[i]+1-ht[i]+bi[i-1];
    for(int i=1;i<=q;i++){
        ask[i].bx=lower_bound(bi+1,bi+1+n,ask[i].x)-bi,
        ask[i].by=lower_bound(bi+1,bi+1+n,ask[i].y)-bi;
        if(~ask[i].lim){
            int tempx=rk[ask[i].backx],tempy=rk[ask[i].backy];
            if(tempx>tempy)swap(tempx,tempy);
            int t=min(ask[i].lim,tempx<tempy?LCP(tempx+1,tempy):N);
            printf("%lld
",ask[i].ans+1LL*t*t);
        }
        else puts("-1");
    }
}