HDU 5442 Favorite Donut 后缀数组

题意：给你一个字符串（环形），求字典序最大的起始位置。

按字典序、起始位置、方向的优先级顺序依次判断。

方法：倍长字符串长度，构建后缀数组。

注意：顺时针时在字符串末尾补0(无穷小）;逆时针时在字符串末尾补'z'+1（无穷大）.

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <set>
#include <algorithm>
#include <math.h>
using namespace std;
typedef long long ll;
const int N=1e5;
char s[N];
int sa[N],height[N],Rank[N],tmp[N],top[N];
void makesa(int n){
	int i,j,len,na;
	na=(n<256?256:n);
	memset(top,0,na*sizeof(int));
	for(i=0;i<n;i++)top[Rank[i]=s[i]&0xff]++;
	for(i=1;i<na;i++)top[i]+=top[i-1];
	for(i=0;i<n;i++)sa[--top[Rank[i]]]=i;
	for(len=1;len<n;len<<=1){
		for(i=0;i<n;i++){
			j=sa[i]-len;if(j<0)j+=n;
			tmp[top[Rank[j]]++]=j;
		}
		sa[tmp[top[0]=0]]=j=0;
		for(i=1;i<n;i++){
			if(Rank[tmp[i]]!=Rank[tmp[i-1]]||
				Rank[tmp[i]+len]!=Rank[tmp[i-1]+len])
				top[++j]=i;
			sa[tmp[i]]=j;
		}
		memcpy(Rank,sa,n*sizeof(int));
		memcpy(sa,tmp,n*sizeof(int));
		if(j>=n-1)break;
	}
	//for(int i=0;i<n;i++)printf("%c",s[i]);
	//printf(",,");
	//for(int i=0;i<n;i++)printf("(%d,%d),",i,Rank[i]);
	//printf("\n");
}
//void lcp(int n){
//	int i,j,k;
//	for(j=Rank[height[i=k=0]=0];i<n-1;i++,k++)
//		while(k>=0&s[i]!=s[sa[j-1]+k])
//			height[j]=(k--),j=Rank[sa[j]+1];
//}
char a[N];
int solve(int x,int y,int n){
	int px,py;
	px=x;py=y;
	y+=n;
	for(int i=0;i<n;i++){
		if(a[x]>a[y])return 0;
		if(a[x]<a[y])return 1;
		x++;y--;
	}
	if(px<py)return 0;
	if(py<px)return 1;
	return 0;
}
void work(){
	int T;
	scanf("%d",&T);
	while(T--){
		int n;
		scanf("%d",&n);
		scanf("%s",s);
		for(int i=0;i<n;i++){
			s[i+n]=s[i];
			a[i+n]=a[i]=s[i];
		}
		s[n+n]=0;
		makesa(n+n+1);
		//lcp(n+n);
		int x,y;
		x=0;
		for(int i=1;i<n;i++){
			if(Rank[i]>Rank[x])
				x=i;
		}
		for(int i=0;i<n;i++){
			s[i]=s[i+n]=a[n-1-i];
		}
		s[n+n]='z'+1;
		makesa(n+n+1);
		//lcp(n+n);
		y=n-1;
		for(int i=n-2;i>=0;i--){
			if(Rank[i]>Rank[y])
				y=i;
		}
		y=n-1-y;
		//printf("(%d,%d)",x+1,y+1);
		if(solve(x,y,n)==0)printf("%d 0\n",x+1);
		else printf("%d 1\n",y+1);
	}
}
int main(){
    freopen("data.in","r",stdin);
    work();
    return 0;
}