F - 程序设计：掎角之势 计蒜客 - A2226 （求三角形四心坐标的方法/求内切圆、外接圆面积的方法）

在《亮剑》中，李云龙发动的平安战役不仅把晋西北搅成一锅粥，还成为了第二次世界大战的转折点。

晋绥军 358 团楚云飞得到探报，有大批日军向平安县城移动，于是他准备帮帮场子。友军在平安附近有三个团的兵力，分别是李云龙的独立团，丁伟的新一团，孔捷的新二团，他们三人互成掎角之势。

作为黄埔军校毕业生的楚云飞，对自己的数学非常自信。在研究了友军的防区分布后，想知道李云龙、丁伟、孔捷这三人指挥部位置所形成的三角形的内切圆和外接圆的面积。因为他们的位置是不断变化的，每次重新算一遍很麻烦，想请你编程帮忙来计算。

输入格式

第一行一个整数 T，表示数据的组数。
接下来T 行，每行 6 个整数 x1,y1,x2,y2,x3,y3分别表示李云龙、丁伟、孔捷这三人指挥部的坐标。

输出格式

对于每组数据，如果不能组成三角形，输出"NO SOLUTION"，否则输出两个空格分隔的实数，分别表示内切圆和外接圆的面积。相对误差或者绝对误差在 10^−6 范围内就认为是正确的。

数据范围

对于 40% 的数据：T≤100x1=y2=x3=y3=0，0≤y1,x2≤100。
对于 100% 的数据：−10000≤x1,y1,x2,y2,x3,y3≤10000。

Sample Input

4
0 3 4 0 0 0
0 0 10 10 -10 -10
3 3 5 3 3 5
0 0 -5 -10 5 -10

Sample Output

3.1415926536 19.6349540849
NO SOLUTION
1.0780241689 6.2831853072
29.9995403716 122.7184630309

求三角形内切圆、外接圆的半径方法：

设三角形三边为 a、b、c

半周长 p=(a+b+c)/2

三角形面积 S=√[p(p-a)(p-b)(p-c)] ……（海伦公式）

 内切圆半径 r = S/p=√[(p-a)(p-b)(p-c)/p]= ½√[(-a+b+c)(a-b+c)(a+b-c)/(a+b+c)]

外接圆半径 R= abc/(4S)= ¼ abc/√[p(p-a)(p-b)(p-c)]= abc/√[(a+b+c)(-a+b+c)(a-b+c)(a+b-c)]

R、r、S 关系rR = S/p * abc/(4S) = abc/[2(a+b+c)]

答案：

#include <iostream>
#include<bits/stdc++.h>
#define ll long long
#define PI 3.1415926536
const int N = 1e5 + 10;
using namespace std;

double f(double x1,double y1,double x2,double y2)  //求长度（两点间距离公式）
{
    double len=sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
    return len;
}

double neixin(double len1,double len2,double len3) //求内心，得到内接圆的半径  （海伦公式）
{
    double p=(len1+len2+len3)/2;
    double s=sqrt(p*(p-len1)*(p-len2)*(p-len3));
    double r=s/p;
    return r;
}

double waixin(double x1,double y1,double x2,double y2,double x3,double y3)  //求外心，得到外接圆的半径
{
    double xx,yy;
    double a1=x2-x1,b1=y2-y1,c1=(a1*a1+b1*b1)/2;
    double a2=x3-x1,b2=y3-y1,c2=(a2*a2+b2*b2)/2;
    double d=a1*b2-a2*b1;
    xx=x1+(c1*b2-c2*b1)/d;
    yy=y1+(a1*c2-a2*c1)/d;
    double R=f(x1,y1,xx,yy);
    return R;
}

int main()
{
    ios::sync_with_stdio(0);
    int t;
    cin>>t;
    while(t--)
    {
        double x1,y1;
        double x2,y2;
        double x3,y3;
        cin>>x1>>y1>>x2>>y2>>x3>>y3;
        double len1=f(x1,y1,x2,y2); //求边长
        double len2=f(x2,y2,x3,y3);
        double len3=f(x3,y3,x1,y1);
        double xp=(x1-x2)*(y1-y3)*1.0;//判断是否组成三角形
        double yp=(x1-x3)*(y1-y2)*1.0;
        if(xp!=yp)
        {
            double r,R;
            r=neixin(len1,len2,len3);
            R=waixin(x1,y1,x2,y2,x3,y3);
            double s,S;
            s=PI*r*r;  //求面积
            S=PI*R*R;
            printf("%.15f %.15f\n",s,S);
        }
        else
            printf("NO SOLUTION\n");
    }
    return 0;
}

求内心的坐标（x，y）的方法

void neixin(double x1,double y1,double x2,double y2,double x3,double y3,double len1,double len2,double len3) //求内心
{
    double x,y;
    double r;
    x=(x1*len1+x2*len2+x3*len3)/(len1+len2+len3);
    y=(y1*len1+y2*len2+y3*len3)/(len1+len2+len3);
}

求外心的坐标（xx，yy）的方法

void waixin(double x1,double y1,double x2,double y2,double x3,double y3)  //求外心
{
    double xx,yy;
    double a1=x2-x1,b1=y2-y1,c1=(a1*a1+b1*b1)/2;
    double a2=x3-x1,b2=y3-y1,c2=(a2*a2+b2*b2)/2;
    double d=a1*b2-a2*b1;
    xx=x1+(c1*b2-c2*b1)/d;
    yy=y1+(a1*c2-a2*c1)/d;
}

求重心的坐标（x，y）的方法

void zhongxin(double x1,double y1,double x2,double y2,double x3,double y3)  //求重心
{
    double x,y;
    x=(x1+x2+x3)/3.0;
    y=(y1+y2+y3)/3.0;

}

求垂心的坐标（cx，cy）的方法

void chuixin(double x1,double y1,double x2,double y2,double x3,double y3)
{
    double x,y;  //求重心
    x=(x1+x2+x3)/3.0;
    y=(y1+y2+y3)/3.0;
    double xx,yy;  //求外心
    double a1=x2-x1,b1=y2-y1,c1=(a1*a1+b1*b1)/2;
    double a2=x3-x1,b2=y3-y1,c2=(a2*a2+b2*b2)/2;
    double d=a1*b2-a2*b1;
    xx=x1+(c1*b2-c2*b1)/d;
    yy=y1+(a1*c2-a2*c1)/d;
    double xc,yc;
    xc=x*3.0-xx*2.0;
    yc=y*3.0-yy*2.0;
}