codeforce_exercise_r35

P1522 [USACO2.4]牛的旅行 Cow Tours ：

问题描述

题目简述

Farmer John 的农场里有很多牧区。有的路径连接一些特定的牧区。一片所有连通的牧区称为一个牧场。但是就目前而言，你能看到至少有两个牧区通过任何路径都不连通。这样，Farmer John 就有多个牧场了。

John 想在牧场里添加恰好一条路径。对这条路径有以下限制：

一个牧场的直径就是牧场中最远的两个牧区的距离（本题中所提到的所有距离指的都是最短的距离）。

请编程找出一条连接两个不同牧场的路径，使得连上这条路径后，这个更大的新牧场有最小的直径。输出在所有牧场中最小的可能的直径。

注意，如果两条路径中途相交，我们不认为它们是连通的。只有两条路径在同一个牧区相交，我们才认为它们是连通的。

输入/输出格式

输入格式：
第一行一个整数 N（1≤N≤150），表示牧区数。

接下来 N行，每行两个整数 X,Y（0≤X,Y≤100000），表示 N 个牧区的坐标。注意每个牧区的坐标都是不一样的。

接下来 N 行，每行 N 个数字，代表邻接矩阵。

输出格式：
只有一行，包括一个实数，表示所求直径。数字保留六位小数。

只需要打到小数点后六位即可，不要做任何特别的四舍五入处理。

样例

输入样例：
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
输出样例：
22.071068

问题分析

解题思路

不知道是不是因为翻译的问题，读题的时候感觉非常的别扭。但实际上是求解最大直径的最小值。这里不同的一点是可以加一条边，这样处理后，直径其实还是很好算的，比如加入边（i,j），此时的直径为：i点的直径+weight(i,j)+j点的直径。有了上面的式子其实整个题就不难解决了，加上题目中给出N最大为150，可以使用floyd算法求出所有结点之间的最短路径。并求出每个结点到其它点的最大距离，这个就是该结点的直径。之后枚举任意两个不连通的结点，按照上式计算加边后的直径，将求出的结果的最小值与之前求出的结点直径的最大值比较，其较大者即为结果。

参考代码

include <iostream>
#include <cmath>
#include <vector>
#include <cstring>

using namespace std;

const int inf=0x3f3f3f3f;

class point
{
public:
	int p_x;
	int p_y;
    point(int xx,int yy)
	{
		p_x=xx;
		p_y=yy;
	} 
};

vector<point> points;
double dist[160][160];
double long_dist[160];
int n;

void init()
{
	memset(long_dist,0,sizeof(long_dist));
	for(int i=1;i<160;i++)
	{
		for(int j=1;j<160;j++)
		{
			dist[i][j]=inf;
		}
		dist[i][i]=0;
	}
	points.clear();
	point sp(-1,-1);
	points.push_back(sp);
}

void addedge(int i,int j)
{
	point p1=points[i];
	point p2=points[j];
	dist[i][j]=sqrt((p1.p_x-p2.p_x)*(p1.p_x-p2.p_x)+(p1.p_y-p2.p_y)*(p1.p_y-p2.p_y));
}

double dis(int i,int j)
{
	point p1=points[i];
	point p2=points[j];
	return sqrt((p1.p_x-p2.p_x)*(p1.p_x-p2.p_x)+(p1.p_y-p2.p_y)*(p1.p_y-p2.p_y));
}

void floyd()
{
	for(int k=1;k<=n;k++)
	{
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=n;j++)
			{
				dist[i][j]=min(dist[i][j],dist[i][k]+dist[k][j]);
			}
		}
	}
}

int main()
{
	init();
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		int x,y;
		scanf("%d %d",&x,&y);
		point tp(x,y);
		points.push_back(tp);
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			int have_edge;
			scanf("%1d",&have_edge);
			if(have_edge) addedge(i,j);
		}
	}
	floyd();
	double ans1=0;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if(dist[i][j]!=inf) 
			{
				ans1=max(ans1,dist[i][j]);
				long_dist[i]=max(long_dist[i],dist[i][j]);
			}
		}
	}
	double ans2=inf;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if(dist[i][j]==inf) ans2=min(ans2,long_dist[i]+dis(i,j)+long_dist[j]);
		}
	}
	printf("%.6f",max(ans1,ans2));
	return 0;
}

心得体会

稍微有点绕的题。一开始的时候WA了一个点，是因为没有和未加边之前的最大值比较导致的。这一点说实话，我也没太整明白。。。理论上ans2应该是大于ans1的，但是不知道为啥出错了。还得再想想。