女神教你字符串——导字符串（回文）

Link：http://acdream.info/contest?cid=1320#problem-F

女神教你字符串——导字符串

Time Limit:  2000/1000MS (Java/Others)     Memory Limit:  128000/64000KB (Java/Others)

Submit  Status

Problem Description

正如大家知道的，女神喜欢字符串，而在字符串中，女神最喜欢回文字符串，但是不是所有的字符串都是回文字符串，但是有一些字符串可以进行“求导”来变成回文字符串。

字符串中只包含小写字母。
求导过程如下，C++：

string dif(const string x)
{
    if(x.length()<=1)
        return "";
    string res="";
    for(int i=1;i<x.length();++i)
        res+=abs(x[i]-x[i-1])+'a';
    return res;
}

C：

void dif(char*x,char*res)//注意有可能会溢出
{
    if(x[0]==0||x[1]==0)
    {
        res[0]=0;
        return;
    }
    int len=1;
    for(int i=1;x[i];++i,++len)
        res[i-1]=abs(x[i]-x[i-1])+'a';
    res[len-1]=0;
}

例如"aa"的导字符串是“a"，”aab“的导字符串是"ab"，"aacfwssg"的导字符串是"acdream"。

那么给定一个字符串，请判断在它各阶导字符串中，最长的回文子串是多长？

二阶导字符串即为导字符串的导字符串。

n阶导字符串即为n-1阶导数的导字符串。

Input

多组数据，每组数据包括一个字符串s(1<=|s|<=1000)

Output

对于每组数据，输出一个整数

Sample Input

abcd
abcba
acdream

Sample Output

3
5
3

Hint

样例一，求一次导字符串后为”aaa“,最长回文子串就是本身，所以长度为3

样例二，本身就是回文串，因此，输出本身长度即可

样例三，acdream->cbonem->bnbji->mmib->aeh->ed->b

最长回文串长度分别为1,1,3,2,1,1,1,因此输出3

Submit

AC code：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
#define LL long long
#define MAXN 1000010
using namespace std;
//char s1[MAXN],s2[MAXN];
string s,s1,s2,s3,s4;
int p[MAXN];
int len,ans,cnt;
void pre()
{
	int i,j,k;
	len=s1.length();
	s2[0]='$';
	s2[1]='#';
	for(i=0;i<len;i++)
	{
		s2[i*2+2]=s1[i];
		s2[i*2+3]='#';
	}
	len=len*2+2;
	s2[len]=0;
}
int Manacher()
{
	int i;
	int mx=0;
	int id;
	ans=0;
	for(i=1;i<len;i++)
	{
		if(mx>i)
			p[i]=min(p[2*id-i],mx-i);
		else
			p[i]=1;
		for(;s2[i+p[i]]==s2[i-p[i]];p[i]++);
		ans=max(ans,p[i]);
		if(p[i]+i>mx)
		{
			mx=p[i]+i;
			id=i;
		}
	}
	return ans-1;
}
string dif(const string x)
{
    if(x.length()<=1)
        return "";
    string res="";
    for(int i=1;i<x.length();++i)
        res+=abs(x[i]-x[i-1])+'a';
    return res;
}
int main()
{
	int mm,fg;
	while(cin>>s1)
	{
		int len=s1.length();
		cnt=0;
	//	fg=0;
		mm=0;
		while(1)
		{
			pre();
			cnt=Manacher();
			if(len==cnt||cnt==len-1)
			{
				//fg=1;
				break;
			}
			mm=max(mm,cnt);
			s1=dif(s1);
		    len=s1.length();
		    if(len==1)
		    	break;
		}
		printf("%d\n",max(mm,cnt));
	}
	return 0;
}