排名为i和j（i＜j）的最长前缀长度=min(height[i+1..j])

对于排名为i的后缀，想要求得f[i]=∑lcp(i,j)f[i]=∑lcp(i,j)(i＜j)我们只用维护好(i,j]区间的height最小值就好，而且如果遇到某个j，height[j]比height[i]小，那么j及其后的后缀对答案的贡献就是f[j]了，j之前的后缀一共是j-i个，对答案的贡献就是height[i]*(j-i)；反之如果height[j]>=height[i]，那么height[j]以后对答案没有任何贡献（因为有比它小的height[i]存在），直接排除它，也就是说对于height的使用是存在单调性的，使用单调栈就好（一开始我还不怎么会单调栈，蛋疼了好久）
注意：
1.起初对栈底放入len+1，使得栈不为空，从而计算各个值
2.对于原式中lcp以外的东西，我们可以把它化成（n是字符串长度）
(n(n+1)(2n+1)6−n(n+1)2)∗32(\frac{n(n+1)(2n+1)}6-\frac{n(n+1)}2)*\frac32
代码：

#include<bits/stdc++.h>
#define M 500004
#define LL long long
using namespace std;
char s[M];
int w[M],cnt[M],sa[M],rank[M],tmp[M],id[M],height[M];
LL ans,f[M];
stack<int>S;
void SA(int len,int up)
{int *rk=rank,p=0,*t=tmp,d=1;for (int i=0;i<len;i++) cnt[rk[i]=w[i]]++;for (int i=1;i<up;i++) cnt[i]+=cnt[i-1];for (int i=len-1;i>=0;i--) sa[--cnt[rk[i]]]=i;for (;;){for (int i=len-d;i<len;i++) id[p++]=i;for (int i=0;i<len;i++)if (sa[i]>=d) id[p++]=sa[i]-d;for (int i=0;i<up;i++) cnt[i]=0;for (int i=0;i<len;i++) cnt[t[i]=rk[id[i]]]++;for (int i=1;i<up;i++) cnt[i]+=cnt[i-1];for (int i=len-1;i>=0;i--) sa[--cnt[t[i]]]=id[i];swap(t,rk);p=1;rk[sa[0]]=0;for (int i=0;i<len-1;i++)if (sa[i]+d<len&&sa[i+1]+d<len&&t[sa[i]]==t[sa[i+1]]&&t[sa[i]+d]==t[sa[i+1]+d])rk[sa[i+1]]=p-1;elserk[sa[i+1]]=p++;if (p==len) break;d<<=1;up=p;p=0;}
}
void Height(int len)
{for (int i=1;i<=len;i++) rank[sa[i]]=i;int k=0,x;for (int i=0;i<len;i++){k=max(k-1,0);x=sa[rank[i]-1];while (w[i+k]==w[x+k]) k++;height[rank[i]]=k;}
}
main()
{scanf("%s",s);int len=strlen(s);ans=((LL)len*(len+1)*(len*2+1)/6-(LL)len*(len+1)/2)*3/2;for (int i=0;i<len;i++) w[i]=s[i]-'a'+1;SA(len+1,28);Height(len);S.push(len+1);for (int i=len;i>=1;i--){while(height[S.top()]>height[i]) S.pop();f[i]=(LL)height[i]*(S.top()-i)+f[S.top()];ans-=f[i]<<1;S.push(i);}printf("%lld",ans);
}

【BZOJ3238】差异，后缀数组+单调栈维护height相关推荐

[BZOJ3238][AHOI2013]差异 [后缀数组+单调栈]
题目地址 - GO-> 题目大意: 给定一个长度为 nn 的字符串SS,令TiTi表示它从第ii个字符开始的后缀,求以下这个式子的值: ∑1≤i<j≤nlen(Ti)+len(Tj)−2× ...
【BZOJ3879】SvT，后缀数组+单调栈维护sum
Time:2016.08.15 Author:xiaoyimi 转载注明出处谢谢如果有不明白的地方可以在下面评论问我传送门思路: 建立后缀数组求出Height 如果说对每次询问暴力求LCP累加的 ...
[Ahoi2013]差异[后缀数组+单调栈]
链接解题思路:很明显前面∑1<=i<j<=nlen(Ti)+len(Tj)\sum_{1<=i<j<=n}len(T_i)+len(T_j)∑1<=i< ...
BZOJ 3238: [Ahoi2013]差异 [后缀数组单调栈]
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 2326 Solved: 1054 [Submit][Stat ...
[BZOJ 3238] [AHOI 2013] 差异【后缀数组 + 单调栈】
题目链接:BZOJ - 3238 题目分析显然,这道题就是求任意两个后缀之间的LCP的和,这与后缀数组的联系十分明显. 求出后缀数组后,求出字典序相邻两个后缀的LCP,即 Height 数组. 那么 ...
POJ - 3415 Common Substrings(后缀数组+单调栈)
题目链接:点击查看题目大意:给出两个字符串,再给出一个k,问两个字符串中长度大于等于k的公共子串有多少个(种类可重复) 题目分析:因为涉及到了子串问题,先用后缀数组跑出height数组来,接下来如果 ...
[bzoj3238]差异（后缀数组+单调栈）
显然我们可以先把len(Ti)+len(Tj)的值先算出来,再把LCP减去.所有len(Ti)+len(Tj)的值为n*(n-1)*(n+1)/2,这个随便在纸上画一画就可以算出来的. 接下来问题就是 ...
bzoj 3238: [Ahoi2013]差异（后缀数组+单调栈）
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 3443 Solved: 1562 [Submit][Stat ...
BZOJ3879: SvT【后缀数组+单调栈】
Description (我并不想告诉你题目名字是什么鬼) 有一个长度为n的仅包含小写字母的字符串S,下标范围为[1,n]. 现在有若干组询问,对于每一个询问,我们给出若干个后缀(以其在S中出现的起始 ...

【BZOJ3238】差异，后缀数组+单调栈维护height

排名为i和j（i＜j）的最长前缀长度=min(height[i+1..j])

【BZOJ3238】差异，后缀数组+单调栈维护height相关推荐

最新文章

热门文章