ZOJ 3203 Light Bulb

在LJOJ上说这是三分,我就懒得去推公式了。学习三分法第一题。十分方便。

//By Richard
#include <cstdio>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <ctime>
#define rep(x,y,z) for (int x=(y);(x)<=(z);(x)++)
#define per(x,y,z) for (int x=(y);(x)>=(z);(x)--)
#define log2(x) (31-__builtin_clz(x))
#define mod (int)(1e9+7)
#define inf 0x3f3f3f3f
#define cls(x) memset(x,0,sizeof(x))
#ifdef DEBUG
#define debugdo(X) X
#define debugndo(X)
#define debugout(X) cout<<(#X)<<"="<<(X)<<endl
#else
#define debugdo(X)
#define debugndo(X) X
#define debugout(X)
#endif // debug
#ifdef ONLINE_JUDGE
#define debugdo(X)
#define debugndo(X)
#define debugout(X)
#endif
#define putarray(x,n) rep(iiii,1,n) printf("%d ",x[iiii])
#define mp make_pair
using namespace std;
typedef pair<int,int> pairs;
typedef long long LL;
/read3.0
template <typename T>
inline void read(T &x){char ch;x=0;bool flag=false;ch=getchar();while (ch>'9'||ch<'0') {ch=getchar();if (ch=='-') flag=true;}while ((ch<='9'&&ch>='0')){x=x*10+ch-'0';ch=getchar();}if (flag) x*=-1;}
template <typename T>
inline void read(T &x,T &y){read(x);read(y);}
/variables&functions


int T;
double H,h,D;
const double eps=0.0000001;
// #define calc(x) (x+((h/(x+(h*D-H*x)/(H-h)))*(h*D-H*x)/(H-h)))
// #define calc(x) (x+((h*D-H*x)/(D-x)))
inline double calc(double x)
{
	double temp=(h*D-H*x)/(D-x);
	if (temp>0) return x+temp;
	else return (D-x)/(H-h);
}
int main()
{
	read(T);
	while (T--)
	{
		scanf("%lf%lf%lf",&H,&h,&D);
		double l=0,r=D;
		while (r-l>eps)
		{
			double ll=(r-l)/3+l,rr=(r-l)*2/3+l;
			double ansll=calc(ll),ansrr=calc(rr);
			if (ansll>ansrr) r=rr;
			else l=ll;
		}
		printf("%.3lf\n",calc(l));
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值