题目:a-z都表示一个0~25的整数,一个长度为len的回文串的价值由前(len+1)/2个字符组成的26进制整数表示,
给出一个字符串,询问第k小的回文子串的价值
思路:保存一下节点的价值
代码:
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<iostream>
#include<algorithm>
#include<ctime>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<list>
#include<numeric>
using namespace std;
#define LL long long
#define ULL unsigned long long
#define INF 0x3f3f3f3f3f3f3f3f
#define mm(a,b) memset(a,b,sizeof(a))
#define PP puts("*********************");
template<class T> T f_abs(T a){ return a > 0 ? a : -a; }
template<class T> T gcd(T a, T b){ return b ? gcd(b, a%b) : a; }
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
// 0x3f3f3f3f3f3f3f3f
const LL MOD=777777777;
const int MAXN = 100005, SIZE = 26;
LL score[20][MAXN],pow26[MAXN],num[20][26];
LL ask[20];
char str[MAXN];
int n,m;
struct Palindromic_Tree {
int next[MAXN][SIZE];//next指针
int fail[MAXN];//fail指针
LL cnt[MAXN];
int len[MAXN];//len[i]表示节点i表示的回文串的长度
int S[MAXN];//存放添加的字符 S[MAXN << 1]
int last;//指向上一个字符所在的节点,方便下一次add
int n;//字符数组指针 //lbd, rbd
int p;//节点指针
int NewNode(int l) {//新建节点
for(int i = 0; i < SIZE; ++i)
next[p][i] = 0;
cnt[p] = 0;
len[p] = l;
return p++;
}
void Init(/*int x*/) {//初始化
p = n = 0;
NewNode(0); NewNode(-1);
last = 0;
S[n] = -1;
fail[0] = 1;
mm(score[0],0);
mm(score[1],0);
}
int GetFail(int x/*, int d*/) {//和KMP一样,失配后找一个尽量最长的
while(S[n - len[x] - 1] != S[n]) x = fail[x];
return x;
}
void Add(int c/*, int d*/) {//d=0表示在首部添加字符,d=1表示在尾部添加字符
S[++n] = c;
int cur = GetFail(last);//通过上一个回文串找这个回文串的匹配位置
if(!next[cur][c]) {//如果这个回文串没有出现过,说明出现了一个新的本质不同的回文串
int now = NewNode(len[cur] + 2);//新建节点
fail[now] = next[GetFail(fail[cur]/*, d*/)][c];//和AC自动机一样建立fail指针,以便失配后跳转
next[cur][c] = now;
for(int i=0;i<m;i++)
score[i][now]=(score[i][cur]+pow26[(len[now]-1)/2]*num[i][c]%MOD)%MOD;
}
last = next[cur][c];
cnt[last]++;
}
void Count() {
for(int i = p - 1; i >= 0; --i)
cnt[fail[i]] += cnt[i];
}
}tree;
int ans[MAXN],now;
bool cmp(int x,int y){
return score[now][x]<score[now][y];
}
void solve(){
int sz=tree.p;
for(int i=0;i<m;i++){
now=i;
for(int j=2;j<sz;j++)
ans[j]=j;
sort(ans+2,ans+sz,cmp);
for(int j=2;j<sz;j++){
if(ask[i]>tree.cnt[ans[j]])
ask[i]-=tree.cnt[ans[j]];
else{
printf("%lld\n",score[i][ans[j]]);
break;
}
}
}
}
int main(){
int T;
pow26[0]=1;
for(int i=1;i<MAXN;i++)
pow26[i]=pow26[i-1]*26%MOD;
scanf("%d",&T);
while(T--){
tree.Init();
scanf("%d%d",&n,&m);
scanf("%s",str);
for(int i=0;i<m;i++){
scanf("%lld",&ask[i]);
for(int j=0;j<26;j++)
scanf("%lld",&num[i][j]);
}
for(int i=0;i<n;i++)
tree.Add(str[i]-'a');
tree.Count();
solve();
printf("\n");
}
return 0;
}