1. 程式人生 > >51 Nod 1629 B君的圓錐

51 Nod 1629 B君的圓錐

基準時間限制:1 秒 空間限制:131072 KB 分值: 10 難度:2級演算法題

 收藏

 關注

B君要用一個表面積為S的圓錐將白山雲包起來。

B君希望包住的白山雲體積儘量大,B君想知道體積最大可以是多少。

注意圓錐的表面積包括底面和側面。

Input

一行一個整數,表示表面積S。(1 <= S <= 10^9)

Output

一行一個實數,表示體積。

Input示例

8

Output示例

1.504506

思路:用三分法求最大值即可。


#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
#include <cmath>
#define PI (3.1415926535)
//#define PI acos(-1.0)
#define EPS (1e-7)
using namespace std;
double S; //圓錐的表面積
double f(double r) //已知圓錐底面半徑r,計算當前圓錐的體積
{     double R=(S-PI*r*r)/PI/r; //側面展開的扇形半徑
      double h=sqrt(R*R-r*r); //
      return PI*r*r*h/3;
}
double three_devide(double low,double up)
{
    double m1,m2;
    while(up-low>=EPS)
    {
        m1=low+(up-low)/3;
        m2=up-(up-low)/3;
        if(f(m1)<=f(m2))
            low=m1;
        else
            up=m2;
    }
    return (m1+m2)/2;
}

int main()
{
#ifndef ONLINE_JUDGE
    //freopen("in.txt","r",stdin);
#endif // ONLINE_JUDGE
    scanf("%lf",&S);
    double ans=f(three_devide(0,S+100));
    printf("%f\n",ans);
    return 0;
}