Manacher算法常数时间内求任意子串是否回文

先上代码

void Manacher(string s, int Len[], int len)
{
    char tmp[2*len+4];
    tmp[0] = '@';
    for(int i = 1; i <= 2*len; i += 2)
    {
        tmp[i] = '#';
        tmp[i+1] = s[i/2];
    }
    tmp[2*len+1] = '#';
    tmp[2*len+2] = '$';

    int mx = 0, po = 0;
    for(int i = 1; i <= 2*len+1; i++)
    {
        if(mx > i)
        {
            Len[i] = min(mx - i, Len[2*po - i]);
        }
        else
        {
            Len[i] = 1;
        }
        while(tmp[i - Len[i]] == tmp[i + Len[i]])
        {
            Len[i]++;
        }
        if(i + Len[i] > mx)
        {
            mx = i + Len[i];
            po = i;
        }
    }
}

bool isPalindrom(int Len[], int left, int right)
{
    int middle = left + right + 2;
    return 2*Len[middle] - 1 >= 2*(right - left) ? true : false;
}

解释

首先我们要明白Manacher算法的原理(如果不太明白的可以点击此处学习)。

如果我们我们要判断 s[L : R] 是否为回文串，首先，我们需要找到 tmp 中与 s[L : R] 对应的子串位置，记为 tmp[L’ : R’]，其中 L’ = 2 * L + 2，R’ = 2 * R + 2。然后求出 L’ 和 R’ 的中点 middle = (L’ + R’) / 2，判断 2 * Len[middle] - 1 是否大于等于 R’ - L’，如果大于等于，s[L : R] 是回文串，否则不是。

这样判断的原因：2 * Len[middle] - 1 记录的是以 tmp[middle] 为中心的最长回文子串的长度，tmp[L’ : R’] 也是以 tmp[middle] 为中心的子串。如果 2 * Len[middle] - 1 大于等于 R’ - L’，tmp[L’ : R’] 一定是回文串，那么 s[L : R] 一定是回文串，否则一定不是回文串。

适用场景

需要频繁地判断任意子串是否回文的题目