WA on #7

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int INF=1e15;
const int N=150+5;
double f[N][N],ma[N];
int x[N],y[N],n,v[N];
bool a[N][N];
queue<int>q;
double dis(int i,int j){
	return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}
signed main(){
	scanf("%lld",&n);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++) f[i][j]=INF;
		f[i][i]=0;
	}
	for(int i=1;i<=n;i++) scanf("%lld%lld",&x[i],&y[i]);
	char c;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			scanf("%c",&c);
			if(c=='0') a[i][j]=0;
			else if(c=='1') a[i][j]=1,f[i][j]=dis(i,j);
			else j--;
		}
	}
	int cnt=0;
	for(int i=1;i<=n;i++){
		if(v[i]!=0) continue;
		q.push(i);
		v[i]=++cnt;
		while(!q.empty()){
			int k=q.front();
			q.pop();
			for(int j=1;j<=n;j++){
				if(a[j][k]&&!v[j]){
					q.push(j);
					v[j]=cnt;
				}
			}
		}
	}
	for(int k=1;k<=n;k++) for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
	for(int i=1;i<=n;i++){
		ma[i]=0;
		for(int j=1;j<=n;j++) if(f[i][j]!=INF) ma[i]=max(ma[i],f[i][j]);
	}
	double ans=INF;
	for(int i=1;i<=n;i++){
		for(int j=i+1;j<=n;j++){
			if(v[i]!=v[j]) ans=min(ans,dis(i,j)+ma[i]+ma[j]);
		}
	}
	printf("%.6lf",ans);
	return 0;
}