二分装水问题-算法笔记

#include<stdio.h>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<fstream>
using namespace std;
const double PI = acos(-1);
const double eps = 1e-5;
double R, r;
double f(double h, double R){
	double alpha = 2 *  acos((R - h) / R);//扇形所对应的角 
	double s = alpha * R * R / 2; // 扇形面积 
	double L = 2 * sqrt(R * R - (R - h) * (R - h));
	double s1 = s - L * (R - h) / 2;//扇形减去三角形 
	double s2 = PI * R * R / 2;
	return s1 / s2; 
}
double cal_h(){
	double left = 0, right = R, mid;
	while(right - left > eps){
		mid = (left + right) / 2;
		if(f(mid, R) > r){
			right = mid;
		}else{
			left = mid;
		}
		
	}
	return mid;
} 
int main() {
	cin >> R >> r;
	cout << cal_h();
	return 0;
}