#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
using namespace std;
using ll=long long;
using ld=long double;

void IO(){
  ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  cout<<fixed<<setprecision(10);
}

ll n,K;
vector<ld> x,y;

ld d1(ll i,ll j){
  ld n2=(x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]);
  return sqrt(n2);
}

ld d2(ll i){
  ld n2=x[i]*x[i]+y[i]*y[i];
  return sqrt(n2);
}

int main(){
  IO();
  cin>>n>>K;
  x.resize(n);
  y.resize(n);
  for(ll i=0;i<n;i++){
    cin>>x[i]>>y[i];
  }
  vector<vector<vector<ld>>> dp(1<<n,vector<vector<ld>>(K+1,vector<ld>(n,1e18)));
  for(ll i=0;i<n;i++){
    dp[1<<i][K-1][i]=d2(i);
  }
  for(ll b=1;b<(1<<n);b++){
    for(ll i=0;i<n;i++){
      for(ll j=0;j<n;j++){
        if(b!=0&&!(b&(1<<j))){
          continue;
        }
        if((b&(1<<i))==0){
          if(i!=j){
            for(ll k=0;k<=K;k++){
              if(k>0){
                dp[b|(1<<i)][k-1][i]=min(dp[b|(1<<i)][k-1][i],dp[b][k][j]+d1(j,i));
              }
              dp[b|(1<<i)][K][i]=min(dp[b|(1<<i)][K][i],dp[b][k][j]+d2(j)+d2(i));
            }
          }
        }
      }
    }
  }
  ld ans=1e18;
  for(ll i=0;i<n;i++){
    for(ll j=0;j<=K;j++){
      ans=min(ans,dp[(1<<n)-1][j][i]+d2(i));
    }
  }
  cout<<ans<<endl;
}