后缀数组学习

接触到后缀数组也有一定的时间了，但是从来没有写过代码　主要是不会
然而后缀数组对于字符串处理来说确实是一大利器，所以不得不学（以前碰到字符串处理自动使用Fail树）
经过zyz大神日复一日的教♂导，蒟蒻终于对后缀数组有了一小点的了解

还是看代码吧，有注释

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;
char s[N];//原始字符串 
int x[N];//第一关键字 
int y[N];//第二关键字 
int sa[N];//排第i的是哪个后缀 
int c[N];//临时数组存储小于等于该字符的后缀个数 
int rank[N];//第i个后缀的排名 
int height[N];//和前一个后缀的最长公共前缀 
int len;
inline void makesa()
{
    int m=26;
    for(int i=0;i<m;i++) c[i]=0;
    for(int i=0;i<len;i++) c[x[i]=s[i]-'a']++;
    for(int i=1;i<m;i++) c[i]+=c[i-1];//得到小于等于该字符的后缀数量 
    for(int i=len-1;i>=0;i--) sa[--c[x[i]]]=i;//sa中从0~len-1标号,所以先减1;由于前缀和从前往后加入,所以再从后往前提取 

    for(int k=1;k<=len;k<<=1)
    {
        int p=0;
        for(int i=len-k;i<len;i++) y[p++]=i;
        for(int i=0;i<len;i++) if(sa[i]>=k) y[p++]=sa[i]-k;

        for(int i=0;i<m;i++) c[i]=0;
        for(int i=0;i<len;i++) c[x[y[i]]]++;
        for(int i=1;i<m;i++) c[i]+=c[i-1];
        for(int i=len-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];

        swap(x,y);
        m=1;x[sa[0]]=0;
        for(int i=1;i<len;i++)
          x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
        if(p>=len) break;
        m=p;
    }
}
inline void makeheight()
{
    for(int i=0;i<len;i++) rank[sa[i]]=i;
    height[0]=0;
    int k=0;
    for(int i=0;i<len;i++)
      if(rank[i])
      {
        int j=sa[rank[i]-1];
        if(k) k--;
        while(s[i+k]==s[j+k]) k++;
        height[rank[i]]=k;
      }
}