洛谷 P1522 牛的旅行 Cow Tours

//P1522 牛的旅行 Cow Tours
//首先说明一点，翻译得很不好，还是看原版的比较好。
//http://blog.csdn.net/supersnow0622/article/details/9763979此文写得不错，摘抄如下：
//分析：首先要连在一起修路的这两个牧区，之前是不连通的。其次就是这两个牧区修路之后，
//相隔最远的两个牧区的距离最小。其实，只要枚举任意不连通的两个牧区在计算这种情况下的最远的距离就好了，
//然后在这些最远的距离中找到距离最小（题目要求半径最小）的一个就是所求了。
//关键是在两个不连通的牧区修路之后求牧场的半径呢？假设已将one牧区种的D和two牧区中的F相连，
//想一下，这个牧场的半径是不是等于与D相距最远的点到D的距离+D与F的距离+与F相距最远的点的距离呢？
//这样的求法是满足题目所说的牧场的半径是相距最远的两个点的距离。
//所以思路很快就出来了：
//1.用floyed计算任意两点间的距离
//2.计算每个点与其他点的最远距离
//3.计算任意两个不连通的点修路之后牧场的半径。
//http://www.xuebuyuan.com/1457314.html此文给出了英文原版题目，读下来，配合中文意思，弄懂题目了。
//http://www.nocow.cn/index.php/USACO/cowtour根据该文，弄懂了，为什么
//那我们只需要计算新连成的区块的直径，然后与原来的所有区块的直径放在一起，取最大直径就可以了
//突然对为什么取最大直径，又心存疑虑，想了3个小时，想不通，暂时搁置吧。
//老外的题，与中国的题，确实存在着文化差异。感叹。
//编码的过程中，真是服了，洛谷的样例输入数据有问题。查找了一通，发现没有问题，数据格式本就是如此。折腾了一场。
#include <stdio.h>
#include <string.h>
#include <math.h>//洛谷中，该头文件，必须以C++形式提交。
#define INF 1000000
int x[150+10],y[150+10];
double a[150+10][150+10],d[150+10],minD=INF;
char s[200];
double dis(int x1,int y1,int x2,int y2){
    return sqrt((double)(x1-x2)*(x1-x2)+(double)(y1-y2)*(y1-y2));//别忘了强转,否者容易溢出100000*100000
}
int main(){
    int i,j,k,n,w;
    memset(d,0,sizeof(d));
    scanf("%d",&n);
    for(i=1;i<=n;i++)
        scanf("%d%d",&x[i],&y[i]);
    for(i=1;i<=n;i++){
        scanf("%s",s+1);
        for(j=1;j<=n;j++){
            if(i!=j){
                if(s[j]=='0')
                    a[i][j]=a[i][j]=INF;
                else{
                    a[i][j]=a[j][i]=dis(x[i],y[i],x[j],y[j]);
                }
            }else
                a[i][j]=0;
        }
    }
    for(k=1;k<=n;k++)//floyd算法
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
                if(a[i][j]>a[i][k]+a[k][j])
                    a[i][j]=a[i][k]+a[k][j];
    for(i=1;i<=n;i++)//连通区块里找最大值
        for(j=1;j<=n;j++)
            if(a[i][j]<INF&&a[i][j]>d[i])
                d[i]=a[i][j];                
    for(i=1;i<=n;i++)//连接两个牧区，找最小值。
        for(j=1;j<=n;j++)
            if(a[i][j]==INF&&minD>d[i]+d[j]+dis(x[i],y[i],x[j],y[j]))
                minD=d[i]+d[j]+dis(x[i],y[i],x[j],y[j]);
    for(i=1;i<=n;i++)//此处想不明白，为什么找最大值。 等水平高了再回看此题。2017-6-20 21:23
        if(minD<d[i])
            minD=d[i];
    printf("%.6lf\n",minD);//1此处写成printf("%d\n",minD);低级错误
    return 0;
}