【SSL 1119】【USACO 2.4 】牛的旅行【最短路 FLoyed】

Description

农民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

分析&说明:

这道题也是用Floyed算法做的,题目给的数据是150*150,三重for肯定可以过。思路:Floyd求一遍每个点的最长直径,进而求出整个农场的最长直径。再依次看看哪些点之间没有道路,一个一个地尝试即可。

善意的提醒

某些变量、函数记得用double型
极小值、极大值的设置
Floyd中间转移点的循环必须放在最外层

不多说,直接CODE:

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
double f[1001][1001],ma[1001],maxx=0,ans,maxn=2147483647;
int n,x,y,m,s,e,a[10001][3],hyy;
string hy[1001];
int main() {
    cin>>n;
    for(int i=1;i<=n;i++) cin>>a[i][1]>>a[i][2];
    for(int i=1;i<=n;i++){
        cin>>hy[i];
    }
    for(int i=1;i<=n;i++){
        for(int j=0;j<=n;j++){
            hyy=hy[i][j-1]-'0';  //注意转数字
            if(hyy) f[i][j]=sqrt((a[i][1]-a[j][1])*(a[i][1]-a[j][1])+(a[i][2]-a[j][2])*(a[i][2]-a[j][2]));
            else if(i!=j) f[i][j]=2147483647;
        }
    }
    for(int k=1;k<=n;k++)                    
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(f[i][k]+f[k][j]<f[i][j])
                    f[i][j]=f[i][k]+f[k][j];
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(f[i][j]!=2147483647){  //如果他们之间有道路 
                ma[i]=max(f[i][j],ma[i]);//判断是否为i的最长直径 
                maxx=max(maxx,ma[i]);//再判断它是不是农场的最大直径 
            }
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(f[i][j]==2147483647){
                maxn=min(ma[i]+sqrt((a[i][1]-a[j][1])*(a[i][1]-a[j][1])+(a[i][2]-a[j][2])*(a[i][2]-a[j][2]))+ma[j],maxn);  //依次尝试,试出之后农场的最大直径 
            }
        }
    }
    ans=max(maxx,maxn); //由于农场之后的直径有可能不是最大的,所以要再次比较 
    printf("%.6f",ans);
    return 0;
}
发布了37 篇原创文章 · 获赞 24 · 访问量 897

猜你喜欢

转载自blog.csdn.net/dgssl_xhy/article/details/103743018