POJ 2217：Secretary（后缀数组）

题目大意：求两个字符串的公共子串。

分析：

模板题，将两个字符串接起来用不会出现的字符分割，然后求分属两个字符串的相邻后缀lcp的最大值即可。

代码：

program work;
type
  arr=array[..]of longint;
var
  sa,rank,b,tmp,lcp:arr;
  n,i,m,l,u,ans,t:longint;
  s,s1,s2:ansistring;
  ch:char;
function compare(i,j,k:longint):longint;
var ri,rj:longint;
begin
  if rank[i]<>rank[j] then exit(ord(rank[i]<rank[j]))
   else
     begin
       if i+k<=n then ri:=rank[i+k] else ri:=-;
       if j+k<=n then rj:=rank[j+k] else rj:=-;
       exit(ord(ri<rj));
     end;
end;
procedure qsort(l,h,k:longint; var a:arr);
var i,j,t,m:longint;
begin i:=l; j:=h;
  m:=a[(i+j) div ];
 repeat
while compare(a[i],m,k)= do inc(i);
while compare(m,a[j],k)= do dec(j);
if i<=j then
 begin   t:=a[i]; a[i]:=a[j]; a[j]:=t;
inc(i); dec(j);  end;  until i>j;
 if i<h then qsort(i,h,k,a); if j>l then qsort(l,j,k,a);  end;
procedure work_sa(s:ansistring; var sa:arr);
var k,i:longint;
begin
  for i:= to n do begin sa[i]:=i;  rank[i]:=ord(s[i]); end;
  k:=;
  while k<=n do
   begin
     qsort(,n,k,sa);
     tmp[sa[]]:=;
     for i:= to n do
       tmp[sa[i]]:=tmp[sa[i-]]+compare(sa[i-],sa[i],k);
     rank:=tmp;
     k:=k*;
   end;
end;
function max(x,y:longint):longint;
begin
   if x>y then max:=x else max:=y;
end;
procedure work_lcp(s:ansistring;var sa,lcp:arr);
var i,j,h:longint;
begin
  for i:= to n do rank[sa[i]]:=i;
  h:=; lcp[]:=; sa[]:=;rank[]:=;
  for i:= to n do
   begin
     j:=sa[rank[i]-];
     if h> then dec(h);
     while (i+h<=n)and(j+h<=n) do
      begin
        if s[i+h]<>s[j+h] then break; inc(h);
      end;
     lcp[rank[i]-]:=h;
   end;
end;
begin
  readln(u);
  for l:= to u do
   begin
     readln(s1); readln(s2);
     s:=s1+'@'+s2;
     t:=length(s1);
     n:=length(s); work_sa(s,sa);
     work_lcp(s,sa,lcp); ans:=;
    for i:= to n- do begin
      if ((sa[i]<=t)and(sa[i+]>=t+))or((sa[i+]<=t)and(sa[i]>=t+)) then ans:=max(ans,lcp[i]);
    end;
     writeln('Nejdelsi spolecny retezec ma delku ',ans,'.');
  end;
end.