关于SA中基数排序的偷懒写法

因为我并没有学习过标准的SA写法，以前写的一直是按照自己瞎yy的n*log^2来做的

然而被卡了好多次

做法就是记录rk[]和rk2[]作为第一和第二关键字进行sort//这样做的话第一关键字有序  基数排序常数巨大

luogu模板n*log^2   70分

改成基数（巨大常数版）排序  9s AC

于是自己yy出一个不明复杂度的算法。。复杂度大概是n*log*log（n/log）

实测5s  AC

算法记录sa[i]表示排名第i的位置是哪

先给sa按照rk排个序和普通的基排没有任何区别

接下来我们考虑 

每次rk[sa[i]]都是有序的，而rk2[sa[i]]则是无序的

于是按照rk分段，每次只对rk相等的一段排序，这样每次排序的长度可以认为是n/log的？

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<ctime>
#include<cstdlib>
#include<algorithm>
#include<queue>
#define For(i,j,k)	for(int i=j;i<=k;++i)
#define Dow(i,j,k)	for(int i=k;i>=j;--i)
#define ll long long
#define inf 1e9
using namespace std;
inline int read()
{
	int t=0,f=1;char c=getchar();
	while(!isdigit(c))	{if(c=='-')	f=-1;c=getchar();}
	while(isdigit(c))	t=t*10+c-'0',c=getchar();
	return t*f;
}
inline void write(int x){if(x>=10)	write(x/10);putchar(x%10+'0');}
inline void writeln(int x){write(x);puts("");}
inline void write_p(int x){write(x);putchar(' ');}
int sa[2000001],n,rk[2000001],rk2[2000001],tmp[2000001];
char s[2000001];
inline bool cmp(int x,int y)
{
	return rk[x]==rk[y]?rk2[x]<rk2[y]:rk[x]<rk[y];
}
inline void resort(int p[])
{
    int las=1;
    For(i,2,n)
        if(rk[p[i]]!=rk[p[i-1]])//处理las~i-1
        {
            sort(p+las,p+i,cmp);
            las=i;
        }
    sort(p+las,p+n+1,cmp);
}
int main()
{
	scanf("%s",s+1);n=strlen(s+1);
	For(i,1,n)	rk[i]=s[i],sa[i]=i;
	sort(sa+1,sa+n+1,cmp);
	for(int len=1;len<=n;len*=2)
	{
		For(i,1,n)	rk2[i]=i+len<=n?rk[i+len]:0;
		resort(sa);
		int tot=1;
		tmp[sa[1]]=1;
		For(i,2,n)	if(rk[sa[i]]==rk[sa[i-1]]&&rk2[sa[i]]==rk2[sa[i-1]])	tmp[sa[i]]=tot;else	tmp[sa[i]]=++tot;
		swap(rk,tmp);
	}
	For(i,1,n)	write_p(sa[i]);
}