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<cstring>
#include<iostream>
#define M 100003
using namespace std;
char s1[M],s2[M];
int l1,l2,len,id[M<<1],sa[M<<1],ranks[M<<1],tmp[M<<1],w[M<<1],height[M<<1],cnt[M<<1];
void SA(int len,int up)
{
    int *rk=ranks,*t=tmp,d=1,p=0;
    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[w[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>=0) 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(rk,t);
        p=1;
        rk[sa[0]]=0;
        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
                rk[sa[i+1]]=p++;
        if (p==len) break;
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX