POJ 2728 JZYZOJ 1636 分数规划 最小生成树 二分 prim

http://172.20.6.3/Problem_Show.asp?id=1636

复习了prim，分数规划大概就是把一个求最小值或最大值的分式移项变成一个可二分求解的式子。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
const int maxn=1100;
const double eps=0.0001;
int n;
long long c[maxn][3]={};
double e[maxn][maxn]={};
double t[maxn][maxn]={};
double w[maxn][maxn]={},k[maxn]={};
bool vis[maxn]={};
inline long long mabs(long long x){ return x>0?x:-x; }
inline long long sqr(long long x){ return x*x; }
bool ke(double v){
    memset(vis,0,sizeof(vis));
    for(int i=1;i<=n;i++){
        for(int j=i;j<=n;j++){
            w[i][j]=w[j][i]=t[i][j]-v*e[i][j];
        }
    }vis[1]=1;double fla=0;
    for(int i=1;i<=n;i++){
        k[i]=w[1][i];
    }
    for(int i=1;i<n;++i){
        int u=-1;
        for(int j=1;j<=n;++j){
            if(!vis[j])
                if(u==-1||k[j]<k[u])u=j;
        }
        fla+=k[u];vis[u]=1;
        for(int j=1;j<=n;++j){
            if(!vis[j])
                if(w[u][j]<k[j])k[j]=w[u][j];
        }
    }
    return fla<=0;
}
double doit(double l,double r){
    while(r-l>eps){
        double mid=(l+r)/2;
        if(ke(mid))r=mid;
        else l=mid+eps;
    }
    return l;
}
int main(){
    //freopen("wtf.in","r",stdin);
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%lld%lld%lld",&c[i][0],&c[i][1],&c[i][2]);
    }
    for(int i=1;i<=n;i++){
        for(int j=i;j<=n;j++){
            e[i][j]=sqrt((double)(sqr(c[i][0]-c[j][0])+sqr(c[i][1]-c[j][1])));
            t[i][j]=(double)mabs(c[i][2]-c[j][2]);
            e[j][i]=e[i][j];t[j][i]=t[i][j];
        }
    }printf("%.3f",doit(1,100));
    return 0;
}