1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#define M 100005
using namespace std;
int len;
char s[M];
int w[M],sa[M],ranks[M],tmp[M],id[M],cnt[M],height[M],ST[17][M],sum[M];
void SA(int len,int up)
{
    int d=1,p=0,*rk=ranks,*t=tmp;
    for (int i=0;i<len;++i) cnt[rk[i]=w[i]]++;
    for (int i=1;i<up;++i) cnt[i]+=cnt[i-1];
    for (int i=len-1;i>=0;--i) sa[--cnt[rk[i]]]=i;
    for (;;)
    {
        for (int i=len-d;i<len;++i) id[p++]=i;
        for (int i=0;i<len;++i)
            if (sa[i]>=d) id[p++]=sa[i]-d;
        for (int i=0;i<up;++i) cnt[i]=0;
        for (int i=0;i<len;++i) ++cnt[t[i]=rk[id[i]]];
        for (int i=1;i<up;++i) cnt[i]+=cnt[i-1];
        for (int i=len-1;i>=0;--i) sa[--cnt[t[i]]]=id[i];
        swap(t,rk);
        rk[sa[0]]=0;
        p=1;
        for (int i=0;i<len-1;++i)
            if (sa[i]+d<len&&sa[i+1]+d<len&&t[sa[i]]==t[sa[i+1]]&&t[sa[i]+d]==t[sa[i+1]+d])
                rk[sa[i+1]]=p-1;
            else
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX