#10016. 「一本通 1.2 练习 3」灯泡

【题目描述】

原题来自：ZOJ 3203

相比 wildleopard 的家，他的弟弟 mildleopard 比较穷。他的房子是狭窄的而且在他的房间里面仅有一个灯泡。每天晚上，他徘徊在自己狭小的房子里，思考如何赚更多的钱。有一天，他发现他的影子的长度随着他在灯泡和墙壁之间走到时发生着变化。一个突然的想法出现在脑海里，他想知道他的影子的最大长度。

【输入格式】

输入文件的第一行包含一个整数 T ，表示测试数据的组数。

对于每组测试数据，仅一行，包含三个实数 H，h 和 D，H 表示灯泡的高度，h 表示 mildleopard 的身高，D 表示灯泡和墙的水平距离。

【输出格式】

输出文件共 T 行，每组数据占一行表示影子的最大长度，保留三位小数。

【样例输入】

3
2 1 0.5
2 0.5 3
4 3 4

【样例输出】

1.000
0.750
4.000

【数据范围与提示】

T≤100，10^−2≤H,h,D≤103，10^−2≤H−h。

大致思路：

首先，显而易见，有两种可能。当 影子恰好都在地上 或 影子一段在地上一段在墙上。

一、对于第一种可能，我们设影子到墙的距离为x（实际为L=0的情况），根据三角形相似，我们可得

$$\frac{h}{x}=\frac{H}{D}$$

移项可得

$$x=\frac{Dh}{H}$$

二，对于第二种可能

我们同样设影子到墙的距离为x，在墙上的距离为L，很明显，人到灯的距离为D-x，作平行于D的直线，同样根据三角形相似列出方程，表示出影子长度x+L

$$\frac{x}{D}=\frac{h-L}{H-L}$$

$$x(H-L)=D(h-L)$$

$$xH-xL=Dh-DL$$

$$L=\frac{Dh-Hx}{D-x}$$

因此可得

影子长度 $$L+x=\frac{Dh-Hx}{D-x}+x$$

到这里，基本的思路就已经很清晰了，但为什么要三分，又为什么能三分呢？

证明：

方法：对函数$L=\frac{Dh-Hx}{D-x}$求导，可得：

$$L^{\prime }=f(x{)}^{\prime }=\frac{Hx-HD+Dh-Hx}{(D-x{)}^2}$$

$找到f(x{)}^{\prime }的零点x=a，可发现当x>a时f(x{)}^{\prime }<f(a{)}^{\prime }且当x<a时f(x{)}^{\prime }<f(a{)}^{\prime }，因此此函数为单峰凸函数，满足三分$

代码实现

1. check函数

double check(double x){
   ans1=h*D/H;
   double ans2=x+((D*h-H*x)/(D-x));
   if(x<=ans1)return ans2;
   else return ans1;
}

这里注意一点，因为x为影子在地上的长度，因此我们在第一种情况所得的x即为影子在地上的长度最大值，虽然实际上x只在地上的长度也是一个函数，但是我们只需要判断其最大值即可（想一想为什么） 若x<$x_{max}$，那么这时即有一段影子在墙上的情况，只有这时才可以套用我们情况二的公式，否则返回$x_{max}$

2. 精度提醒！！！ 精度提醒！！精度提醒！！精度提醒！ 精度提醒

double eps=1e-11

本题精度要求较高，务必将精度设置的尽量高

3. 剩余部分基本为三分模板，不做太多解释

node

#include<bits/stdc++.h>
using namespace std;
int n,t;
double l,r,lmid,rmid,D,H,h,ans1;
double check(double x){
	double ans2=x+((D*h-H*x)/(D-x));
	if(x<=ans1)return ans2;
	else return ans1;
}
int main(){
	cin>>t;
	for(int tt=1;tt<=t;tt++){
		cin>>H>>h>>D;
		ans1=h*D/H;
		l=0;r=D;
		double eps=1e-11;
		while(r-l>eps){
			lmid=l+(r-l)/3;
			rmid=r-(r-l)/3;
			if(check(lmid)>=check(rmid))r=rmid;
			else l=lmid;
		}
		printf("%.3lf\n",max(ans1,check(l)));
	}
	return 0;
}