题目:
BZOJ2565
分析:
首先看到回文串,肯定能想到Manacher算法。下文中字符串 s s 是输入的字符串在Manacher算法中添加了字符‘#’后的字符串 (构造方式如下)
string s = "#";
for (int i = 0; i < str.size(); i++)
{
s += str[i];
s += '#';
}
如果用
maxli
m
a
x
l
i
表示以第
i
i
个字符结尾的最长回文串的长度,表示以第
i
i
个字符开头的最长回文串的长度,那么题目中要求的可以转化为在中找一个位置
i
i
,满足是’#’且
maxli+maxri
m
a
x
l
i
+
m
a
x
r
i
最大。在原串
str
s
t
r
中,它是两个长度分别为
maxli−12
m
a
x
l
i
−
1
2
和
maxri−12
m
a
x
r
i
−
1
2
的回文串 (要减掉额外加进去的’#’字符) 。因此,算出
maxl
m
a
x
l
和
maxr
m
a
x
r
后,就可以枚举所有’#’字符来得到答案。
怎么算
maxl
m
a
x
l
和
maxr
m
a
x
r
呢?对于一个位置
pos
p
o
s
,显然以它结尾的最长回文串的中心是一个最小的
i
i
满足 (
pi
p
i
是Manacher中求出的以
i
i
为中心的回文串的“半径”),此时。那么带着单调队列从左往右扫一遍就能算出
maxr
m
a
x
r
,详见代码。同理,从右往左扫一遍可以算出
maxl
m
a
x
l
代码:
我WA一下午,只因为局部变量没初始化……
#include <iostream>
#include <string>
using namespace std;
namespace zyt
{
const int M = 1e5 * 2 + 10;
int p[M];
void manacher(const string &str)
{
string s = "#";
int id = 0, right = 0;
for (int i = 0; i < str.size(); i++)
{
s += str[i];
s += '#';
}
for (int i = 0; i < s.size(); i++)
{
if (i < right)
p[i] = min(p[id * 2 - i], right - i);
else p[i] = 1;
while (i - p[i] >= 0 && i + p[i] < s.size() && s[i - p[i]] == s[i + p[i]])
p[i]++;
if (i + p[i] > right)
right = i + p[i], id = i;
}
}
inline int abs(const int x)
{
return x >= 0 ? x : -x;
}
void mk_max(int *maxx, const int len, const bool flag)
{
static int q[M];
int h = 0, t = 0;
for (int i = 0; i < len; i++)
{
int pos = flag ? i : len - i - 1;
q[t++] = pos;
while (h < t && abs(pos - q[h]) >= p[q[h]])
h++;
maxx[pos] = abs(pos - q[h]) * 2 + 1;
}
}
void work()
{
string s;
static int maxl[M], maxr[M];
ios::sync_with_stdio(false);
cin >> s;
manacher(s);
mk_max(maxl, s.size() * 2 + 1, true);
mk_max(maxr, s.size() * 2 + 1, false);
int ans = 0;
for (int i = 0; i < s.size() * 2 + 1; i += 2)
if (maxl[i] > 1 && maxr[i] > 1)
ans = max(ans, (maxl[i] - 1) / 2 + (maxr[i] - 1) / 2);
cout << ans << endl;
}
}
int main()
{
zyt::work();
return 0;
}