后缀数组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.初始化工作