BZOJ 3238: [Ahoi2013]差异 [后缀自动机]

3238: [Ahoi2013]差异

Time Limit: 20 Sec  Memory Limit: 512 MB
Submit: 2512  Solved: 1140
[Submit][Status][Discuss]

Description

Input

一行，一个字符串S

Output

一行，一个整数，表示所求值

---

后缀数组看这里 http://www.cnblogs.com/candy99/p/6250732.html

反串建SAM然后Parent Tree就是后缀树了

后缀树上两点的LCP就是LCA啊

然后树形DP，没必要显示建树，直接基数排序倒推行了

状态u作为LCA的时候，就是u的不同子树v的|Right|两两相乘的和

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int N=1e6+;
typedef long long ll;
int n,d[N];
char s[N];
struct State{
    int ch[],par,val;
}t[N];
int sz,root,last;
inline int nw(int _){t[++sz].val=_;return sz;}
inline void iniSAM(){sz=;root=last=nw();}
void extend(int c){
    int p=last,np=nw(t[p].val+); d[np]=;
    for(;p&&!t[p].ch[c];p=t[p].par) t[p].ch[c]=np;
    if(!p) t[np].par=root;
    else{
        int q=t[p].ch[c];
        if(t[q].val==t[p].val+) t[np].par=q;
        else{
            int nq=nw(t[p].val+);
            memcpy(t[nq].ch,t[q].ch,sizeof(t[q].ch));
            t[nq].par=t[q].par;
            t[q].par=t[np].par=nq;
            for(;p&&t[p].ch[c]==q;p=t[p].par) t[p].ch[c]=nq;
        }
    }
    last=np;
}
 
int c[N],a[N];
void RadixSort(){
    for(int i=;i<=sz;i++) c[t[i].val]++;
    for(int i=;i<=n;i++) c[i]+=c[i-];
    for(int i=sz;i>=;i--) a[c[t[i].val]--]=i;
}
ll dp(){
    ll _=;
    for(int i=sz;i>=;i--){
        int v=a[i],u=t[v].par;
        _+=(ll)t[u].val*d[u]*d[v];
        d[u]+=d[v];
    }
    return _;
}
void solve(){
    ll ans=(ll)n*(n+)*(n-)/;
    iniSAM();
    for(int i=;i<=n;i++) extend(s[i]-'a');
    RadixSort();
    printf("%lld",ans-*dp());
}
int main(int argc, const char * argv[]) {
    freopen("in","r",stdin);
    scanf("%s",s+);
    n=strlen(s+);
    reverse(s+,s++n);
    solve();
    return ;
}