hdu【最长回文子串——Manacher算法】POJ 3974/HDU 3068

最新推荐文章于 2016-04-13 19:47:09 发布
最新推荐文章于 2016-04-13 19:47:09 发布
阅读量612 收藏
点赞数
分类专栏： ACM--数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u010793761/article/details/10180495
版权
转载
O（n）回文子串算法
    这里，我介绍一下O（n）回文串处理的一种方法。Manacher算法.
 原文地址：
http://zhuhongcheng.wordpress.com/2009/08/02/a-simple-linear-time-algorithm-for-finding-longest-palindrome-sub-string/
     其实原文说得是比较清楚的，只是英文的，我这里写一份中文的吧。
     首先：大家都知道什么叫回文串吧，这个算法要解决的就是一个字符串中最长的回文子串有多长。这个算法可以在O（n）的时间复杂度内既线性时间复杂度的情况下，求出以每个字符为中心的最长回文有多长，
     这个算法有一个很巧妙的地方，它把奇数的回文串和偶数的回文串统一起来考虑了。这一点一直是在做回文串问题中时比较烦的地方。这个算法还有一个很好的地方就是充分利用了字符匹配的特殊性，避免了大量不必要的重复匹配。
     算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符，当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了（见下面的一个例子，回文串长度全为奇数了），然后用一个辅助数组P记录以每个字符为中心的最长回文串的信息。P［id］记录的是以字符str［id］为中心的最长回文串，当以str［id］为第一个字符，这个最长回文串向右延伸了P［id］个字符。
     原串：    w aa bwsw f d 
     新串：   # w # a # a # b # w # s # w # f # d #
 辅助数组P：  1 2 1 2 3 2 1 2 1 2 1 4 1 2 1 2 1 2 1
    这里有一个很好的性质，P［id］-1就是该回文子串在原串中的长度（包括‘#’）。如果这里不是特别清楚，可以自己拿出纸来画一画，自己体会体会。当然这里可能每个人写法不尽相同，不过我想大致思路应该是一样的吧。
    好，我们继续。现在的关键问题就在于怎么在O（n）时间复杂度内求出P数组了。只要把这个P数组求出来，最长回文子串就可以直接扫一遍得出来了。
    由于这个算法是线性从前往后扫的。那么当我们准备求P［i］的时候，i以前的P［j］我们是已经得到了的。我们用mx记在i之前的回文串中，延伸至最右端的位置。同时用id这个变量记下取得这个最优mx时的id值。（注：为了防止字符比较的时候越界，我在这个加了‘#’的字符串之前还加了另一个特殊字符‘$’，故我的新串下标是从1开始的）
 好，到这里，我们可以先贴一份代码了。

[cpp]  view plain copy
  1. void pk()  
  2. {  
  3.     int i;  
  4.     int mx = 0;  
  5.     int id;  
  6.     for(i=1; i<n; i++)  
  7.     {  
  8.         if( mx > i )  
  9.             p[i] = MIN( p[2*id-i], mx-i );          
  10.         else  
  11.             p[i] = 1;  
  12.         for(; str[i+p[i]] == str[i-p[i]]; p[i]++)  
  13.             ;  
  14.         if( p[i] + i > mx )  
  15.         {  
  16.             mx = p[i] + i;  
  17.             id = i;  
  18.         }  
  19.     }  
  20. }  

代码是不是很短啊，而且相当好写。很方便吧，还记得我上面说的这个算法避免了很多不必要的重复匹配吧。这是什么意思呢，其实这就是一句代码。

if (  mx  >  i  )
    p [ i ]   =  MIN (  p [ 2 * id - i ],  mx - i  );

就是当前面比较的最远长度mx>i的时候，P［i］有一个最小值。这个算法的核心思想就在这里，为什么P数组满足这样一个性质呢?
   （下面的部分为图片形式）

 
    看完这个算法，你有可能会觉得这种算法在哪会用到呢？其实回文串后缀数组也可以做。只是复杂度是O（n log n）的，而且一般情况下也不会刻意去卡一个log n的算法。可正好hdu就有这么一题，你用后缀数组写怎么都得T（当然应该是我写得太烂了）。不信的话大家也可以去试试这题。
        http://acm.hdu.edu.cn/showproblem.php?pid=3068O （ n ）回文子串算法
    这里，我介绍一下O（n）回文串处理的一种方法。Manacher算法.
 原文地址：
