STL6----K - Light Bulb

最新推荐文章于 2019-09-18 23:56:31 发布

原创最新推荐文章于 2019-09-18 23:56:31 发布 · 159 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#三分

细节题同时被 3 个专栏收录

18 篇文章

订阅专栏

经典题型

17 篇文章

订阅专栏

三分

1 篇文章

订阅专栏

本文探讨了在ZOJ3203问题中，如何利用高精度计算解决复杂数学问题。通过采用特定制导搜索算法，如三等分法，实现了对函数solve(x)的精确求解。该文详细介绍了算法的实现过程，包括初始化参数、设置搜索区间及迭代求解，最终输出了求解结果。

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

zoj3203

这道题对精度的要求很高。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
const double MAX = 1e4;
const double ESP = 1e-10;
double H, h, D;
double solve(double x)
{
 return D - x + H - (H - h)*D / x;
 //h * (D + x) / H - x + x * H / (D + x);//(h * (D + x) * (D + x) + H * H * x - D * H * x - x * x * H) / (H*(D + x));
 //return h * (D + x) / H - x + x * h / ();//(h * (D + x) * (D + x) + H * H * x - D * H * x - x * x * H) / (H*(D + x));
}
int main()
{
 int t; cin >> t;
 while (t--)
 {
  cin >> H >> h >> D;
  double lt = (H - h)*D / H, rt = D, mid, midmid;
  //cout << "ckjjv" << endl;
  while (rt - lt >= ESP)
  {
   //cout << lt << " " << rt << endl;
   mid = lt + (rt - lt) / 3;
   midmid = rt - (rt - lt) / 3;
   if (solve(mid) >= solve(midmid)) rt = midmid;
   else lt = mid;
  }
  printf("%.3lf\n", solve(lt));
 }
 return 0;
}