#include <bits/stdc++.h>
using namespace std;
#define For(Ti,Ta,Tb) for(int Ti=(Ta);Ti<=(Tb);++Ti)
#define Dec(Ti,Ta,Tb) for(int Ti=(Ta);Ti>=(Tb);--Ti)
#define Debug(...) fprintf(stderr,__VA_ARGS__)
using ll=long long;
const int N=1e6+5,LogN=20;
int n,tote;char s[N];
ll ans;
struct Edge{
	int v,nxt;
}e[N*2];
struct Graph{
	int head;
	struct Iterator{
		int i;
		void operator++(){i=e[i].nxt;}
		bool operator!=(const Iterator &it) const{return i!=it.i;}
		int operator*() const{return e[i].v;}
	};
	Iterator begin(){return {head};}
	Iterator end(){return {0};}
}g[N*2];
void AddEdge(int u,int v){
	e[++tote]={v,g[u].head},g[u].head=tote;
}
struct SAM{
	struct Node{
		int len,link,ed,nxt[26];
	}t[N*2];
	int size,last,pres[N];
	int New(){
		t[size].len=t[size].link=t[size].ed=0;
		memset(t[size].nxt,-1,sizeof(t[size].nxt));
		return size++;
	}
	void Init(){
		size=last=0;
		New(),t[0].link=-1;
	}
	void Extend(char _c){
		int cur=New(),c=_c-'a';
		t[cur].len=t[last].len+1,t[cur].ed=1;
		int p=last;last=pres[t[cur].len]=cur;
		while(~p&&!~t[p].nxt[c]) t[p].nxt[c]=cur,p=t[p].link;
		if(!~p) return;
		int q=t[p].nxt[c];
		if(t[p].len+1==t[q].len) t[cur].link=q;
		else{
			int cl=size++;
			t[cl]={t[p].len+1,t[q].link,0};
			memcpy(t[cl].nxt,t[q].nxt,sizeof(t[cl].nxt));
			while(~p&&t[p].nxt[c]==q) t[p].nxt[c]=cl,p=t[p].link;
			t[q].link=t[cur].link=cl;
		}
	}
	vector<int> qry[N*2];
	pair<int,int> cur[N*2];
	int siz[N*2],curl;
	void GetSize(int u){
		siz[u]=t[u].ed;
		for(int v:g[u]) GetSize(v),siz[u]+=siz[v];
	}
	void DFS(int u){
		cur[++curl]={t[u].len,u};
		for(int v:g[u]) DFS(v);
		for(int len:qry[u]){
			int v=lower_bound(cur+1,cur+curl+1,make_pair(len,0))->second;
			ans=max(ans,1LL*siz[v]*len*len);
		}
		--curl;
	}
	void Build(){
		For(i,1,size-1) AddEdge(t[i].link,i);
		GetSize(0),DFS(0);
	}
}sam;
char t[N*2];int d[N*2],L[N*2],R[N*2],orig[N*2];
int main(){
	ios::sync_with_stdio(false),cin.tie(nullptr);
	cin>>(s+1),n=strlen(s+1);
	sam.Init();
	For(i,1,n) sam.Extend(s[i]);
	int len=0;
	t[0]='@',t[++len]='#';
	For(i,1,n) t[++len]=s[i],orig[len]=i,t[++len]='#';
	For(i,1,len) R[i]=(orig[i]?orig[i]:R[i-1]);
	Dec(i,len,1) L[i]=(orig[i]?orig[i]:L[i+1]);
	For(i,0,25) ans=max<ll>(ans,count(s+1,s+n+1,i+'a'));
	for(int i=1,l=0,r=0;i<=len;++i){
		if(r>=i) d[i]=min(r-i+1,d[l+r-i]);
		else d[i]=1;
		while(t[i-d[i]]==t[i+d[i]]){
			int _l=L[i-d[i]],_r=R[i+d[i]];
			if(_l<=_r&&orig[i-d[i]]) sam.qry[sam.pres[_r]].push_back(_r-_l+1);
			++d[i];
		}
		if(i+d[i]>r) l=i-d[i]+1,r=i+d[i]-1;
	}
	sam.Build();
	cout<<ans<<'\n';
	return 0;
}