HDU 5954 Do not pour out2016ACM-ICPC亚洲区沈阳站

本文介绍了一种使用二分法求解特定几何问题的算法实现,即求解刚好能将水完全倒出的临界面积。通过分析椭圆投影到圆的面积变化规律,结合数学公式推导,最终实现了一个有效的解决方案。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题意不多说了,写的时候脑残了,没看懂题目,以为求最大的面积,。。。。。。。结果求刚好倒出水的临界面积

思路:

参考一下这个博客推导很清晰:http://blog.youkuaiyun.com/danliwoo/article/details/53002695

求出刚好倒出的角度,用二分倾斜角度的方法求出答案

要记得椭圆投影成圆的面积的公式 S圆=S椭圆cosX

还有过不去的时候把精度1E-8改成1E-10 [黑脸][黑脸][黑脸]

#include<bits/stdc++.h>
using namespace std;
const double pi=acos(-1.0);
const double epx=1e-10;
double f(double x)
{
	double y=acos(2*tan(x)-1);
	double a1,a2;
	a1=pi*cos(y)-y*cos(y)+sin(y)-(1.0/3.0)*sin(y)*sin(y)*sin(y);
	y=pi;
	a2=pi*cos(y)-y*cos(y)+sin(y)-(1.0/3.0)*sin(y)*sin(y)*sin(y);
	return (1.0/tan(x))*(a1-a2);
}
double fen(double L)
{
	double l=0.0,r=pi/4.0;
	double mid;
	while(r-l>=epx)
	{
		mid=(l+r)/2.0;
		if(f(mid)>L) r=mid;
		else l=mid;
	}
	return mid;
}
int main()
{
	double a,b,c,x,y;
	int t;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lf",&a);
		if(a==0) printf("%.5f\n",0.0);
	    else if(a>=1.0) {
	    	b=2.0-a;
	    	b=sqrt(1+b*b);
	    	printf("%.5f\n",b*pi);
	    }
	    else
	    {
	        x=fen(a*pi);
	        y=acos(2*tan(x)-1);
	        y=pi-y+cos(y)*sin(y);
	        printf("%.5f\n",y/sin(x));
	    }
	}
	return 0;
}




评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值