ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

PAM练习

2020-03-01 11:06:49  阅读:223  来源: 互联网

标签:cnt const int ll 练习 PAM mod


题:https://www.luogu.com.cn/problem/P1659

题意:问前k大的奇数长度的回文串的长度乘积;

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int M=1e6+6;
const int maxn=26;
const int mod=19930726;
char s[M];
struct node{
    ll val,num;
    bool operator<(const node &b)const{
        return val>b.val;
    }
}b[M];

struct pam{
    int  son[M][maxn],cnt[M],num[M],fail[M],len[M];
    char s[M];
    int tot,last;
    int newnode(int Len){
        for(int i=0;i<maxn;i++)
            son[tot][i]=0;
        cnt[tot]=0;
        num[tot]=0;
        fail[tot]=0;
        len[tot]=Len;
        return tot++;
    }
    void init(){
        s[0]='#';
        tot=0;
        last=0;
        newnode(0);
        newnode(-1);
        fail[0]=1;
    }
    int getfail(int p,int i){
        while(s[i-len[p]-1]!=s[i])
            p=fail[p];
        return p;
    }
    void solve(const char *buf){
        init();
        int n=strlen(buf+1);
        for(int i=1;i<=n;i++){
            s[i]=buf[i]-'a';
            int cur=getfail(last,i);
            if(!son[cur][s[i]]){
                int now=newnode(len[cur]+2);
                fail[now]=son[getfail(fail[cur],i)][s[i]];
                son[cur][s[i]]=now;
                num[now]=num[fail[now]]++;
            }
            cnt[last=son[cur][s[i]]]++;
        }
        for(int i=tot-1;i>=0;i--)
            cnt[fail[i]]+=cnt[i];
    }
}PAM;
ll ksm(ll a,ll b){
    ll t=1ll;
    while(b){
        if(b&1){
            t=(t*a)%mod;
        }
        b>>=1ll;
        a=(a*a)%mod;
    }
    return t;
}
int main(){
    ll n,k;
    scanf("%lld%lld",&n,&k);
    scanf("%s",s+1);
    PAM.solve(s);
    int m=0;
    for(int i=0;i<PAM.tot;i++){
    //    cout<<PAM.len[i]<<"!!"<<PAM.cnt[i]<<endl;
        if(PAM.len[i]&1){
            b[m].val=PAM.len[i];
            b[m].num=PAM.cnt[i];
            m++;
        }
        
    }
    sort(b,b+m);
    ll ans=1ll;
    for(int i=0;i<m;i++){
        ans=(1ll*ans*ksm(b[i].val,min(1ll*b[i].num,k)))%mod;
    //    cout<<ans<<endl;
        k-=b[i].num;
        if(k<=0)
            break;
    }
    printf("%lld\n",ans);
    return 0;
}
View Code

 

标签:cnt,const,int,ll,练习,PAM,mod
来源: https://www.cnblogs.com/starve/p/12388643.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有