hdu 3756 Dome of Circus 三分

来源：互联网发布：数据分析技术发展趋势编辑：程序博客网时间：2024/05/09 06:10

求能包围 n个点的最小体积的圆锥的高和底面半径。

对于一个点而言，通过分析体积公式。求导后，公式是单点的。所以体积公式是先下降后上升（推导和分析好一会）。

而对于n个点， Fi （h）为高为h 时最小体积。

对于高 h，得到最小 V = max( Fi (h)) i 为 1--n

分析的 n 个对钩函数的max 也是对钩函数（这里对钩函数是先下降在上升的函数即可）。

因此对于n个点的最小包围圆锥的体积对于高h 而言就是一个先下降在上升。这样就可以三分。

#include <vector>#include <list>#include <map>#include <set>#include <deque>#include <stack>#include <cstring>#include <bitset>#include <algorithm>#include <functional>#include <numeric>#include <utility>#include <sstream>#include <iostream>#include <iomanip>#include <cstdio>#include <cmath>#include <cstdlib>#include <ctime>#include <assert.h>#include <queue>#define REP(i,n) for(int i=0;i<n;i++)#define TR(i,x) for(typeof(x.begin()) i=x.begin();i!=x.end();i++)#define ALLL(x) x.begin(),x.end()#define SORT(x) sort(ALLL(x))#define CLEAR(x) memset(x,0,sizeof(x))#define FILLL(x,c) memset(x,c,sizeof(x))using namespace std;const double eps = 1e-9;#define LL long long #define pb push_backconst int maxn = 10100;const double PI = 3.1415926;double x[maxn],y[maxn],z[maxn];int n ;double f(double a){    double ret = 0;    for(int i=1;i<=n;i++){         double l = sqrt(x[i]*x[i]+y[i]*y[i]);         double h = z[i];                 double r = h*l/(a-h)+l;        //  cout << l << "  "<< h<< "  "<<r <<endl;         ret  = max(ret , r);    }    return ret;}double s(double a){      double  r = f(a);      // cout << a << "  "<< r<<endl;      return PI*r*r *a/3;}double down = 0 ;void solve(){    //cout << f(3)<<"fasdfas"<<endl;      double left = down + 1e-6;      double right = 1e7;      while(right - left >1e-4){      //    cout << left << "  "<<right <<endl;             double lmid = (left + right)/2;             double rmid = (lmid + right )/2;             //cout << lmid << "**** "<<rmid<<endl;             if(s(lmid)>s(rmid)){                    left = lmid ;             }else{                    right = rmid;             }      }      double h = left ;      double r = f(left);      printf("%.3lf %.3lf\n",h,r);}int main(){   int t ;   cin >>t ;   while(t--){          scanf("%d",&n);          down = 0 ;          for(int i=1;i<=n;i++){                 scanf("%lf%lf%lf",&x[i],&y[i],&z[i]);                 down = max(down , z[i]);          }                  solve();   }    return 0;}