http://zhuhongcheng.wordpress.com/2009/08/02/a-simple-linear-time-algorithm-for-finding-longest-palindrome-sub-string/
     其实原文说得是比较清楚的，只是英文的，我这里写一份中文的吧。
     首先：大家都知道什么叫回文串吧，这个算法要解决的就是一个字符串中最长的回文子串有多长。这个算法可以在O（n）的时间复杂度内既线性时间复杂度的情况下，求出以每个字符为中心的最长回文有多长，
     这个算法有一个很巧妙的地方，它把奇数的回文串和偶数的回文串统一起来考虑了。这一点一直是在做回文串问题中时比较烦的地方。这个算法还有一个很好的地方就是充分利用了字符匹配的特殊性，避免了大量不必要的重复匹配。
     算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符，当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了（见下面的一个例子，回文串长度全为奇数了），然后用一个辅助数组P记录以每个字符为中心的最长回文串的信息。P［id］记录的是以字符str［id］为中心的最长回文串，当以str［id］为第一个字符，这个最长回文串向右延伸了P［id］个字符。
     原串：    w aa bwsw f d 
     新串：   # w # a # a # b # w # s # w # f # d #
 辅助数组P：  1 2 1 2 3 2 1 2 1 2 1 4 1 2 1 2 1 2 1
    这里有一个很好的性质，P［id］-1就是该回文子串在原串中的长度（包括‘#’）。如果这里不是特别清楚，可以自己拿出纸来画一画，自己体会体会。当然这里可能每个人写法不尽相同，不过我想大致思路应该是一样的吧。
    好，我们继续。现在的关键问题就在于怎么在O（n）时间复杂度内求出P数组了。只要把这个P数组求出来，最长回文子串就可以直接扫一遍得出来了。
    由于这个算法是线性从前往后扫的。那么当我们准备求P［i］的时候，i以前的P［j］我们是已经得到了的。我们用mx记在i之前的回文串中，延伸至最右端的位置。同时用id这个变量记下取得这个最优mx时的id值。（注：为了防止字符比较的时候越界，我在这个加了‘#’的字符串之前还加了另一个特殊字符‘$’，故我的新串下标是从1开始的）
 好，到这里，我们可以先贴一份代码了。

[cpp]  view plain copy
  1. void pk()  
  2. {  
  3.     int i;  
  4.     int mx = 0;  
  5.     int id;  
  6.     for(i=1; i<n; i++)  
  7.     {  
  8.         if( mx > i )  
  9.             p[i] = MIN( p[2*id-i], mx-i );          
  10.         else  
  11.             p[i] = 1;  
  12.         for(; str[i+p[i]] == str[i-p[i]]; p[i]++)  
  13.             ;  
  14.         if( p[i] + i > mx )  
  15.         {  
  16.             mx = p[i] + i;  
  17.             id = i;  
  18.         }  
  19.     }  
  20. }  

代码是不是很短啊，而且相当好写。很方便吧，还记得我上面说的这个算法避免了很多不必要的重复匹配吧。这是什么意思呢，其实这就是一句代码。

if (  mx  >  i  )
    p [ i ]   =  MIN (  p [ 2 * id - i ],  mx - i  );

就是当前面比较的最远长度mx>i的时候，P［i］有一个最小值。这个算法的核心思想就在这里，为什么P数组满足这样一个性质呢?
   （下面的部分为图片形式）

 
    看完这个算法，你有可能会觉得这种算法在哪会用到呢？其实回文串后缀数组也可以做。只是复杂度是O（n log n）的，而且一般情况下也不会刻意去卡一个log n的算法。可正好hdu就有这么一题，你用后缀数组写怎么都得T（当然应该是我写得太烂了）。不信的话大家也可以去试试这题。
        http://acm.hdu.edu.cn/showproblem.php?pid=3068 

