BZOJ3238：[AHOI 2013]差异

求一个字符串的∑ ∑ len[i] + len[j] - 2 * lcp(i, j)，其中i，j表示从i，j开始的后缀。

方法一：SA+单调栈，自行yy。

方法二：SAM构造出来，然后每个状态对答案的贡献就是：C(|right_s|,2)*(Max_s-Max_parent_s)。前面使用的变量名含义与CLJ论文里的一致。right集合的大小需要构造后缀树然后从那些叶子节点开始算。

由于可以用马拉车在O(N)时间内枚举所有的回文串，所以在parent树中倍增找到每个回文串的出现次数即可。

跑得巨慢。

 #include<stdio.h>
 #include<string.h>
 #include<stdlib.h>
 #include<algorithm>
 //#include<iostream>
 using namespace std;

 int n;
 #define maxn 600011
 char s[maxn];

 struct samnode
 {
     int ch[],pre;
     int pos;
     samnode() {memset(ch,,sizeof(ch)); pre=;}
 };
 int pp[maxn],val[maxn];
 struct SAM
 {
     samnode a[maxn];
     int root,last,size;
     SAM() {root=; a[].pos=; size=; last=root;}
     int idx(char c) {return c-'a';}
     void insert(char c,int p)
     {
         int id=idx(c),x=++size;
         a[x].pos=p; pp[p]=x; val[x]=;
         int y=last;
         for (;y && !a[y].ch[id];y=a[y].pre) a[y].ch[id]=x;
         last=x;
         if (!y) a[x].pre=root;
         else if (a[a[y].ch[id]].pos==a[y].pos+) a[x].pre=a[y].ch[id];
         else
         {
             int z=a[y].ch[id],w=++size;
             a[w]=a[z]; a[w].pos=a[y].pos+;
             a[z].pre=a[x].pre=w;
             for (;y && a[y].ch[id]==z;y=a[y].pre) a[y].ch[id]=w;
         }
     }
 }sam;

 int list[maxn],fa[maxn][];
 bool cmp(const int &i,const int &j) {return sam.a[i].pos>sam.a[j].pos;}
 void build()
 {
     for (int i=;i<=sam.size;i++) list[i]=i,fa[i][]=sam.a[i].pre;
     sort(list+,list++sam.size,cmp);
     for (int j=;j<=;j++)
         for (int i=;i<=sam.size;i++)
             fa[i][j]=fa[fa[i][j-]][j-];
     for (int i=;i<=sam.size;i++)
         val[fa[list[i]][]]+=val[list[i]];
 }

 #define LL long long
 LL getans(int x,int len)
 {
     int p=pp[x];
     for (int j=;j>=;j--) if (fa[p][j] && sam.a[fa[p][j]].pos>=len) p=fa[p][j];
 //    cout<<x<<' '<<len<<' '<<val[p]<<endl;
     return 1ll*len*val[p];
 }

 int mnc[maxn];
 int main()
 {
     scanf("%s",s+);n=strlen(s+);
     for (int i=;i<=n;i++) sam.insert(s[i],i);
     build();

     for (int i=n;i>=;i--)
     {
         s[(i<<)-]=s[i];
         s[i<<]='#';
     }
     n=n+n-;
 //    cout<<s+1<<endl;
     int id=;LL ans=;
     for (int i=;i<=n;i++)
     {
         if (mnc[id]+id>i) mnc[i]=min(mnc[id*-i],mnc[id]+id-i);
         else
         {
             mnc[i]=;
             if (s[i]!='#') ans=max(ans,getans((i+)>>,));
         }
         while (i+mnc[i]<=n && i-mnc[i]> && s[i+mnc[i]]==s[i-mnc[i]])
         {
             if (s[i+mnc[i]]!='#') ans=max(ans,getans((i+mnc[i]+)>>,mnc[i]+));
             mnc[i]++;
         }
         if (mnc[id]+id<mnc[i]+i) id=i;
 //        cout<<i<<' '<<mnc[i]<<endl;
     }
     printf("%lld\n",ans);
     return ;
 }

方法三：PAM。不会。