后缀数组Pascal代码实现

以前不会想，也不会写。现在会写了，也大概想明白了。后缀数组，真是个麻烦的东西。

题目大意如下：给定一个长度在5000以内的字符串，求出出现过两次或两次以上的最长字串（保证唯一解）。

不多说了，裸的后缀数组：

Program Suffix;//By_Thispoet
Const 
	maxn=10000;
Var
	i,j,k,m,n,p,q,ans,pos,sum			:Longint;
	rank,sa,a,x,y						:Array[1..maxn]of Longint;
	st									:Ansistring;

Procedure Qsort(l,r:Longint);
var i,j,k,p,temp:Longint;
begin
	i:=l;j:=r;k:=x[(i+j)>>1];p:=y[(i+j)>>1];
	repeat
		while (x[i]<k)or((x[i]=k)and(y[i]<p)) do inc(i);
		while (x[j]>k)or((x[j]=k)and(y[j]>p)) do dec(j);
		if i<=j then
			begin
				temp:=x[i];x[i]:=x[j];x[j]:=temp;
				temp:=y[i];y[i]:=y[j];y[j]:=temp;
				temp:=sa[i];sa[i]:=sa[j];sa[j]:=temp;
				inc(i);dec(j);
			end;
	until i>j;
	if i<r then Qsort(i,r);
	if l<j then Qsort(l,j);
end;
	
	
BEGIN
	
	readln(st);
	n:=length(st);
	p:=1;
	for i:=1 to n do rank[i]:=ord(st[i])-64;
	while p<=n-1 do 
		begin
			
			for i:=1 to n do x[i]:=rank[i];
			for i:=1 to n do y[i]:=rank[i+p];
			for i:=1 to n do sa[i]:=i;
			
			Qsort(1,n);
			
			sum:=1;
			rank[sa[1]]:=1;
			for i:=2 to n do
				begin
					if (x[i]<>x[i-1])or(y[i]<>y[i-1]) then inc(sum);
					rank[sa[i]]:=sum;
				end;
			if sum=n then break;
			
			p:=p << 1;
		end;
	
	for i:=2 to n do 
		begin
			
			p:=sa[i];q:=sa[i-1];
			j:=0;
			
			while (p<=n)and(q<=n)and(st[p]=st[q]) do 
				begin
					inc(j);
					inc(p);inc(q);
				end;
			if j>ans then 
				begin
					ans:=j;
					pos:=p-1;
				end;
			
		end;
	
	writeln(ans);
	for i:=pos-ans+1 to pos do write(st[i]);
	writeln;
	
END.