【BZOJ 4199】[Noi2015]品酒大会 后缀自动机构造后缀树+dp

题面：http://uoj.ac/problem/131（BZOJ肿么了qaq）

如果连后缀自动机为什么可以转成后缀树都不知道的，先去切3238差异：http://blog.csdn.net/pbihao/article/details/53905995

我们先构造出后缀自动机，然后转成后缀树，在后缀树上跑dp，只用维护一个最大最小，和size就好了#include<cstdio>
#include<cstring>
#include<iostream>
#define maxn 600020
#define LL long long
using namespace std;
int n,son[maxn][26],fail[maxn],size[maxn],len[maxn],tot=1,rt=1,last=1,cnt=1,head[maxn],val[maxn];
LL ans1[maxn],ans2[maxn],max1[maxn],min1[maxn];
bool is[maxn];
int dep[maxn],w[maxn];
char s[maxn];
struct edge{int v,next;}e[maxn*2];
void adde(int a,int b){e[cnt].v=b,e[cnt].next=head[a];head[a]=cnt++;}

void insert(int x,int v){
int p=last,np=++tot,q,nq;
size[np]=1,len[np]=len[p]+1,is[np]=true;
while(p&&!son[p][x])son[p][x]=np,p=fail[p];
if(p==0)fail[np]=rt;
else{
int q=son[p][x];
if(len[q]==len[p]+1)fail[np]=q;
else{
nq=++tot;len[nq]=len[p]+1;
memcpy(son[nq],son[q],sizeof(son[q]));
fail[nq]=fail[q];
fail[q]=fail[np]=nq;
while(p&&son[p][x]==q)son[p][x]=nq,p=fail[p];
}
}
last=np;w[np]=v;
}

void dfs(int u,int fa){
dep[u]=len[u];
if(is[u])max1[u]=min1[u]=w[u];
else max1[u]=-1e15,min1[u]=1e15;
for(int v,i=head[u];i;i=e[i].next){
dfs(v=e[i].v,u);

if(max1[u]!=-1e15&&max1[v]!=-1e15&&min1[v]!=1e15&&min1[u]!=1e15)
ans1[dep[u]]=max(max(ans1[dep[u]],(LL)max1[u]*max1[v]),(LL)min1[u]*min1[v]);
ans2[dep[u]]+=(LL)size[u]*size[v];

max1[u]=max(max1[u],max1[v]);
min1[u]=min(min1[u],min1[v]);

size[u]+=size[v];
}
}

int main(){
scanf("%d",&n);
scanf("%s",s+1);
for(int i=1;i<=n;i++)scanf("%d",val+i);
for(int i=n;i>=1;i--)insert(s[i]-'a',val[i]);
for(int i=2;i<=tot;i++)
adde(fail[i],i);
for(int i=0;i<=n;i++)ans1[i]=-1e18;
dep[0]=-1;
dfs(1,0);
for(int i=n-1;i>=0;i--)
ans1[i]=max(ans1[i],ans1[i+1]),ans2[i]+=ans2[i+1];
for(int i=0;i<n;i++)
if(ans2[i])printf("%lld %lld\n",ans2[i],ans1[i]);
else puts("0 0");
return 0;
}