BZOJ_4516_[Sdoi2016]生成魔咒_后缀数组+ST表+splay
BZOJ_4516_[Sdoi2016]生成魔咒_后缀数组+ST表+splay
Description
魔咒串由许多魔咒字符组成，魔咒字符可以用数字表示。例如可以将魔咒字符 1、2 拼凑起来形成一个魔咒串 [1,2]。
一个魔咒串 S 的非空字串被称为魔咒串 S 的生成魔咒。
例如 S=[1,2,1] 时，它的生成魔咒有 [1]、[2]、[1,2]、[2,1]、[1,2,1] 五种。S=[1,1,1] 时，它的生成魔咒有 [1]、
[1,1]、[1,1,1] 三种。最初 S 为空串。共进行 n 次操作，每次操作是在 S 的结尾加入一个魔咒字符。每次操作后都
需要求出，当前的魔咒串 S 共有多少种生成魔咒。
Input
第一行一个整数 n。
第二行 n 个数，第 i 个数表示第 i 次操作加入的魔咒字符。
1≤n≤100000。,用来表示魔咒字符的数字 x 满足 1≤x≤10^9
Output
输出 n 行，每行一个数。第 i 行的数表示第 i 次操作后 S 的生成魔咒数量
Sample Input
7
1 2 3 3 3 1 2
1 2 3 3 3 1 2
Sample Output
1
3
6
9
12
17
22
3
6
9
12
17
22
在后面添加字符每次会产生若干个后缀。于是翻转一下相当于每次只添加一个后缀。
然后如何求不同的子串个数呢？
首先对于一个后缀$Suffix(i)$,会贡献$n-i+1$个子串，但这些子串会有一些重复的，于是找$rank[i]$的前驱和后继所在的后缀设为$j$。
那么会有$Lcp(suffix(i),suffix(j))$这么多是重复的。
于是用splay维护前驱后继，然后ST表求两个后缀的LCP。
代码：
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 100050
#define ls ch[p][0]
#define rs ch[p][1]
#define get(x) (ch[f[x]][1]==x)
typedef long long ll;
int n,ws[N],wv[N],wa[N],wb[N],rank[N],sa[N],height[N],r[N];
int f[N],ch[N][2],siz[N],rt,val[N],reimu;
//////////////////////////////////////////////
struct A {
int num,id,v;
}a[N];
bool cmp1(const A &x,const A &y){return x.num<y.num;}
bool cmp2(const A &x,const A &y){return x.id<y.id;}
///////////////////////////////////////////////
void build_suffix_array() {
int i,j,p,*x=wa,*y=wb,*t,m=n;
for(i=0;i<m;i++) ws[i]=0;
for(i=0;i<n;i++) ws[x[i]=r[i]]++;
for(i=1;i<m;i++) ws[i]+=ws[i-1];
for(i=n-1;i>=0;i--) sa[--ws[x[i]]]=i;
for(j=p=1;p<n;j<<=1,m=p) {
for(p=0,i=n-j;i<n;i++) y[p++]=i;
for(i=0;i<n;i++) if(sa[i]-j>=0) y[p++]=sa[i]-j;
for(i=0;i<n;i++) wv[i]=x[y[i]];
for(i=0;i<m;i++) ws[i]=0;
for(i=0;i<n;i++) ws[wv[i]]++;
for(i=1;i<m;i++) ws[i]+=ws[i-1];
for(i=n-1;i>=0;i--) sa[--ws[wv[i]]]=y[i];
for(t=x,x=y,y=t,i=p=1,x[sa[0]]=0;i<n;i++) {
if(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+j]==y[sa[i-1]+j]) x[sa[i]]=p-1;
else x[sa[i]]=p++;
}
}
for(i=1;i<n;i++) rank[sa[i]]=i;
for(i=p=0;i<n-1;height[rank[i++]]=p)
for(p?p--:0,j=sa[rank[i]-1];r[i+p]==r[j+p];p++);
}
/////////////////////////////////////////////
int newnode(int x) {
siz[++reimu]=1; val[reimu]=x; return reimu;
}
void pushup(int p) {
siz[p]=siz[ls]+siz[rs]+1;
}
void rotate(int x) {
int y=f[x],z=f[y],k=get(x);
ch[y][k]=ch[x][!k]; f[ch[y][k]]=y;
ch[x][!k]=y; f[y]=x; f[x]=z;
if(z) ch[z][ch[z][1]==y]=x;
pushup(y); pushup(x);
if(rt==y) rt=x;
}
void splay(int x,int y) {
for(int fa;(fa=f[x])!=y;rotate(x))
if(f[fa]!=y)
rotate(get(x)==get(fa)?fa:x);
}
void insert(int x) {
int p=rt,l,r;
while(p) {
if(val[p]>=x) r=p,p=ls;
else l=p,p=rs;
}
splay(l,0); splay(r,rt);
ch[r][0]=newnode(x);
f[reimu]=r; pushup(r); pushup(l);
}
int pre(int x) {
int p=rt,ans;
while(p) {
if(val[p]>=x) p=ls;
else ans=p,p=rs;
}
return val[ans];
}
int nxt(int x) {
int p=rt,ans;
while(p) {
if(val[p]<=x) p=rs;
else ans=p,p=ls;
}
return val[ans];
}
/////////////////////////////////////////
struct ST {
int f[N][20],L[N];
void init() {
int i,j;
memset(f,0x3f,sizeof(f));
L[1]=0;
for(i=2;i<=n;i++) L[i]=L[i>>1]+1;
for(i=0;i<=n;i++) {
f[i][0]=height[i];
}
for(j=1;(1<<j)<=n;j++) {
for(i=0;i+(1<<j)-1<=n;i++) {
f[i][j]=min(f[i][j-1],f[i+(1<<j-1)][j-1]);
}
}
}
int get_min(int l,int r) {
int len=L[r-l+1];
return min(f[l][len],f[r-(1<<len)+1][len]);
}
}S;
int main() {
scanf("%d",&n);
int i;
for(i=1;i<=n;i++) scanf("%d",&a[n-i+1].num),a[i].id=i;
sort(a+1,a+n+1,cmp1);
int j=0;a[0].num=23333443;
for(i=1;i<=n;i++) {
if(a[i].num!=a[i-1].num) j++;
a[i].v=j;
}
sort(a+1,a+n+1,cmp2);
for(i=0;i<n;i++) {
r[i]=a[i+1].v;
}
r[n++]=0;
build_suffix_array();
/*for(i=0;i<n;i++) printf("%d ",r[i]); puts("");
for(i=0;i<n;i++) printf("%d ",sa[i]); puts("");
for(i=0;i<n;i++) printf("%d ",height[i]); puts("");
for(i=0;i<n;i++) printf("%d ",rank[i]);puts("");*/
ll ans=0;
rt=newnode(-100000000);
ch[rt][1]=newnode(100000000);
f[ch[rt][1]]=rt;
pushup(rt);
S.init();
for(i=n-2;i>=0;i--) {
ans+=n-i-1;
int pr=pre(rank[i]);
int tmp=0;
if(pr>=0) {
tmp=S.get_min(pr+1,rank[i]);
}
int nx=nxt(rank[i]);
if(nx<=n) {
tmp=max(tmp,S.get_min(rank[i]+1,nx));
}
ans-=tmp;
insert(rank[i]);
printf("%lld\n",ans);
}
}
/*
7
1 2 3 3 3 1 2
*/
