bzoj2251

最新推荐文章于 2019-09-07 11:17:00 发布

原创最新推荐文章于 2019-09-07 11:17:00 发布 · 301 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

乱搞专栏收录该内容

4 篇文章

订阅专栏

本文介绍了一种求解后缀数组的算法，时间复杂度为O(n^2+nlogn)，并通过一种称为“乱搞”的方法来寻找最优解。尽管整体方法看似简单粗暴，但针对特定的数据集却能有效地找到答案。

乱搞

前面求后缀数组是O(n^2+nlog n)的，后面求答案是乱搞，貌似可以卡到O(n^3)，不过数据很水

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
typedef long long LL;
inline int read()
{
	int x=0;bool f=0;char c=getchar();
	for (;c<'0'||c>'9';c=getchar()) f=c=='-'?1:0;
	for (;c>='0'&&c<='9';c=getchar()) x=x*10+c-'0';
	return f?-x:x;
}
const int N=3010;
int n,lcp[N][N],sa[N],rk[N];
char str[N];
inline bool cmp(int a,int b) {return str[a+lcp[a][b]]<str[b+lcp[a][b]];}

int main()
{
	n=read();scanf("%s",str+1);
	for (int i=n;i;i--)
		for (int j=i;j;j--)
			lcp[j][i]=lcp[i][j]=(str[i]==str[j])?lcp[i+1][j+1]+1:0;
	for (int i=1;i<=n;i++) sa[i]=i;
	sort(sa+1,sa+n+1,cmp);
	for (int i=1;i<=n;i++) rk[sa[i]]=i;
	for (int i=1;i<=n;i++)
		for (int j=lcp[sa[i-1]][sa[i]]+1;j<=n-sa[i]+1;j++)
		{
			int rec=0;
			for (int k=i;k<=n&&lcp[sa[i]][sa[k]]>=j;k++,rec++);
			if (rec>1) printf("%d\n",rec);
		}
	return 0;
}