【BZOJ 3676】[Apio2014]回文串
【链接】 链接
【题意】
给你一个字符串s.
定义一个子串的出现值为它出现的次数*字符串的长度。
让你求里面的回文子串的最大出现值
|s|
【题解】
马拉车算法里面.
只有在回文往外扩展超过mx的时候,才出现了不同的回文。
只有这个时候需要计算。
快速计算一个子串在字符串中出现了多少次。
可以用二分。
比如子串为s[x..y]
则先获取temp = Rank[x];
然后看看往左height>=y-x+1的能走多远。
再看看往右height>=y-x+1的有多少。
用RMQ看看height在temp,mid(二分的中点)之间的height最小值是多少
显然有单调性.
找到最左和最右,就是出现的次数了。
马拉车算法虽然改造了原串。
但可以记录每个位置原本是属于字符串的哪一个位置的。
【错的次数】
在这里输入错的次数
【反思】
在这里输入反思
【代码】
#include<bits/stdc++.h>
using namespace std;
const int N = 3e5;
const int MAX_CHAR = 255;//每个数字的最大值。
char s[N + 10];//如果是数字,就写成int s[N+10]就好,从0开始存
int Sa[N + 10], T1[N + 10], T2[N + 10], C[N+10];
int Height[N + 10], Rank[N + 10];
void build_Sa(int n, int m) {
int i, *x = T1, *y = T2;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[i] = s[i]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[i]]] = i;
for (int k = 1; k <= n; k <<= 1)
{
int p = 0;
for (i = n - k; i<n; i++) y[p++] = i;
for (i = 0; i<n; i++) if (Sa[i] >= k) y[p++] = Sa[i] - k;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[y[i]]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[y[i]]]] = y[i];
swap(x, y);
p = 1; x[Sa[0]] = 0;
for (i = 1; i<n; i++)
x[Sa[i]] = y[Sa[i - 1]] == y[Sa[i]] && y[Sa[i - 1] + k] == y[Sa[i] + k] ? p - 1 : p++;
if (p >= n) break;
m = p;
}
}
void getHeight(int n)
{
int i, j, k = 0;
for (i = 1; i <= n; i++) Rank[Sa[i]] = i;
for (i = 0; i<n; i++) {
if (k) k--;
j = Sa[Rank[i] - 1];
while (s[i + k] == s[j + k]) k++;
Height[Rank[i]] = k;
}
}
const int MAXL = 19;//log2数组的最大长度
const int INF = 0x3f3f3f3f;//数值绝对值的最大值
struct abc{
int pre2[MAXL+5],need[N+10];
int fmin[N+10][MAXL+5];
void init(int n)
{
pre2[0] = 1;
for (int i = 1;i <= MAXL;i++)
{
pre2[i] = pre2[i-1]<<1;
}
need[1] = 0; need[2] = 1;
int temp = 2;
for (int i = 3; i <= n; i++)//need[i]表示长度为i是2的多少次方,可以理解为[log2i]
if (pre2[temp] == i)
need[i] = need[i - 1] + 1, temp++;
else
need[i] = need[i - 1];
}
void getst(int *a,int n)
{
memset(fmin,INF,sizeof fmin);
for (int i = 1;i <= n;i++)//下标从0开始就改成对应的就好
fmin[i][0] = a[i];
for (int l = 1;pre2[l]<=n;l++)
for (int i = 1;i <= n;i++)
if (i+pre2[l]-1<=n)
fmin[i][l] = min(fmin[i][l-1],fmin[i+pre2[l-1]][l-1]);
}
int getmin(int l,int r)
{
int len = need[r-l+1];
return min(fmin[l][len],fmin[r-pre2[len]+1][len]);
}
}ST;
char S[2*N+100];
int belong[2*N+100],p[2*N+100];
int n;
long long ans = 0;
int lcp(int x,int y)
{
return ST.getmin(x+1,y);
}
void cal(int x,int y)
{
int len = y - x + 1;//子串s[x..y]的长度
int temp = Rank[x];//这个后缀x的排名
int pr = temp;//最右能到哪里
int l = temp+1,r = n;
while (l <= r)
{
int mid = (l+r)>>1;
if (lcp(temp,mid)>=len)
{
pr = mid;
l = mid + 1;
}else
r = mid - 1;
}
int pl = temp;
l = 1,r = temp - 1;
while (l <= r)
{
int mid = (l+r)>>1;
if (lcp(mid,temp)>=len)
{
pl = mid;
r = mid - 1;
}else
l = mid + 1;
}
ans = max(ans,1LL*(pr-pl+1)*len);
}
int main() {
//freopen("F:\\rush.txt","r",stdin);
scanf("%s", s);
n = strlen(s);
s[n] = 0;
build_Sa(n + 1, MAX_CHAR);//注意调用n+1
getHeight(n);
ST.init(n);
ST.getst(Height,n);
int len = 0;
S[len++] = '$';
S[len++] = '#';
for (int i = 0;i < n;i++)
{
S[len++] = s[i];
belong[len-1] = i;
S[len++] = '#';
}
int mx = 0,id = 0;
for (int i = 1;i < len;i++)
{
if (i > mx)
p[i] = 1;
else
p[i] = min(p[2*id-i],mx-i+1);
while (S[i-p[i]]==S[i+p[i]])
{
if (i+p[i] > mx && S[i-p[i]+1]!='#')
{
cal(belong[i-p[i]+1],belong[i+p[i]-1]);
}
p[i]++;
}
if (i + p[i]-1 > mx) mx = i + p[i] - 1,id = i;
}
printf("%lld\n",ans);
return 0;
}
【BZOJ 3676】[Apio2014]回文串的更多相关文章
- BZOJ 3676: [Apio2014]回文串
3676: [Apio2014]回文串 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 2013 Solved: 863[Submit][Status ...
- bzoj 3676: [Apio2014]回文串 回文自动机
3676: [Apio2014]回文串 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 844 Solved: 331[Submit][Status] ...
- 字符串(马拉车算法,后缀数组,稀疏表):BZOJ 3676 [Apio2014]回文串
Description 考虑一个只包含小写拉丁字母的字符串s.我们定义s的一个子串t的“出 现值”为t在s中的出现次数乘以t的长度.请你求出s的所有回文子串中的最 大出现值. Input 输入只有一行 ...
- ●BZOJ 3676 [Apio2014]回文串
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3676 题解: 后缀数组,Manacher,二分 首先有一个结论:一个串的本质不同的回文串的个 ...
- BZOJ 3676 [Apio2014]回文串(回文树)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3676 [题目大意] 考虑一个只包含小写拉丁字母的字符串s. 我们定义s的一个子串t的& ...
- bzoj 3676: [Apio2014]回文串【回文自动机】
回文自动机板子 或者是SAM+manacher+倍增,就是manacher求本质不同回文串(让f++的串),然后在SAM倍增查询对应点出现次数 #include<iostream> #in ...
- BZOJ 3676 [Apio2014]回文串 (后缀自动机+manacher/回文自动机)
题目大意: 给你一个字符串,求其中回文子串的长度*出现次数的最大值 明明是PAM裸题我干嘛要用SAM做 回文子串有一个神奇的性质,一个字符串本质不同的回文子串个数是$O(n)$级别的 用$manach ...
- bzoj 3676 [Apio2014]回文串(Manacher+SAM)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3676 [题意] 给定一个字符串,定义一个串的权值为长度*出现次数,求最大权的回文子串. ...
- BZOJ.3676.[APIO2014]回文串(回文树)
BZOJ 洛谷 很久之前写(抄)过一个Hash+Manacher的做法,当时十分懵逼=-= 然而是道回文树模板题. 回文树教程可以看这里(真的挺妙的). 顺便再放上MilkyWay的笔记~ //351 ...
- bzoj 3676: [Apio2014]回文串【后缀自动机+manacher】
用manacher找出本质不同的回文子串放在SAM上跑 #include<iostream> #include<cstdio> #include<cstring> ...
随机推荐
- GO语言学习(二)Windows 平台下 LiteIDE 的安装和使用
1. 安装 Go 语言并设置环境变量 参考GO语言学习(一) 2. MinGW 的下载和安装 Windows 下的 Go 调试还需要安装 MinGW. 2.1 下载安装工具的安装 最新版本下载安装工具 ...
- HibernateCRUD基础框架(2)-HQL语句构造器(HqlQueryBuilder,HqlUpdateBuilder)
上篇讲述了最基本的实体类,本篇接着讲述HQL语句构造器,包括查询和更新等. 优点:通过面向对象的方式构造HQL语句,更快捷,不需要手动拼接HQL. 缺点:封装可能降低性能,只能支持常用的和较为简单的H ...
- IOS总结
1.Difference between shallow copy and deep copy? 浅复制和深复制的区别? 答案:浅层复制:只复制指向对象的指针,而不复制引用对象本身. 深层复制:复制引 ...
- Altium Designer如何重命名文件
- HDU 2147kiki's game
KIKI和zz一起玩跳棋游戏,KIKI先.跳棋棋盘有n行m列.在顶行的最右侧位置放上一枚硬币.每次每个人可以把硬币移动到左边,下边或是左下边的空格中.最后不能移动硬币的那个人将输掉比赛. P点:即必败 ...
- eclipse编译器错误、警告设置
颜色配置步骤:Window->Preferences->General->Editors->Text Editors->Annotations
- [Angular2 Router] Auxiliary Routes bit by bit
Article Github Auxiliary Routes is is little bit hard to understand. Here is demo, link You can see ...
- 【Codeforces Round #185 (Div. 2) B】Archer
[链接] 链接 [题意] 在这里输入题意 [题解] 概率水题. 枚举它是第几轮成功的. 直到满足精度就好 [错的次数] 1 [反思] long double最让人安心. [代码] #include & ...
- 关于JavaScript对象概念的总结
原文 https://www.jianshu.com/p/88213b499c4b 大纲 前言 1.对象的相关概念 2.对象的创建(简单创建) 3.对象的属性 3.1.数据属性 3.2.访问器属性 4 ...
- 互联网支付的"二清"介绍
所谓"二清"公司,是针对"一清"机构而言. 的"一清"和"二清"一般属于POS收单业务范畴,"一清" ...