#include<bits/stdc++.h>
using namespace std;
typedef long double ld;
const int N=1e6+99;
int n;
ld d;
struct node{
    ld x,y;
}p[N];
ld dist(node A,node B){
    return sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y));
}
ld length(node A){
    return sqrt(A.x*A.x+A.y*A.y);
}
void rotate(node&A,ld alpha){
    ld len=length(A);
    if(len==0) return;
    ld c=A.x/len,s=A.y/len;
    ld cc=c*cos(alpha)+s*sin(alpha);
    ld ss=s*cos(alpha)-c*sin(alpha);
    A=(node){cc*len,ss*len};
}
bool cmp(node A,node B){
    return A.x<B.x||(A.x==B.x&&A.y<B.y);
}
bool cmpp(int A,int B){
    return p[A].y<p[B].y;
}
int q[N],cnt;
void merge(int l,int r){
    if(l==r) return;
    int mid=((l+r)>>1);
    merge(l,mid); merge(mid+1,r);
    cnt=0;
    for(int i=l;i<=r;i++)
        if(fabs(p[i].x-p[mid].x)<d/2)
            q[++cnt]=i;
    sort(q+1,q+cnt+1,cmpp);
    for(int i=1;i<=cnt;i++)
        for(int j=i+1;j<=cnt&&p[q[j]].y-p[q[i]].y<d/2;j++)
            for(int k=j+1;k<=cnt&&p[q[k]].y-p[q[i]].y<d/2;k++)
                d=min(d,dist(p[q[i]],p[q[j]]) + dist(p[q[j]],p[q[k]]) + dist(p[q[k]],p[q[i]]));
    return;
}
int main(){
//    freopen("triangle.in","r",stdin);
//    freopen("triangle.out","w",stdout);
    std::ios::sync_with_stdio(false);
    cin.tie(0);
    int T; cin>>T;
    int tt=0;
    while(T--){
        tt++;
        cin>>n;
        int alpha=abs(rand());
        for(int i=1;i<=n;i++){
            cin>>p[i].x>>p[i].y;
            rotate(p[i],alpha);
        }
        sort(p+1,p+n+1,cmp);
        d=dist(p[1],p[2])+dist(p[2],p[3])+dist(p[3],p[1]);
        merge(1,n);
        cout<<"Case #"<<tt<<": ";
        cout<<fixed<<setprecision(6)<<d<<'\n';
    }
    return 0;
}