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

51Nod - 1629 B君的圓錐

ace true href space input inpu 體積 line pri

51Nod - 1629 B君的圓錐

B君要用一個表面積為S的圓錐將白山雲包起來。 B君希望包住的白山雲體積盡量大,B君想知道體積最大可以是多少。 註意圓錐的表面積包括底面和側面。 Input
一行一個整數,表示表面積S。(1 <= S <= 10^9)
Output
一行一個實數,表示體積。
Input示例
8
Output示例
1.504506

題解:

    想使用三分法來解決,但是在9個test中只通過了第一個。(不知道出錯在哪裏了)

    後轉 直接公式法解決。

#include <iostream> 
#include <cstdio> 
#include <cmath> 
using namespace std; 
const double PI = 3.14159265358979; 

inline double fun(double x, const int &n){
	return PI*x*x*sqrt(n*n/(PI*PI*x*x) - 2*n/PI)/3;   
}

double solve(double l, double r, const int &n){
	double tmp_l, tmp_r; 
	while( (r-l) >= 1e-8 ){
		tmp_l = (r-l)/3.0 + l; 
		tmp_r = (r-l)*2.0/3.0 + l; 
		if( fun(tmp_r, n) > fun(tmp_l, n) ){
			l = tmp_l; 
		}else{
			r = tmp_r; 
		}
	}
	return fun(l + (r-l)/2.0, n); 
}

int main(){

	int n; 
	double l, r, ans; 
	while(scanf("%d", &n) != EOF){
// 		l = 0.0000001; 
// 		r = sqrt(n*1.0 / PI); 
// 		ans = solve(l, r, n); 
		ans = n*1.0*sqrt(n*1.0/(72.0*PI)); 
		printf("%.7f\n", ans );
	}
	return 0; 
}

  

51Nod - 1629 B君的圓錐