1. 程式人生 > 其它 >[學習筆記]Manacher

[學習筆記]Manacher

實現

為解決迴文串長度奇偶性之間的差距,在相鄰兩個字元之間和串頭串尾插入\('\#'\);
為避免越界,在字串前加\('\$'\),串尾有\('\setminus0'\).
\(r[i]\)表示以第\(i\)個字元為對稱軸的迴文串的迴文半徑,顯然現在迴文串半徑為\(l\)的迴文串在原串的迴文串長度為\(l-1\).
\(mx\)表示在以\(1\sim\;i-1\)為中心的迴文串中,串尾下標的最大值\(+1\);\(id\)表示\(mx\)值所對應的迴文串中心.

核心程式碼:

if(mx>i) r[i]=min(r[(id<<1)-i],mx-i);

\(id\times2-i\)

\(i\)關於\(id\)對稱,\(a[id-r[id]+1...id]=a[id...id+r[id]-1]\),所以兩者之間取最小值(超出部分匹配不一定成功).
時間複雜度:\(O(n)\).

#define N 200005
using namespace std;
int r[N],m,n,mx,id,ans;
char a[N];
inline int manacher(){
	for(int i=n;i;--i){
		a[i<<1]=a[i];
		a[i<<1|1]='#';
	}
	n=n<<1|1;mx=id=0;
	a[0]='$';a[1]=a[n+1]='#';
	for(int i=1;i<=n;++i){
		r[i]=i<mx?min(r[(id<<1)-i],mx-i):1;
		while(a[i+r[i]]==a[i-r[i]]) ++r[i];
		if(i+r[i]>mx) mx=i+r[i],id=i;
		ans=max(ans,r[i]-1);
	}
	return ans;
}

推薦

http://blog.csdn.net/ggggiqnypgjg/article/details/6645824/