[清华集训2015]灯泡(浙江大学ZOJ 3203 Light Bulb）

Time Limit: 1 Second      Memory Limit: 32768 KB

Compared to wildleopard's wealthiness, his brother mildleopard is rather poor. His house is narrow and he has only one light bulb in his house. Every night, he is wandering in his incommodious house, thinking of how to earn more money. One day, he found that the length of his shadow was changing from time to time while walking between the light bulb and the wall of his house. A sudden thought ran through his mind and he wanted to know the maximum length of his shadow.

Input

The first line of the input contains an integer T (T <= 100), indicating the number of cases.

Each test case contains three real numbers H, h and D in one line. H is the height of the light bulb while h is the height of mildleopard. D is distance between the light bulb and the wall. All numbers are in range from 10^-2 to 10³, both inclusive, and H - h >= 10^-2.

Output

For each test case, output the maximum length of mildleopard's shadow in one line, accurate up to three decimal places..

Sample Input

3
2 1 0.5
2 0.5 3
4 3 4

Sample Output

1.000
0.750
4.000

题目连接

题目大意：

他的房子狭窄，他家里只有一个灯泡。天天晚上，他都在他不知名的房子里徘徊，想着如何赚更多的钱。有一天，他发现他的影子长度在灯泡和房子的墙壁之间行走时不时变化。忽然想到了他的思绪，他想知道他的影子的最大长度。

输入的第一行包含整数T（T <= 100），表示个案数。

每一个测试用例在一行中 包含三个实数H，h和D. H是灯泡的高度，而h是轻度高度的高度。 D是灯泡和墙壁之间的距离。全部数字的范围均为10 ^-2至10 ³，包括二者，H - h > = 10 ^-2。

对于每一个测试用例，在一行中输出mildleopard阴影的最大长度，精确到三位小数。

算法分析

参考：

https://blog.csdn.net/Mrx_Nh/article/details/52745348    分析理解应该是对的，可是代码有误，估计是没有讨论影子不落到墙上的状况。

https://blog.csdn.net/qq_38944163/article/details/80870928  代码是正确的。

这个题要分两种状况讨论：影子没有落到墙上；影子落到墙上。

影子没有落到墙上的状况以下图所示：

影子落到墙上的状况以下图所示

 

算影子的长度能够用类似三角形，相信各位上过初中的朋友都会。。

那么咱们发现这是一个单峰函数，要求峰值，这里咱们能够用一个很强大的东西，三分法。

具体看代码吧

 1 #include<stdio.h>
 2 double f(double x,double H,double h,double D)
 3 {
 4     double tanA,Y;
 5     tanA=(H-h)/x;//这里你们能够画个图帮助理解
 6     Y=H/tanA;//下面那条的长（在没有墙的状况下）
 7     if(Y<=D) return Y-x;//若是没到墙，就直接返回地上影子的长度。
 8     Y=Y-D;
 9     return D-x+Y*tanA;//D-X是地上影子的长度,Y*tanA是墙上影子的长度.
10 }
11 int main()
12 {
13     freopen("p3382.in","r",stdin);
14     int i,T;
15     double H,h,D;
16     double L,R,m1,m2,fm1,fm2;
17     scanf("%d",&T);
18     for(i=0;i<T;i++)
19     {
20         scanf("%lf%lf%lf",&H,&h,&D);
21         L=0;R=D;
22         while(L+1e-11<R)
23         {
24             m1=L+(R-L)/3; m2=R-(R-L)/3;
25             fm1=f(m1,H,h,D);  fm2=f(m2,H,h,D);
26             if(fm1<fm2) L=m1;
27             else R=m2;
28         }
29         printf("%.3lf\n",f(L,H,h,D));
30     }
31     return 0;
32 }

 

为什么是一个单峰，能够考虑看一下第一个参考的博客。

理解时参考下图：