#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
int n,m;
double l[21],a[21],b[21];
int data[21];
double res;

void dfs(int v,int nm){
	if(v==n){
		double val=0.0;
		for(int i=0;i<nm;i++){
			for(int j=i+1;j<nm;j++){
				val+=pow(l[data[i]]-l[data[j]],2)+pow(a[data[i]]-a[data[j]],2)+pow(b[data[i]]-b[data[j]],2);
			}
		}
		res=max(res,val);
		return;
	}
	dfs(v+1,nm);
	if(nm<m){
		data[nm]=v;
		dfs(v+1,nm+1);
	}
}

int main(void){
	scanf("%d %d",&n,&m);
	for(int i=0;i<n;i++){
		scanf("%lf %lf %lf",&l[i],&a[i],&b[i]);
	}
	dfs(0,0);
	printf("%.20f\n",res);
	return 0;
}