【GDOI2017模拟9.21】妮厨的愤怒

来源：互联网发布：安卓游戏开发知乎编辑：程序博客网时间：2024/04/29 00:11

题目大意：

给出一个长度为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();}

阅读全文

0 0