原文地址:http://www.cnblogs.com/GXZlegend/p/6802558.html


题目描述

顺序和逆序读起来完全一样的串叫做回文串。比如acbca是回文串,而abc不是(abc的顺序为“abc”,逆序为“cba”,不相同)。
输入长度为n的串S,求S的最长双回文子串T,即可将T分为两部分X,Y,(|X|,|Y|≥1)且X和Y都是回文串。

输入

一行由小写英文字母组成的字符串S。

输出

一行一个整数,表示最长双回文子串的长度。

样例输入

baacaabbacabb

样例输出

12


题解

个人并不是擅长单调数据结构,于是用了Manacher后写了树状数组,稍微慢了点。

先求出以每个字符为中心的最大回文半径。

然后想办法求以某个字符结尾的最大回文半径,显然s[i]=i-min{j}+1,其中j≤i且j+p[j]-1≥i。

可以用树状数组维护j+p[j]-1≥i的最小的j,其中需要反过来减一下,因为树状数组是向下查询。

这样求出以某个字符结尾的最大回文半径,并得到最大回文长度。

之后同理求出以某个字符开头的最大回文长度。

最后一个向后一个向前加起来求一下就好了。

时间复杂度O(nlogn),存在更优秀的O(n)算法,参考CQzhangyu's blog

#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 2000010
using namespace std;
int p[N] , f[N] , s1[N] , s2[N] , n;
char str[N] , tmp[N];
void update1(int x , int a)
{
int i;
for(i = x ; i <= n ; i += i & -i) f[i] = min(f[i] , a);
}
void update2(int x , int a)
{
int i;
for(i = x ; i <= n ; i += i & -i) f[i] = max(f[i] , a);
}
int query1(int x)
{
int i , ans = 0x7f7f7f7f;
for(i = x ; i ; i -= i & -i) ans = min(ans , f[i]);
return ans;
}
int query2(int x)
{
int i , ans = 0;
for(i = x ; i ; i -= i & -i) ans = max(ans , f[i]);
return ans;
}
int main()
{
int i , mx = 0 , last = 0 , ans = 0;
scanf("%s" , str + 1) , n = strlen(str + 1);
tmp[0] = '0';
for(i = 1 ; i <= n ; i ++ ) tmp[i * 2 - 1] = '#' , tmp[i * 2] = str[i];
n = n * 2 + 1 , tmp[n] = '#';
for(i = 1 ; i <= n ; i ++ )
{
if(mx >= i) p[i] = min(p[last * 2 - i] , mx - i + 1);
else p[i] = 1;
while(tmp[i - p[i]] == tmp[i + p[i]]) p[i] ++ ;
if(mx < i + p[i] - 1) mx = i + p[i] - 1 , last = i;
}
memset(f , 0x7f , sizeof(f));
for(i = 1 ; i <= n ; i ++ )
update1(n - (i + p[i] - 1) + 1 , i) , s1[i] = 2 * i - 2 * query1(n - i + 1) + 1;
memset(f , 0 , sizeof(f));
for(i = n ; i >= 1 ; i -- )
update2(i - p[i] + 1 , i) , s2[i] = 2 * query2(i) - 2 * i + 1;
for(i = 2 ; i <= n ; i ++ ) ans = max(ans , (s1[i - 1] + s2[i]) / 2);
printf("%d\n" , ans);
return 0;
}

