POj 3581 看起来不像后缀数组的后缀数组

先离散化数据
离散化：当数据只与它们之间的相对大小有关，而与具体是多少无关时，可以进行离散化。
思路 1：输入时反转输入，因为题目要求
2：找到最小的 i 满足 sa[i] >1 把这个后缀输出
3：翻倍剩余的字符，翻倍是因为 如 2 0 1 0 1这样的串会因为长度的问题导致字典序排序失败
4：之后对翻倍的那一段数组进行一次后缀数组
5: 求出 第一个小于第一次的边界值并且不为0的下标，之后输出都到 边界位置
6：之后把剩余的前半部分输出
7： 具体看代码

#include <stdio.h>
#include <algorithm>
#include <string.h>
#define maxs 200020
using namespace std;
struct node
{
    int id;
    int ns;
}nd[maxs];

bool nodecmp(node a,node b)
{
    if(a.ns!=b.ns)
        return a.ns<b.ns;
    return a.id<b.id;
}

int wa[maxs],wb[maxs],wv[maxs],ws[maxs];

int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b]&&r[a+l]==r[b+l];
}

void getsa(int *r,int *sa,int n,int m)
{
    int i,j,p,*x=wa,*y=wb;
    for(i=0;i<m;i++) ws[i]=0;
    for(i=0;i<n;i++) ws[x[i]=r[i]]++;
    for(i=1;i<m;i++) ws[i]+=ws[i-1];
    for(i=n-1;i>=0;i--) sa[--ws[x[i]]]=i;
    for(j=1,p=1;p<n;j*=2,m=p)
    {
        for(p=0,i=n-j;i<n;i++) y[p++]=i;
        for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
        for(i=0;i<n;i++) wv[i]=x[y[i]];
        for(i=0;i<m;i++) ws[i]=0;
        for(i=0;i<n;i++) ws[wv[i]]++;
        for(i=1;i<m;i++) ws[i]+=ws[i-1];
        for(i=n-1;i>=0;i--) sa[--ws[wv[i]]]=y[i];
        for(swap(x,y),p=1,x[sa[0]]=0,i=1;i<n;i++)
            x[sa[i]]=cmp(y,sa[i],sa[i-1],j)? p-1:p++;
    }
}

int s[maxs],sa[maxs];

int main()
{
    int n,i,j;
    scanf("%d",&n);
    for(i=n-1;i>=0;i--)
    {
        scanf("%d",&nd[i].ns);
        nd[i].id=i;
    }
    sort(nd,nd+n,nodecmp);

    for(i=0;i<n;i++)
    {
        if(i&&nd[i].ns==nd[i-1].ns)
            s[nd[i].id]=s[nd[i-1].id];
        else s[nd[i].id]=i+1;
    }
    s[n]=0;
    getsa(s,sa,n+1,n+10);
    //printf("******************\n");
    for(i=1;i<=n&&sa[i]<2;i++);
    int first=sa[i];
    for(i=first;i<n;i++)
        printf("%d\n",nd[s[i]-1].ns);
    for(i=0;i<first;i++)
        s[i+first]=s[i];
    first*=2;
    s[first]=0;
    getsa(s,sa,first+1,n+10);

    for(i=1;i<first&&(sa[i]==0||sa[i]>=first/2);i++);
    int second=sa[i];
    for(i=second;i<first/2;i++)
        printf("%d\n",nd[s[i]-1].ns);

    for(i=0;i<second;i++)
        printf("%d\n",nd[s[i]-1].ns);

    return 0;
}