用 二分+哈希 求后缀数组

　　个人感觉后缀数组的板子太难背了，听了小火车讲二分+哈希可以实现求后缀数组，貌似很好理解，代码如下。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
using namespace std;
typedef unsigned long long ULL;
const ULL BASE=127;
const int maxn=6000000;
ULL base[maxn],hash[maxn];
int len;
char s[maxn];
inline ULL get_hash(int l,int r){
    return hash[r]-hash[l-1]*base[r-l+1];
}
struct S{
    int pos;
}suf[maxn];
inline int find(int l,int r,int f1,int f2){
    if(l+1>=r){
        if(get_hash(f1,f1+r-1)==get_hash(f2,f2+r-1)) return r;
        else return l;
    }
    int mid=(l+r)>>1;
    if(get_hash(f1,f1+mid-1)==get_hash(f2,f2+mid-1)) return find(mid,r,f1,f2);
    else return find(l,mid-1,f1,f2);
}
inline bool cmp(const S & x,const S & y){
    int len1=len-x.pos+1,len2=len-y.pos+1;
    int maxx=find(0,min(len1,len2)+1,x.pos,y.pos);
    if(maxx==min(len1,len2)+1) return x.pos>y.pos;
    else return s[x.pos+maxx]<s[y.pos+maxx];
}
int sa[maxn],rank[maxn],height[maxn];
int main(){
    freopen("makedata.out","r",stdin);
    freopen("mine.out","w",stdout);
    scanf("%s",s+1); len=strlen(s+1);
    base[0]=1; for(int i=1;i<=len;i++) base[i]=base[i-1]*BASE;
    for(int i=1;i<=len;i++) hash[i]=hash[i-1]*BASE+s[i]-'a'+1;
    for(int i=1;i<=len;i++) suf[i].pos=i;
    sort(suf+1,suf+len+1,cmp);
    for(int i=1;i<=len;i++){
        rank[suf[i].pos]=i;
        sa[i]=suf[i].pos;
    }
    for(int i=2;i<=len;i++){
        int len1=len-sa[i-1]+1,len2=len-sa[i]+1;
        height[i]=find(0,min(len1,len2),sa[i-1],sa[i]);
    }
    for(int i=1;i<=len;i++){
        printf("%d %d
",sa[i],height[i]);
    }
    return 0;
}