[BZOJ3238] [AHOI2013] 差异

3238: [Ahoi2013]差异

Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 1968 Solved: 896
[Submit][Status][Discuss]

Description

Input

一行，一个字符串S

Output

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

Sample Input

cacao

Sample Output

HINT

2<=N<=500000,S由小写英文字母组成

Source

[Submit][Status][Discuss]

感觉窝会求后缀数组然后完全不会用……

好伐我还是写了SAM2333333。

显然题意只要我们求出∑LCP(i,j)(i!=j)就好

首先我们跑一遍SAM求出后缀树，然后拓扑一下就能求出所有right集合的大小。

显然两个后缀的最长公共前缀长就是他们在后缀树上的lca的深度。

我们可以考虑对每个深度进行分段，即点x对答案的贡献只包含其在后缀树上的父节点的深度以上且在自身深度以下的答案。

显然这样做的答案就是∑deep[i]*∑right[i]-right[son[i]]

又因为right[i]=∑right[son[i]]

那么每个点的贡献就是deep[i]*right[i]*right[i-1]了

#include"bits/stdc++.h"
using namespace std;
typedef long long ll;const int N=1000005,A=26;
struct SAM{char s[N];int t,l,q[N],last,a[N];int fail[N],e[N][A],v[N],d[N];SAM(){t=last=1;}void extend(int c){int np=++t,p=last;d[np]=d[p]+1,v[last=np]=1;for(;!e[p][c]&&p;p=fail[p])e[p][c]=np;if(!p)fail[np]=1;else {int q=e[p][c],nq;if(d[q]!=d[p]+1){nq=++t;d[nq]=d[p]+1;memcpy(e[nq],e[q],sizeof(int[A]));fail[nq]=fail[q];fail[np]=fail[q]=nq;for(;e[p][c]==q;p=fail[p])e[p][c]=nq;} else fail[np]=q;}}void read(){scanf("%s",s);l=strlen(s);}void build(){for(int i=0;i<l;i++)extend(s[i]-'a');}ll getans(){for(int i=1;i<=t;i++)a[d[i]]++;for(int i=1;i<=l;i++)a[i]+=a[i-1];for(int i=t;i;i--)q[a[d[i]]--]=i;for(int i=t;i;i--)v[fail[q[i]]]+=v[q[i]];v[1]=v[0]=0;ll ret=(ll)(l-1)*l*(l+1)/2;for(int i=2;i<=t;i++)ret-=(ll)v[i]*(v[i]-1)*(d[i]-d[fail[i]]);return ret;}
} sam;int main(){sam.read();sam.build();printf("%lld\n",sam.getans());return 0;
}

[BZOJ3238] [AHOI2013] 差异 - 后缀自动机相关推荐

bzoj3238 [Ahoi2013]差异后缀自动机
这个题虽然看起来比较麻烦,但其实有一点搞懂了就好做了每一个后缀的前缀的lcp一定是这个串里的子串, 对于全部子串,如果它不是后缀,那一定可以找到这个串对应的后缀, 而这样的子串和后缀子串是没有区别的 ...
bzoj3238 [Ahoi2013]差异后缀自动机
题意:给出一个串,求其中任意两个字串的lcp的总和. 我们可以对于这个串建一颗后缀自动机,实际上,他的parent边树就是一颗后缀树,我们在后缀树上统计答案,设f表示right集合的大小,可以理解为后 ...
[BZOJ3238][AHOI2013]差异 [后缀数组+单调栈]
题目地址 - GO-> 题目大意: 给定一个长度为 nn 的字符串SS,令TiTi表示它从第ii个字符开始的后缀,求以下这个式子的值: ∑1≤i<j≤nlen(Ti)+len(Tj)−2× ...
【bzoj3238】差异后缀自动机
AC通道:http://www.lydsy.com/JudgeOnline/problem.php?id=3238 [题解] 题上所给公式可以化为其中: =1×(n-1)+2×(n-2)+3×(n- ...
BZOJ3238[Ahoi2013]差异
BZOJ3238[Ahoi2013]差异题目描述 n<=500000n<=500000n,都是小写字母输入一行,一个字符串S 输出一行,一个整数,表示所求值 Solution 公式 ...
bzoj3238: [Ahoi2013]差异（后缀自动机）
传送门后缀自动机好题. 题意: 做法:samsamsam 废话考虑翻转字串,这样后缀的最长公共前缀等于前缀的最长公共后缀. 然后想到parentparentparent树上面两个串的最长公共后缀跟 ...
BZOJ3238: [Ahoi2013]差异
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 4840 Solved: 2298 [Submit][Stat ...
BZOJ 3238: [Ahoi2013]差异 [后缀数组单调栈]
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 2326 Solved: 1054 [Submit][Stat ...
BZOJ 3238 差异 [后缀自动机]
题目一个长度为n的字符串S,令TiT_i表示它从第i个字符开始的后缀.求 ∑1≤i<j≤nlen(Ti)+len(Tj)−2∗lcp(Ti,Tj) \sum_{1\le i 其中len(a)表 ...

[BZOJ3238] [AHOI2013] 差异 - 后缀自动机