#include <map>  
#include <set>  
#include <list>  
#include <queue>  
#include <deque>  
#include <stack>  
#include <string>  
#include <time.h>  
#include <cstdio>  
#include <math.h>  
#include <iomanip>  
#include <cstdlib>  
#include <limits.h>  
#include <string.h>  
#include <iostream>  
#include <fstream>  
#include <algorithm>  
using namespace std;  
//HDU 3068  
#define LL long long  
#define MIN INT_MIN  
#define MAX INT_MAX  
#define PI acos(-1.0)  
#define FRE freopen("input.txt","r",stdin)  
#define FF freopen("output.txt","w",stdout)  
#define N 210010  
char str[N];  
int p[N];  
char s[N];  
void pk (int n) {  
    int i,j;  
    int mx = 0;  
    int id;  
    for (i = 1; i < n; i++) {  
        if (mx > i) {  
            p[i] = min(p[2*id - i], mx - i);  
        } else p[i] = 1;  
        for (; str[i - p[i]] == str[i + p[i]]; p[i]++);  
        if (p[i] + i > mx) {  
            mx = p[i] + i;  
            id = i;  
        }  
    }  
}  
void gao(int n) {  
    int i,j;  
    int mx = 0;  
    for (i = 1; i < n; i++) {  
        if (p[i] > mx) {  
            mx = p[i];  
        }  
    }  
    printf("%d\n",mx-1);  
}  
int main () {  
    while (scanf("%s",s+1) != -1) {  
        int i,j;  
        int len = strlen(s+1);  
        str[0] = '@';  
        str[1] = '#';  
        for (i = 2, j = 1; j <= len ; j++) {  
            str[i++] = s[j];  
            str[i++] = '#';  
        }  
        str[i] = '\0';  
        //cout<<str<<endl;  
        pk(i);  
        gao(i);  
    }  
    return 0;  
}  

[cpp] view plaincopy
#include <map>  
#include <set>  
#include <list>  
#include <queue>  
#include <deque>  
#include <stack>  
#include <string>  
#include <time.h>  
#include <cstdio>  
#include <math.h>  
#include <iomanip>  
#include <cstdlib>  
#include <limits.h>  
#include <string.h>  
#include <iostream>  
#include <fstream>  
#include <algorithm>  
using namespace std;  
//POJ 3974  
#define LL long long  
#define MIN INT_MIN  
#define MAX INT_MAX  
#define PI acos(-1.0)  
#define FRE freopen("input.txt","r",stdin)  
#define FF freopen("output.txt","w",stdout)  
#define N 2000020  
char str[N];  
int p[N];  
char s[N/2];  
void pk (int n) {  
    int i,j;  
    int mx = 0;  
    int id;  
    for (i = 1; i < n; i++) {  
        if (mx > i) {  
            p[i] = min(p[2*id - i], mx - i);  
        } else p[i] = 1;  
        for (; str[i - p[i]] == str[i + p[i]]; p[i]++);  
        if (p[i] + i > mx) {  
            mx = p[i] + i;  
            id = i;  
        }  
    }  
}  
void gao(int n) {  
    int i,j;  
    int mx = 0;  
    for (i = 1; i < n; i++) {  
        if (p[i] > mx) {  
            mx = p[i];  
        }  
    }  
    printf("%d\n",mx-1);  
}  
int main () {  
    int ca = 1;  
    while (scanf("%s",s+1)) {  
        if (strcmp(s+1,"END") == 0) break;  
        int i,j;  
        int len = strlen(s+1);  
        str[0] = '@';  
        str[1] = '#';  
        for (i = 2, j = 1; j <= len ; j++) {  
            str[i++] = s[j];  
            str[i++] = '#';  
        }  
        str[i] = '\0';  
        printf("Case %d: ",ca++);  
        pk(i);  
        gao(i);  
    }  
    return 0;  
}




