R=range
J=lambda:map(int,input().split())
N,K=J()
v=[list(J())for i in R(N)]+[[0,0]]
M=N+1
w=[[9**20]*M for i in R(M)]
for i in R(M):
	for j in R(i,M):w[i][j]=w[j][i]=((v[i][0]-v[j][0])**2+(v[i][1]-v[j][1])**2)**0.5
S=R(1<<N)
D=[[[9**30]*(K+1)for i in R(M)]for p in S]
D[0][N][K]=0
for p in S:
	d=D[p]
	for i in R(N):
		if(p>>i)&1:
			q=p^(1<<i);e=D[q];d[i][K-1]=e[N][K]+w[N][i]
			for j in R(N):
				if(q>>j)&1:
					for k in R(K):d[i][k]=min(d[i][k],e[j][k+1]+w[j][i])
	for j in R(N):
		if(p>>j)&1:
			for k in R(K):d[N][K]=min(d[N][K],d[j][k]+w[j][N])
print(D[(1<<N)-1][N][K])