E - 拉拉队排练 HYSBZ - 2160

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

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

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

Sample Input
5 3
ababa
Sample Output
45
【样例说明】
和谐小群体女生所拿牌子上写的字母从左到右按照女生个数降序排序后为ababa, aba, aba, bab, a, a, a,

题解：回文自动机+快速幂，我们要统计连续的女生有奇数个，且是回文，我们把len数组是奇数的取下来。然后再用cnt就能知道长度多少的出现了几次，用快速幂做。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const int maxn=1e6+1000;
const int mod=19930726;
int m;
ll k;
ll powquick(ll a,ll b)
{
  ll res=1;
  while(b)
  {
    if(b&1) res=res*a%mod;
    a=a*a%mod;
    b>>=1;
  }
  return res;
  cout<<res<<endl;
}
struct PAM
{
  int len[maxn],fail[maxn];
  int ch[maxn][26],cnt[maxn];
  int num[maxn],S[maxn];
  int last,p,n;
  int nowcode(int x)
  {
    for(int i=0;i<26;i++) ch[p][i]=0;
    cnt[p]=0;num[p]=0;;len[p]=x;
    return p++;
  }
  void init()
  {
    last=p=n=0;
    nowcode(0);
    nowcode(-1);
    S[0]=-1;
    fail[0]=1;
  }
  int get_fail(int x)
  {
    while(S[n-len[x]-1]!=S[n]) x=fail[x];
    return x;
  }
  void add(int c)
  {
    c=c-'a';
    S[++n]=c;
    int cur=get_fail(last);
    if(!ch[cur][c])
    {
      int now=nowcode(len[cur]+2);
      fail[now]=ch[get_fail(fail[cur])][c];
      ch[cur][c]=now;
      num[now]=num[fail[cur]]+1;
    }
    last=ch[cur][c];
    cnt[last]++;
  }
  void count()
  {
    for(int i=p-1;i>=2;i--) cnt[fail[i]]+=cnt[i];
  }
}T;
struct node
{
  int val,num;
  friend bool operator <(const node &a,const node &b)
  {
    return a.val>b.val;
  }
}p[maxn];
char s[maxn];
int main()
{
  T.init();
  scanf("%d %lld",&m,&k);
  scanf("%s",s);
  int n=strlen(s);
  for(int i=0;i<n;i++) {
    T.add(s[i]);
  };
  T.count();
  int tot=0;
  for(int i=2;i<=T.p;i++) if(T.len[i]&1) p[++tot]=node{T.len[i],T.cnt[i]};
  sort(p+1,p+1+tot);
//  for(int i=1;i<=tot;i++) cout<<p[i].val<<"   "<<p[i].num<<endl;
  ll ans=1;
  for(int i=1;i<=tot;i++)
  {
    if(k>=p[i].num)
    {
      k-=p[i].num;
    //  cout<<1ll*p[i].val<<"    "<<1ll*p[i].num<<endl;
      ans=ans*powquick(1ll*p[i].val,1ll*p[i].num)%mod;
    }
    else{
      ans=ans*powquick(1ll*p[i].val,k)%mod;
      k=0;
    }
  }
  if(k!=0) printf("-1\n");
  else printf("%lld\n",ans);
}