BZOJ 2342: [Shoi2011]双倍回文 (Manacher+set)

题意

如果一个回文串的前一半和后一半也是回文串，那就是一个双倍回文，求最长的双倍回文的长度。

思路

先manacher预处理出来所有p[i]是必须要的。
然后我选择了枚举所有p[i]，也就是枚举出双倍回文的大串的中点，然后判断前半部分的中点的p[j]是不是等于二分之一的p[i]，然后WA了。因为p[i]代表的是i这个点的最长回文半径像这组数据：ccddddcc就是不能通过的。
然后窝们就要好好yy一下了。。。首先我们能发现，不管是长串还是半串，中心都是在manacher处理后的’#’上面的。然后我们假设j为前半串的中心，i为长串的中心，那么j肯定是大于i-p[i]/2的并且j+p[j]肯定是大于等于i的，这个如果不满足这两个条件是会出现长串不是回文串的情况或者半串不是回文串的情况。
所以基于这两个条件我们可以枚举i然后把j放到set里面每次upper_bound一下找到满足条件的最大ans就可以了。

代码

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;
#define LL long long
#define Lowbit(x) ((x)&(-x))
#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1|1
#define MP(a, b) make_pair(a, b)
const int INF = 0x3f3f3f3f;
const int MOD = 1000000007;
const int maxn = 1e6 + 10;
const double eps = 1e-8;
const double PI = acos(-1.0);
typedef pair<int, int> pii;
char s[maxn], str[maxn];
int p[maxn], q[maxn];

void Manacher(int len)
{
    int mx = 0, id = 0;
    for (int i = 0; i < len; i++)
    {
        if (mx > i) p[i] = min(p[2*id-i], mx - i);
        else p[i] = 1;
        while (str[i-p[i]] == str[i+p[i]]) p[i]++;
        if (p[i] + i > mx)
            mx = i + p[i], id = i;
    }
}

int init()
{
    int len = strlen(s);
    int l = 0;
    str[l++] = '$';
    str[l++] = '#';
    for (int i = 0; i < len; i++)
    {
        str[l++] = s[i];
        str[l++] = '#';
    }
    str[l] = 0;
    return l;
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int len;
    scanf("%d%s", &len, s);
    int l = init();
    Manacher(l);
    set<int> st;
    for (int i = 1; i < l; i += 2)
        st.insert(i);
    int ans = 0;
    for (int i = 1; i < l; i += 2)
    {
        int left = i - p[i] / 2;
        if (left % 2 == 0) left++;
        while (left != i)
        {
            if (left + p[left] >= i)
            {
                ans = max(ans, (i - left) * 2);
                break;
            }
            else
                st.erase(left), left = *st.upper_bound(left);

        }
    }
    printf("%d\n", ans);
    return 0;
}