Manacher回文串(板子总结)

板子1（以hdu3068为例）

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=110005;
char w[maxn],t[maxn<<1];
int R[maxn<<1],tot;
//R[i]: 以i为中心的包含i和井号的最长回文半径 减1为回文串长 
int init(char w[],char t[])//原串:w 新串:t  
{
	tot=0;
	memset(t,0,sizeof t);
	memset(R,0,sizeof R);
	int len=strlen(w);
	t[0]='@';
	for(int i=0;i<len;++i)
	{
		t[++tot]='#';
		t[++tot]=w[i];
	}
	t[++tot]='#';
	t[++tot]='*';
	return 2*len+1;
	//新串长度 为#3#2#3#的长度 不包括@和*
	//新串下标为 [1,2*len+1] 
}
int manacher(int len)
{
	int r=0,ans=0,pos=0;
	for(int i=1;i<=len;++i)
	{
		if(r>i)R[i]=min(r-i,R[2*pos-i]);
		else R[i]=1;
		while(t[i-R[i]]==t[i+R[i]])R[i]++;
		if(R[i]+i>r)
		{
		 pos=i;//最右位置对应的回文中心 
		 r=R[i]+i;//最右位置 
	    }
		ans=max(ans,R[i]);
	}
	return ans-1;
}
int main()
{
	while(~scanf("%s",w))
    {
    	int a=init(w,t);
    	printf("%d\n",manacher(a));
    }
	return 0;
}

板子2（福州大学模板）

增加了能判断原串中[L,R]对应子串是否为回文串功能，

这对调下标万年调不对的我来说，简直是福音

2020年3月20日，调了一下下标，增加了index函数

#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define rep(i, a, b) for(int i=(a); i<(b); i++)
#define per(i, a, b) for(int i=(b)-1; i>=(a); i--)
#define sz(a) (int)a.size()
#define de(a) cout << #a << " = " << a << endl
#define dd(a) cout << #a << " = " << a << " "
#define all(a) a.begin(), a.end()
#define pw(x) (1ll<<(x))
#define endl "\n"
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef double db;
#define mem(a,x) memset(a,x,sizeof(a))

const int N=3e5+5;
char s[N];
int pa[N<<1],n,ret[N];

void Manacher(char *s,int n,int *pa){
	//s[0]对应pa[0]
	//s[1]对应pa[2]
	//s[i]对应pa[2*i]
	//奇数位置为原#号位置 此处pa半径是不考虑#的 但考虑中心 如aba对于b的回文半径为2 
    pa[0] = 1;
    for(int i=1,j=0;i<(n<<1)-1;++i){
        int p = i >> 1 , q = i - p , r = ((j + 1)>>1) + pa[j] - 1;
        pa[i] = r < q ? 0 : min(r - q + 1 , pa[(j<<1) - i]);
        while(0 <= p - pa[i] && q + pa[i] < n && s[p - pa[i]] == s[q + pa[i]])
            pa[i]++;
        if(q + pa[i] - 1 > r) j = i;
    } 
}

inline bool check(int L,int R) {//询问[L,R]是否为回文串 0<=L<=R<n
    int mid=L+R>>1;
    //实质是 回文半径 应为 串长/2向上取整 aba对应2 bb对应1 
    if ((L&1)==(R&1)) return pa[mid<<1]>=R-L+2>>1;//长度为奇的串 a#b#a 0-2 (R-L+2)/2==2 即ab的长度 
    return pa[mid<<1|1]>=R-L+1>>1;//长度为偶的串 b#b 0-1 (R-L+1)/2==1 即b的长度 (不考虑#)
}

inline void index(int j,int &L,int &R) {//对于pa[j]确定的回文中心，对应的回文串在原串的[L,R]区间 
	int len=pa[j]<<1;
	if(!(j&1))len--;//字母计入回文半径故要减1
	L=(j>>1)-pa[j]+1;
	R=L+len-1;
}

inline void Put(string &s,vi &V) {
    V.clear();
    rep(i,0,sz(s)) V.pb(s[i]-'a');
}

int main() {
    
    while (scanf("%s",s)!=EOF) {
        n=strlen(s);
        Manacher(s,n,pa);
    }
    
    return 0;
}