[BZOJ2160]拉拉队排练

Description

艾利斯顿商学院篮球队要参加一年一度的市篮球比赛了。拉拉队是篮球比赛的一个看点，好的拉拉队往往能帮助球队增加士气，赢得最终的比赛。所以作为拉拉队队长的楚雨荨同学知道，帮助篮球队训练好拉拉队有多么的重要。拉拉队的选拔工作已经结束，在雨荨和校长的挑选下，n位集优秀的身材、舞技于一体的美女从众多报名的女生中脱颖而出。这些女生将随着篮球队的小伙子们一起，和对手抗衡，为艾利斯顿篮球队加油助威。一个阳光明媚的早晨，雨荨带领拉拉队的队员们开始了排练。n个女生从左到右排成一行，每个人手中都举了一个写有26个小写字母中的某一个的牌子，在比赛的时候挥舞，为小伙子们呐喊、加油。雨荨发现，如果连续的一段女生，有奇数个，并且他们手中的牌子所写的字母，从左到右和从右到左读起来一样，那么这一段女生就被称作和谐小群体。现在雨荨想找出所有和谐小群体，并且按照女生的个数降序排序之后，前K个和谐小群体的女生个数的乘积是多少。由于答案可能很大，雨荨只要你告诉她，答案除以19930726的余数是多少就行了。

Input

输入为标准输入。第一行为两个正整数n和K，代表的东西在题目描述中已经叙述。接下来一行为n个字符，代表从左到右女生拿的牌子上写的字母。

Output

输出为标准输出。输出一个整数，代表题目描述中所写的乘积除以19930726的余数，如果总的和谐小群体个数小于K，输出一个整数-1。

Sample Input
5 3
ababa

Sample Output
45

先跑遍Manacher，然后记录各个长度的回文串的出现次数。由于长度大的回文串肯定包含长度小的回文串，因此我们继续处理一下，就可以贪心了

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline int read(){
	int x=0,f=1;char ch=getchar();
	for (;ch<'0'||ch>'9';ch=getchar())	if (ch=='-')    f=-1;
	for (;ch>='0'&&ch<='9';ch=getchar())	x=(x<<1)+(x<<3)+ch-'0';
	return x*f;
}
inline void print(int x){
	if (x>=10)     print(x/10);
	putchar(x%10+'0');
}
const int N=1e6,P=19930726;
char s[N*2+10];
int p[N*2+10],cnt[N*2+10];
int mlt(int a,int b){
	int res=1;
	for (;b;b>>=1,a=1ll*a*a%P)	if (b&1)	res=1ll*res*a%P;
	return res;
}
int main(){
	int len=read();
	ll k;
	scanf("%lld%s",&k,s+1);
	for (int i=len;i;i--)	s[i<<1]=s[i],s[i<<1|1]='&';
	len=len<<1|1;
	s[0]='#',s[1]='&',s[len+1]='^';
	int Max=0,ID=0,ans=1;
	for (int i=1;i<=len;i++){
		p[i]=Max>i?min(p[ID*2-i],Max-i):1;
		while (s[i+p[i]]==s[i-p[i]])	p[i]++;
		if (Max<p[i]+i)	Max=p[ID=i]+i;
	}
	for (int i=2;i<len;i+=2)	cnt[p[i]/2]++;
	for (int i=len/2;i;i--)		cnt[i]+=cnt[i+1];
	for (int i=len/2;i;i--){
		if (k>cnt[i])	k-=cnt[i],ans=1ll*ans*mlt(i*2-1,cnt[i])%P;
		else{
			ans=1ll*ans*mlt(i*2-1,k)%P;
			break;
		}
	}
	printf("%d\n",ans);
	return 0;
}