ACM学习感悟——ACdream字符串专场F（manacher）

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

字符串中只包含小写字母。
求导过程如下，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阶导数的导字符串。

abcd
abcba
acdream

3
5
3

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

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

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

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

起始我写这个只是想记得模板罢了，好吧，一道manacher算法的裸题，虾米那直接上代，有一点注意的就是p数组要开的大一点，不然会越界，因为预处理会将字符扩大两倍。

AC代码：

//
Created by
CQUWEL
//
Copyright (c) 2015年 CQUWEL. All rights reserved.
//
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <queue>
#include <stack>
#include <set>
#include <algorithm>
#include <numeric>
#include <functional>
#define INF 0x3f3f3f3f
#define cir(i,a,b)
for (int i=a;i<=b;i++)
#define CIR(j,a,b)
for (int j=a;j>=b;j--)
#define CLR(x) memset(x,0,sizeof(x))
typedef long long
ll;
using namespace std;
string s,tt;
int p[1010*3];
int nn;
string dif(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;
}
void preprocess(string &x){
s[0]='@';s[1]='#';nn=2;
for (int i=0;x[i];i++){
s[nn++]=x[i];
s[nn++]='#';
//
cout << "s"<<endl;
}
s[nn]='';
}
int manacher(string &x){
preprocess(x);
int mxl=0,id=0;
for (int i=1;i<nn;i++){
if (mxl>i) p[i]=min(mxl-i,p[2*id-i]);
else p[i]=1;
while (s[i-p[i]]==s[i+p[i]]) p[i]++;
if (mxl<p[i]+i) mxl=p[i]+i,id=i;
}
int ret=0;
for (int i=1;i<nn;i++) ret=max(ret,p[i]);
return ret-1;
//这里要注意！！！！
}
int main()
{
while (getline(cin,tt)){
int ans=0;
//
cout << dif(tt);
//
preprocess(tt);
//
for (int i=0;s[i];i++) cout << s[i];
int l=tt.size();
if (l==0){
cout << "0n";
continue;
}
//这里只是放置有空串，不过一般题目是不会给空串的。
while (ans<tt.size()){
ans=max(ans,manacher(tt));
//
cout << tt << "||" << ans <<endl;
tt=dif(tt);
if (tt.size()<=ans) break;
}
cout << ans << endl;
}
return 0;
}