Description:
顺序和逆序读起来完全一样的串叫做回文串。比如acbca是回文串,而abc不是(abc的顺序为“abc”,逆序为“cba”,不相同)。
输入长度为n的串S,求S的最长双回文子串T,即可将T分为两部分X,Y,(|X|,|Y|≥1)且X和Y都是回文串。
对于100%的数据,
2≤|S|≤105
。
题解:
manacher的裸题。
先用马拉车求出r数组。
显然可以先枚举X部分的回文中心i,对于每一个i,目标就是找到最大的j,满足(j>i+r[i]-1)and(j - r[j]<=i+r[i]-1)。
倒着枚举i,用个树状数组维护一下即可。
Code:
#include<cstdio>
#include<cstring>
#define fo(i, x, y) for(int i = x; i <= y; i ++)
#define fd(i, x, y) for(int i = x; i >= y; i --)
#define min(a, b) ((a) < (b) ? (a) : (b))
#define max(a, b) ((a) > (b) ? (a) : (b))
#define low(x) ((x) & -(x))
using namespace std;
const int N = 100005 * 2;
int n, r[N], f[N];
char s[N];
void Build() {
int p = 0, mx = 0;
fo(i, 1, n) {
r[i] = mx > i ? min(r[p * 2 - i], mx - i + 1) : 1;
while (i > r[i] && s[i - r[i]] == s[i + r[i]]) r[i] ++;
if(i + r[i] - 1 > mx) p = i, mx = i + r[i] - 1;
}
}
void add(int x, int y) {
while(x <= n) {
f[x] = max(f[x], y);
x += low(x);
}
}
int find(int x) {
int s = 0;
while(x) {
s = max(s, f[x]);
x -= low(x);
}
return s;
}
int main() {
freopen("a.in", "r", stdin);
scanf("%s", s + 1);
n = strlen(s + 1);
fd(i, n, 1) s[i * 2] = s[i];
fo(i, 1, n) s[i * 2 - 1] = '@';
n = n * 2 + 1; s[n] = '@';
Build();
int ans = 0;
fd(i, n, 1) {
int s = i + r[i] > n ? find(i + r[i] - 1) : find(i + r[i]);
if(s > 0) ans = max(ans, s - i);
add(i - r[i] + 1, i);
}
printf("%d", ans);
}