# BZOJ_2160_拉拉队排练_manacher

BZOJ_2160_拉拉队排练_manacher

5 3
ababa

## Sample Output

45

manacher过程中求极长回文子串然后区间加差分一下。

#include <cstdio>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long ll;
#define mod 19930726
#define N 1000050
int c[N],p[N<<1],a[N<<1],n,d[N],e[N];
char w[N<<1];
ll K;
ll qp(ll x,ll y) {
ll re=1; for(;y;y>>=1ll,x=x*x%mod) if(y&1ll) re=re*x%mod; return re;
}
int main() {
scanf("%d%lld",&n,&K);
int i;
scanf("%s",w+1);
for(i=1;i<=n;i++) a[2*i-1]='$',a[2*i]=w[i]; n=n<<1|1; a[n]='$';
int mx=0,lst;
for(i=1;i<=n;i++) {
if(i<=mx) p[i]=min(p[2*lst-i],mx-i+1);
else p[i]=1;
while(i-p[i]>=1&&i+p[i]<=n&&a[i-p[i]]==a[i+p[i]]) p[i]++;
if(mx<i+p[i]-1) {mx=i+p[i]-1; lst=i;}
int lim=(i+p[i]-2)/2-(i-p[i]+2)/2+1;
if(i%2==0) {
d[1]++; d[(lim+1)/2+1]--;
}else {
e[1]++; e[lim/2+1]--;
}
}
n>>=1;
ll sum=0;
for(c[1]=sum=d[1],i=3;i<=n;i+=2) {
c[i]=c[i-2]+d[(i+1)>>1]; sum+=c[i];
}
for(i=2;i<=n;i+=2) {
c[i]=c[i-2]+e[i>>1];
}
if(sum<K) {
puts("-1"); return 0;
}
ll ans=1;
i=n; if(n%2==0) i--;
for(;i>=1;i-=2) {
// printf("%d\n",c[i]);
if(K<=c[i]) {
printf("%lld\n",ans*qp(i,K)%mod); return 0;
}
K-=c[i]; ans=ans*qp(i,c[i])%mod;
}
}


posted @ 2018-06-02 20:07  fcwww  阅读(185)  评论(0编辑  收藏  举报