算法专题
前言
| 所讲例题 |
|---|
| 洛谷 3805【模板】manacher算法 |
| JZOJ 2682 洛谷 4555 最长双回文串 |
| JZOJ 1950 洛谷 1659 拉拉队排练 |
洛谷 3805【模板】manacher算法
题目
找到一个字符串中的最长回文子串
分析
首先字符串的长度为
，显然只能用
的时间完成这道题目，首先朴素的方法就是
，也就是枚举中间点向外扩展，但是这样的坏处在哪里呢，就是这样扩展很有可能找不到最优解而浪费了时间，那应该怎么做呢?
那就可以引进一种神奇的算法，manacher
那么这是什么神奇的东西呢，首先对于长度奇偶性可以在字符串中间填充分隔符，这样字符串就一定是奇数，为了避免越界，还要在开头加一个分隔符。
那要记录三样东西，首先是
表示中点为
时的最长回文串的半径，然后那么显而易见答案就是
，那问题是答案怎么算，这是一个大问题
再首先，我们要记录两个东西，
表示当前找到的最长回文串的中点，
表示该回文串的右边界。
首先对于每一个
，若
说明还是有希望的，否则就只能像纯模拟一样
当然不管怎么样还是得加上这句话
那
是为什么，显然可知是更新最长回文子串
但问题是
如何优化呢，那就是
的情况了
那么首先设
，也就是对于
，
的对称点，可以通过
得到。
那么分类讨论，首先先配两张图(蓝色所示为j的答案范围，浅绿所示为i的答案范围，区间为mid的答案范围）
- 的答案范围超出 的答案范围,那么如果 像 那样求答案，就违背了 的答案，故假设不成立，所以
- 的答案范围小于 的答案范围,那么如果 像 那样求答案，就违背了 的答案，故假设不成立，所以
- 的答案范围等于 的答案范围，则两者均可
综上所述，就可以得到上面的代码
代码
#include <cstdio>
#include <cctype>
#define rr register
using namespace std;
char s[22000015]; int n,ans,p[22000015];
inline signed max(int a,int b){return (a>b)?a:b;}
inline signed min(int a,int b){return (a<b)?a:b;}
inline signed manacher(char *s,int len){
rr int maxlen=1,mx=0,mid=0;
for (rr int i=1;i<len;++i){
p[i]=mx>i?min(p[(mid<<1)-i],mx-i):1;
while (s[i-p[i]]==s[i+p[i]]) ++p[i];
if (mx<i+p[i]) mx=i+p[i],mid=i;
maxlen=max(maxlen,p[i]-1);
}
return maxlen;
}
signed main(){
rr char c=getchar(); s[0]='!';
while (!isalpha(c)) c=getchar();
while (isalpha(c)) s[++n]='|',s[++n]=c,c=getchar();
s[++n]='|'; ans=manacher(s,n);
return !printf("%d",ans);
}
洛谷 4555 最长双回文串
题目
输入长度为 的串 ，求 的最长双回文子串 ,即可将 分为两部分 ， ， 且 和 都是回文串。
分析
那么在manacher的同时，可以记录最长回文串的最右端的长度 和最左端的长度 ,然后用O(n)的时间更新，也就是 为什么呢，因为首先最长回文串内的回文子串答案是没有更新的，那为什么要更新呢，因为最长双回文串的两个回文串不一定是最长的，那么统计 即可
代码
#include <cstdio>
#include <cctype>
#define rr register
using namespace std;
char s[200015]; int n,ans,p[200015],l[200015],r[200015];
inline signed max(int a,int b){return (a>b)?a:b;}
inline signed min(int a,int b){return (a<b)?a:b;}
inline void manacher(char *s,int len){
rr int mx=0,mid=0;
for (rr int i=1;i<len;++i){
p[i]=mx>i?min(p[(mid<<1)-i],mx-i):1;
while (s[i-p[i]]==s[i+p[i]]) ++p[i];
if (mx<i+p[i]) mx=i+p[i],mid=i;
l[i+p[i]-1]=max(l[i+p[i]-1],p[i]-1);
r[i-p[i]+1]=max(r[i-p[i]+1],p[i]-1);
}
for (rr int i=1;i<=len;i+=2) r[i]=max(r[i],r[i-2]-2);
for (rr int i=len;i>0;i-=2) l[i]=max(l[i],l[i+2]-2);
for (rr int i=1;i<=len;i+=2) if (l[i]&&r[i]) ans=max(ans,l[i]+r[i]);
}
signed main(){
rr char c=getchar(); s[0]='!';
while (!isalpha(c)) c=getchar();
while (isalpha(c)) s[++n]='|',s[++n]=c,c=getchar();
s[++n]='|'; manacher(s,n);
return !printf("%d",ans);
}
洛谷 1659 拉拉队排练
题目
问前 个最长奇回文串长度的乘积
分析
在manacher的基础上开个桶统计长度个数,当然是一个前缀和，因为每个最长回文串其实包含着一些回文子串，然后一定要记得快速幂
代码
#include <cstdio>
#include <cctype>
#define rr register
using namespace std;
const int mod=19930726; typedef long long ll;
char s[2000015]; int n,p[2000015],b[2000015];
inline signed max(int a,int b){return (a>b)?a:b;}
inline signed min(int a,int b){return (a<b)?a:b;}
inline void manacher(char *s,int len){
rr int mx=0,mid=0;
for (rr int i=1;i<len;++i){
p[i]=mx>i?min(p[(mid<<1)-i],mx-i):1;
while (s[i-p[i]]==s[i+p[i]]) ++p[i];
if (mx<i+p[i]) mx=i+p[i],mid=i;
if (!(p[i]&1)) ++b[p[i]-1];
}
}
inline ll ksm(ll x,ll y){
rr ll ans=1;
for (;y;y>>=1,x=(x*x)%mod)
if (y&1) ans=(ans*x)%mod;
return ans;
}
signed main(){
rr ll res=1,k;
scanf("%*d%lld",&k);
rr char c=getchar(); s[0]='!';
while (!isalpha(c)) c=getchar();
while (isalpha(c)) s[++n]='|',s[++n]=c,c=getchar();
s[++n]='|'; manacher(s,n);
for (rr int i=n,sum=0;i>0;i-=2){
sum+=b[i];
if (k>=sum){
res=(res*(ksm(i,sum)))%mod;
k-=sum;
}else{
res=(res*(ksm(i,k)))%mod;
k=0; break;
}
}
if (k>0) res=-1;
return !printf("%lld",res);
}