后缀数组dc3算法模版（待补）

模版：

const int maxn = ;

#define F(x) ((x)/3+((x)%3==1?0:tb))
#define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)
int wa[maxn],wb[maxn],wv[maxn],wss[maxn];
char s[maxn];
int r[maxn],sa[maxn];
int c0(int *r,int a,int b)
{
    return r[a]==r[b]&&r[a+]==r[b+]&&r[a+]==r[b+];
}
int c12(int k,int *r,int a,int b)
{
    if(k==) return r[a]<r[b]||r[a]==r[b]&&c12(,r,a+,b+);
    else return r[a]<r[b]||r[a]==r[b]&&wv[a+]<wv[b+];
}
void sort(int *r,int *a,int *b,int n,int m)
{
    int i;
    for(i=;i<n;i++) wv[i]=r[a[i]];
    for(i=;i<m;i++) wss[i]=;
    for(i=;i<n;i++) wss[wv[i]]++;
    for(i=;i<m;i++) wss[i]+=wss[i-];
    for(i=n-;i>=;i--) b[--wss[wv[i]]]=a[i];
    return;
}
void dc3(int *r,int *sa,int n,int m)
{
    int i,j,*rn=r+n,*san=sa+n,ta=,tb=(n+)/,tbc=,p;
    r[n]=r[n+]=;
    for(i=;i<n;i++) if(i%!=) wa[tbc++]=i;
    sort(r+,wa,wb,tbc,m);
    sort(r+,wb,wa,tbc,m);
    sort(r,wa,wb,tbc,m);
    for(p=,rn[F(wb[])]=,i=;i<tbc;i++)
        rn[F(wb[i])]=c0(r,wb[i-],wb[i])?p-:p++;
    if(p<tbc) dc3(rn,san,tbc,p);
    else for(i=;i<tbc;i++) san[rn[i]]=i;
    for(i=;i<tbc;i++) if(san[i]<tb) wb[ta++]=san[i]*;
    if(n%==) wb[ta++]=n-;
    sort(r,wb,wa,ta,m);
    for(i=;i<tbc;i++) wv[wb[i]=G(san[i])]=i;
    for(i=,j=,p=;i<ta && j<tbc;p++)
        sa[p]=c12(wb[j]%,r,wa[i],wb[j])?wa[i++]:wb[j++];
    for(;i<ta;p++) sa[p]=wa[i++];
    for(;j<tbc;p++) sa[p]=wb[j++];
    return;
}
int rank[maxn],height[maxn];
void calheight(int *r,int *sa,int n)
{
    int i,j,k=;
    for(i=;i<=n;i++) rank[sa[i]]=i;
    for(i=;i<n;height[rank[i++]]=k)
        for(k?k--:,j=sa[rank[i]-];r[i+k]==r[j+k];k++);
    return;
}
int RMQ[maxn];
int mm[maxn];
int best[][maxn];
void initRMQ(int n)
{
    int i,j,a,b;
    for(mm[]=-,i=;i<=n;i++)
        mm[i]=((i&(i-))==)?mm[i-]+:mm[i-];
    for(i=;i<=n;i++)   best[][i]=i;
    for(i=;i<=mm[n];i++)
        for(j=;j<=n+-(<<i);j++)
    {
        a=best[i-][j];
        b=best[i-][j+(<<(i-))];
        if(RMQ[a]<RMQ[b]) best[i][j]=a;
        else best[i][j]=b;
    }
    return;
}
int askRMQ(int a,int b)
{
    int t;
    t=mm[b-a+]; b-=(<<t)-;
    a=best[t][a];b=best[t][b];
    return RMQ[a]<RMQ[b]?a:b;
}
int lcp(int a,int b)
{
    int t;
    a=rank[a];b=rank[b];
    if(a>b) {t=a;a=b;b=t;}
    return height[askRMQ(a+,b)];
}

等下问问大洲rmq的一些问题要补回这里。

1.使用方法

2.初始化工作