ZOJ - 3806 ：Incircle and Circumcircle

最新推荐文章于 2019-06-23 21:05:26 发布

Mitsuha_

最新推荐文章于 2019-06-23 21:05:26 发布

阅读量362

点赞数 1

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

本文链接：https://blog.csdn.net/Mitsuha_/article/details/77685099

版权

计算几何专栏收录该内容

43 篇文章 0 订阅

订阅专栏

A triangle is one the basic shapes in geometry. It's a polygon with three vertices and three sides which are line segments. A triangle with vertices A, B, C is denoted ΔABC. And its three sides, BC, CA, AB are often denoted a, b and c.

The incircle of a triangle is the largest circle contained in the triangle, it is tangent to the three sides. The center of the incircle is called the triangle's incenter and the radius of the incircle is called inradius, which is denoted r.

The circumcircle of a triangle is the circle which passes through all its three vertices. The center of the this circle is called circumcenter and its radius is called circumradius, which is denoted R.

It's very easy to calculate r and R after knowing a, b and c. Now you are given rand R, can you calculate a valid triple (a, b, c)?

Input

There are multiple cases. Each case has two positive integers r and R in one line. (r and R are less than 10⁵)

Ouput

For each case, print three floating numbers a, b and c in one line if it's possible to get r and R. If there are no possible tuples, print "NO Solution!".

The judge program uses your a, b and c to calculate the inradius and circumradius. Only the relative error of your inradius and circumradius less than 10^-8 will be accepted.

Sample Input

1 2
2 5
9 9

Sample Ouput

3.464101615137754587 3.464101615137754587 3.464101615137754587
6 8 10
NO Solution!

思路：http://www.cnblogs.com/Silence-AC/p/3935014.html

#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<iostream>
using namespace std;
int main()
{
    double r,R;
    while(cin>>r>>R)
    {
        if(R<2*r){puts("NO Solution!");continue;}
        double down=0,up=sqrt(3)*R,mid,a;
        while(up-down>1e-10) //二分枚举底边长
        {
            mid=(up+down)/2;
            a=sqrt((sqrt(R*R-(mid*mid)/4)+R)*(sqrt(R*R-(mid*mid)/4)+R)+(mid*mid)/4); //等腰三角形的腰长
            double p=(a+a+mid)/2;

            //三角形内切圆半径公式 r=2*S/(a+b+c) (S为三角形面积，a，b，c分别是三角形三条边)
            if(2*sqrt(p*(p-a)*(p-a)*(p-mid))/(a+a+mid)-r>1e-12)
            {
                up=mid;
            }
            else down=mid;
        }
        printf("%0.16lf %0.16lf %0.16lf\n",a,a,mid);
    }
    return 0;
}