HDU 3756 Dome of Circus

不会做，参见别人的程序：

/* 
底面为xy平面和轴为z轴的圆锥，给定一些点，使得圆锥覆盖所有点并且体积最小 
点都可以投射到xz平面，问题转换为确定一条直线（交x，z与正半轴）使得与x的截距r 
和与z轴的截距h满足h*r*r最小。 
三分，对于确定的h可以找到最佳的r，并且h*r*r的曲线必定只有一个极小值 
*/  
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
struct po  
{  
    double x,y;  
}p[10005];  
const double eps=1e-9;  
double Y;  
int n;  
double makeR(double h)  
{  
    double R=0;  
    for(int i=0;i<n;i++)  
    {   
       if(R*(h-p[i].y)<h*p[i].x)  
       R=h*p[i].x/(h-p[i].y);  
    }  
    return R;  
}  
void solve()  
{  
    double L=Y,R=1<<30,tmp;  
    while(R-L>eps)  
    {  
         tmp=(R-L)/3.0;  
        double mid1=L+tmp;  
        double mid2=L+2.0*tmp;  
        double R1=makeR(mid1);  
        double R2=makeR(mid2);    
        if(R1*R1*mid1>R2*R2*mid2)  
        {  
            L=mid1;  
        }  
        else  
        {  
            R=mid2;  
        }  
    }  
    printf("%.3lf %.3lf
",L+eps,makeR(L)+eps);  
}  
int main()  
{  
    int i,j,k;  
    int ca;  
    scanf("%d",&ca);  
    while(ca--)  
    {  
        scanf("%d",&n);  
        double tx,ty,Left=0;  
        Y=0;  
        for(i=0;i<n;i++)  
        {  
            scanf("%lf%lf%lf",&tx,&ty,&p[i].y);  
            p[i].x=sqrt(tx*tx+ty*ty);  
            Y=max(p[i].y,Y);  
        }  
        solve();  
    }  
    return 0;  
}