usaco Training Cow Tours

农民John 的农场里有很多牧区.有的路径连接一些特定的牧区.一片所有连通的牧区称为一个牧场.
但是就目前而言,你能看到有两个牧区不连通.这样,农民John 就有多个牧区了. John 想在农
场里添加一条路径(注意,恰好一条).对这条路径有以下限制:
一个牧场的直径就是牧场中最远的两个牧区的距离(本题中所提到的所有距离指的都是最短的距
离).考虑如下的有5 个牧区的牧场,牧区用“*”表示,路径用直线表示.每一个牧区都有自己的坐

标:

这个牧场的直径大约是12.07106, 最远的两个牧区是A 和E,它们之间的最短路径是A-B-E.
这里是另一个牧场:

这两个牧场都在John 的农场上.John 将会在两个牧场中各选一个牧区,然后用一条路径连起来,使得连通后这个新的更大的牧场有最小的直径.
注意,如果两条路径中途相交,我们不认为它们是连通的.只有两条路径在同一个牧区相交,我们才认为它们是连通的.
输入文件包括牧区、它们各自的坐标,还有一个如下的对称邻接矩阵:
A B C D E F G H
A 0 1 0 0 0 0 0 0
B 1 0 1 1 1 0 0 0
C 0 1 0 0 1 0 0 0
D 0 1 0 0 1 0 0 0
E 0 1 1 1 0 0 0 0
F 0 0 0 0 0 0 1 0
G 0 0 0 0 0 1 0 1
H 0 0 0 0 0 0 1 0
输入包括两个不连通的牧区.
请编程找出一条连接这两个不同牧场的路径,使得连上这条路径后,这个更大的新牧场有最小的直径.

Input

第1 行: 一个整数N (1 <= N <= 150), 表示牧区数
第2 到N+1 行: 每行两个整数X,Y (0 <= X ,Y<= 100000), 表示N 个牧区的坐标.注意每个 牧区
的坐标都是不一样的.
第N+2 行到第2*N+1 行: 每行包括N 个数字(0 或1) 表示如上文描述的对称邻接矩阵.

Output

只有一行,包括一个实数,表示所求答案.数字保留六位小数.

Sample Input

8

10 10

15 10

20 10

15 15

20 15

30 15

25 10

30 10

01000000

10111000

01001000

01001000

01110000

00000010

00000101

00000010

Sample Output

22.071068

这道题主要是floyd，但需要预处理，因为牛只有2拨，所以用mia[i]存最小值即可n^2查询,floyd n^3,可以ac。

以下为代码：

#include<stdio.h>
#include<string.h>
#include<math.h>
int a[151][151];
double f[151][151];
double x[151],y[151];
double mia[151];
double distance(int i,int j)
{
    return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}
double max(double x,double y)
{
    if(x>y)
        return x;
    return y;
}
double min(double x,double y)
{
    if(x>y)
        return y;
    return x;
}
int main()
{
    int n,i,j,k;
	double ans=1e12;
    scanf("%d",&n);
    for(i=1;i<=n;i++)
        scanf("%lf%lf",&x[i],&y[i]);
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            scanf("%1d",&a[i][j]);
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
        {
        	if(a[i][j]==1)
        		f[i][j]=distance(i,j); 
			else
				f[i][j]=1e12;      	
        }
    for(k=1;k<=n;k++)
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
            	if(i!=j&&i!=k&&j!=k&&f[i][k]<1e12&&f[k][j]<1e12)
                	f[i][j]=min(f[i][k]+f[k][j],f[i][j]);
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            if(f[i][j]>mia[i]&&f[i][j]<1e12)
                mia[i]=f[i][j];
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            if(f[i][j]==1e12&&i!=j)
            {
            	double tmp=distance(i,j);
                ans=min(ans,mia[i]+mia[j]+tmp);
            }
    for(i=1;i<=n;i++)
    {
    	if(mia[i]>ans)
    		ans=mia[i];
    }
    printf("%.6lf",ans);                
}