首页 > 其他 > 详细

最长回文串长度——manacher算法

时间:2014-02-19 13:55:42      阅读:331      评论:0      收藏:0      [点我收藏+]

网上有不少对manacher算法的解释,但总会存在某些漏洞,在这里给一个整理过的版本,并用Java给予实现


首先,在字符串srad[i]表示第i个字符的回文半径rad[i]尽可能大且满足

s[i-rad[i],i-1] = s[i+1,i+rad[i]]

很明显求出了所有的rad,就求出了所有的长度为奇数的回文子串

至于偶数的怎么求,最后再讲


假设现在求出了rad[1..i-1],现在要求后面的rad并且通过前面的操作得知了当前字符irad值至少为j,现在通过试图扩大j来扫描求出了rad[i]。再假设现在有个指针k,1循环到rad[i],试图通过某些手段来求出[i+1,i+rad[i]]rad

根据定义黑色的部分是一个回文子串两段红色的区间全等

因为之前已经求出了rad[i-k],所以直接用它.3种情况

(1)rad[i]-k < rad[i-k]

bubuko.com,布布扣

如图,rad[i-k]的范围为青色因为黑色的部分是回文的且青色的部分超过了黑色的部分所以rad[i+k]肯定至少为rad[i]-k,即橙色的部分那橙色以外的部分就不是了吗这是肯定的因为如果橙色以外的部分也是回文的那么根据青色和红色部分的关系可以证明黑色部分再往外延伸一点也是一个回文子串这肯定不可能,因此rad[i+k] = rad[i]-k


(2)rad[i]-k > rad[i-k]

bubuko.com,布布扣

如图,rad[i-k]的范围为青色因为黑色的部分是回文的且青色的部分在黑色的部分里面根据定义很容易得出:rad[i+k] = rad[i-k]

 

根据上面两种情况可以得出结论rad[i]-k! = rad[i-k]的时候,rad[i+k] = min(rad[i]-k,rad[i-k])。


(3)rad[i]-k = rad[i-k]

bubuko.com,布布扣

如图通过和第一种情况对比之后会发现因为青色的部分没有超出黑色的部分所以即使橙色的部分全等也无法像第一种情况一样引出矛盾因此橙色的部分是有可能全等的但是根据已知的信息我们不知道橙色的部分是多长因此就把i指针移到i+k的位置,j=rad[i-k](因为它的rad值至少为rad[i-k]),等下次循环的时候再做了

整个算法就这样

至于时间复杂度为什么是O(n)


上文还留有一个问题就是这样只能算出奇数长度的回文子串偶数的就不行怎么办呢?

有一种直接但比较笨的方法,就是做两遍(因为两个程序是差不多的只是rad值的意义和一些下标变了而已)。但是写两个差不多的程序是很痛苦的而且容易错所以一种比较好的方法就是在原来的串中每两个字符之间加入一个特殊字符再做

:aabbaca把它变成a#a#b#b#a#c#a,这样的话无论原来的回文子串长度是偶数还是奇数现在都变成奇数了


/**
 * 对称字符串的最大长度(字符串)<br/>
 * 题目:输入一个字符串,输出该字符串中对称的子字符串的最大长度。<br/>
 * 比如输入字符串“google”,由于该字符串里最长的对称子字符串是“goog”,因此输出4。
 * 
 * @author Joeson
 * 
 */
public class Algorithm073
{
	public static void main(String[] args)
	{
		String s = "acbcaacbcd";
		char[] str = init(s);
		print(str);
		manacher(str);

	}

	/**
	 * 返回例如 #a#c#b#c#a#a#c#b#c#d#形式的字符串数组
	 * 
	 * @param s
	 * @return
	 */
	public static char[] init(String s)
	{
		char[] str = new char[s.length() * 2 + 1];

		int i = 0;
		for (; i < s.length(); i++)
		{
			str[2 * i] = ‘#‘;
			str[2 * i + 1] = s.charAt(i);
		}
		str[2 * i] = ‘#‘;

		return str;
	}

	/**
	 * 
	 * 
	 * @param str
	 */
	public static void manacher(char[] s)
	{
		int rad[] = new int[s.length];

		int i = 1, j = 0, k;

		// 记录最长的回文串的长度
		int maxLen = 0;
		while (i < s.length)
		{
			// 扫描得出rad值
			while (i - j - 1 > -1 && i + j + 1 < s.length
					&& s[i - j - 1] == s[i + j + 1])
				j++;
			rad[i] = j;

			maxLen = maxLen > j ? maxLen : j;

			k = 1;
			while (k <= rad[i] && rad[i - k] != rad[i] - k)
			{
				rad[i + k] = Math.min(rad[i - k], rad[i] - k);
				k++;
			}
			i = i + k;
			j = Math.max(j - k, 0);
		}

		print(rad);
		System.out.println("最长回文串长度: " + maxLen);
	}

	private static void print(char[] rad)
	{
		for (int i = 0; i < rad.length; i++)
		{
			System.out.print(rad[i] + " ");
		}

		System.out.println();
	}

	private static void print(int[] rad)
	{
		for (int i = 0; i < rad.length; i++)
		{
			System.out.print(rad[i] + " ");
		}

		System.out.println();
	}
}


结果:

bubuko.com,布布扣



最长回文串长度——manacher算法

原文:http://blog.csdn.net/chenxuegui1234/article/details/19418907

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!