[bzoj3238]差异（后缀数组+单调栈）

显然我们可以先把len(Ti)+len(Tj)的值先算出来，再把LCP减去。所有len(Ti)+len(Tj)的值为n*(n-1)*(n+1)/2，这个随便在纸上画一画就可以算出来的。

接下来问题就是如何把LCP减去。我们先用后缀数组把height求出来，当有一段区间l~r，height[i]为height[l]~height[r]中的最小值，那么随便取rk[l]~rk[r]中的两个后缀，他们的LCP则都是height[i]，这个很好理解吧。那么l~r这个区间里有(l-i+1)*(r-i+1)对后缀，所以我们最后的答案就要减去2*height[i]*(l-i+1)*(r-i+1)【1≤i≤n】。

然后就是如何求出每一个i的l~r了，暴力枚举+RMQ显然不行，那我们就用一个单调栈，栈里存着i前面height值比height[i]小的height值的编号，记为j，如果height[j]比height[i]大那么就弹出，那么这段区间的左端点则为栈顶的j+1，右端点同理。这样我们就可以求出每个height的l和r了。

奇丑无比的代码如下：

var
  s:ansistring;
  i:longint;
  n,m,l,r,ans,top:int64;
  rk,trk,sa,tsa,sum,h,ll,rr,st:array[..]of int64;

procedure suffix;
var
  i,j,p:longint;
begin
  for i:= to n do begin trk[i]:=ord(s[i]);inc(sum[trk[i]]);end;
  for i:= to  do inc(sum[i],sum[i-]);
  for i:=n downto  do begin sa[sum[trk[i]]]:=i;dec(sum[trk[i]]);end;
  rk[sa[]]:=;p:=;
  for i:= to n do begin if trk[sa[i]]<>trk[sa[i-]] then inc(p);rk[sa[i]]:=p;end;
  m:=p;j:=;
  while m<n do
  begin
    move(rk,trk,sizeof(rk));fillchar(sum,sizeof(sum),);p:=;
    for i:=n-j+ to n do begin inc(p);tsa[p]:=i;end;
    for i:= to n do if sa[i]>j then begin inc(p);tsa[p]:=sa[i]-j;end;
    for i:= to n do begin rk[i]:=trk[tsa[i]];inc(sum[rk[i]]);end;
    for i:= to n do inc(sum[i],sum[i-]);
    for i:=n downto  do begin sa[sum[rk[i]]]:=tsa[i];dec(sum[rk[i]]);end;
    rk[sa[]]:=;p:=;
    for i:= to n do
    begin
      if (trk[sa[i]]<>trk[sa[i-]])or(trk[sa[i]+j]<>trk[sa[i-]+j])then inc(p);
      rk[sa[i]]:=p;
    end;
    m:=p;j:=j*;
  end;
  h[]:=;p:=;
  for i:= to n do
  begin
    if rk[i]= then continue;
    j:=sa[rk[i]-];
    while s[i+p]=s[j+p] do inc(p);
    h[rk[i]]:=p;
    if p> then dec(p);
  end;
end;

begin
  readln(s);
  n:=length(s);
  s:=s+' ';
  suffix;
  ans:=n*(n-)*(n+)div ;
  h[]:=-maxlongint;
  for i:= to n do
  begin
    while h[i]<=h[st[top]] do dec(top);
    if st[top]= then ll[i]:=
    else ll[i]:=st[top]+;
    inc(top);
    st[top]:=i;
  end;
  h[n+]:=-maxlongint;top:=;st[]:=n+;
  for i:=n downto  do
  begin
    while h[i]<h[st[top]] do dec(top);
    if st[top]=n+ then rr[i]:=n
    else rr[i]:=st[top]-;
    inc(top);
    st[top]:=i;
  end;
  for i:= to n do
  ans:=ans-*(i-ll[i]+)*(rr[i]-i+)*h[i];
  writeln(ans);
end.