[bzoj2160]拉拉队排练

神犇曰：“这道题不就是manacher裸题么，随便切一切吧”

我：“manacher是个啥。。。我只会hash啊。。”

看完manacher教程后：嗯？这都是个卵啊。。我去写hash吧。。

于是hash的pow数组写挂了。。。一直在WA。。。

终于A掉了。。

正着hash一遍然后反着hash一遍，最后暴力二分就ok。。就酱紫。。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <string>
#include <vector>
using namespace std;

const int MOD=19930726;
const int N=2000000;


typedef unsigned long long int ull;

int n,f[N];
char str[N];
ull h1[N],h2[N],p[N]={1};
long long ans=1,k,v[N];

ull mul(ull a,ull b){
    ull t=1;
    while(b){
        if(b&1)t=(t*a)%MOD;
        a=(a*a)%MOD;
        b>>=1;
    }
    return t;
}

int main(){
    scanf("%d%lld",&n,&k);
    scanf("%s",str+1);
    for(int i=1;i<=n;i++)p[i]=p[i-1]*10007;
    for(int i=1;i<=n;i++)h1[i]=h1[i-1]+str[i]*p[i-1];
    for(int i=n;i>=1;i--)h2[i]=h2[i+1]+str[i]*p[n-i];

    for(int i=1;i<=n;i++){
        f[i]=1;
        for(int l=1,r=min(i-1,n-i);l<=r;){
            int mid=(l+r)>>1;
            /*
                s1: i-mid~i-1
                s2: i+1~i+mid
            */
            if((h1[i-1]-h1[i-mid-1])*p[n-i-mid]==(h2[i+1]-h2[i+mid+1])*p[i-mid-1]){
                l=mid+1;
                f[i]=mid*2+1;
            }else{
                r=mid-1;
            }
        }
    }

    for(int i=1;i<=n;i++)v[f[i]]++;
    for(int i=n;i>=1;i--){
        (ans*=mul(i,min(v[i],k)))%=MOD;
        k-=v[i];
        if(i!=1)v[i-2]+=v[i];
        if(k<=0)break;
    }
    printf("%lld\n",ans);
}