[国家集训队]最长双回文串

嘟嘟嘟

这道题,我的大致思想是先用manacher求出所有回文串,然后用刚好拼接在一起的两个回文串的总长的最大值更新答案。

manacher是O(n)的,但是暴力的枚举回文中心能达到O(n2),所以得想办法优化枚举。

令pre[i]表示离 i最远且左半部分包含 i 的回文串的回文中心的位置,suf[i]表示离 i 最远且右半部分包含 i 的回文串的回文中心的位置。有人会问,为什么是包含,而不是到 i 结束的回文串?因为最优解对于这个回文串可能取不完,比如说abaab,单纯看aba的话,只能取到ab + a = abaa,但实际上应该是a + baab = abaab。维护好这两个数组后,ans = max(suf[i] - pre[i])。

这两个数组的维护跟manacher有点像:对于pre,建一个指针 j,然后扫维护好的回文半径数组,每一次把[j, i + p[i])的pre都更新成 i。suf的维护反过来即可。

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cmath>
 4 #include<algorithm>
 5 #include<cstring>
 6 #include<cstdlib>
 7 #include<cctype>
 8 #include<vector>
 9 #include<stack>
10 #include<queue>
11 using namespace std;
12 #define enter puts("") 
13 #define space putchar(' ')
14 #define Mem(a, x) memset(a, x, sizeof(a))
15 #define rg register
16 typedef long long ll;
17 typedef double db;
18 const int INF = 0x3f3f3f3f;
19 const db eps = 1e-8;
20 const int maxn = 1e5 + 5;
21 inline ll read()
22 {
23     ll ans = 0;
24     char ch = getchar(), last = ' ';
25     while(!isdigit(ch)) {last = ch; ch = getchar();}
26     while(isdigit(ch)) {ans = ans * 10 + ch - '0'; ch = getchar();}
27     if(last == '-') ans = -ans;
28     return ans;
29 }
30 inline void write(ll x)
31 {
32     if(x < 0) x = -x, putchar('-');
33     if(x >= 10) write(x / 10);
34     putchar(x % 10 + '0');
35 }
36 
37 int n;
38 char s[maxn], t[maxn << 1];
39 int p[maxn << 1];
40 
41 void init()
42 {
43   n = strlen(s);
44   t[0] = '@';
45   for(int i = 0; i < n; ++i) t[i << 1 | 1] = '#', t[(i << 1) + 2] = s[i];
46   n = (n + 1) << 1;
47   t[n - 1] = '#'; t[n] = '$';
48 }
49 int pre[maxn << 2], suf[maxn << 2];
50 void manacher()
51 {
52   int mx = 0, id;
53   for(int i = 1; i < n; ++i)
54     {
55       if(i < mx) p[i] = min(p[(id << 1) - i], mx - i);
56       else p[i] = 1;
57       while(t[i - p[i]] == t[i + p[i]]) p[i]++;
58       if(i + p[i] > mx) mx = i + p[i], id = i; 
59     }
60 }
61 
62 int main()
63 {
64   scanf("%s", s);
65   init(); manacher();
66   for(int i = 1, j = 1; i < n; ++i)
67     for(; j < i + p[i]; ++j) pre[j] = i;
68   for(int i = n - 1, j = n - 1; i; --i)
69     for(; j > i - p[i]; --j) suf[j] = i;
70   int ans = 0;
71   for(int i = 1; i < n; ++i) ans = max(ans, suf[i] - pre[i]);
72   write(ans); enter;
73   return 0;
74 }
View Code
原文地址:https://www.cnblogs.com/mrclr/p/9772084.html