POJ3581 后缀数组

http://poj.org/problem?id=3581

这题说是给了N个数字组成的序列A1 A2 ..An 其中A1 大于其他的数字 , 现在要把序列分成三段并将每段分别反转求最小字典序

以后还是老老实实用基数排序 用sort wa 了 一晚

#include <iostream>
#include <algorithm>
#include <string.h>
#include <cstdio>
#include <vector>
using namespace std;
/* run this program using the console pauser or add your own getch, system("pause") or input loop */
const int maxn =;
struct SuffixArry
{
    int c[maxn*],t[maxn*],t1[maxn*],sa[maxn*];
    void build(int n,int m,int *S)
    {
         int *x=t,*y=t1,i;
         for( i=; i<m; i++)c[i]=;
         for( i=; i<n; i++)c[x[i]=S[i] ]++;
         for( i=; i<m; i++)c[i]+=c[i-];
         for( i=n-; i>=; i--)sa[ --c[ x[i] ] ]=i;
         int p;
         for(int k=; k<=n; k<<= )
         {
            p=;
            for( i = n-k; i<n; i++)y[p++]=i;
            for( i=; i<n; i++) if(sa[i]>=k)y[p++]=sa[i]-k;
            for( i=; i<m; i++)c[ i ] = ;
            for( i= ; i<n; i++ )c[ x[ y[i] ] ]++;
            for( i=; i<m; i++) c[i]+=c[i-];
            for( i=n-; i>=; i--) sa[--c[ x[ y[i] ] ] ]=y[i];
            swap(x,y);
            p=; x[ sa[  ] ] = ;
            for(int i=; i<n; i++)
             x[ sa[i] ] = y[ sa[i] ]== y[ sa[i-] ]&& y[ sa[i]+k ]==y[sa[i-]+k]?p-:p++;
            if(p>=n)break;
            m=p;
         }
    }
    void clear(){
        memset(t,-,sizeof(t));
        memset(t1,-,sizeof(t1));
    }
}T;
int A[maxn],B[maxn],C[maxn],rev[maxn*];
void revercopy(int *a, int *b, int len)
{
     for(int i=; i<len; i++)
     {
          b[i]=a[len--i];
     }
}
void rever(int *a, int len)
{
      for(int i=; i<len/; i++ )
      {
           int t= a[i];
           a[i]=a[len--i];
           a[len--i]=t;
      }
}
int main(int argc, char *argv[]) {

     int n;
     scanf("%d",&n);
     for(int i=; i<n; i++)
        {
          scanf("%d",&A[i]);
          B[i]=A[i];
        }
     sort(B,B+n);
     int L=unique(B,B+n)-B;
     for(int i=; i<n; i++)
         {
            C[i]=lower_bound(B,B+L,A[i])-B;
         }
     int p1;
     T.clear();
     revercopy(C,rev,n);
     T.build(n,L,rev);

      for(int i=; i<n; i++)
      {
           p1 = n-T.sa[i];
           if(p1>&&n-p1>=)break;
      }

      if(p1<||n-p1<)while(true){};
      int m= n-p1;
      revercopy(C+p1,rev,m); //for(int i=0; i<m; i++)printf("%d ",rev[i]);
      revercopy(C+p1,rev+m,m);//for(int i=0; i<m*2; i++)printf("%d ",rev[i]);
      T.build(m*,L,rev);
      int p2;
      for(int i=; i<m*; i++)
      {
          p2 = p1 + m - T.sa[ i ];
          if(p2>p1&&n-p2>=)break;
      }
      if(p2<=p1||n-p2<)while(true){
      };
      rever(A,p1);
      rever(A+p1,p2-p1);
      rever(A+p2,n-p2);
      for(int i=; i<n; i++) printf("%d\n",A[i]);
    return ;
}