BZOJ 2160: 拉拉队排练 manacher

2160: 拉拉队排练

Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 1673 Solved: 664
[Submit][Status][Discuss]
Description

艾利斯顿商学院篮球队要参加一年一度的市篮球比赛了。拉拉队是篮球比赛的一个看点,好的拉拉队往往能帮助球队增加士气,赢得最终的比赛。所以作为拉拉队队长的楚雨荨同学知道,帮助篮球队训练好拉拉队有多么的重要。拉拉队的选拔工作已经结束,在雨荨和校长的挑选下,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个测试点,数据范围满足: 这里写图片描述

题解:不难看出是求回文串,先跑一遍manacher,然后把回文串长度加入一个数组中(类似于桶排的方法),这样我们就可以很容易地知道每种长度的回文串有多少个。然后从大到小枚举,如果没到k的话,就用快速幂把这个长度的所有都乘起来,然后把这些的长度-2,加入下一个桶中,(一个长度为7的回文串一定有一个长度为5的回文串),继续判断。如果到了k个,就只乘k次就好了,然后输出答案。
这道题其实并不难,只是有点坑而已……k就要用long long,以及模数。某大佬还往我身上泼脏水,我的快速幂明明没问题!

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;

const int mod = 19930726;
const int N = 1000010;

char s[N],s1[N<<1];
int cnt=0,val[N<<1];
ll len[N];
int n;
ll k;
int num[N];

inline int Min(int a,int b){
    return a<b?a:b;
}

void manacher(){
    for(register int i=0;i<n;i++)
        s1[++cnt]='#',s1[++cnt]=s[i];
    s1[0]='+',s1[++cnt]='#';
    int mx=0,id;
    for(register int i=1;i<=cnt;i++){
        if(i<=mx) val[i]=Min(mx-i+1,val[2*id-i]);
        else val[i]=1;
        while(s1[val[i]+i]==s1[i-val[i]]) val[i]++;
        if(val[i]+i-1>mx) mx=val[i]+i-1,id=i;
    }
}

ll mpow(ll a,ll b){
    ll rt=1;
    for(rt;b;b>>=1,a=a*a%mod)
        if(b&1) rt=rt*a%mod;
    return rt;
}

ll ans=1;

int main(){
    scanf("%d%lld",&n,&k);
    scanf("%s",s);
    manacher(); 
    for(register int i=2;i<=cnt;i+=2) ++num[val[i]-1];
    if(!(n&1)) n-=1;
    int i;
    for(i=n;i>=1;i-=2){
        if(num[i]>=k) {(ans*=mpow(i,k))%=mod;break;}
        else ans=(ans*=mpow(i,num[i]))%=mod;num[i-2]+=num[i];k-=num[i];
    }
    if(i) printf("%lld",ans);
    else printf("-1");
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值