1. 程式人生 > >[BZOJ4516][Sdoi2016]生成魔咒(字尾陣列+連結串列)

[BZOJ4516][Sdoi2016]生成魔咒(字尾陣列+連結串列)

題目

傳送門

題解

這道題還是比較好的;
要求出每一個字首本質不同的字尾的個數,那麼我們可以把原序列倒過來,然後實際上就是對於每一個字尾求與其它字尾不重複的字首個數,也即是字尾長度減去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; }

總結

連結串列的使用比較巧妙;