分式规划问题（二分）2

TZOJ：8532

题目大意是给定n个点的二维坐标，每个点都有一个高度h值，将这n个点通过n-1条边相连，并求和每条边相连的两个点之间的水平距离X，和两点的高度差H，最终要使得H/X最小。

首先，分式规划问题的思路一般都是二分，所以先把二分模板写出来，之后就是判断在当前二分答案是否可行，我们可以将每个点之间的权值当成他们之间的距离，那么我们只需要用prime最小生成树来生成一颗n-1条边的树，就可以得到总距离的最小值，从而判断该二分答案是否可行

#include<bits/stdc++.h>
using namespace std;
const int N=1e3+10;
struct d
{
	int x;
	int y;
	int h;
}s[N];
double a[N][N];
int n;

double han(double x,double y,double x1,double y1)
{
	double p=x-x1,q=y-y1;
	return sqrt(q*q+p*p);
}
int check(double x)
{
	for(int i=1;i<=n;i++)a[i][i]=1e12;
	for(int i=1;i<=n;i++)
		for(int j=i+1;j<=n;j++)a[i][j]=a[j][i]=abs(s[i].h-s[j].h)-x*han(s[i].x,s[i].y,s[j].x,s[j].y);
	double sum=0,bj[N],dis[N];
	memset(bj,0,sizeof(bj));
	for(int i=1;i<=n;i++)dis[i]=1e12;
	dis[1]=0;
	for(int i=1;i<=n;i++)
	{
		double Min=1e12;
		int t;
		for(int j=1;j<=n;j++)
		{
			if(!bj[j]&&dis[j]<Min)
			{
				Min=dis[j];
				t=j;
			}
		}
		bj[t]=1;
		sum+=dis[t];
		for(int j=1;j<=n;j++)
			if(!bj[j]&&dis[j]>a[t][j])dis[j]=a[t][j];
	}
//	cout<<sum<<endl;
	return sum<=0;
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	while(cin>>n&&n)
	{
		for(int i=1;i<=n;i++)cin>>s[i].x>>s[i].y>>s[i].h;
		double l=0,r=1e7;
		while(r-l>0.00001)
		{
		
			double mid=(l+r)/2;	
			if(check(mid))r=mid;
			else l=mid;
		}
		cout<<setiosflags(ios::fixed);
		cout<<setprecision(3)<<r<<endl;
	}
}

要注意当前二分权值为（两村庄之间的高度差）-x*（两村庄之间的水平距离）。

如果在水平距离为0的时候，其高度差最大值为1e7，也就是说最大的结果也会是1e7，最小的结果就是0。

在初始化到各个村庄的最短距离dis时，dis中能出现的最大值是1e7*1e4，所以初始化要为1e12。不过实际好像1e6也能过，有可能是数据比较水，也可能是这种情况应该是不会出现的。不过赋值到数据最大值总归是好的