[BZOJ4516][Sdoi2016]生成魔咒(字尾陣列+連結串列)
阿新 • • 發佈:2019-02-09
題目
題解
這道題還是比較好的;
要求出每一個字首本質不同的字尾的個數,那麼我們可以把原序列倒過來,然後實際上就是對於每一個字尾求與其它字尾不重複的字首個數,也即是字尾長度減去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;
}
總結
連結串列的使用比較巧妙;