蒟蒻求助 50分

#include<cmath>
#include<cstdio>
#include<iostream>
using namespace std;
int n;
const int inf=0x7f;
double dis[200][200],d[200],t1,t2=inf;
struct node{
    int x,y;
}a[200];
int chw(int i,int j){
    return sqrt((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y));
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d %d",&a[i].x,&a[i].y);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++){
            int b;
            scanf("%1d",&b);
            if(b)
                dis[i][j]=chw(i,j);
            else if(i!=j)
                dis[i][j]=inf;
        }
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(dis[i][j]>dis[i][k]+dis[k][j])
                    dis[i][j]=dis[i][k]+dis[k][j];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++){
            if(dis[i][j]!=inf)
                d[i]=max(dis[i][j],d[i]);
            t1=max(d[i],t1);
        }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(dis[i][j]==inf)
                t2=min(t2,chw(i,j)+d[i]+d[j]);
    printf("%.6lf",max(t1,t2));
    return 0;
}