【HDOJ 5371】 Hotaru's problem

【HDOJ 5371】 Hotaru’s problem

Manacher算法+穷举/set

Manacher算法一好文:http://blog.csdn.net/yzl_rex/article/details/7908259

套一个Manacher算出回文半径数组p之后 有两种方法
穷举法:
枚举-1的点(根据题意只必为偶数回文) 找在该点回文半径内与其相隔最远 并且回文半径等于他俩距离(即两点为中心的回文串相同) 的点 记录找到时的距离 不断枚举找最大值即为最大回文串长 串长/2*3即为答案

set法:
记录-1的位置 根据每个-1对应的回文半径排序 从大到小枚举加入set 每加入一个位置i后 查找大于i-p[i]的第一个数(lower_bound) 和小于等于i+p[i]的第一个数(–upper_bound) 该位置与这两个位置的距离即为各自组成的回文串长度的二分之一 找到最大的*3即为答案

代码如下:

//穷举
#include <iostream>
#include <cstdio>
#include <cstdlib>

using namespace std;

int nwm[200002];
int p[200002];

void Manacher(int n)//马拉车模板
{
    nwm[0] = -2;
    int i;
    for(i = 0; i < n; i++)
    {
        nwm[i*2+1] = -1;
        scanf("%d",&nwm[i*2+2]);
    }
    nwm[n*2+1] = -1;
    nwm[n*2+2] = -3;
    int maxid = 0,id;
    n = n*2+2;
    for(i = 2; i < n; ++i)
    {
        if(maxid > i) p[i] = min(p[id*2-i],maxid-i);
        else p[i] = 1;

        while(nwm[i+p[i]] == nwm[i-p[i]]) p[i]++;

        if(p[i]+i > maxid)
        {
            maxid = p[i]+i;
            id = i;
        }
    }
}

int main()
{
    int t,n,z = 0,i,mm,x;
    scanf("%d",&t);
    while(t--)
    {
        mm = 0;
        scanf("%d",&n);
        Manacher(n);

        for(i = 3; i +4< n*2+2; i+=2)//枚举-1
        {
            if(p[i]-1 > mm)
            {
                x = p[i]-1;//记录-1为中心的最大回文串长
                while(x > mm && p[i+x] < x)//枚举找满足题意的第二回文中心
                    x--;
                mm = max(mm,x);
            }
        }

        printf("Case #%d: %d\n",++z,mm/2*3);
    }
    return 0;
}

//set方法
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <set>
#include <algorithm>

using namespace std;

int nwm[200002];
int p[200002];
int id[200002],tp;
set <int> s;

void Manacher(int n)
{
    nwm[0] = -2;
    int i;
    for(i = 0; i < n; i++)
    {
        nwm[i*2+1] = -1;
        id[tp++] = i*2+1;//将-1位置记录
        scanf("%d",&nwm[i*2+2]);
    }
    nwm[n*2+1] = -1;
    nwm[n*2+2] = -3;
    int maxid = 0,id;
    n = n*2+2;
    for(i = 2; i < n; ++i)
    {
        if(maxid > i) p[i] = min(p[id*2-i],maxid-i);
        else p[i] = 1;

        while(nwm[i+p[i]] == nwm[i-p[i]]) p[i]++;

        if(p[i]+i > maxid)
        {
            maxid = p[i]+i;
            id = i;
        }
    }
}

bool cmp(int a,int b)//按-1对应回文距离由大到小排序
{
    return p[a] > p[b];
}

int main()
{
    int t,n,z = 0,i,mm,x,l,r;

    scanf("%d",&t);
    while(t--)
    {
        s.clear();
        tp = 0;
        mm = 0;
        scanf("%d",&n);
        Manacher(n);
        sort(id,id+tp,cmp);
        for(i = 0; i < tp; ++i)//回文距离从大到小枚举位置
        {
            s.insert(id[i]);
            l = *s.lower_bound(id[i]-p[id[i]]);//找左边第一个>=id[i]-p[id[i]]的位置
            r = *(--s.upper_bound(id[i]+p[id[i]]));//找右边第一个<=id[i]+p[id[i]]的位置
            mm = max(mm,max(id[i]-l,r-id[i]));
        }
        printf("Case #%d: %d\n",++z,mm/2*3);
    }

    return 0;
}