bzoj 2160（manacher+差分）

传送门

问题：对于给出的长度为n的字符串如果其回文串的数量比K少则输出-1, 否则输出所有回文串中长度为奇数的最长的前K个回文串的长度的乘积, 结果对于19930726取模输出。

由于只考虑奇数长度的回文子串，所以不用构造新串，直接在原串头尾添加”+”,”-“做标记后跑manacher，求出回文半径p[i]后用差分前缀和思想统计各长度回文串个数，最后乘的时候套个快速幂即可。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e6+4;
const ll mod=19930726;
char b[maxn];
int p[maxn];
int n;
ll k,sum[maxn];
inline void manacher() {
int id=0,mx=0;
for (register int i=1;i<=n;++i) {
p[i]=(mx>i)?min(mx-i,p[(id<<1)-i]):1;
while (b[i+p[i]]==b[i-p[i]]) ++p[i];
if (mx<i+p[i]) mx=i+p[i],id=i;
++sum[0],--sum[p[i]+1];
}
}
ll fpow(ll a,ll b,ll mod) {
ll ans=1;
while (b) {
if (b&1) ans=ans*a%mod;
a=a*a%mod;
b>>=1;
}
return ans;
}
int main() {
//  freopen("bzoj 2160.in","r",stdin);
memset(sum,0,sizeof(sum));
scanf("%d%lld",&n,&k);
scanf("%s",b+1);
b[0]='+',b[n+1]='-';
manacher();
for (register int i=1;i<=n;++i)
sum[i]+=sum[i-1];
ll cnt=0,res=1;
for (register int i=n;i>0;--i) {
if (sum[i]+cnt<k) res=res*fpow((i<<1)-1,sum[i],mod)%mod,cnt+=sum[i];
else {res=res*fpow((i<<1)-1,k-cnt,mod)%mod,cnt=k;break;}
}
if (cnt^k) puts("-1");
else printf("%lld\n",res);
return 0;
}