BZOJ 2565 最长双回文串 Hash+二分

23 篇文章 0 订阅
31 篇文章 0 订阅

题目大意:给定一个字符串,求一个最长的子串,该字串可以分解为两个回文子串

傻逼的我又忘了Manacher怎么写了= = 无奈Hash+二分吧

首先将字符串用分隔符倍增,然后求出以每个点为中心的最长回文半径

然后考虑两个回文串怎么合并成一个


我们发现图中以i为中心的回文串和以j为中心的回文串合并后长度恰好为(j-i)*2

能合并的前提是以两个点为中心的回文串有交点

那么对于每个j我们要求出有交点的最左侧的i

维护一个后缀min随便搞搞就可以了

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 200200
#define BASE 131
using namespace std;
typedef unsigned long long ll;
int n,ans;
int f[M],min_pos[M];
char s[M];
ll hash1[M],hash2[M],power[M];
bool Judge(int mid,int len)
{
	ll _hash1=hash1[mid+len-1]-hash1[mid-1]*power[len];
	ll _hash2=hash2[mid-len+1]-hash2[mid+1]*power[len];
	return _hash1==_hash2;
}
int Bisecition(int x)
{
	int l=1,r=min(x,n-x+1);
	while(l+1<r)
	{
		int mid=l+r>>1;
		if( Judge(x,mid) )
			l=mid;
		else
			r=mid;
	}
	return Judge(x,r)?r:l;
}
int main()
{
	int i;

	static char str[M];
	scanf("%s",str+1);
	for(s[1]='#',i=1;str[i];i++)
		s[i<<1]=str[i],s[i<<1|1]='#';
	n=strlen(s+1);

	for(power[0]=1,i=1;i<=n;i++)
		power[i]=power[i-1]*BASE;
	for(i=1;i<=n;i++)
		hash1[i]=hash1[i-1]*BASE+s[i];
	for(i=n;i;i--)
		hash2[i]=hash2[i+1]*BASE+s[i];

	memset(min_pos,0x3f,sizeof min_pos);
	for(i=1;i<=n;i++)
	{
		f[i]=Bisecition(i);
		min_pos[i+f[i]-1]=min(min_pos[i+f[i]-1],i);
	}
	for(i=n-1;i;i--)
		min_pos[i]=min(min_pos[i],min_pos[i+1]);
	for(i=1;i<=n;i++)
	{
		int temp=min_pos[i-f[i]+1];
		ans=max(ans,i-temp<<1);
	}
	cout<<(ans>>1)<<endl;
	return 0;
}


评论 4 您还未登录,请先 登录 后发表或查看评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
©️2022 CSDN 皮肤主题:大白 设计师:CSDN官方博客 返回首页

打赏作者

PoPoQQQ

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值