牛的旅行 Cow Tours

P1522 牛的旅行 Cow Tours
题目概述：有的路径连接一些特定的牧区。一片所有连通的牧区称为一个牧场。但是就目前而言，你能看到至少有两个牧区通过任何路径都不连通。在牧场里添加一条路径(注意，恰好一条)。对这条路径有以下限制：一个牧场的直径就是牧场中最远的两个牧区的距离(本题中所提到的所有距离指的都是最短的距离)。请编程找出一条连接两个不同牧场的路径，使得连上这条路径后，这个更大的新牧场有最小的直径。输出在所有牧场中最小的可能的直径。
数据规模：N (1 <= N <= 150)，X,Y (0 <= X ,Y<= 100000)
思路：这个题首先要预处理距离，用勾股定理，然后考虑到N的值比较小（事实上小于300就可以了），我们考虑用O(N^3)的Floyd算法，然后如果暴力枚举一次就是O(N^2)*O(N^3)的时间，就算N再小也不能用这种算法，所以考虑是求一个最值，所以用最大值来记录，这样时间变为O(N^2)，就可以过了。
代码：

#include<iostream>#include<algorithm>#include<cstdio>#include<cstring>#include<cmath>using namespace std;int x[155],y[155];double a[155][155];int temp;int i,j,k,n,t,m;double maxx,minn=10000000000;int b[155];double c[155][155];int rr(){    int ans=0,f=1;    char ch=getchar();    while(ch<'0'||ch>'9')    {        if(ch=='-')            f=-1;        ch=getchar();    }    while(ch>='0'&&ch<='9')    {        ans*=10;        ans+=ch-'0';        ch=getchar();    }    return ans*f;}double dis(int xx,int yy){    return sqrt((x[xx]-x[yy])*(x[xx]-x[yy])    +(y[xx]-y[yy])*(y[xx]-y[yy]));}void solve(){    for(int i=1;i<=n;i++)    for(int j=1;j<=n;j++)    {        for(int k=1;k<=n;k++)        {        if(a[i][j]>a[i][k]+a[k][j])        a[i][j]=a[i][k]+a[k][j];        }    }    for(int i=1;i<=n;i++)    for(int j=1;j<=n;j++)    {        for(int k=1;k<=n;k++)        {        if(a[i][j]>a[i][k]+a[k][j])        a[i][j]=a[i][k]+a[k][j];        }    }    for(int i=1;i<=n;i++)    for(int j=1;j<=n;j++)    if(a[i][j]!=1000000000)    maxx=max(maxx,a[i][j]);    for(int i=1;i<=n;i++)    for(int j=1;j<=n;j++)    a[i][j]=c[i][j];}int main(){    n=rr();    for(i=1;i<=n;i++)    for(j=1;j<=n;j++)    a[i][j]=1000000000;    char ch;    for(i=1;i<=n;i++)    x[i]=rr(),y[i]=rr();    for(i=1;i<=n;i++)    for(j=1;j<=n;j++)    {        cin>>ch;        if(ch=='1')        a[i][j]=a[j][i]=dis(i,j);    }    for(i=1;i<=n;i++)    a[i][i]=0;    for(i=1;i<=n;i++)    for(j=1;j<=n;j++)    for(k=1;k<=n;k++)    {        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++)    c[i][j]=a[j][i];    for(int i=1;i<=n;i++)    for(int j=1;j<i;j++)    {        if(a[i][j]==1000000000)        {            maxx=0;            a[i][j]=a[j][i]=dis(i,j);            solve();            minn=min(minn,maxx);            a[i][j]=a[j][i]=1000000000;        }    }    printf("%.6f",minn);    return 0;}/**/