poj 3974 
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int LEN=2200050;
int p[LEN];
char str[LEN], dstr[LEN];
int getstring(int n)
{
	int i, mx=0, size=0;
	for(i=1; i<n; i++)
	{
		if( mx>i )
			p[i]=min(p[2*size-i], mx-i);
		else
			p[i]=1;
		for(; dstr[i+p[i]]==dstr[i-p[i]]; p[i]++);
		if( i+p[i]>mx )
		{
			mx=p[i]+i;
			size=i;
		}
	}
	mx=-1;
	for(i=1; i<n; i++)
	{
		if( mx<p[i] )
		{
			mx=p[i];
		}
	}
	return mx;
}
int main()
{
	int i, j, len, cases=1;
	while( scanf("%s", str+1)!=EOF )
	{
		if( strcmp(str+1, "END")==0 )
			break;
		len=strlen(str+1);
		dstr[0]='@';
		dstr[1]='#';
		for(i=2, j=1; j<=len; dstr[i++]=str[j++], dstr[i++]='#');
		dstr[i]='\0';
		int ans=getstring(i);
		printf("Case %d: %d\n",cases++, ans-1);
	}
	return 0;
}





u010793761
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
algorithm:记录我的在线裁判（hdu，poj等）的算法问题代码
05-17
算法从在线裁判（hdu，poj等）记录我的算法问题代码。
【最长回文子串——Manacher算法】POJ 3974/HDU 3068
leolin_的专栏
11-12 971
转载： O（n）回文子串算法     这里，我介绍一下O（n）回文串处理的一种方法。Manacher算法. 原文地址： http://zhuhongcheng.wordpress.com/2009/08/02/a-simple-linear-time-algorithm-for-finding-longest-palindrome-sub-string/     其实原文
POJ 3974 HDU 3068 最长回文子串 Manacher 算法
lentty1452的专栏
07-22 334
今天看扩展KMP的时候无意看到这个算法的，这个算法能在线性时间内算出最长回文串的长度，真是强大啊，花了好长时间看懂原理，看懂网上的代码，顺便A了两道题 HAPPY   给几个好的链接http://blog.sina.com.cn/s/blog_70811e1a01014esn.html http://bbs.dlut.edu.cn/bbstcon.php?board=Competition&
回文子串算法——Manacher学习
theArcticOcean
04-13 616
首次学习manacher算法，感谢分享优秀博文的博主： http://blog.csdn.net/xingyeyongheng/article/details/9310555 读完后，我对其有了大体的认识。经过manacher算法处理后的字符串，其内容完全改变。比如字符串212： 设p[i]是i位置处向外拓展的半径。 after the change， it’s type conte
字符串 专题训练 · 最长回文子串之Manacher算法
山外小楼夜听雨
08-16 591
在字符串问题中，有一类经典的问题是求字符串中的最长回文子串，而解决这类问题的算法也很多，例如后缀树或者分治+拓展KMP。但是后缀树的极为复杂，没有实用性；分治+拓展KMP的复杂度达到了O(NlogN)，并不能算是非常高效。对此，我们可以学习一种专门的算法：Manacher算法。 Manacher算法的核心可以用一句话来概括：在回文子串中找回文子串，带上这个思想应该会更好理解。 可以参看星夜永恒
北大杭电acm题解（详细）
06-19
北大杭电acm题解（详细）北大杭电acm题解（详细）北大杭电acm题解（详细）北大杭电acm题解（详细）北大杭电acm题解（详细）北大杭电acm题解（详细）
poj题目3034
10-23
poj上的题目3034，利用动态规划的思想解答。
leetcode和oj-ACM:算法和编码练习。试题来自POJ、HDUOJ和Leetcode
06-30
leetcode 和 oj ACM Practice for algorithm and coding. The exam questions come from POJ, HDU OJ and Leetcode. 母函数 数学 排序
leetcode和hdoj-Algorithm-Packages:算法包
07-06
主要用来记录算法刷题记录和一些模板 文件结构 leetcode 存放leetcode题目和周赛 atcoder 用于存放参与和vp的atcoder比赛 codeforces 用于存放参与和vp的cf比赛，比赛文件夹以比赛序号和div描述，平时练习的题目放在...
求解x=a^b(mod m)的几种方法
热门推荐
janvinal的专栏
03-29 1万+
求解x=a^b(mod m) Posted on 2012 年 11 月 6 日 by NARUTOACM — 14 Comments ↓ 本文致力于解决如下问题：求解x≡ab(mod m)，其中a,b,m都是正整数。 如果b足够小，则可直接用逐次平方法求解，如果你不知道逐次平方法，可以先看这里。所以这里假设b足够大（这不是说是一个64位整数，而是可以上百上千位的一个数），大到逐次平
输出1/n的小数位，是循环小数输出一个循环节
janvinal的专栏
10-22 1893
/* 从这里充分体现了C++输入输出和C的输入输出 在数大的时候速度快很多。 */ #include #include #include #include #include #include using namespace std; const int LEN=100010; int main() { int k, n, T, vis[LEN], digit[LEN]; cin
poj2769之寻找最小完全剩余系
janvinal的专栏
12-11 1636
题目解析：该题所求即为寻找到最小的数 m 满足 m 模G个数两两不同余，即为找到最小的 m 使得这些数在 m 的最小完全剩余系中。 #include #include #include using namespace std; #define LENN 1000010 #define LENV 100001 int num[LENN], vis[LENV]; int main(int argc,
POJ 2447 破解RSA(经典公钥算法)
janvinal的专栏
10-07 1402
周五刚好在俞研的网络安全课上学了RSA，回来想实现以下，由于以前数论方面的积累还算比较深厚，很快就过了这一题，呵呵:-) 总结一下吧，这题可以说是数论部分的一个大综合题，因为它将算法导论上数论这部分的知识点全部包含了进来，包括gcd,扩展gcd,模线性方程，a^b mod c(还是比较难的那种，相关题目可以看一下FOZ上面的2道题),miller-rabin素数测试，pollard_rho质因数
数论专题学习阶段总结
janvinal的专栏
08-23 984
转载：http://www.cnblogs.com/kuangbin/archive/2012/08/28/2661066.html 从放暑假前周sir给我讲了一个用polya计数法和burnside定理做的题目（pku2409）后，突然觉得组合数学挺有意思，然后从那时起到现在几乎都在做这类的题目。   做到现在感觉这类题目的一些基本知识点都差不多有所了解了，水
hdu2179 求 PI，精确到小数点后1500位，麦金公式，java大数
janvinal的专栏
03-12 922
麦金公式 #include #include int main() { double s; int b,x,n,c,i,j,d,l,a[1500]; while(scanf("%d",&n)!=EOF&&n!=0) { x=5*n; for(s=0,n=1; n<=5000; n++) {
N！最末非零位数值
janvinal的专栏
11-10 882
给定一个数N（N 在N很小的情况下，有一个简便的方法：求出1到N之间每个数的2的因子数和5的因子数，记为F(2)和F(5)，显然F(2) >= F(5)。由于在末尾只有2和5相乘才能产生0，如果我们把2和5抛去，那么肯定不会有0，这样就可以一边乘一边模10，防止溢出。剩下的一堆2和5如何处理呢？因为2肯定比5多，因此最末位肯定是偶数（0的阶乘和1的阶乘除外）。而一个偶数不停地乘2，最末位的规律
ACM数论题集
janvinal的专栏
10-22 837
Volume 1 http://acm.hdu.edu.cn/showproblem.php?pid=1005 http://acm.hdu.edu.cn/showproblem.php?pid=1014 http://acm.hdu.edu.cn/showproblem.php?pid=1019 Volume 2 http://acm.hdu.edu.cn/showproblem.
hdu 2207 子网掩码转换确定
janvinal的专栏
04-14 805
不知道你是否注意过你的电脑上下面的网络信息： “连接状态： 地址类型：          通过DHCP指派   IP地址：          60.176.38.23   子网掩码：          255.255.252.0   默认网关：          60.176.36.1   首选DNS服务器：      210.32.32.1   备用DNS服务器：      210
请给出hdu5113的算法思路
最新发布
05-19
算法思路： 首先我们可以计算出数组的平均值 $avg=\dfrac{a_1+a_2+...+a_n}{n}$，如果平均值不是整数，那么无论如何都不能通过操作将所有元素变为相等的。如果平均值是整数，那么我们可以枚举将所有元素变成 $avg,...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值