[模板] manacher(教程)

还是不会马拉车啊.今天又学了一遍,在这里讲一下.

其实就是一个很妙的思路,就是设置一个辅助的数组len,记录每个点的最大对称长度,然后再存一个mx记录最大的对称子串的右端点.先开二倍数组,然后一点点扩大.有两种情况,一种i比mx大,这种只好暴力匹配了.还有一种就是i<mx,然后比较mx-i和len[2*po - i](po就是mx的对称中心),假如mx>i,那么说明无法再扩展,len[i] = len[2 * po - i];否则还需要暴力匹配.假如不好想,自己画个图就行了.

这里还有一个更负责的教程:

https://www.cnblogs.com/z360/p/6375514.html

大家可以对照看一看.

代码:

#include<iostream>
#include<cstdio>
#include<cmath>
#include<ctime>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
#define duke(i,a,n) for(register int i = a;i <= n;i++)
#define lv(i,a,n) for(register int i = a;i >= n;i--)
#define clean(a) memset(a,0,sizeof(a))
const int INF =  << ;
typedef long long ll;
typedef double db;
template <class T>
void read(T &x)
{
    char c;
    bool op = ;
    while(c = getchar(), c < '' || c > '')
        if(c == '-') op = ;
    x = c - '';
    while(c = getchar(), c >= '' && c <= '')
        x = x *  + c - '';
    if(op) x = -x;
}
template <class T>
void write(T x)
{
    if(x < ) putchar('-'), x = -x;
    if(x >= ) write(x / );
    putchar('' + x % );
}
char s1[];
char s[];
int len[];
int l = ;
int po = ;
int maxr = ;
int main()
{
    scanf("%s",s1 + );
    l = strlen(s1 + );
    s[] = '@';
    duke(i,, * l + )
    {
        if(i %  == )
        s[i] = '*';
        else
        s[i] = s1[i / ];
    }
    /*duke(i,0,2 * l + 1)
    printf("%c ",s[i]);
    puts("");*/
    int mx = -,ans = ,po = ;
    duke(i,, * l + )
    {
        if(mx > i)
        len[i] = min(mx - i,len[ * po - i]);
        else
        len[i] = ;
        while(s[i - len[i]] == s[i + len[i]])
        len[i]++;
        if(len[i] + i > mx)
        {
            mx = len[i] + i;
            po = i;
        }
        ans = max(ans,len[i] - );
    }
    printf("%d\n",ans);
    return ;
}