bzoj 3238: [Ahoi2013]差异（后缀数组+单调栈）

3238: [Ahoi2013]差异

Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 3443 Solved: 1562
[Submit][Status][Discuss]

Description

Input

一行，一个字符串S

Output

一行，一个整数，表示所求值

Sample Input

cacao

Sample Output

后缀数组height数组有个很好的性质：

Lcp(i, j) = min(height[k], k∈[i+1, j] )

这样在算∑lcp(i, j)时，就可以用单调栈计算每个最小值的贡献了

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define LL long long
typedef struct
{int l;int r;
}Point;
Point s[500010];
char str[500010];
int q, n, st[500010], a[500010], Rank[500010], Temp[500010], cot[500050], Sa[500010], height[500010];
int main(void)
{LL ans;int i, m, p, k, top;scanf("%s", str+1);n = strlen(str+1);for(i=1;i<=n;i++)a[i] = str[i];m = 128;for(i=1;i<=n;i++){Rank[i] = a[i];cot[a[i]]++;}for(i=1;i<=m;i++)cot[i] += cot[i-1];for(i=n;i>=1;i--)Sa[cot[a[i]]--] = i;for(k=1;k<=n;k*=2){p = 0;for(i=n-k+1;i<=n;i++)Temp[++p] = i;for(i=1;i<=n;i++){if(Sa[i]>k)Temp[++p] = Sa[i]-k;}memset(cot, 0, sizeof(cot));for(i=1;i<=n;i++)cot[Rank[i]]++;for(i=1;i<=m;i++)cot[i] += cot[i-1];for(i=n;i>=1;i--)Sa[cot[Rank[Temp[i]]]--] = Temp[i];swap(Rank, Temp);p = 1;Rank[Sa[1]] = 1;for(i=2;i<=n;i++){if(Temp[Sa[i-1]]==Temp[Sa[i]] && Temp[Sa[i-1]+k]==Temp[Sa[i]+k])Rank[Sa[i]] = p;elseRank[Sa[i]] = ++p;}m = p;}k = 0;for(i=1;i<=n;i++){k = max(k-1, 0);m = Sa[Rank[i]-1];while(a[i+k]==a[m+k])k++;height[Rank[i]] = k;}ans = 0;for(i=n;i>=1;i--)ans += (LL)i*(n-1);top = 0;for(i=1;i<=n;i++){while(top>=1 && height[st[top]]>=height[i]){s[st[top]].r = i-1;top--;}s[i].l = st[top];st[++top] = i;}while(top){s[st[top]].r = n;top--;}for(i=1;i<=n;i++)ans -= (LL)2*(s[i].r-i+1)*(i-s[i].l)*height[i];printf("%lld\n", ans);return 0;
}
/*
aaaaaa
*/

bzoj 3238: [Ahoi2013]差异（后缀数组+单调栈）相关推荐

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