关于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
#include
#include
#include
#include
#include
#include
#include
#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]

你可能感兴趣的:(算法,其他)