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 longconst 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\n",ask[i].ans+1LL*t*t);        }        else puts("-1");    }}