C++ :

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
double x[151],y[151],maxint=1e12,m[151],minx,r,temp,f[151][151];

double dist(int i,int j)
{
	return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}

int main()
{
	int i,j,n,k;char c;
	cin>>n;
	for(int i=1;i<=n;i++) cin>>x[i]>>y[i];
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	{
		cin>>c;
		if(c=='1') f[i][j]=dist(i,j);
		else f[i][j]=maxint;
	}
	for(int k=1;k<=n;k++)
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	{
		if(i!=j&&i!=k&&k!=j)//本步勿忘 
		if(f[i][k]<maxint-1&&f[k][j]<maxint-1)
		if(f[i][j]>f[i][k]+f[k][j])
		f[i][j]=f[i][k]+f[k][j];
	}
	 memset(m,0,sizeof(m));
	 for(int i=1;i<=n;i++)
	 for(int j=1;j<=n;j++)
	 if(f[i][j]<maxint-1&&m[i]<f[i][j]) m[i]=f[i][j];
	 minx=1e20;
	 for(int i=1;i<=n;i++)
	 for(int j=1;j<=n;j++)
	 {
	 	if(i!=j&&f[i][j]>maxint-1)
	 	{
	 		temp=dist(i,j);
	 		if(minx>m[i]+m[j]+temp) minx=m[i]+m[j]+temp;//使总距离最小（新直径） 
		 }
	  } 
	  r=0;
	  for(int i=1;i<=n;i++) if(m[i]>minx) minx=m[i];//不要忘记其他原先路程也存在，故他们的最小值也要扫一遍。 
	  printf("%.6lf",minx);
	  return 0;
}