[JSOI2016]扭动的回文串

题目

非常板子了

看到求什么最长的回文，我们就想到枚举回文中心的方法

首先对于这个回文串只包含在一个串内的情况，我们随便一搞就可以了，大概\(Manacher\)一下就没有了

对于那种扭动的回文串，我们枚举回文中心，求一下回文半径，我们发现其必须先在一个串内扩展一个最长回文半径的长度，再去另外一个串扩展才是最优的

于是我们对这个两个串正反建\(SA\)求一下\(lcp\)就没了

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define re register
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
const int maxn=4e5+10;
char A[100005],B[100005],S[maxn];
int sa[maxn],rk[maxn],tax[maxn],tp[maxn],het[maxn];
int lg[maxn],St[maxn][20];
int ans=0,L,n,m;
int a[200005][2],b[200005][2];
inline void qsort() {
	for(re int i=0;i<=m;i++) tax[i]=0;
	for(re int i=1;i<=L;i++) tax[rk[i]]++;
	for(re int i=1;i<=m;i++) tax[i]+=tax[i-1];
	for(re int i=L;i;--i) sa[tax[rk[tp[i]]]--]=tp[i];
}
inline int ask(int l,int r) {
	int k=lg[r-l+1];
	return min(St[l][k],St[r-(1<<k)+1][k]);
}
inline int LCP(int i,int j) {
	if(rk[i]>rk[j]) std::swap(i,j);
	return ask(rk[i]+1,rk[j]);
}
int main() {
	scanf("%d",&n);
	scanf("%s",A+1);scanf("%s",B+1);
	for(re int i=1;i<=n;i++) S[i]=A[i],a[i][0]=i;
	L=n+1;S[L]='$';
	for(re int i=1;i<=n;i++)S[++L]=B[i],b[i][0]=L;
	S[++L]='#';
	for(re int i=n;i;--i) S[++L]=A[i],a[i][1]=L;
	S[++L]='&';
	for(re int i=n;i;--i) S[++L]=B[i],b[i][1]=L;
	for(re int i=1;i<=L;i++) rk[i]=S[i],tp[i]=i;
	m=255;qsort();
	for(re int w=1,p=0;p<L;w<<=1,m=p) {
		p=0;
		for(re int i=1;i<=w;i++) tp[++p]=L-w+i;
		for(re int i=1;i<=L;i++) if(sa[i]>w) tp[++p]=sa[i]-w;
		qsort();
		for(re int i=1;i<=L;i++) std::swap(rk[i],tp[i]);
		rk[sa[1]]=p=1;
		for(re int i=2;i<=L;i++)
			rk[sa[i]]=(tp[sa[i]]==tp[sa[i-1]]&&tp[sa[i]+w]==tp[sa[i-1]+w])?p:++p;
	}
	int k=0;
	for(re int i=1;i<=L;i++) {
		if(k) --k;
		int j=sa[rk[i]-1];
		while(S[i+k]==S[j+k]) ++k;
		het[rk[i]]=k;
	}
	for(re int i=2;i<=L;i++) lg[i]=lg[i>>1]+1;
	for(re int i=1;i<=L;i++) St[i][0]=het[i];
	for(re int j=1;j<=lg[L];j++)
		for(re int i=1;i+(1<<j)-1<=L;i++)
			St[i][j]=min(St[i][j-1],St[i+(1<<(j-1))][j-1]);
	for(re int i=1;i<n;i++) {
		int now=LCP(a[i+1][0],a[i][1]);
		ans=max(2*now,ans);
		if(i-now>=1) {
			int t=LCP(a[i-now][1],b[i+now][0]);
			ans=max(ans,2*(now+t));
		}
	}
	for(re int i=1;i<n;i++) {
		int now=LCP(b[i+1][0],b[i][1]);
		ans=max(2*now,ans);
		if(i+now+1<=n) {
			int t=LCP(b[i+1+now][0],a[i-now+1][1]);
			ans=max(ans,2*(now+t));
		}
	}
	for(re int i=2;i<n;i++) {
		int now=LCP(a[i+1][0],a[i-1][1]);
		ans=max(2*now+1,ans);
		if(i-now-1>=1) {
			int t=LCP(a[i-now-1][1],b[i+now][0]);
			ans=max(ans,2*(now+t)+1);
		}
	}
	for(re int i=2;i<n;i++) {
		int now=LCP(b[i+1][0],b[i-1][1]);
		ans=max(2*now+1,ans);
		if(i+now+1<=n) {
			int t=LCP(b[i+now+1][0],a[i-now][1]);
			ans=max(ans,2*(now+t)+1);
		}
	}
	printf("%d\n",ans);
	return 0;
}