k-substrings

题目

给定一个长度为n的字符串S
我们设S的k-子串是S[k…n-k+1]，设字符串t是字符串T的“奇正确前后缀“当且仅当满足以下条件：

t长度为奇数
|t|<|T|
t是T的border（既是前缀又是后缀）
对于k=1…n/2上取整，求S的k-子串的最长“奇正确前后缀“长度。无解输出-1
2≤n≤1,000,000

代码

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <iostream>
#include <cstdlib>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define mo 1000000009
#define LL long long
#define N 1000005
using namespace std;
char ch[N];
int n,f[N];
LL h[N],K=119,cf[N];
bool pd(int p,int q,int le)
{
	int l=p-le+1,r=p+le-1,x=q-le+1,y=q+le-1;
	return ((h[r]-h[l-1]*cf[2*le-1]%mo+mo)%mo==(h[y]-h[x-1]*cf[2*le-1]%mo+mo)%mo);
}
int main()
{
	cin>>n;
	scanf("\n");
	cf[0]=1;
	fo(i,1,n) 
	{
		ch[i]=getchar();
		h[i]=(h[i-1]*K%mo+ch[i]-'a')%mo;
		cf[i]=cf[i-1]*K%mo;
	}
	int n1=(n+1)>>1;
	fo(i,1,n1)
	{
		if(n-i+1==i) break;
		int l=1,r=i;
		while(l+1<r)
		{
			int mid=(l+r)>>1;
			if(pd(i,n-i+1,mid)) l=mid;
			else r=mid;
		}
		if(pd(i,n-i+1,r)) f[i-r+1]=max(f[i-r+1],2*r-1);
		else if(pd(i,n-i+1,l)) f[i-l+1]=max(f[i-l+1],2*l-1);
	}
	fo(i,1,n1)
	{
		f[i]=max(f[i],f[i-1]-2);
		if(f[i]==0) printf("-1 ");
		else printf("%d ",f[i]);
	}
}