马拉车manacher
目的:线性查找一个串的最长回文子串
时间复杂度:O(n)
len[i]表示以i为中心的回文串的半径,mx即为当前计算回文串最右边字符的最大值,p是中心点mid,mx-i和2*p-1关于p对称
https://blog.csdn.net/csdn_kou/article/details/82917937
hdu3068,板子题,求最长回文长度。
#include<bits/stdc++.h>
using namespace std;
const int maxn=;
int t,len[maxn*];
char S[maxn*],T[maxn*],s[maxn*]; int init(char *str)
{
int n=strlen(str);
for(int i=,j=;i<=*n;j++,i+=)
{
s[i]='#';
s[i+]=str[j];
}
s[]='$';
s[*n+]='#';
s[*n+]='@';
s[*n+]='\n';
return *n+;
}
void manacher(int n)
{
int mx=,p=;
for(int i=;i<=n;i++)
{
if(mx>i)len[i]=min(mx-i,len[*p-i]);
else len[i]=;
while(s[i-len[i]]==s[i+len[i]])len[i]++;
if(len[i]+i>mx)mx=len[i]+i,p=i;
}
} int main()
{
while(scanf("%s",S)!=EOF)
{
int Len=strlen(S),n=init(S);
for(int i=;i<=n;i++)len[i]=;
manacher(n); int ans=;
for(int i=;i<=n;i++)
{
ans=max(ans,len[i]-);
}
printf("%d\n",ans);
}
return ;
}
2019徐州G colorful string,求所有回文子串的value之和,一个串的value为串中字母种类,dfs预处理了第i位前一个a-z的位置,复杂度26*n。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=3e5+;
int t,len[maxn*];
char S[maxn*],T[maxn*],s[maxn*]; int init(char *str){
int n=strlen(str);
for(int i=,j=;i<=*n;j++,i+=){
s[i]='#';
s[i+]=str[j];
}
s[]='$';
s[*n+]='#';
s[*n+]='@';
s[*n+]='\n';
return *n+;
} void manacher(int n)
{
int mx=,p=;
for(int i=;i<=n;i++){
if(mx>i) len[i]=min(mx-i,len[*p-i]);
else len[i]=;
while(s[i-len[i]]==s[i+len[i]]) len[i]++;
if(len[i]+i>mx) mx=len[i]+i,p=i;
}
} int dp[maxn*][],place[];
int main()
{
scanf("%s",S);
int Len=strlen(S),n=init(S);
for(int i=;i<=n;i++)len[i]=;
manacher(n); ll ans=;
int k=;
for(int i=;i<;i++)place[i]=-;
for(int i=;i<=n;i++)
{
if(i%==)place[S[k++]-'a']=i;
for(int j=;j<;j++)dp[i][j]=place[j];
} for(int i=;i<=n;i++)
{
for(int j=;j<;j++)
{
if(i-dp[i][j]<len[i])
{
ans+=1ll*(len[i]-(i-dp[i][j]))/;
}
}
}
printf("%lld\n",ans); return ;
}
hdu3613,一个串割成两个串,如果是回文串则val为所有字母val之和,否则为零。
字母的val题目给出,求使总串的val最高的割法的val值。
#include<bits/stdc++.h>
using namespace std;
const int maxn=5e5+;
int t,len[maxn<<];
char S[maxn<<],T[maxn<<],s[maxn<<]; int init(char *str)
{
int n=strlen(str);
for(int i=,j=;i<=*n;j++,i+=)
{
s[i]='#';
s[i+]=str[j];
}
s[]='$';
s[*n+]='#';
s[*n+]='@';
s[*n+]='\n';
return *n+;
}
void manacher(int n)
{
int mx=,p=;
for(int i=;i<=n;i++)
{
if(mx>i)len[i]=min(mx-i,len[*p-i]);
else len[i]=;
while(s[i-len[i]]==s[i+len[i]])len[i]++;
if(len[i]+i>mx)mx=len[i]+i,p=i;
}
} int main()
{
int T;
scanf("%d",&T);
while(T--)
{
int val[]={};
for(int i=;i<;i++)scanf("%d",&val[i]);
cin>>S;
int Len=strlen(S),n=init(S);
for(int i=;i<=n;i++)len[i]=;
manacher(n); int lg=,rg=,ans=;
for(int i=;i<Len;i++)rg+=val[S[i]-'a'];
for(int i=;i<=n;i++)
{
if(i%==)
{
lg+=val[s[i]-'a'];
rg-=val[s[i]-'a'];
}
else
{
int l=,r=*Len+,tmplg=,tmprg=;
if(i!=l&&i!=r)
{
int mid1=(l+i)>>,mid2=(i+r)>>;
if(len[mid1]==mid1-l+)tmplg=lg;
if(len[mid2]==r-mid2+)tmprg=rg;
ans=max(ans,tmplg+tmprg);
}
}
}
printf("%d\n",ans);
}
return ;
}
hdu3294,求最早出现的最长回文串的l、r区间,并且转换后输出
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+;
int t,len[maxn<<];
char S[maxn<<],T[maxn<<],s[maxn<<];
int init(char *str)
{
int n=strlen(str);
for(int i=,j=;i<=*n;j++,i+=)
{
s[i]='#';
s[i+]=str[j];
}
s[]='$';
s[*n+]='#';
s[*n+]='@';
s[*n+]='\n';
return *n+;
}
void manacher(int n)
{
int mx=,p=;
for(int i=;i<=n;i++)
{
if(mx>i)len[i]=min(mx-i,len[*p-i]);
else len[i]=;
while(s[i-len[i]]==s[i+len[i]])len[i]++;
if(len[i]+i>mx)mx=len[i]+i,p=i;
}
} int main()
{
char ch;
while(scanf("%c %s",&ch,S)!=EOF)
{
int Len=strlen(S),n=init(S);
for(int i=;i<=n;i++)len[i]=;
manacher(n); int ans=,l=,r=;
for(int i=;i<=n;i++)
{
if(len[i]->ans)
{
ans=len[i]-;
l=(i-len[i]+)/;
r=(i+len[i]-)/;
}
}
if(ans<)printf("No solution!\n");
else
{
printf("%d %d\n",l,r);
for(int i=l;i<=r;i++)
{
S[i]=S[i]-ch+'a';
if(S[i]>'z')S[i]=S[i]-'z'+'a'-;
if(S[i]<'a')S[i]=S[i]+'z'-'a'+;
}
S[r+]='\0';
printf("%s\n",S+l);
}
getchar();
}
return ;
}
hdu4513,最长递增(不降)回文串,manacher里只需要改个while
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+;
int t,len[maxn*];
int s[maxn*],T[maxn*],S[maxn*]; int init(int *str,int n)
{
for(int i=,j=;i<=*n;j++,i+=)
{
s[i]=-;
s[i+]=str[j];
}
s[]=-;
s[*n+]=-;
s[*n+]=-;
return *n+;
}
void manacher(int n)
{
int mx=,p=;
for(int i=;i<=n;i++)
{
if(mx>i)len[i]=min(mx-i,len[*p-i]);
else len[i]=;
while(s[i-len[i]]==s[i+len[i]]&&s[i-len[i]]<=s[i-len[i]+])len[i]++;
if(len[i]+i>mx)mx=len[i]+i,p=i;
}
} int main()
{
int T;
scanf("%d",&T);
while(T--)
{
int n;
scanf("%d",&n);
for(int i=;i<n;i++)
scanf("%d",&S[i]);
int Len=n,len2=init(S,n);
for(int i=;i<=n;i++)len[i]=;
manacher(len2); int ans=;
for(int i=;i<=len2;i++)
{
ans=max(ans,len[i]-);
}
printf("%d\n",ans);
}
return ;
}
...
马拉车manacher的更多相关文章
- 51nod 1595 回文度 | 马拉车Manacher DP
51nod 1595 回文度 题目描述 如果长度为n的字符串是一个回文串,同时它的长度为floor(n/2)的前缀和后缀是K-1度回文串,则这个字符串被称为K度回文串.根据定义,任何字符串(即使是空字 ...
- HDU - 3068 最长回文(马拉车Manacher)题解
思路:马拉车裸题,我们用一个p[i]数组代表以i为中心的最大回文半径.这里用了一个小技巧,如果一个串是aaaa这样的,那我们插入不相干的字符使它成为#a#a#a#a#,这样无论这个串是奇数还是偶数都会 ...
- hdu3068-最长回文-马拉车(Manacher)算法
http://acm.hdu.edu.cn/showproblem.php?pid=3068 脑子转个弯总算看懂马拉车算法了.记录一下思路和模板. 马拉车算法是在O(n)的时间内求出最大回文子串. 一 ...
- 马拉车——Manacher一篇看上去很靠谱的理解(代码显然易懂)
由于回文分为偶回文(比如 bccb)和奇回文(比如 bcacb),而在处理奇偶问题上会比较繁琐,所以这里我们使用一个技巧,在字符间插入一个字符(前提这个字符未出现在串里).举个例子:s="a ...
- 一点总结-关于debug比赛
上午的题目是: 1. main里面定义的变量必须手动初始化,使用memset或者其他,函数外或者函数内,会进行初始化为0. 2. 最长回文子串的马拉车manacher算法,不会写! 3. 数字三角形d ...
- ACM-ICPC 2018 南京网络赛
题目顺序:A C E G I J L A. An Olympian Math Problem 打表,找规律,发现答案为n-1 C. GDY 题意: m张卡片,标号1-13: n个玩家,标号1-n:每个 ...
- ACM模板_axiomofchoice
目录 语法 c++ java 动态规划 多重背包 最长不下降子序列 计算几何 向量(结构体) 平面集合基本操作 二维凸包 旋转卡壳 最大空矩形 | 扫描法 平面最近点对 | 分治 最小圆覆盖 | 随机 ...
- [菜b]Isaunoya 的一些学习笔记…[保持咕咕咕]
fread/fwrite标记永久化 分块 树链剖分 莫比乌斯反演 斜率优化/单调队列 kruskal重构树 回滚莫队 可持久化线段树/trie树 Link-Cut-Tree dsu on tree F ...
- Manacher's Algorithm 马拉车算法
这个马拉车算法Manacher‘s Algorithm是用来查找一个字符串的最长回文子串的线性方法,由一个叫Manacher的人在1975年发明的,这个方法的最大贡献是在于将时间复杂度提升到了线性,这 ...
随机推荐
- ThreadLocal<T> 源码解析
在activeJDBC框架内部的实现中看到了 ThreadLocal 这个类,记录下了每个线程独有的连接 private static final ThreadLocal<HashMap< ...
- [LC]643题 Maximum Average Subarray I(子数组最大平均数 I)
①英文题目 Given an array consisting of n integers, find the contiguous subarray of given length k that h ...
- hdu 2255 奔小康赚大钱 (KM)
奔小康赚大钱Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submi ...
- setState何时同步,何时异步,为什么?
setState何时同步,何时异步,为什么 答案:在React库控制时,异步:否则同步. 示例代码如下: constructor(props){ super(porps); this.state = ...
- NIO流的学习以及Buffer的相关操作
NIO的使用 一).什么叫NIO? 定义:是一套新的Java I/O标准, 在java1.4中被纳入JDK中. 二).NIO的实现方法 NIO是基于块的, 以块为基本单位处理数据. 标准的I/O是基于 ...
- vue监听移动设备的返回事件
在公共方法文件common.js中实现一个存储当前历史记录的方法 common.js // 存储当前历史记录点,实现控制手机物理返回键的按钮事件 var pushHistory = function ...
- vim的查找功能
vim是一款强大的编辑器. 在vim下要查找字符串: 一,全匹配: 1,从上往下查找,比如“string” : /string 2,从下往上查找,比如“string” : ?string 二,模糊 ...
- 如何通过swoole加速laravel的问题?
这篇文章主要介绍了关于如何使用swoole加速laravel,有着一定的参考价值,现在分享给大家,有需要的朋友可以参考一下 再来复习一下吧,导致 php 慢的各种因素中解析性语言的特性可以说是罪魁祸首 ...
- 预训练语言模型整理(ELMo/GPT/BERT...)
目录 简介 预训练任务简介 自回归语言模型 自编码语言模型 预训练模型的简介与对比 ELMo 细节 ELMo的下游使用 GPT/GPT2 GPT 细节 微调 GPT2 优缺点 BERT BERT的预训 ...
- Nginx热部署 平滑升级 日志切割
1.重载 修改nginx配置文件之后,在不影响服务的前提下想加载最新的配置,就可以重载配置即可. 操作如下: 1)修改nginx配置文件 2)nginx -t 检查nginx文件语法是否有误 ...