BZOJ 3676 UOJ 103 APIO 2014 后缀自动机 Manacher
来源:互联网 发布:淘宝哪里申请换货 编辑:程序博客网 时间:2024/04/29 06:17
对于所有的本质不同的回文子串,在SAM中查找对应节点,由于我们已知回文子串的右端点,因此可以在Parent树中快速定位我们要的状态。
回文自动机很厉害
所以我肯定学不会
#include <cstdio>#include <cstring>#include <algorithm>#define FOR(i,j,k) for(i=j;i<=k;++i)using namespace std;const int N = 610000, rt = 1;int last = 1, cnt = 1;long long ans = 0;int ch[N][26], ma[N], fa[N][19];int dep[N], pos[N], p[N], b[N], c[N], v[N], n;char a[N];void add(char c) { int np = ++cnt, p = last; last = np; ma[np] = ma[p] + 1; v[np] = 1; pos[ma[np]] = np; while (p && !ch[p][c]) ch[p][c] = np, p = fa[p][0]; if (!p) fa[np][0] = rt; else { int q = ch[p][c]; if (ma[q] == ma[p] + 1) fa[np][0] = q; else { int nq = ++cnt; ma[nq] = ma[p] + 1; memcpy(ch[nq], ch[q], sizeof ch[q]); fa[nq][0] = fa[q][0]; fa[q][0] = fa[np][0] = nq; while (p && ch[p][c] == q) ch[p][c] = nq, p = fa[p][0]; } }}void preprocess() { int i, j, x; FOR(i,1,cnt) ++b[ma[i]]; FOR(i,1,n) b[i] += b[i - 1]; for(i=cnt;i;--i) c[b[ma[i]]--] = i; for(i=cnt;i;--i) v[fa[c[i]][0]] += v[c[i]]; FOR(i,1,cnt) { x = c[i]; dep[x] = dep[fa[x][0]] + 1; for (j = 1; (1 << j) <= dep[x]; ++j) fa[x][j] = fa[fa[x][j - 1]][j - 1]; }}void query(int l, int r) { int p = pos[r], i, x; for (i = 18; i >= 0; --i) { x = fa[p][i]; if (ma[x] >= r - l + 1) p = x; } ans = max(ans, 1ll * v[p] * (r - l + 1));}void solve() { int mx, i, id; a[0] = '+'; a[n + 1] = '-'; for (i = 1, mx = 0; i <= n; ++i) { if (mx > i) p[i] = min(mx - i, p[2 * id - i - 1]); else p[i] = 0; while (a[i + p[i] + 1] == a[i - p[i]]) { ++p[i]; query(i - p[i] + 1, i + p[i]); } if (i + p[i] > mx) mx = i + p[i], id = i; } for (i = 1, mx = 0; i <= n; ++i) { if (mx > i) p[i] = min(mx - i - 1, p[2 * id - i]); else { p[i] = 1; query(i - p[i] + 1, i + p[i] - 1); } while (a[i - p[i]] == a[i + p[i]]) { ++p[i]; query(i - p[i] + 1, i + p[i] - 1); } if (i + p[i] > mx) mx = i + p[i], id = i; }}int main() { scanf("%s", a + 1); for (int i = 1; a[i]; ++i) add(a[i] - 'a'); n = ma[last]; preprocess(); solve(); printf("%lld", ans); return 0;}
为了URAL 2059,还是得回来打PAM QwQ
#include <cstdio>#include <algorithm>using namespace std;const int N = 300005;typedef long long ll;struct PAM { int cnt, last; char str[N]; int ch[N][26], fa[N], ma[N], sz[N]; PAM(): cnt(1) { fa[0] = fa[1] = 1; ma[1] = -1; } void extend(int c, int n) { int p = last; while (str[n - ma[p] - 1] != str[n]) p = fa[p]; if (!ch[p][c]) { int np = ++cnt, k = fa[p]; ma[np] = ma[p] + 2; while (str[n - ma[k] - 1] != str[n]) k = fa[k]; fa[np] = ch[k][c]; ch[p][c] = np; } sz[last = ch[p][c]]++; } void construct() { for (int i = 1; str[i]; ++i) extend(str[i] - 'a', i); for (int i = cnt; i; --i) sz[fa[i]] += sz[i]; } ll solve() { ll ans = 0; for (int i = cnt; i; --i) ans = max(ans, (ll)sz[i] * ma[i]); return ans; }} pam;int main() { scanf("%s", pam.str + 1); pam.construct(); printf("%lld\n", pam.solve()); return 0;}
3676: [Apio2014]回文串
Time Limit: 20 Sec Memory Limit: 128 MB
Submit: 1487 Solved: 623
[Submit][Status][Discuss]
Description
考虑一个只包含小写拉丁字母的字符串s。我们定义s的一个子串t的“出
现值”为t在s中的出现次数乘以t的长度。请你求出s的所有回文子串中的最
大出现值。
Input
输入只有一行,为一个只包含小写字母(a -z)的非空字符串s。
Output
输出一个整数,为逝查回文子串的最大出现值。
Sample Input
【样例输入l】
abacaba
【样例输入2]
www
Sample Output
【样例输出l】
7
【样例输出2]
4
HINT
一个串是回文的,当且仅当它从左到右读和从右到左读完全一样。
在第一个样例中,回文子串有7个:a,b,c,aba,aca,bacab,abacaba,其中:
● a出现4次,其出现值为4:1:1=4
● b出现2次,其出现值为2:1:1=2
● c出现1次,其出现值为l:1:l=l
● aba出现2次,其出现值为2:1:3=6
● aca出现1次,其出现值为1=1:3=3
●bacab出现1次,其出现值为1:1:5=5
● abacaba出现1次,其出现值为1:1:7=7
故最大回文子串出现值为7。
【数据规模与评分】
数据满足1≤字符串长度≤300000。
0 0
- BZOJ 3676 UOJ 103 APIO 2014 后缀自动机 Manacher
- [manacher 后缀自动机 || 回文自动机] BZOJ 3676 [Apio2014]回文串
- NKOJ 2844 (APIO 2014)回文串(Manacher+后缀自动机+倍增/回文树)
- [BZOJ 3676][APIO 2014]回文串
- [bzoj 3676][uoj #103]【APIO2014】Palindromes回文串 后缀数组+manachar
- bzoj 3998(后缀自动机)
- bzoj 2946(后缀自动机)
- BZOJ 4516 后缀自动机
- BZOJ 3676 & UOJ 103 [Apio2014]回文串
- BZOJ 2882 工艺 后缀自动机
- BZOJ 2555 SubString 后缀自动机
- 【后缀自动机】 BZOJ 2555 SubString
- BZOJ 2555 SubString 后缀自动机
- 【后缀自动机】[BZOJ 2806]Cheat
- bzoj 2555(linkcuttree+后缀自动机)
- BZOJ 3238 差异 [后缀自动机]
- manacher 后缀数组 AC自动机 回文自动机 知识点讲解 课件
- BZOJ3676: [Apio2014]回文串(后缀自动机+manacher/回文自动机)
- Java网络编程-初识Socket
- Python实现MapReduce
- 支付宝钱包系统架构内部剖析(架构图)
- iOS开发之xib 加载xib视图需要用到的方法
- Activity重新创建之recreate
- BZOJ 3676 UOJ 103 APIO 2014 后缀自动机 Manacher
- FOJ 2214 Knapsack problem 第六届福建省大学生程序设计竞赛 C 01背包DP变种
- SQLite数据库笔记
- dstat默认输出说明
- Eclipse使用入门教程
- 什么是库?(linux平台下的理解)
- 正则表达式在iOS中的运用
- 修改JAVA_HOME无效,java版本保持不变的问题解决
- 简单的Timer