拉拉队排练 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, b, b，前三个长度的乘积为。

Hint

 

总共20个测试点，数据范围满足：

题目大意：将回文串按照长度从大到小选择k个出来,并求出长度乘积。

解题思路：求出相同长度的回文子串的个数保存在数组中，对其进行遍历，根据题目要求求出乘积即可。因为后面的k的数据量很大，所以在进行乘积的时候需要用快速幂优化。

/*
@Author: Top_Spirit
@Language: C++
*/
#include <bits/stdc++.h>
using namespace std ;
typedef unsigned long long ull ;
typedef long long ll ;
const int Maxn = 1e6 + 10 ;
const int INF = 0x3f3f3f3f ;
const double PI = acos(-1.0) ;
const ull seed = 133 ;
const int Mod = 19930726 ;

ll n, k ;

struct palindromic_tree{
    int Next[Maxn][26] ;
    int fail[Maxn] ;
    ll cnt[Maxn] ;
    int num[Maxn] ;
    int len[Maxn] ;
    int s[Maxn] ;
    int last ;
    int n, p ;
    int newNode (int k){
        for (int i = 0; i < 26; i++) Next[p][i] = 0 ;
        cnt[p] = 0 ;
        num[p] = 0 ;
        len[p] = k ;
        return p++ ;
    }
    void init(){
        p = 0 ;
        newNode(0) ;
        newNode(-1) ;
        last = 0 ;
        n = 0 ;
        s[n] = -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 -= 'a' ;
        s[++n] = c ;
        int cur = get_fail(last) ;
        if (!Next[cur][c]){
            int Now = newNode(len[cur] + 2) ;
            fail[Now] = Next[get_fail(fail[cur])][c] ;
            Next[cur][c] = Now ;
            num[Now] = num[fail[Now]] + 1 ;
        }
        last = Next[cur][c] ;
        cnt[last]++ ;
    }
    void Count(){
        for (int i = p - 1; i >= 0; i--){
            cnt[fail[i]] += cnt[i] ;
        }
    }
}Tree;

char str[Maxn] ;
ll lenNum[Maxn] ;

ll Pow (ll a, ll b){
    ll res = 1ll ;
    while (b) {
        if (b & 1) res = res * a % Mod ;
        a = a * a % Mod ;
        b >>= 1 ;
    }
    return res ;
}

int main (){
    cin >> n >> k ;
    cin >> str ;
//    cout << str << endl ;
    Tree.init() ;
    for (int i = 0; i < n; i++) Tree.add(str[i]) ;
    Tree.Count() ;
    for (int i = 2; i < Tree.p; i++){
        lenNum[Tree.len[i]] += Tree.cnt[i] ;
    }
    ll tmp = 0 ;
    ll ans = 1ll ;
    for (int i = n; i > 0; i--){
        if (!(i & 1)) continue ;
        if (!lenNum[i]) continue ;
        if (tmp + lenNum[i] >= k){
            ll sub = k - tmp ;
            ans = ans * Pow(i, sub) % Mod ;
            tmp = k ;
            break ;
        }
        ans = ans * Pow(i, lenNum[i]) % Mod ;
        tmp += lenNum[i] ;
    }
    if (tmp < k) cout << -1 << endl ;
    else cout << ans << endl ;
    return 0 ;
}