bzoj2342SHOI2011双倍回文

Posted AaronPolaris

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj2342SHOI2011双倍回文相关的知识,希望对你有一定的参考价值。

2342: [Shoi2011]双倍回文

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1798  Solved: 670
[Submit][Status][Discuss]

Description

技术分享

Input

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

Output

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

Sample Input

16
ggabaabaabaaball

Sample Output

12

HINT

N<=500000




manacher+set,思路好题

首先用manacher可以求出以i和i+1中间为对称轴,最长回文串能扩增的长度p[i]。

然后4*(y-x)能更新答案,当且仅当y≤x+p[x]/2且y-p[y]≤x。

按i-p[i]将所有点排序,依次插入set中,查询x+p[x]/2的前驱更新答案,时间复杂度O(n*logn)。




#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<set>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define maxn 1000005
using namespace std;
int m,n,ans,f[maxn],p[maxn],q[maxn];
char s[maxn],a[maxn];
set<int> st;
set<int>::iterator it;
inline bool cmp(int x,int y)
{
	return x-p[x]<y-p[y];
}
void manacher()
{
	int mx=0,id=0;
	F(i,1,m)
	{
		if (mx>i) f[i]=min(mx-i,f[id*2-i]);
		else f[i]=0;
		while (a[i+f[i]]==a[i-f[i]]) f[i]++;
		if (i+f[i]>mx) mx=i+f[i],id=i;
	}
}
int main()
{
	scanf("%d%s",&n,s+1);
	m=n<<1|1;a[1]='#';a[0]='*';a[m+1]='&';
	F(i,1,n) a[i<<1]=s[i],a[i<<1|1]='#';
	manacher();
	F(i,1,n) p[i]=(f[i<<1|1]-1)/2;
	F(i,1,n) q[i]=i;
	sort(q+1,q+n+1,cmp);
	int t=1;
	F(i,1,n)
	{
		while (t<=n&&q[t]-p[q[t]]<=i) st.insert(q[t++]);
		it=st.upper_bound(i+p[i]/2);
		if (it!=st.begin())  it--,ans=max(ans,(*it-i)*4);
	}
	printf("%d\n",ans);
	return 0;
}













以上是关于bzoj2342SHOI2011双倍回文的主要内容,如果未能解决你的问题,请参考以下文章

[BZOJ2342][Shoi2011]双倍回文

bzoj2342 [Shoi2011]双倍回文 (manacher)

bzoj2342 [Shoi2011]双倍回文

bzoj 2342: [Shoi2011]双倍回文

[BZOJ2342][SHOI2011]双倍回文

BZOJ 2342: [Shoi2011]双倍回文