acm-马拉车算法学习笔记

引言

本文主要介绍马拉车算法的原理及相关应用，并给出经典习题

原理

马拉车(manacher)算法主要用于解决回文串问题，它能够$O(n)$求出一个字符串每个位置的$h$数组值，$h[i]$代表的是从字符串的第$i$个位置开始向右延伸能够到达的最大长度使得$s[i-h[i]+1\sim i+h[i]-1]$是一个回文串。

首先对字符串预处理，即再每个字符间加上一个不在字符集中的字符(比如#)。

考虑维护一个变量$now$，改变量记录当前使得$i+h[i]$最大的$i$的值，也就是说$now+h[now]$是已求解的位置中值最大的。那么考虑$i$关于$now$对称的位置$j$，满足 h [ i ] ≥ m i n { h [ j ] , n o w + h [ n o w ] − i } h[i]\ge min\{h[j],now+h[now]-i\} h[i]≥min{h[j],now+h[now]−i}，这个很容易理解，不理解的话看这里。然后每次更新完$h[i]$的时候再更新一下$now$即可，对于每个位置而言的回文串长度其实就是$h[i]-1$。

代码

这里以LuoGu P3805 【模板】manacher算法为例给出一份参考代码。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 2e7;
char s[maxn],ss[maxn<<1];
int h[maxn<<1];
int manacher(char *s,int len){
	int now=0,ans=0;
	h[0]=1;
	for(int i=1;i<=len;++i){
		h[i]=min(h[(now<<1)-i],now+h[now]-i);
		while(i+h[i]<=len && i-h[i]>=1 && s[i+h[i]]==s[i-h[i]])h[i]++;
		if(i+h[i]>now+h[now])now=i;
		ans=max(ans,h[i]-1);
	}
	return ans;
}
int del(char *src,char *tar,int len){
	int new_len=0;
	for(int i=1;i<=len;++i){
		tar[++new_len]='#';
		tar[++new_len]=src[i];
	}
	tar[++new_len]='#';
	return new_len;
}

int main(){
	scanf("%s",s+1);
	int len=del(s,ss,strlen(s+1));
	printf("%d\n",manacher(ss,len));
}

习题

例题一
题目来源：LuoGuP4555 [国家集训队]最长双回文串

题面：

题解：首先用马拉车处理出$h[]$数组，然后考虑枚举双回文串的分界点所在的位置，这个位置恰好就是一个#符号，因此枚举每个#所在的位置，分别向左边延伸和向右边延伸看能够延伸多长的回文串。

不过向左延伸和向右边延伸回文串的最大长度可以提前预处理，利用h数组可以处理出$l[i],r[i]$两个数组，分别代表从$i$点向左边延伸能够达到的最远点$l[i]$满足$s[l[i]\sim i]$是一个回文串，以及从$i$点向右延伸能够到达的最远点$r[i]$满足$s[i\sim r[i]]$是一个回文串。

如何处理$l[i],h[i]$呢？以$l[i]$为例，我们维护一个队列，储存回文串中心点，并且要求该回文串右端点不小于当前枚举的点$i$，按照中心点位置从小到大摆放，注意当某个中心点对应的回文串的右端点小于$i$时要及时pop出队列。这样每次就取队头来更新$l[i]$即可，对应的是队头的回文串的左端点。
代码：

#include <bits/stdc++.h>
using namespace std;

const int maxn = 2e5+10;
char s[maxn],ss[maxn];
int h[maxn],l[maxn],r[maxn];

void manacher(char *s,int len){
	int now=0;
	h[0]=1;
	for(int i=1;i<=len;++i){
		h[i]=min(now+h[now]-i,h[(now<<1)-i]);
		while(i+h[i]<=len && i-h[i]>=1 && s[i+h[i]]==s[i-h[i]])h[i]++;
		if(now+h[now]<i+h[i])now=i;
	}
}
int deal(char *src,char *tar,int len){
	int ans=0;
	for(int i=1;i<=len;++i){
		tar[++ans]='#';
		tar[++ans]=src[i];
	}
	tar[++ans]='#';
	return ans;
}
queue<int>q; 
int main(){
	scanf("%s",s+1);
	int len=deal(s,ss,strlen(s+1));
	manacher(ss,len);
	ss[len+1]='\0';
	for(int i=1;i<=len;++i){
		q.push(i);
		while(q.front()+h[q.front()]-1<i)q.pop();
		l[i]=(q.front()<<1)-i;
	}
	while(!q.empty())q.pop();
	for(int i=len;i>=1;--i){
		q.push(i);
		while(q.front()-h[q.front()]+1>i)q.pop();
		r[i]=(q.front()<<1)-i;
	}
	int ans=0;
	for(int i=1;i<=len;++i){
		if(ss[i]=='#' && (i-l[i]) && (r[i]-i)){
			ans=max(ans,((i-l[i])>>1)+((r[i]-i)>>1));
		}
	}
	printf("%d\n",ans);
} 

例题二
题目来源：SPOJ7586 Number of Palindromes

题面：


题解：本题让求解所有回文子串个数，板题就不多解释了。
代码：

#include <bits/stdc++.h>
using namespace std;

const int maxn = 10000;
typedef long long ll;

char s[maxn],ss[maxn];
int h[maxn];
int init(char *src,char *tar,int len){
	int ans=0;
	for(int i=1;i<=len;++i){
		tar[++ans]='#';
		tar[++ans]=src[i];
	}
	tar[++ans]='#';
	return ans;
}

void manacher(char *s,int len){
	int now=0;h[0]=1;
	for(int i=1;i<=len;++i){
		h[i]=min(now?h[(now<<1)-i]:0,now+h[now]-i);
		while(i+h[i]<=len && i-h[i]>=1 && s[i-h[i]]==s[i+h[i]])h[i]++;
		if(i+h[i]>now+h[now])now=i;
	}
} 

int main(){
	scanf("%s",s+1);
	int len=strlen(s+1);
	len=init(s,ss,len);
	manacher(ss,len);
	ss[len+1]='\0';
	ll ans=0;
	for(int i=1;i<=len;++i){
		ans+=h[i]/2;
	}
	printf("%lld\n",ans);
}