hdu 5442 (后缀数组)

稍微学习了下第一次用后缀数组- - ， 强行凑出答案 ， 感觉现在最大的问题是很多算法都不知道 ，导致有的题一点头绪都没有（就像本题）。  /*推荐 《后缀数组——处理字符串的有力工具》——罗穗骞

后缀数组sa，  sa[ i ] = a表示字符串从第a个开始到结尾的字典序排序为i

本题是绕成了一个环，所以我将字符串重复一遍再用后缀数组模板跑一次。

顺序的话，取最后一个即可

反序的话，因为相同情况下，要求取得尽可能小。

hight求的是i和i-1的最长公共前缀，如果串i-1整个是公共最长前缀，即说明他们的字典序本该是一样的，只是在双倍的情况下i多算一部分，所以有先取i-1

//思路应该没问题，代码是在不停wa中莫名改出来的- -

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<string>
#include<map>
#include<set>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
typedef long long ll;

const int maxn=200000+10;

int wa[maxn],wb[maxn],wv[maxn],wss[maxn],rak[maxn];
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void da(int *r,int *sa,int n,int m)
{
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0; i<m; i++) wss[i]=0;
    for(i=0; i<n; i++) wss[x[i]=r[i]]++;
    for(i=1; i<m; i++) wss[i]+=wss[i-1];
    for(i=n-1; i>=0; i--) sa[--wss[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++) wss[i]=0;
        for(i=0; i<n; i++) wss[wv[i]]++;
        for(i=1; i<m; i++) wss[i]+=wss[i-1];
        for(i=n-1; i>=0; i--) sa[--wss[wv[i]]]=y[i];
        for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1; i<n; i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
    }
    return ;
}

char s[2*maxn],rev[2*maxn];
int sa[maxn],r[maxn],height[maxn];
int n;
int p0,p1;
void build_height(int n)
{
    int i,k=0;
    for(i=0; i<n; i++)rak[sa[i]]=i;
    for(i=0; i<n; i++)
    {
        if(k)k--;
        int j=sa[rak[i]-1];
        while(s[i+k]==s[j+k])k++;
        height[rak[i]]=k;
    }
}

void solve(int x)
{
    if(!x)
    {
        int t = 0;
        memset(sa,0,sizeof(sa));
        for(int i = 0; i < n; i++)
            r[t++]=s[i]-'a'+1;
        r[t++] = 0;

        da(r,sa,t,30);
        p0 = sa[t-1]+1;

    }
    else
    {
        int t = 0;
        memset(sa,0,sizeof(sa));
        for(int i = 0; i < n; i++)
            r[t++]=rev[i]-'a'+1;

        r[t++] = 0;
        da(r,sa,t,30);
        build_height(t);
        int k = t-2;
        while(k>=4&&height[k]&&height[k]==n-sa[k-1])k--;
        p1 = sa[k];

    }
}


int cmp(char*s,char*t,int p,int q)
{
    for(int i=0; i<n; i++)
        if(s[(p+i)%n]!=t[(q+i)%n])
            return s[(p+i)%n]>t[(q+i)%n];
    return 2;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        scanf("%s",s);
        for(int i = 0; i < n; i++)
            s[n+i] = s[i];
        n = 2*n;
        for(int i=0; i<n; i++)
            rev[i]=s[n-1-i];


        solve(0);
        solve(1);
        int ans,dir;
        int tmp = p1;
       // printf("%d %d
",p0,tmp);
        p1 = n-p1;
        n/=2;
        int d=cmp(s,rev,p0-1,tmp);


        //printf("%d %d
",p0,tmp);

        if(d==2)
        {
            if(p0<=p1)ans=p0,dir=0; 
            else ans=p1,dir=1;
        }
        else if(d>0)ans=p0,dir=0;
        else ans=p1,dir=1;
        printf("%d %d
",ans,dir);
    }
    return 0;
}