ZOJ-3806

最新推荐文章于 2018-08-01 20:00:53 发布

心之痕

最新推荐文章于 2018-08-01 20:00:53 发布

阅读量404

点赞数

分类专栏： ZOJ 计算几何文章标签： ZOJ

本文链接：https://blog.csdn.net/xhldtc/article/details/38820377

版权

ZOJ 同时被 2 个专栏收录

546 篇文章 3 订阅

订阅专栏

计算几何

8 篇文章 0 订阅

订阅专栏

几何基础没怎么做过。。太弱太弱了，在基本的精度问题上WA了好久。。首先二分的方法没写对，开始是以内切圆的半径误差为条件来跳出循环，这样一直TLE，看了看别人的代码是以二分左右边界误差小于eps后跳出循环的，改了之后还是WA，然后读输入的时候要以double来读，你妹的题目明明说的是正整数啊，用int读有什么问题么。。判断那个无解的条件时也要加上double的误差控制，坑死。。

#include<cstdio>
#include<cmath>

using namespace std;

namespace
{
	double inner_radius(double a, double b, double c)
	{
		double p = (a + b + c) / 2;
		double s = sqrt(p * (p - a) * (p - b) * (p - c));
		return s / p;
	}
}

int main()
{
	double r, R;
	const double eps = 1e-10;
	while (scanf("%lf %lf", &r, &R) != EOF)
	{
		if (R + eps < 2 * r)
		{
			puts("NO Solution!");
			continue;
		}
		double left = 0, right = sqrt(3.0) * R, center, len, rr, temp;
		while (right - left > eps)
		{
			center = (left + right) / 2;
			temp = sqrt(R * R - center * center / 4) + R;
			len = sqrt(temp * temp + center * center / 4);
			rr = inner_radius(len, len, center);
			if (rr < r)
				left = center;
			else
				right = center;
		}
		printf("%.18lf %.18lf %.18lf\n", center, len, len);
	}
	return 0;
}