#include using namespace std; typedef signed long long ll; #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define FORR2(x,y,arr) for(auto& [x,y]:arr) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) template bool chmax(T &a, const T &b) { if(a bool chmin(T &a, const T &b) { if(a>b){a=b;return 1;}return 0;} //------------------------------------------------------- using VT = string; struct RollingHash { static const ll mo0=1000000021,mo1=1000000009; static ll mul0,mul1; static const ll add0=1000010007, add1=1003333331; static vector pmo[2]; VT s; int l; vector hash_[2]; void init(VT s) { this->s=s; l=s.size(); int i,j; hash_[0]=hash_[1]=vector(1,0); if(!mul0) mul0=10009+(((ll)&mul0+time(NULL))>>5)%1259,mul1=10007+(time(NULL)+((ll)&mul1)>>5)%2257; if(pmo[0].empty()) pmo[0].push_back(1),pmo[1].push_back(1); FOR(i,l) hash_[0].push_back((hash_[0].back()*mul0+add0+s[i])%mo0); FOR(i,l) hash_[1].push_back((hash_[1].back()*mul1+add1+s[i])%mo1); } ll hash(int l,int r) { // s[l..r] if(l>r) return 0; while(pmo[0].size()>32,Lb=(L<<32)>>32,Ra=R>>32,Rb=(R<<32)>>32; return (((Ra + La*pmo[0][RL])%mo0)<<32)|((Rb + Lb*pmo[1][RL])%mo1); } }; vector RollingHash::pmo[2]; ll RollingHash::mul0,RollingHash::mul1; pair,pair > manacher(string S) { int L=S.size(),i,j,k; vector rad(2*L-1,0); for(i=j=0;i<2*L-1;i+=k,j=max(j-k,0)) { while(i-j>=0 && i+j+1<=2*L-1 && S[(i-j)/2]==S[(i+j+1)/2]) j++; rad[i]=j; for(k=1;i-k>=0 && rad[i]-k>=0 && rad[i-k]!=rad[i]-k&&i+k> R[202020]; unordered_map num[202020]; RollingHash rh; void solve() { int i,j,k,l,r,x,y; string s; cin>>S; N=S.size(); rh.init(S); vector V=manacher(S).first; FOR(i,V.size()) if(V[i]>0) { int l,r; l=(i-(V[i]-1))/2; r=(i+(V[i]-1))/2; auto h=rh.hash(l,r); R[V[i]][h]={l,r}; num[V[i]][h]++; } for(i=N;i>=3;i--) { FORR2(a,b,R[i]) { auto h=rh.hash(b.first+1,b.second-1); num[i-2][h]+=num[i][a]; R[i-2][h]={b.first+1,b.second-1}; } } ll ma=0; FOR(i,N+1){ FORR2(a,b,R[i]) { num[i][a]*=i; int l=b.first; int r=b.second-1; while(r>=l) { auto h=rh.hash(l,r); if(R[r-l+1].count(h)) { num[i][a]+=num[r-l+1][h]; break; } r--; } ma=max(ma,num[i][a]); } } cout<