BZOJ 2342: [Shoi2011]双倍回文 马拉车算法/并查集

2342: [Shoi2011]双倍回文

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1123  Solved: 408

题目连接

http://www.lydsy.com/JudgeOnline/problem.php?id=2342

Description

输入分为两行，第一行为一个整数n，表示字符串的长度，第二行有n个连续的小写的英文字符，表示字符串的内容。

Input

输入分为两行，第一行为一个整数，表示字符串的长度，第二行有个连续的小写的英文字符，表示字符串的内容。

Output

输出文件只有一行，即：输入数据中字符串的最长双倍回文子串的长度，如果双倍回文子串不存在，则输出0。

Sample Input

16
ggabaabaabaaball

Sample Output

12

HINT

 

题解：

用传说中的马拉车算法，算出所有的回文中心，以及P数组，然后我们就枚举i-(i-p[i])/2位置到第i位置的p[j]是否大于i

然后如果大于就说明这是一个双倍回文

但是，这样子会T，那么我们肿么优化呢？

那就献祭出并查集吧~！

♪(^∇^*)，100ms就跑过去啦

 

代码：

 

//qscqesze

#include <cstdio>

#include <cmath>

#include <cstring>

#include <ctime>

#include <iostream>

#include <algorithm>

#include <set>

#include <vector>

#include <sstream>

#include <queue>

#include <typeinfo>

#include <fstream>

#include <map>

typedef long long ll;

using namespace std;

//freopen("D.in","r",stdin);

//freopen("D.out","w",stdout);

#define sspeed ios_base::sync_with_stdio(0);cin.tie(0)

#define maxn 500005

#define mod 10007

#define eps 1e-9

const int inf=0x7fffffff;   //无限大

/*



*/

//**************************************************************************************

char s[maxn];

char str[maxn*2];

int p[maxn*2];

int fa[maxn*2];

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}

    return x*f;

}

int l;

int manacher(char s[],int l)

{

   int i,j,k,ans=0;

   for(i=1;i<=l;++i)str[i<<1]=s[i],str[(i<<1)+1]='#';

   str[1]='#';str[l*2+1]='#';str[0]='&';str[l*2+2]='$';

   l=l*2+1;j=0;

   for(i=1;i<=l;)

   {

       while(str[i-j-1]==str[i+j+1])++j;

       p[i]=j;if(j>ans)ans=j;

       for(k=1;k<=j&&p[i]-k!=p[i-k];++k)p[i+k]=min(p[i-k],p[i]-k);

       i+=k;j=max(j-k,0);

   }

   return ans;

}

int get(int p)

{

    if(fa[p]==p)

        return p;

    fa[p]=get(fa[p]);

    return fa[p];

}

int main()

{

    l=read();

    scanf("%s",s+1);

    manacher(s,l);

    l=l*2+1;

    for(int i=1;i<=l;i++)

    {

        if(str[i]=='#')

            fa[i]=i;

        else

            fa[i]=i+1;

    }

    int j;

    int ans=0;

    for(int i=3;i<l;i+=2)

    {

        j=get(max(i-p[i]/2,1));

        for(;j<i&&j+p[j]<i;fa[j]=get(j+1),j=fa[j]);

        if(j<i)

            if((i-j)*2>ans)

            ans=(i-j)*2;

    }

    cout<<ans<<endl;

}