[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;
}

总结

链表的使用比较巧妙;

猜你喜欢

转载自blog.csdn.net/A_Comme_Amour/article/details/79987498