【TJOI & HEOI 2016】【BZOJ 4556】【JZOJ 4614】 游戏

Description

$Time Limits:2000ms$

Analysis

题解法是离线的啊
可是此题可以在线做
先把SA弄出来。
然后对于每个询问，二分答案，转化为判定性问题。
然后，再二分出与suffix(c)的LCP>=mid的左右边界。
这个的判定，根据那个lcp的性质，两个后缀的lcp是连续一段height数组的min
RMQ就好了。
问题变成在左右边界中是否能找到[a,b-mid+1]开始的。
这个的话，主席树就好了嘛
$O(mlog_2^2n)$
注：此题本人被卡，发现log2函数极慢，所以预处理log2，再卡卡常，强行1100ms+碾过去了。

Code

#include<cstdio>
#include<cmath>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,b,a) for(int i=b;i>=a;i--)
using namespace std;
const int N=100010,M=N*18;
int n,m,a,b,c,d,x[N],y[N],ws[N],wv[N],SA[N],rank[N],height[N];
int _2[N],lg2[N],rmq[N][20],root[N],lch[M],rch[M],sum[M];
char ch,s[N];
bool p;
int read()
{
    int t=0,p=1;
    for(ch=getchar();ch<'0' || ch>'9';ch=getchar())
        if(ch=='-') p=-1;
    for(;'0'<=ch && ch<='9';ch=getchar()) t=t*10+ch-'0';
    return t*p;
}
bool cmp(int i,int j,int l)
{
    return y[i]==y[j] && y[i+l]==y[j+l];
}
void DA()
{
    int mx=0;
    fo(i,1,n) mx=max(mx,wv[i]=x[i]=s[i]-'a'+1);
    fill(ws,ws+mx+1,0);
    fo(i,1,n) ws[wv[i]]++;
    fo(i,1,mx) ws[i]+=ws[i-1];
    fd(i,n,1) SA[ws[wv[i]]--]=i;
    for(int len=1;len<=n;len*=2)
    {
        int m=0;
        fo(i,n-len+1,n) y[++m]=i;
        fo(i,1,n)
            if(SA[i]>len) y[++m]=SA[i]-len;
        mx=0;
        fo(i,1,n) mx=max(mx,wv[i]=x[y[i]]);
        fill(ws,ws+mx+1,0);
        fo(i,1,n) ws[wv[i]]++;
        fo(i,1,mx) ws[i]+=ws[i-1];
        fd(i,n,1) SA[ws[wv[i]]--]=y[i];
        fo(i,1,n) y[i]=x[i],x[i]=0;
        x[SA[m=1]]=1;
        fo(i,2,n) x[SA[i]]=cmp(SA[i-1],SA[i],len)?m:++m;
        if(m==n) break;
    }
    fo(i,1,n) rank[SA[i]]=i;
}
void getheight()
{
    int len=0;
    fo(i,1,n)
    {
        if(len) len--;
        if(rank[i]==1) continue;
        int j=SA[rank[i]-1];
        while(i+len<=n && j+len<=n && s[i+len]==s[j+len]) len++;
        height[rank[i]]=len;
    }
}
void getrmq()
{
    fo(i,2,n) rmq[i][0]=height[i];
    fo(j,1,lg2[n])
        fo(i,2,n-_2[j]+1) rmq[i][j]=min(rmq[i][j-1],rmq[i+_2[j-1]][j-1]);
}
int minh(int x,int y)
{
    int l=lg2[y-x+1];
    int a=rmq[x][l],b=rmq[y-_2[l]+1][l];
    return a<b?a:b;
}
void build(int u,int &v,int l,int r,int k)
{
    sum[v=++m]=sum[u]+1;
    lch[v]=lch[u],rch[v]=rch[u];
    if(l==r) return;
    int mid=(l+r)>>1;
    if(k<=mid) build(lch[u],lch[v],l,mid,k);
    else build(rch[u],rch[v],mid+1,r,k);
}
void query(int u,int v,int l,int r,int x,int y)
{
    if(sum[v]<=sum[u] || p) return;
    if(l==x && r==y)
    {
        p=1;
        return;
    }
    int mid=(l+r)>>1;
    if(y<=mid) query(lch[u],lch[v],l,mid,x,y);
    else
    if(x>mid) query(rch[u],rch[v],mid+1,r,x,y);
    else
    query(lch[u],lch[v],l,mid,x,mid),query(rch[u],rch[v],mid+1,r,mid+1,y);
}
int binary(int l,int r,int k,bool bz)
{
    int mid;
    while(l<r)
    {
        mid=(l+r)>>1;
        if(!bz)
            if(minh(mid+1,rank[c])>=k) r=mid;
            else l=mid+1;
        else
            if(minh(rank[c]+1,mid+1)>=k) l=mid+1;
            else r=mid;
    }
    return l;
}
bool check(int k)
{
    int l=binary(1,rank[c],k,0),r=binary(rank[c],n,k,1);
    p=0;
    query(root[l-1],root[r],1,n,a,min(n,b-k+1));
    if(p) return 1;
    return 0;
}
int main()
{
    _2[0]=1;
    fo(i,1,20) _2[i]=_2[i-1]*2;
    int _;
    n=read(),_=read();
    fo(i,1,n) lg2[i]=log2(i);
    scanf("%s",s+1);
    DA();
    getheight();
    getrmq();
    fo(i,1,n) build(root[i-1],root[i],1,n,SA[i]);
    while(_--)
    {
        a=read(),b=read(),c=read(),d=read();
        int l=0,r=min(b-a+1,d-c+1),mid;
        while(l<r)
        {
            mid=(l+r)>>1;
            if(check(mid)) l=mid+1;
            else r=mid;
        }
        while(!check(l)) l--;
        printf("%d\n",l);
    }
    return 0;
}