【HDU3068】最长回文-O(n)复杂度的Manacher算法

测试地址（题目是中文，自己看就好了）：最长回文

做法：求字符串内的最长回文子串也是一个老生常谈的问题了，不过大多数OIer以前学的都是O(n^2)的做法（穷举或者分治），而Manacher算法就是一个优秀的，能用O(n)的复杂度解决这个问题的算法。关于Manacher算法的教程网上有很多，我这里贴一个我觉得写得最好的，让我一看就看得懂的：Manacher算法。所以这道题就可以说是Manacher算法的模板题了。

以下是本人代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int po,len[220010],lens;
char s[220010],ts[110010];

void init()
{
  lens=strlen(ts);
  s[0]='@';
  for(int i=1;i<=lens;i++)
  {
    s[2*i-1]='#';
	s[2*i]=ts[i-1];
  }
  s[2*lens+1]='#';
  s[2*lens+2]='$';
  s[2*lens+3]=0;
  lens=2*lens+1;
}

void Manacher()
{
  len[1]=1;po=1;
  for(int i=2;i<=lens;i++)
  {
    if (i>po+len[po]-1)
	{
	  int j=1;
	  while(s[i-j]==s[i+j]) j++;
	  len[i]=j;
	  po=i;
	}
	else
	{
	  int k=2*po-i;
	  if (i+len[k]<po+len[po]) len[i]=len[k];
	  else
	  {
	    int j=po+len[po]-i;
		while(s[i-j]==s[i+j]) j++;
		len[i]=j;
		po=i;
	  }
	}
  }
}

int main()
{
  while(scanf("%s",ts)!=EOF)
  {
    init();
    Manacher();
    int ans=0;
    for(int i=1;i<=lens;i++)
      ans=max(ans,len[i]);
    printf("%d\n",ans-1);
  }
  
  return 0;
}