#include<bits/stdc++.h>
typedef long long ll;
typedef long double ld;
using namespace std;
const int N=16,M=160,INF=0x3f3f3f3f;
inline int max(int x,int y){return x>y?x:y;}
inline int min(int x,int y){return x<y?x:y;}
inline void swap(int &x,int &y){x^=y^=x^=y;}
int n;
ld ans=INF,g[N][N],f[1<<12][N];
struct node{
	ld x,y;
	ld mdis(node a){
		return fabs(x-a.x)+fabs(y-a.y);
	}
	ld dist(node a){
		return sqrt((x-a.x)*(x-a.x)+(y-a.y)*(y-a.y));
	}
	node rotate(ld sin,ld cos){//clockwise
		return (node){x*cos+y*sin,-x*sin+y*cos};
	}
}p[N],w[N];
ld dp(){
	ld res=INF;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			g[i-1][j-1]=i==j?INF:w[i].mdis(w[j]);
	for(int S=0;S<1<<n;S++)
		for(int i=0;i<n;i++)f[S][i]=INF;
	for(int i=0;i<n;i++)f[1<<i][i]=0;
	for(int S=1;S<1<<n;S++)
		for(int i=0;i<n;i++)if(!((S>>i)&1))
			for(int j=0;j<n;j++)if((S>>j)&1)
				f[S|(1<<i)][i]=min(f[S|(1<<i)][i],f[S][j]+g[i][j]);
	for(int i=0;i<n;i++)
		res=min(res,f[(1<<n)-1][i]);
	return res;
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%Lf%Lf",&p[i].x,&p[i].y);
	for(int i=1;i<=n;i++)
		for(int j=i+1;j<=n;j++){
			ld r=p[i].dist(p[j]);
			ld x=p[i].x-p[j].x;
			ld y=p[i].y-p[j].y;
			for(int k=1;k<=n;k++)
				w[k]=p[k].rotate(y/r,x/r);
			ans=std::min(ans,dp());
		}
	printf("%.15Lf",ans);
	return 0;
}