BZOJ2160: 拉拉队排练

题目大意：给定一个字符串，将这个串中所有的奇数长度回文串挑出来按长度降序排序，求前k个长度的乘积

manacher裸题啊....

先用manacher，求出每个节点的P值，然后把他们推到一个（队列？数组？）里，然后从大到小一个一个的乘起来就好了，注意长的字符串乘完之后长度减二还是一个合法字符串

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1000010
using namespace std;
int n;
long long k;
char c[N],s[N<<1];
int p[N<<1];
int A[N];
void manacher()
{
	int i,j;
	for(i=1;i<=n;i++)
	{
		s[2*i]=c[i];
		s[2*i-1]='#';
	}
	s[2*n+1]='#';
	s[0]='$';s[2*n+2]='&';
	int maxn=0,maxb=0;
	for(i=1;i<=2*n;i++)
	{
		if(i<maxn) p[i]=min(maxn-i,p[2*maxb-i]);
		else p[i]=0;
		while(s[i+p[i]+1]==s[i-p[i]-1]) p[i]++;
		if(i+p[i]>maxn)
		{
			maxn=i+p[i];
			maxb=i;
		}
	}
	for(i=2;i<=2*n;i+=2)
	A[p[i]]++;
}
int mod=19930726;
long long ksm(long long d,long long c)
{
	long long ret=1;
	while(c)
	{
		if(c&1) ret=ret*d%mod;
		d=d*d%mod;c/=2;
	}
	return ret;
}
int main()
{
	scanf("%d%lld",&n,&k);
	scanf("%s",c+1);
	manacher();
	long long ans=1,i;
	for(i=n;i>=1;i--)
	{
		if(A[i]>=k)
		{
			(ans*=ksm(i,k))%=mod;
			break;
		}
		else
		{
			(ans*=ksm(i,A[i]))%=mod;
			k-=A[i];
			A[i-2]+=A[i];
		}
	}
	if(i==0) puts("-1");
	else printf("%lld",ans);
}