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=[[0]*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:
	for i in R(N):
		if(p>>i)&1:
			for j in R(N):
				if(p>>j)&1<1:
					for k in R(K):D[p|(1<<j)][j][k]=min(D[p|(1<<j)][j][k],D[p][i][k+1]+w[i][j])
			for k in R(K):D[p][N][K]=min(D[p][N][K],D[p][i][k]+w[i][N])
		else:D[p|(1<<i)][i][K-1]=D[p][N][K]+w[N][i]
print(D[(1<<N)-1][N][K])