HYSBZ - 2160 (回文树）

传送门

题意：给定一个字符串，找出该串中长度为奇数的最长$k$个子回文串，并且求出它们的乘积

题解: 由于$k$很大，所以不能简单的$sort$

用回文树处理出每个节点的回文串长度，并且记录下该节点下该回文串长度的串有多少个，然后快速幂模拟一下。
注意细节！

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAXN = 1e6+5;
const int N = 27;
const int mod = 19930726;

char s1[MAXN];

ll quick(ll a,ll b){
    ll res = 1;
    while(b){
        if(b&1) res = res*a%mod;
        a = a*a%mod;
        b >>= 1;
    }
    return res;
}
struct Palindromic_Tree {
	int next[MAXN][N] ;//next指针，next指针和字典树类似，指向的串为当前串两端加上同一个字符构成
	int fail[MAXN] ;//fail指针，失配后跳转到fail指针指向的节点
	int cnt[MAXN] ; //表示节点i表示的本质不同的串的个数（建树时求出的不是完全的，最后count()函数跑一遍以后才是正确的）
	int num[MAXN] ; //表示以节点i表示的最长回文串的最右端点为回文串结尾的回文串个数。
	int len[MAXN] ;//len[i]表示节点i表示的回文串的长度
	int S[MAXN] ;//存放添加的字符
	int last ;//指向上一个字符所在的节点，方便下一次add
	int n ;//字符数组指针
	int p ;//节点指针
	ll sum[MAXN];
    ll tol;
	int newnode ( int l ) {//新建节点
		for ( int i = 0 ; i < N ; ++ i ) next[p][i] = 0 ;
		cnt[p] = 0 ;
		num[p] = 0 ;
		len[p] = l ;
		return p ++ ;
	}

	void init () {//初始化
		p = tol = 0 ;
		newnode (  0 ) ;
		newnode ( -1 ) ;
		last = 0 ;
		n = 0 ;
		S[n] = -1 ;//开头放一个字符集中没有的字符，减少特判
		fail[0] = 1 ;
	}

	int get_fail ( int x ) {//和KMP一样，失配后找一个尽量最长的
		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] ;//和AC自动机一样建立fail指针，以便失配后跳转
			next[cur][c] = now ;
			num[now] = num[fail[now]] + 1 ;
		}
		last = next[cur][c] ;
		cnt[last] ++ ;
        //return len[last];
	}

	void count () {
		//父亲累加儿子的cnt，因为如果fail[v]=u，则u一定是v的子回文串！
        for(int i = p-1; i >= 0; i--){
            cnt[fail[i]] += cnt[i] ;
            cnt[fail[i]] %= mod;
        }
	}
	void solve(ll k){
        for(int i = 2; i <= p-1; i++){
                //if(len[i] <= 0) continue;　　！！
                if(len[i]%2==0) continue;
                sum[len[i]] = (sum[len[i]] + cnt[i])%mod;
                tol += cnt[i];
        }
        if(tol < k){
            puts("-1");
        }else{
            ll ans = 1;
            for(int i = n; i >= 1; i--){ //枚举回文串长度
                if(i%2==0) continue;
                if(sum[i] <= k){
                    ans = (ans*quick(i,sum[i]))%mod;
                    k -= sum[i];
                }else{
                    ans = (ans*quick(i,k))%mod;
                    k = 0;
                }
                if(!k) break;
            }
            printf("%lld\n",ans);
        }

	}

} pat;



int n,tmp[MAXN];
ll k;
int main(){
    scanf("%d%lld",&n,&k);
    scanf("%s",s1);
    pat.init();
    for(int i = 0; i < n; i++){
        pat.add(s1[i]);  //x1表示以i为结尾的最大回文串长度
    }
    pat.count();
    pat.solve(k);
    return 0;
}