HDU - 3068 最长回文(manacher算法)
题意:给出一个只由小写英文字符a,b,c...y,z组成的字符串S,求S中最长回文串的长度.
分析:
manacher算法:
1、将字符串中每个字符的两边都插入一个特殊字符。(此操作的目的是,将字符串长度统一变成奇数,道理很容易想---奇数+偶数=奇数or偶数+奇数=奇数)
eg:abba--->#a#b#b#a#
此外,为了便于处理边界,可在字符串开始处再加另外一个特殊字符。
eg:#a#b#b#a#--->$#a#b#b#a#
2、数组 P[i] 来记录以字符S[i]为中心的最长回文子串向左/右扩张的长度,显然,P[i]-1正好是原字符串中以字符S[i]为中心的回文串长度。
3、计算P[i]。
id---已知的 右边界最大 的回文子串ss的中心
ma---id+P[id],子串ss的右边界
j--- i 关于 id 的对称点,因此j=id * 2 - i
(1)如果ma > i,则
(i)若p[j] < ma - i, 则p[i] = p[j]。
原因:因为以id为中心是个大的回文串,若已知 在位置id左侧的回文串内 有 某点j ,且以点j为中心有个长度为p[j](p[j]是小于ma-i的)的回文串,那么 以点i(i是j关于id的对称点,i在id右侧且因为对称,所以一定在回文串内)为中心一定有一个回文串,且该回文串长度为p[j]。
(ii)若p[j] >= ma - i, 则p[i] = ma - i。
原因:因为以id为中心的大的回文串右边界只延伸到ma,根据对称性,虽然p[j]大于ma-i,但是,对于以点i为中心的回文串,只能知道右边界为ma-i的范围内一定是个回文串,再往右情况未知,只能左右延伸依次判断。
(2)如果ma <= i,只能左右延伸依次判断。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<cmath>
#include<iostream>
#include<sstream>
#include<iterator>
#include<algorithm>
#include<string>
#include<vector>
#include<set>
#include<map>
#include<stack>
#include<deque>
#include<queue>
#include<list>
#define lowbit(x) (x & (-x))
const double eps = 1e-8;
inline int dcmp(double a, double b){
if(fabs(a - b) < eps) return 0;
return a > b ? 1 : -1;
}
typedef long long LL;
typedef unsigned long long ULL;
const int INT_INF = 0x3f3f3f3f;
const int INT_M_INF = 0x7f7f7f7f;
const LL LL_INF = 0x3f3f3f3f3f3f3f3f;
const LL LL_M_INF = 0x7f7f7f7f7f7f7f7f;
const int dr[] = {0, 0, -1, 1, -1, -1, 1, 1};
const int dc[] = {-1, 1, 0, 0, -1, 1, -1, 1};
const int MOD = 1e9 + 7;
const double pi = acos(-1.0);
const int MAXN = 110000 + 10;
const int MAXT = 10000 + 10;
using namespace std;
char s[MAXN];
int p[MAXN << 1];
char t[MAXN << 1];
int len;
int cnt;
void init(){
memset(t, 0, sizeof t);
cnt = 0;
t[cnt++] = '$';
for(int i = 0; i < len; ++i){
t[cnt++] = '#';
t[cnt++] = s[i];
}
t[cnt++] = '#';
}
int manacher(){
int ma = 0, id = 0, ans = 1;
for(int i = 1; i < cnt; ++i){
p[i] = ma > i ? min(p[id * 2 - i], ma - i) : 1;
while(t[i + p[i]] == t[i - p[i]]) ++p[i];
if(i + p[i] > ma){
ma = i + p[i];
id = i;
}
ans = max(ans, p[i]);
}
return ans - 1;
}
int main(){
while(scanf("%s", s) == 1){
len = strlen(s);
init();
printf("%d\n", manacher());
}
return 0;
}
HDU - 3068 最长回文(manacher算法)的更多相关文章
- hdu 3068 最长回文 manacher算法(视频)
感悟: 首先我要Orz一下qsc,我在网上很难找到关于acm的教学视频,但偶然发现了这个,感觉做的很好,链接:戳戳戳 感觉这种花费自己时间去教别人的人真的很伟大. manacher算法把所有的回文都变 ...
- HDU 3068 最长回文 manacher 算法,基本上是O(n)复杂度
下面有别人的比较详细的解题报告: http://wenku.baidu.com/view/3031d2d3360cba1aa811da42.html 下面贴我的代码,注释在代码中: #include ...
- HDU 3068 最长回文 Manacher算法
Manacher算法是个解决Palindrome问题的O(n)算法,能够说是个超级算法了,秒杀其它一切Palindrome解决方式,包含复杂的后缀数组. 网上非常多解释,最好的解析文章当然是Leetc ...
- HDU - 3068 最长回文manacher马拉车算法
# a # b # b # a # 当我们遇到回判断最长回文字符串问题的时候,若果用暴力的方法来做,就是在字符串中间添加 #,然后遍历每一个字符,找到最长的回文字符串.那么马拉车算法就是在这个基础上进 ...
- hdu 3068 最长回文 manacher
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3068 给出一个只由小写英文字符a,b,c...y,z组成的字符串S,求S中最长回文串的长度.回文就是正 ...
- hdu 3068 最长回文(manacher&最长回文子串)
最长回文 Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submi ...
- hdu 3068 最长回文(manacher入门)
最长回文 Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submis ...
- HDU - 3068 最长回文(manacher)
HDU - 3068 最长回文 Time Limit: 2000MS Memory Limit: 32768KB 64bit IO Format: %I64d & %I64u Subm ...
- hdu 3068 最长回文(manachar求最长回文子串)
题目连接:hdu 3068 最长回文 解题思路:通过manachar算法求最长回文子串,如果用遍历的话绝对超时. #include <stdio.h> #include <strin ...
- hdu_3068 最长回文(Manacher算法)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3068 最长回文 Time Limit: 4000/2000 MS (Java/Others) M ...
随机推荐
- redis api-String
- Mac如何升级自带的vim
brew install vim --with-lua --with-override-system-vi brew install macvim --with-lua --with-override ...
- 范数(norm)
[范数定义] 非负实值函数(非线性) 1)非负性: || a || >= 0 2)齐次性: || ka || = |k| ||a|| 3)三角不等式: || a + b || <= || ...
- 一 Hibernate入门
Hibernate环境搭建 Hibernate的API Hibernate的CRUD EE三层结构: web层 业务逻辑层 持久层 jdbc,DBUTils,Hibernate Hib ...
- 设置zabbix (3.4.2)添加监控项,触发器,让CPU使用超过85%就报警:
zabbix (3.4.2)添加监控项,触发器,让CPU使用超过85%就报警: zabbix自带模板有一个 Template OS Linux模板.这个模板有监控CPU的监控项,如果没有添加一个监控项 ...
- MQTT 协议学习:006-订阅主题 与 对应报文(SUBSCRIBE、SUBACK、UNSUBSCRIBE、UNSUBACK)
背景 之前我们提到了怎么发布消息对应的报文:现在我们来看,订阅一个主题的报文是怎么样的. SUBSCRIBE - 订阅主题 客户端向服务端发送SUBSCRIBE报文用于创建一个或多个订阅.每个订阅注册 ...
- 51nod 1378:夹克老爷的愤怒 很好玩的一道树状dp
1378 夹克老爷的愤怒 基准时间限制:1 秒 空间限制:131072 KB 分值: 80 难度:5级算法题 收藏 取消关注 夹克老爷逢三抽一之后,由于采用了新师爷的策略,乡民们叫苦不堪,开始组织 ...
- ESX/ESXi 主机上的每个插槽中安装了多少内存
要确定在 ESX/ESXi 主机上的每个插槽中安装了多少内存,请执行以下操作: 1. 启动ssh服务 2. 登陆esxi主机查看 使用 SSH 客户端登录主机,以 root 用户身份运行以下命令之一: ...
- oracle 查看表中有多少字段
select count(*) from user_tab_columns t where t.TABLE_NAME='WPM_CHECK_ORDER'
- window 命令行telnet 不能用问题
如图 解决办法 打开window控制面板,启用或关闭window功能,勾选telnet选项.