BZOJ2803: [Poi2012]Prefixuffix

题目大意：给一个字符串，求一个最大长度L≤n/2,使得前L个字符与后L个字符循环同构

首先如果头和尾长度为L的子串循环同构，就一定能找到一个合理的X，使得原字符串的X前缀与X后缀匹配并且砍掉他们之后，剩下的L-X前缀与L-X后缀匹配

大概意思就是这样：

ab abbabab abba ab

所以不妨设f[I]表示原字符串砍掉i前缀和i后缀之后，首尾最大匹配长度

这样我们就只需要把满足原字符串i前缀和i后缀匹配的f[i]+i取个max即可

然后怎么求f[i]呢？

首先我们在纸上画一画会发现f[i]≥f[i-1]-2  =>   f[i]≤f[i+1]+2

所以我们就可以从f[n/2]开始倒着递推求出f数组

据说卡hash？反正我写了个双hash

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 1000010
using namespace std;
char s[N];
long long bas[2]={233,707},mod=1e9+7;
long long cf[N][2],h[N][2];
bool check(long long x,long long y,long long l)
{
	long long x4,x5,y4,y5;
	x4=((h[x+l-1][0]-h[x-1][0]*cf[l][0])%mod+mod)%mod;
	x5=((h[x+l-1][1]-h[x-1][1]*cf[l][1])%mod+mod)%mod;
	y4=((h[y+l-1][0]-h[y-1][0]*cf[l][0])%mod+mod)%mod;
	y5=((h[y+l-1][1]-h[y-1][1]*cf[l][1])%mod+mod)%mod;
	if(x4!=y4||x5!=y5) return false;
	return true;
}
long long f[N];
int main()
{
	long long n,i,j;
	scanf("%lld%s",&n,s+1);
	cf[0][0]=cf[0][1]=1;
	for(i=1;i<=n;i++)
	for(j=0;j<2;j++)
	{
		cf[i][j]=cf[i-1][j]*bas[j]%mod;
		h[i][j]=(h[i-1][j]*bas[j]+s[i]-96)%mod;
	}
	for(i=n/2;i>=1;i--)
	{
		f[i]=min(f[i+1]+2,(n-2*i)/2);
		while(f[i]&&!check(i+1,n-i-f[i]+1,f[i])) f[i]--;
	}
	long long ans=0;
	for(i=1;i<=n/2;i++)
	{
		if(!check(1,n-i+1,i)) continue;
		ans=max(i+f[i],ans);
	}
	printf("%lld",ans);
}