Description

題庫鏈接

記一個字符串為 \(X\) ，它的倒置為 \(X^R\) 。現在給你一個長度為 \(n\) 的字符串 \(S\) ，詢問其最長的形同 \(XX^RXX^R\) 的子串長度為多少。

\(1\leq n\leq 500000\)

Solution

顯然求回文之類的東西還是要 \(manacher\) 的。

假設我們求出了 \(manacher\) 的 \(len\) 數組，考慮怎麽去用這個 \(len\) 。

值得肯定的是我們可以用 \(len\) 來刻畫“雙倍回文”。對於雙倍回文的中點 \(x\) ，考慮記其第三段為 \((x,y-1)\) 。

顯然我們可以去枚舉 \(x\)

我們可以用不等式 \(\begin{cases} \begin{aligned} x &\geq y-len_y \\ y &\leq x+\frac{len_x}{2} \end{aligned} \end{cases}\) 來刻畫一對 \((x,y)\) 。

現在就變成了枚舉 \(x\) 然後去找一個滿足上述情況對應的最大的 \(y\) 。

排序+ \(set\) 維護即可。

Code

//It is made by Awson on 2018.3.16
#include <bits/stdc++.h>
#define LL long long
 

#define dob complex<double>
#define Abs(a) ((a) < 0 ? (-(a)) : (a))
#define Max(a, b) ((a) > (b) ? (a) : (b))
#define Min(a, b) ((a) < (b) ? (a) : (b))
#define Swap(a, b) ((a) ^= (b), (b) ^= (a), (a) ^= (b))
#define writeln(x) (write(x), putchar('\n'))
#define lowbit(x) ((x)&(-(x)))
 

using namespace std;
const int N = 500000;
void read(int &x) {
    char ch; bool flag = 0;
    for (ch = getchar(); !isdigit(ch) && ((flag |= (ch == '-')) || 1); ch = getchar());
    for (x = 0; isdigit(ch); x = (x<<1)+(x<<3)+ch-48, ch = getchar());
    x *= 1-2*flag;
}
void print(int x) {if (x > 9) print(x/10); putchar(x%10+48); }
void write(int x) {if (x < 0) putchar('-'); print(Abs(x)); }

int n, loc, len[(N<<1)+5], f[N+5];
char s[N+5], ch[(N<<1)+5];
struct tt {
    int id, x;
    bool operator < (const tt &b) const {return id-x < b.id-b.x; }
}a[N+5];
set<int>S;

void work() {
    read(n); scanf("%s", s+1);
    for (int i = 1; i <= n; i++) ch[++loc] = '#', ch[++loc] = s[i];
    ch[++loc] = '#', ch[++loc] = '$';
    int mx = 0, po = 0;
    for (int i = 1; i <= loc; i++) {
    if (mx > i) len[i] = Min(len[(po<<1)-i], mx-i);
    else len[i] = 1;
    while (ch[i+len[i]] == ch[i-len[i]]) ++len[i];
    if (i+len[i] > mx) mx = i+len[i], po = i;
    }
    for (int i = 1; i <= n; i++) f[i] = a[i].x = ((len[(i<<1)+1]-1)>>1), a[i].id = i;
    sort(a+1, a+n+1); loc = 1; int ans = 0;
    for (int i = 1; i <= n; i++) {
    while (loc <= n && a[loc].id-a[loc].x <= i) S.insert(a[loc].id), ++loc;
    set<int>::iterator it = S.upper_bound(i+(f[i]>>1));
    if (it == S.begin()) continue;
    int t = *(--it)-i; ans = Max(ans, t);
    }
    writeln(ans<<2);
}
int main() {
    work(); return 0;
}

[SHOI 2011]雙倍回文