【bzoj2565】最长双回文串 Manacher+树状数组的更多相关文章

  1. BZOJ2565:最长双回文串(Manacher)

    Description 顺序和逆序读起来完全一样的串叫做回文串.比如acbca是回文串,而abc不是(abc的顺序为“abc”,逆序为“cba”,不相同). 输入长度为n的串S,求S的最长双回文子串T ...

  2. BZOJ2565最长双回文串——manacher

    题目描述 顺序和逆序读起来完全一样的串叫做回文串.比如acbca是回文串,而abc不是(abc的顺序为“abc”,逆序为“cba”,不相同).输入长度为n的串S,求S的最长双回文子串T,即可将T分为两 ...

  3. BZOJ2565 最长双回文串 【Manacher】

    BZOJ2565 最长双回文串 Description 顺序和逆序读起来完全一样的串叫做回文串.比如acbca是回文串,而abc不是(abc的顺序为"abc",逆序为"c ...

  4. 2021.12.10 P5041 [HAOI2009]求回文串(树状数组求逆序对)

    2021.12.10 P5041 [HAOI2009]求回文串(树状数组求逆序对) https://www.luogu.com.cn/problem/P5041 题意: 给一个字符串 \(S\) ,每 ...

  5. 【BZOJ2565】最长双回文串 Manacher

    [BZOJ2565]最长双回文串 Description 顺序和逆序读起来完全一样的串叫做回文串.比如acbca是回文串,而abc不是(abc的顺序为“abc”,逆序为“cba”,不相同).输入长度为 ...

  6. BZOJ 2565: 最长双回文串 [Manacher]

    2565: 最长双回文串 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1842  Solved: 935[Submit][Status][Discu ...

  7. BZOJ.2565.[国家集训队]最长双回文串(Manacher/回文树)

    BZOJ 洛谷 求给定串的最长双回文串. \(n\leq10^5\). Manacher: 记\(R_i\)表示以\(i\)位置为结尾的最长回文串长度,\(L_i\)表示以\(i\)开头的最长回文串长 ...

  8. bzoj 2565: 最长双回文串 manacher算法

    2565: 最长双回文串 Time Limit: 20 Sec  Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeOnline/problem. ...

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

    ---题面--- 题解: 首先有一个直观的想法,如果我们可以求出对于位置i的最长后缀回文串和最长前缀回文串,那么我们枚举分界点然后合并前缀和后缀不就可以得到答案了么? 所以我们的目标就是求出这两个数列 ...

随机推荐

  1. VMware运行时“内部错误”的解决方法

    解决方法:打开虚拟机实体目录,如下:发现有两个虚拟机配置文件,一个文件大小为4KB,另一个为空.现在虚拟机默认使用为空的配置文件了. 将大小为空的虚拟机配置文件删除掉,然后将另一个配置文件重名命. 接 ...

  2. ES6笔记04-class的基本语法

    JavaScript 语言中,生成实例对象的传统方法是通过构造函数. ES6 提供了更接近传统语言的写法,引入了 Class(类)这个概念,作为对象的模板.通过class关键字,可以定义类. clas ...

  3. 多通道CNN

    在读Convolutional Neural Networks for Sentence Classification 这个文章的时候,它在论文中提出一种模型变种就是 CNN-multichannel ...

  4. 26-dotnet watch run 和attach到进程调试

    1-打开vscode, 按下Ctrl+`,打开命令行窗口 创建一个donet core mvc项目 2-打开刚刚创建的文件夹 3-输入 dotnet run 访问网站 4 -F5键即可调试 5-更改代 ...

  5. 使wlr写cnblog的博客-2 设置cnblog帐号

    ref:http://www.cnblogs.com/liuxianan/archive/2013/04/13/3018732.html   使用: 打开Windows Live Writer,第一次 ...

  6. 十三、MySQL之IDE工具介绍及数据备份

    一.IDE工具介绍 生产环境还是推荐使用mysql命令行,但为了方便我们测试,可以使用IDE工具 下载链接:https://pan.baidu.com/s/1bpo5mqj 二.MySQL数据备份 # ...

  7. 2-安装linux7

    1.操作系统简介 操作系统: 桌面操作系统 redhat fedora slackware ubuntu debian suse linux centos 服务器操作系统 linux redhat s ...

  8. 阿里巴巴Java开发规约Eclipse插件安装及使用

    技术交流群:233513714 插件安装 环境:JDK1.8,Eclipse4+.有同学遇到过这样的情况,安装插件重启后,发现没有对应的菜单项,从日志上也看不到相关的异常信息,最后把JDK从1.6升级 ...

  9. laravel5.5任务调度

    目录 1. 定义调度 1.1 使用Closure 1.2 Artisan 命令调度 1.3 队列任务调度 1.4 Shell 命令调度 1.5 调度频率设置 1.6 闭包测试限制 1.7 避免任务重复 ...

  10. Android保持屏幕常亮唤醒状态

    第一步:  首先添加权限: <uses-permission android:name="android.permission.WAKE_LOCK"></uses ...