因为我并没有学习过标准的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]);
}