BZOJ 3238 [Ahoi2013]差异

后缀自动机+后缀树

有一种不难想到的后缀数组的做法，详见黄学长题解

然而我做这题的目的主要是熟悉一下后缀自动机和后缀树- -

一个有趣的性质：反串的后缀自动机的parent树就是正串的后缀树，其树上的边权长度=parent树上儿子的len-父亲的len

构出后缀树之后就是一个小DP啦。

#include<cstdio>
#include<cstring>
#define N 500010
#define S 28
using namespace std;
namespace ziqian
{typedef long long ll;char s[N];ll ans = 0;int last[N*2], ecnt, siz[N*2]; struct edge{int next, to;}e[N<<1];void addedge(int a, int b){e[++ecnt] = (edge){last[a], b}; last[a] = ecnt;}struct SAM{SAM *par, *ch[S];int len;}mem[N*2], *tot, *null, *tai, *root;SAM *newSAM(int l){SAM *r = ++tot;*r = *null; r->len = l;return r;}void init_SAM(){tot = mem;null = ++tot;null -> par = null;for(int i = 0; i < S; i++) null->ch[i] = null;null->len = 0;root = newSAM(0); tai = root;}void extend(int v){SAM *p = tai, *np = newSAM(p->len + 1); tai = np; siz[np - mem] = 1; for(; p != null && p->ch[v] == null; p = p->par) p->ch[v] = np;if(p == null) np -> par = root;else{SAM *q = p->ch[v];if(p->len + 1 == q->len) np->par = q;else{SAM *nq = newSAM(p->len + 1);*nq = *q;nq -> len = p->len + 1;q->par = np->par = nq;for(; p != null && p->ch[v] == q; p = p->par) p->ch[v] = nq; }}}void dp(int x){for(int i = last[x]; i; i = e[i].next){int y = e[i].to;dp(y);ans -= 2ll * siz[x] * siz[y] * (mem+x)->len;siz[x] += siz[y];}}void main(){scanf("%s",s+1);int n = strlen(s+1);init_SAM();for(int i = n; i; i--) extend(s[i] - 'a');for(int i = 1; i <= n; i++) ans += (ll) i * (n-1);for(SAM *cur = tot; cur != mem; cur--)addedge(cur->par - mem, cur - mem);dp(root - mem);printf("%lld\n",ans);}
}
int main()
{ziqian::main();
}

BZOJ 3238 [Ahoi2013]差异相关推荐

BZOJ 3238: [Ahoi2013]差异 [后缀数组单调栈]
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 2326 Solved: 1054 [Submit][Stat ...
bzoj 3238: [Ahoi2013]差异（后缀数组+单调栈）
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 3443 Solved: 1562 [Submit][Stat ...
@bzoj - 3238@ [Ahoi2013]差异
目录 @description@ @solution@ @accepted code@ @details@ @description@ 给定一个长度为 n 的字符串 S,令 Ti 表示它从第 i 个字 ...
bzoj 3238: [Ahoi2013]差异
Description Input 一行,一个字符串S Output 一行,一个整数,表示所求值 Sample Input cacao Sample Output 54 HINT 2<=N< ...
3238: [Ahoi2013]差异
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 2106 Solved: 953 [Submit][Statu ...
【BZOJ】 3238: [Ahoi2013]差异
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3238 求:$${\sum _{i=1}^{n-1}\sum _{j=i+1}^{n}len ...
【BZOJ】3238: [Ahoi2013]差异
[题意]给定长度为n的小写字母字符串,令Ti表示以i开头的后缀,求Σ[Ti+Tj-2*lcp(Ti,Tj)],1<=i<j<=n. [算法]后缀自动机 [题解]Σ(Ti+Tj)只与n ...
【bzoj 3238】差异
传送门~ 解题思路求后缀数组以后搞个单调栈. 因为height数组的意义让这个单调栈很烦,注意要理清思路再写. 代码: #include<algorithm> #include<i ...
BZOJ3238: [Ahoi2013]差异
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 4840 Solved: 2298 [Submit][Stat ...

BZOJ 3238 [Ahoi2013]差异

BZOJ 3238 [Ahoi2013]差异相关推荐

最新文章

热门文章