hdu5371 Hotaru's problem (Manacher)

/*像我这种菜鸡果然还是要写写博客，不然有的题目真是过去就过去了的状态。自己在硬盘里的代码也不太好管理，等以后写个软件来搞就好了 -  - */

题意:

给出一个串，找出满足以下条件的最长子串：

1. the first part is the same as the thrid part,

第一个部分和第三个部分一样。
2. the first part and the second part are symmetrical.

第一个部分和第二个部分对称。

思路：

这道题其实就是马拉车算法的简单应用

从题意出发  我们发现其实他想说的是 ：

（设把s的一个子串分成三个子串   第i个子串为si）

|s1|==|s2|==|s3|

s1+s2是回文串

s2+s3是回文串

所以我们用马拉车(就爱说这个 -   -）可以知道其实就是求最长的一个子串满足：

设:s1和s2中间的分隔符为中心的回文串的长度p1

     s2和s3中间的分隔符为中心的回文串的长度p2

min(p1,p2)-1>=|s2|

一开始是对于每个分隔符暴力找  若找到则记录长度 

但是TLE了 -  -

所以  对j的循环做了一些优化 若小于所记录的就跳出

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cassert>
using namespace std;

const int MAXN = 1e6+10;

int p[MAXN];

//适用于正负整数
template <class T>
inline bool scan_d(T &ret) {
   char c; int sgn;
   if(c=getchar(),c==EOF) return 0; //EOF
   while(c!='-'&&(c<'0'||c>'9')) c=getchar();
   sgn=(c=='-')?-1:1;
   ret=(c=='-')?0:(c-'0');
   while(c=getchar(),c>='0'&&c<='9') ret=ret*10+(c-'0');
   ret*=sgn;
   return 1;
}

int Proc(int pszIn[],int pszOut[],int len)
{
    int nLen=1;
    pszOut[0]=0;
    int i=0;
    while(i<len)
    {
        pszOut[nLen++]=-1;
        pszOut[nLen++]=pszIn[i];
        i++;
    }
    pszOut[nLen++]=-1;
    pszOut[nLen]=-2;
    return nLen;
}

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

int check(int len)
{
    long long cnt=0;
    int ret=0;
    for(int i=1;i<=len;i+=2)
    {
        if(p[i]>2)
            for(int j=(i+min(((p[i]>>1)<<1),((len>>1)/3)*2));j>i&&ret<(j-i)+1;j-=2,cnt++)
                if(min(p[i],p[j])>=(j-i)+1)
                    ret=(j-i)+1;
    }
    return ret;
}

int in[MAXN],out[MAXN<<2];
int main()
{
    int T,cas=1;
    scan_d(T);
    while(T--)
    {
        int n;
        scan_d(n);
        for(int i=0;i<n;i++)
            scan_d(in[i]);
        n=Proc(in,out,n);
        Manacher(p,out,n);
        printf("Case #%d: %d\n",cas++,(check(n)/2)*3);
    }
    return 0;
}