BZOJ4516SDOI2016生成魔咒 [SAM]
Posted BearChild
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ4516SDOI2016生成魔咒 [SAM]相关的知识,希望对你有一定的参考价值。
生成魔咒
Time Limit: 10 Sec Memory Limit: 128 MB[Submit][Status][Discuss]
Description
魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示。例如可以将魔咒字符 1、2 拼凑起来形成一个魔咒串 [1,2]。
一个魔咒串 S 的非空字串被称为魔咒串 S 的生成魔咒。
例如 S=[1,2,1] 时,它的生成魔咒有 [1]、[2]、[1,2]、[2,1]、[1,2,1] 五种。
S=[1,1,1] 时,它的生成魔咒有 [1]、[1,1]、[1,1,1] 三种。
最初 S 为空串。共进行 n 次操作,每次操作是在 S 的结尾加入一个魔咒字符。
每次操作后都需要求出,当前的魔咒串 S 共有多少种生成魔咒。
Input
第一行一个整数 n。
第二行 n 个数,第 i 个数表示第 i 次操作加入的魔咒字符。
Output
输出 n 行,每行一个数。第 i 行的数表示第 i 次操作后 S 的生成魔咒数量
Sample Input
7
1 2 3 3 3 1 2
1 2 3 3 3 1 2
Sample Output
1
3
6
9
12
17
22
3
6
9
12
17
22
HINT
1≤n≤100000
Main idea
询问在加入每一个字符后当前有多少个本质不同的子串。
Solution
直接用SAM,根据SAM的性质,每次增多的子串个数就是len[New] - len[fa[New]]。
Code
1 #include<iostream> 2 #include<string> 3 #include<algorithm> 4 #include<cstdio> 5 #include<cstring> 6 #include<cstdlib> 7 #include<cmath> 8 #include<map> 9 using namespace std; 10 typedef long long s64; 11 12 const int ONE = 400005; 13 const int INF = 2147483640; 14 15 int n,x; 16 s64 Ans; 17 18 int get() 19 { 20 int res=1,Q=1;char c; 21 while( (c=getchar())<48 || c>57 ) 22 if(c==‘-‘)Q=-1; 23 res=c-48; 24 while( (c=getchar())>=48 && c<=57 ) 25 res=res*10+c-48; 26 return res*Q; 27 } 28 29 struct SAM 30 { 31 map <int, int> a[ONE]; 32 int len[ONE], fa[ONE]; 33 int last, cnt; 34 SAM() {last = cnt = 1;} 35 void Add(int c) 36 { 37 int x = last, New = last = ++cnt; 38 len[New] = len[x] + 1; 39 while(x && !a[x][c]) a[x][c] = New, x = fa[x]; 40 if(!x) {fa[New] = 1; Ans += len[New] - len[fa[New]]; return;} 41 42 int q = a[x][c]; 43 if(len[x] + 1 == len[q]) fa[New] = q; 44 else 45 { 46 int Nq = ++cnt; len[Nq] = len[x] + 1; 47 a[Nq] = a[q]; 48 fa[Nq] = fa[q]; 49 fa[New] = fa[q] = Nq; 50 while(a[x][c] == q) a[x][c] = Nq, x = fa[x]; 51 } 52 Ans += len[New] - len[fa[New]]; 53 } 54 }S; 55 56 int main() 57 { 58 n = get(); 59 while(n--) 60 { 61 x = get(); 62 S.Add(x); 63 printf("%lld\n", Ans); 64 } 65 66 }
以上是关于BZOJ4516SDOI2016生成魔咒 [SAM]的主要内容,如果未能解决你的问题,请参考以下文章