题目
题解
这道题还是比较好的;
要求出每一个前缀本质不同的后缀的个数,那么我们可以把原序列倒过来,然后实际上就是对于每一个后缀求与其它后缀不重复的前缀个数,也即是后缀长度减去height值;
求出某一个后缀对答案的贡献之后,他不应该停留在元序列中对后续答案的求解产生影响,所以应该把它删除;
在实现方式上,可以使用链表,与平衡树的操作有些类似
代码
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define ll long long
const int maxn=1e6+10;
const int inf=1e9;
int n,m=200,x[maxn],y[maxn],c[maxn],sa[maxn],rnk[maxn],height[maxn];
int s[maxn];
void build_sa()
{
for (int i=0; i<m; i++) c[i]=0;
for (int i=0; i<n; i++) c[x[i]=s[i]]++;
for (int i=1; i<m; i++) c[i]+=c[i-1];
for (int i=n-1; i>=0; i--) sa[--c[x[i]]]=i;
for (int k=1; k<=n; k<<=1)
{
int p=0;
for (int i=n-k; i<n; i++) y[p++]=i;
for (int i=0; i<n; i++) if (sa[i]>=k) y[p++]=sa[i]-k;
for (int i=0; i<m; i++) c[i]=0;
for (int i=0; i<n; i++) c[x[i]]++;
for (int i=1; i<m; i++) c[i]+=c[i-1];
for (int i=n-1; i>=0; i--) sa[--c[x[y[i]]]]=y[i];
swap(x,y);//又忘记了swap
p=1; x[sa[0]]=0;
for (int i=0; i<n; i++)
x[sa[i]] = y[sa[i-1]]==y[sa[i]] && ((sa[i-1]+k>=n?-1:y[sa[i-1]+k])==(sa[i]+k>=n?-1:y[sa[i]+k])) ?p-1:p++;
if (p>n) break;
m=p;
}
}
void build_height()
{
for (int i=0; i<n; i++) rnk[sa[i]]=i;
int k=0; height[0]=0;
for (int i=0; i<n; i++)
{
if (!rnk[i]) continue;
if (k) k--;
int j=sa[rnk[i]-1];
while (i+k<n && j+k<n && s[i+k]==s[j+k]) k++;
height[rnk[i]]=k;
}
}
void debug(int *A,int len) {for (int i=0; i<len; i++) printf("%d ",A[i]);}
int tmp[maxn],a[maxn],pre[maxn],next[maxn];
ll ans[maxn];
int main()
{
scanf("%d",&n);
for (int i=0; i<n; i++) scanf("%d",&a[i]);
for (int i=0; i<n; i++) tmp[i]=a[i];
sort(tmp,tmp+n);
int nn=unique(tmp,tmp+n)-tmp;
for (int i=0; i<n; i++)
a[i]=lower_bound(tmp,tmp+nn,a[i])-tmp+1;
for (int i=0; i<n; i++) s[i]=a[n-i-1];
build_sa();
build_height();
for (int i=0; i<n-1; i++) next[i]=i+1;
for (int i=1; i<n; i++) pre[i]=i-1;
for (int i=0; i<n; i++)
{
int rk=rnk[i];
int now=n-i-max(height[rk],height[next[rk]]);
ans[i]=(ll)now;
height[next[rk]]=min(height[next[rk]],height[rk]);
height[rk]=0;
if (rk) next[pre[rk]]=next[rk];
pre[next[rk]]=pre[rk];//相当于在平衡树中删除了一个元素
}
for (int i=n-1; i>=0; i--) ans[i]+=ans[i+1];
for (int i=n-1; i>=0; i--) printf("%lld\n",ans[i]);
return 0;
}
总结
链表的使用比较巧妙;