HYSBZ 2342 最长双倍回文串 Manacher

本文介绍了一种使用Manacher算法寻找字符串中最大双倍回文子串的方法。通过预处理得到每个位置上的最长回文半径,再通过枚举中心的方式找到满足条件的最大双倍回文串。

题目链接

题意: 双倍回文串 为:一个回文串是一个更小的回文串的两倍

首先我们用Manacher将字符串扫一遍 得到dis数组
由题意这个回文串一定是4的倍数

枚举这个串的中心i 和 小回文串的半径k dis[i] - 1 >= 4 * k
判断是否存在这样的小回文串
如果存在 则小回文串的中心 j = i + (4 * k + 1) / 2
如果dis[j] - 1 >= 2*k 那么存在这样的双倍回文串 长度为4 * k

因为要求最大的长度 所以我们可以在枚举完i之后 从最大的k ((dis[i] - 1) / 4) 开始向下枚举 找到一个K 就更新ans 推出循环

#include <cstdio>
#include <iostream>
#include <vector>
#include <cstring>
#include <algorithm>
#define sf scanf
#define pf printf
using namespace std;
const int maxn = 500000  * 2 + 50;
int dis[maxn];
char str1[maxn],str2[maxn];

int get_dis(){
    int len = strlen(str1);
    str2[0] = '$';
    char* str_a = str2 + 1;
    for(int i = 0;i <= len;++i){
        str_a[i * 2] = '#';
        str_a[i * 2 + 1] = str1[i];
    }
    int id = 0, mx = 1,len2 = strlen(str2);
    int ans = 0;
    for(int i = 1;i < len2;++i){
        if(mx > i){
            dis[i] = (dis[id * 2 - i] < (mx - i) ? dis[2 * id - i] : (mx - i));
        }
        else dis[i] = 1;
        while(str2[i - dis[i]] == str2[i + dis[i]]) dis[i]++;
        if(i + dis[i] > mx){
        mx = i + dis[i];
            id = i;
        }
    }
    return len2;
}
int main(){
    int ca = 0,N,ans;
    while(~sf("%d %s",&N,str1) ){
        N = get_dis();
        ans = 0;
        for(int i = 1;i < N;i += 2){
            for(int k = (dis[i] - 1) / 4;k;--k){
                int n_dis = 4 * k + 1;
                int j = n_dis / 2 + i;
                if(dis[j] - 1 >= 2*k){
                    ans = max(ans,k * 4);
                    break;
                }
            }
        }
        pf("%d\n",ans);
    }
    return 0;
}


# P4287 [SHOI2011] 双倍回文 ## 题目描述 记字符串 $w$ 的倒置为 $w^{\mathsf R}$。例如$\tt (abcd)^{\mathsf R}=dcba$,$\tt (abba)^{\mathsf R}=abba$。 对字符串 $x$,如果 $x$ 满足 $x^{\mathsf R}=x$,则称之为回文。例如 $\tt abba$ 是一个回文,而 $\tt abed$ 不是。 如果 $x$ 能够写成 $ww^{\mathsf R} ww^{\mathsf R}$ 形式,则称它是一个“双倍回文”。换句话说,若要 $x$ 是双倍回文,它的长度必须是 $4$ 的倍数,而且 $x$,$x$ 的前半部分,$x$ 的后半部分都要是回文。例如 $\tt abbaabba$ 是一个双倍回文,而 $\tt abaaba$ 不是,因为它的长度不是 $4$ 的倍数。 - $x$ 的子串是指在 $x$ 中连续的一段字符所组成的字符串。例如 $\tt be$ 是 $\tt abed$ 的子串,而 $\tt ac$ 不是。 - $x$ 的回文子串,就是指满足回文性质的 $x$ 的子串。 - $x$ 的双倍回文子串,就是指满足双倍回文性质的 $x$ 的子串。 你的任务是,对于给定的字符串,计算它的最长双倍回文子串的长度。 ## 输入格式 输入分为两行。 第一行为一个整数,表示字符串的长度。 第二行有个连续的小写的英文字符,表示字符串的内容。 ## 输出格式 输出文件只有一行即输入数据中字符串的最长双倍回文子串的长度,如果双倍回文子串不存在,则输出 $0$。 ## 输入输出样例 #1 ### 输入 #1 ``` 16 ggabaabaabaaball ``` ### 输出 #1 ``` 12 ``` ## 说明/提示 ### 数据范围及约定 对于全部数据,$1\le N \le 500000$。 为什么WA了 ```cpp #include <iostream> using namespace std; const int MAX_N = 500050; int n, d[MAX_N * 2], ans; char s[MAX_N], t[MAX_N * 2] = "?"; int main() { cin >> n >> s + 1; for (int i = 1; i <= n; i++) t[2 * i] = s[i], t[2 * i + 1] = '#'; for (int i = 1, l, r = 0; i <= 2 * n + 1; i++) { if (i < r) d[i] = min(d[r - i + l], r - i + 1); for (; t[i - d[i]] == t[i + d[i]]; d[i]++); if (i + d[i] - 1 > r) { for (int j = max(r + 1, i + 4); j <= i + d[i] - 1; j++) if ((j - i) % 4 == 0 && d[i - (j - i) / 2] > (j - i) / 2) ans = max(ans, j - i); l = i - d[i] + 1, r = i + d[i] - 1; } } cout << ans << '\n'; return 0; } ```
最新发布
08-09
评论 1
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值