蓝桥杯 2014A组决赛第五题 重复模式 （后缀数组）

作为 drd 的好朋友，技术男 atm 在 drd 生日时送给他一个超长字符串 S 。atm 要 drd 在其中找出一个最长的字符串 T ，使得 T 在 S 中至少出现了两次，而他想说的秘密就藏在 T 中。
由于字符串实在是太长了，drd 总是找不到合适的 T 。于是 drd 请你帮他找到这个 T 的长度。
【输入格式】
一行。一个字符串，即题目中说的S 。
【输出格式】
一行。一个整数，表示最长的 T 的长度。
【样例输入】
ababa
【样例输出】
3
「数据范围」
对于 30% 的数据，S长度 <= 100
对于 60% 的数据，S长度 <= 8000
对于 100% 的数据，S长度 <= 500000
资源约定：
峰值内存消耗 < 256M
CPU消耗 < 1000ms

#include<iostream>
#include<cstring>
#include<algorithm>
#define max 1000000
using namespace std;
char str[max],*p[max];
 int Rank[max], sa[max],wa[max],wb[max],c[max],height[max];
void DA(int n,int m)
{
	int *x = wa, *y = wb;
	int i, j, p;
	for (i = 0;i < n;i++)       c[x[i] = str[i] - 'A'] += 1;
	for (j = 1;j < m;j++)       c[j] += c[j - 1];
	for (i = n - 1;i >= 0;i--)  sa[--c[x[i]]] = i;    //第一次排序
	for (p = 1, j = 1, i = 0;p <n;j = j * 2)
	{
		for (p = 0, i = n - 1;i >=n-j;i--)
	        y[p++] = i;
		for (i = 0;i < n;i++) if (sa[i] >= j) y[p++] = sa[i] - j;
		memset(c, 0, sizeof(int)*m);
		for ( i = 0;i<n;++i) c[x[y[i]]]++;
		for ( i = 1;i<m;++i) c[i] += c[i - 1];
		for ( i = 0;i<n;++i) sa[--c[x[y[i]]]] = y[i];
		swap(x, y);
		x[sa[0]] = 0; p = 1;
		for ( i = 1;i<n;++i)
			x[sa[i]] = (y[sa[i]] == y[sa[i - 1]] && y[sa[i] + j] == y[sa[i - 1] + j]) ? p -1 : p++;
		m = p;
	}
}
void calheight(char *r, int n)
{
	int i, j, k = 0;
	for (i = 1;i <= n;i++) Rank[sa[i]] = i;
	for (i = 0;i < n;height[Rank[i++]] = k)
		for (k?--k:0,j = sa[Rank[i] - 1];i+k<n&&j+k<n&&r[i + k] == r[j + k];k++); //求h[i] = height[rank[i]];
}
	int main()
	{
		freopen("in.txt", "r", stdin);
		int i = 0;
		int ans = 0;
		while (scanf("%c", &str[i++]), str[i - 1] != '\0');
		int z = strlen(str);
		DA(z, 59);
		calheight(str, z);
		for (i = 0;i < z;i++)
		{
			int num = 0;
			ans = ans < height[i] ? height[i]: ans;
			
		}
		cout << z << endl;
		cout << ans;
		return 0;
	}