#include <bits/stdc++.h>
#define int long long
using namespace std;
const int N=5e4+100;
int n,k;
char a[N];
struct node{
	int len,link,min,max;
	int ch[26];
	void clr(){ memset(this,0,sizeof(node)); }
};
node t[N<<1];
int lst,nc;
void init_sam(){
	t[0].clr();
	t[0].link=-1;
	lst=nc=0;
}
void ext(int c){
	c-='a';
	int now=++nc;
	t[now].clr();
	t[now].len=t[lst].len+1;
	int p=lst;
	while(~p&&!t[p].ch[c]){
		t[p].ch[c]=now;
		p=t[p].link;
	}
	if(p==-1){
		t[now].link=0;
	}else{
		int q=t[p].ch[c];
		if(t[q].len==t[p].len+1){
			t[now].link=q;
		}else{
			int cln=++nc;
			t[cln]=t[q];
			t[cln].len=t[p].len+1;
			while(~p&&t[p].ch[c]==q){
				t[p].ch[c]=cln;
				p=t[p].link;
			}
			t[q].link=t[now].link=cln;
		}
	}
	lst=now;
}
vector<int> glink[N<<1];
int lin[N],con[N];
int qval(int l,int r){
	init_sam();
	for(int i=l;i<=r;i++){
		ext(a[i]);
		t[lst].min=t[lst].max=i;
		lin[i]=0;
		con[i]=0;
	}
	for(int i=0;i<=nc;i++){
		glink[i].clear();
	}
	for(int i=1;i<=nc;i++){
		glink[t[i].link].push_back(i);
	}
	function<void(int)> dfs=[&](int u){
		if(!t[u].min){
			t[u].min=1e9;
			t[u].max=0;
		}
		for(int v:glink[u]){
			dfs(v);
			t[u].min=min(t[u].min,t[v].min);
			t[u].max=max(t[u].max,t[v].max);
		}
	}; dfs(0);
//	for(int i=0;i<=nc;i++){
//		printf("node #%lld\n",i);
//		printf("  len: 	%lld\n",t[i].len);
//		printf("  link:	%lld\n",t[i].link);
//		printf("  rng: 	[%lld %lld]\n",t[i].min,t[i].max);
//		for(int j=0;j<26;j++) if(t[i].ch[j]){
//			printf("  %c -> %lld\n",j+'a',t[i].ch[j]);
//		}
//		puts("");
//	}
	// note:
	// (IF EXIST) [min,max-r-1] += r-l+1
	// (IF EXIST) [MAX(min,max-r),max-l] += (max-l) - i + 1
	for(int i=1;i<=nc;i++){
		int minn=t[i].min,maxn=t[i].max;
		int L=t[t[i].link].len+1,R=t[i].len;
//		printf("%d %d %d %d\n",minn,maxn,L,R);
		int l1=minn,r1=maxn-R-1,lin1=0,con1=R-L+1;
		int l2=max(minn,maxn-R),r2=maxn-L,lin2=-1,con2=maxn-L+1;
//		printf("[%d %d] += %dx%+d\n",l1,r1,lin1,con1);
//		printf("[%d %d] += %dx%+d\n",l2,r2,lin2,con2);
		if(l1<=r1){
//			printf("[%d %d] += %dx%+d\n",l1,r1,lin1,con1);
			lin[l1]+=lin1;
			lin[r1+1]-=lin1;
			con[l1]+=con1;
			con[r1+1]-=con1;
		}
		if(l2<=r2){
//			printf("[%d %d] += %dx%+d\n",l2,r2,lin2,con2);
			lin[l2]+=lin2;
			lin[r2+1]-=lin2;
			con[l2]+=con2;
			con[r2+1]-=con2;
		}
	}
	int ans=0;
	for(int i=l;i<r;i++){
		if(i>l){
			lin[i]+=lin[i-1];
			con[i]+=con[i-1];
		}
		int cur=lin[i]*i+con[i];
		ans+=cur*cur;
	}
	return ans;
}
bool chk(int x){
	int cur=1,cnt=0;
	while(true){
		int len=1;
		auto nxt=[&](int l){return min(2*l,n-cur+1);};
		while(len<n-cur+1&&qval(cur,cur+nxt(len)-1)<=x) len=nxt(len);
		int l=len,r=nxt(len),mid,ans=len;
		while(l<=r){
			mid=(l+r)>>1;
			if(qval(cur,cur+mid-1)<=x){
				ans=mid;
				l=mid+1;
			}else{
				r=mid-1;
			}
		}
		++cnt;
		if(cnt>k) return false;
		cur=cur+ans;
		if(cur==n+1){
			return true;
		}
	}
}
signed main(){
	scanf("%lld%lld%s",&n,&k,a+1);
	int l=0,r=1e18,mid,ans=0;
	while(l<=r){
		mid=(l+r)>>1;
		if(chk(mid)){
			ans=mid;
			r=mid-1;
		}else{
			l=mid+1;
		}
	}
	printf("%lld\n",ans);
}