pku 2728 最优比率生成树

题目描述：

给你n个水库，告诉你每个水库的xy坐标（米）和这个水库的花费。

如果两个水库需要连接，连接的费用是这两个水库的高度差的绝对值。

让你连接一些水库，让所有的水库都能到达1号水库并且只有一条通路（就是一棵树）。

然你使每米的花费最少（就是总花费 比 总长度）

解题报告：实际上就求一棵最小比率生成树，黑书上介绍的挺详细的，就是在代码实现上有点困难

有俩种方法实现，一种是二分法，不过在求上下界方面有点难，并且慢很多，一种是用迭代法，标准的01分数规划，不断的逼近迭代最优解；

关于构造式子以及用二分法实现，这个博客上解答的很详细

http://hi.baidu.com/zzningxp/blog/item/b2d1b4ec1f8bbc2262d09fc9.html

#include<stdio.h>
#include<math.h>
#include<string>
#define   Maxvalue 99999999
int n,x[1000],y[1000],h[1000];
double cost[1000][1000],dist[1000][1000],dis[1000];
double prim(double p);
bool flag[1000];
int main()
{
    int i,j;
    double a,b;
    while(scanf("%d",&n)==1&&n)
    {
          for (i=0;i<n;i++)
                scanf("%d%d%d",&x[i],&y[i],&h[i]);
          for (i=0;i<n;i++)
          {
                for (j=i+1;j<n;j++)
                {
                      b=(x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]);
                      dist[i][j]=sqrt(b);
                      dist[j][i]=dist[i][j];
                      cost[i][j]=h[i]-h[j];
                      if (cost[i][j]<0)
                            cost[i][j]=-cost[i][j];
                      cost[j][i]=cost[i][j];
                }
          }
          a=0;
          while (1)
          {
                b=prim(a);
                if (fabs(b-a)<0.0001)
                      break;
                a=b;
          }
          printf("%.3lf\n",b);
    }
    return 0;
}
double prim(double p)
{
    int i,j,k,pre[1000]={0};
    double mincost,v,totcost=0,totdist=0;
    dis[0]=0;
	memset(flag,false,sizeof(flag));
    for (j=1;j<n;j++)
          dis[j]=cost[j][0]-p*dist[j][0];
    for (i=1;i<n;i++)
    {
          mincost=Maxvalue;
          for (j=1;j<n;j++)
          {
                if (flag[j]==false && mincost>dis[j])
                {
                      mincost=dis[j];
                      k=j;
                }
          }
          flag[k]=true;
          totcost+=cost[pre[k]][k];
          totdist+=dist[pre[k]][k];
          for (j=0;j<n;j++)
          {
                if (flag[j]==false)
                {
                      if ((v=cost[k][j]-p*dist[k][j])<dis[j])
                      {
                            dis[j]=v;
                            pre[j]=k;
                      }
                }
          }
    }
    return totcost/totdist;
}