HDU 4080 Stammering Aliens（后缀数组）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4080

这个题目知道后缀数组的估计都会做，但是想AC估计要wa个一两下才能AC，其实这个题目也没什么陷阱，就是后缀数组的基本处理方法

二分枚举，但是这个要输出最后出现的位置其实还是比较简单的，但是要注意一个细节，在用height数组比较的时候可能会只计算后面的那个

的位置，忽略的第一个的位置，所以造成结果不对，这个错误还不怎么好测！其实保险一点多加一句代码就OK了！

if(height[i] >= len) 

temp++,re_pos=MAX(re_pos,sa[i-1]),re_pos=MAX(re_pos,sa[i]);

总之还是后缀数组的基本处理方法，注意点细节就可以了，不过最近几次写的后缀数组都是用的模板！

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <math.h>
#include <string.h>
using namespace std;
#define MAX(a,b) (a>b?a:b)
#define maxn 110000
#define ws ws1
int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
int cmp(int *r,int a,int b,int l)
{return r[a]==r[b]&&r[a+l]==r[b+l];}
void da(char *r,int *sa,int n,int m){
int i,j,p,*x=wa,*y=wb,*t;
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(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;
}
int rank[maxn],height[maxn];
void calheight(char *r,int *sa,int n){
int i,j,k=0;
for(i=1;i<=n;i++) rank[sa[i]]=i;//在计算height的时候顺便就把rank计算出来了，反正也要用
for(i=0;i<n;height[rank[i++]]=k)
for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
return;
}
char rec[maxn];
int sa[maxn],n,ans,num;
int is_ok(int len){//如果失败返回0，成功返回最后一个位置
    int temp=1,pos=-1,re_pos=-1;
    for(int i=2;i<=n;i++)
    if(height[i] >= len)
    temp++,re_pos=MAX(re_pos,sa[i-1]),re_pos=MAX(re_pos,sa[i]);
    else{
        if(temp >= num)
        pos=MAX(pos,re_pos);//因为不排除pos和re_pos都是符合条件的，那么如果都符合条件那么就肯定要取最大的
        temp=1,re_pos=-1;
    }
    if(temp >= num)
    pos=MAX(pos,re_pos);
    if(pos==-1)
    return 0;
    else
    return pos;
}
int find_ans(){
    int l=1,r=n-num+1,mid,pos=-1,temp;
    ans=-1;
    while(l<=r){
        mid=(l+r)>>1;
        temp=is_ok(mid);
        if(temp!=0) pos=temp,ans=mid,l=mid+1;
        else r=mid-1;
    }
    if(pos!=-1) printf("%d %d\n",ans,pos);
    else printf("none\n");
    return 0;
}
int main(){
    int i,j,k;
    while(scanf("%d",&num),num){
        scanf("%s",rec);
        n=strlen(rec);
        da(rec,sa,n+1,200);//这里是n+1 因为看这个函数里面是 < n 的
        calheight(rec,sa,n);//注意这里面是 n 了因为看函数里面是 <=n 的，所以这里要注意
        if(num==1){
            printf("%d %d\n",n,0);
            continue;
        }
        find_ans();
    }
    return 0;
}