【GDOI2017模拟9.21】妮厨的愤怒

题目大意：

给出一个长度为n的字符串S，Q个询问，每次给出l,r，求S的第l位到第r位的最长回文子串。

(1<=n,Q<=10^5)

题解：

记得以前做过一道很像的题。

是YJL的题，二维RMQ，对于每一个查询也是二分答案。

不说了，首先添加特殊字符，manacher求出r。

对于每一个询问，二分答案，于是我们得到了可能的回文中心的范围，RMQ查询区间最大值，判一下即可。

Code：

#include<cmath>
#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))
using namespace std;

const int N = 100005, M = 17;

char s[N * 2]; int r[N * 2];
int n, Q, x, y, f[M + 1][N * 2], a2[M + 1];

void Init() {
a2[0] = 1; fo(i, 1, M) a2[i] = a2[i - 1] * 2;
scanf("%d", &n); scanf("%s", s + 1);
fd(i, n, 1) s[i * 2] = s[i];
fo(i, 1, n) s[i * 2 - 1] = '@'; s[2 * n + 1] = '@';
n = 2 * n + 1;
}

void Build() {
int mx = 0, p = 0;
fo(i, 1, n) {
r[i] = mx >= i ? min(mx - i + 1, r[2 * p - i]) : 1;
while(i != r[i] && s[i - r[i]] == s[i + r[i]]) r[i] ++;
if(i + r[i] - 1 > mx) mx = i + r[i] - 1, p = i;
}
}

void Build_b() {
fo(i, 1, n) f[0][i] = r[i];
fo(i, 1, M) fo(j, 1, n) f[i][j] = max(f[i - 1][j], f[i - 1][min(j + a2[i - 1], n)]);
}

int Max(int x, int y) {
int l = log2(y - x + 1);
return max(f[l][x], f[l][y - a2[l] + 1]);
}

void End() {
for(scanf("%d", &Q); Q; Q --)    {
scanf("%d %d", &x, &y); x ++; y ++;
x = x * 2 - 1; y = y * 2 + 1;
int ans = 0;
for(int l = 1, r = (y - x) / 2 + 1; l <= r; ) {
int m = (l + r) / 2;
if(Max(x + m - 1, y - m + 1) >= m)
ans = m, l = m + 1; else r = m - 1;
}
printf("%d\n", ans - 1);
}
}

int main() {
Init();
Build();
Build_b();
End();
}