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


    嘟嘟嘟

    这道题,我的大致思想是先用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
  • 相关阅读:
    IDE有毒
    Netbeans 8.2关于PHP的新特性
    什么是人格
    谁该赋予一款产品灵魂?
    自从升级到macOS后,整个人都不好了
    公司不是大家庭
    性能各个指标分析
    Sqlserver2012 alwayson部署攻略
    初探Backbone
    SQL Server AlwaysOn架构及原理
  • 原文地址:https://www.cnblogs.com/mrclr/p/9772084.html
Copyright © 2020-2023  润新知