cf 1c(计算几何，海伦，三角函数，精度，外接圆)

Ancient Berland Circus

原题地址
vj地址

思路：

题目给了一个三角形，一开始没什么思路，就先连边吧。

1. 把三个点先连起来。（求出两两之间的距离）
2. 三角形先求下外接圆(一个正多边形的三个顶点在外接圆上，那么这个三角形的外接圆就是该正多边形的外接圆)
3. 有了外接圆，可还是没思路，先求个半径吧
   1)中学知识，p=（a+b+c）/ 2,
   2)s=sqrt(p*(p-a)(p-b)(p-c)) (海伦公式求面积)
   3）a / sinA = b / sinB = c / sinC = 2R
   4）s= (absinC) / 2
   联立上式可得 R=（abc）/ 4s。
4. 有了外接圆半径又陷入了沉思，该怎么做呢，
   。。。。。。。。。。。。
   。。。。。。。。。。。。
   求这个外接圆里的正多边形的area，可以先求每一小块，三角形开始。把正n边形分成n个三角形。
   eg四边形–四个三角形
   eg六边形–六个三角形

---

那每一块的三角形area怎么求呢
。。。。。。。。。。。。。。。。。。。。。。。。。

（盗用洛谷的一张图）
要求出每一小块三角形的area，即要求出每一块三角形在 2pi 中的角度t。
（因为s=r *r *sint / 2）
根据图，题目给出的三角形，可以求出每一条边的长度a，b，c，
a这条弦对应的圆心角，b这条弦所对应的圆心角，c这条弦所对应的圆心角
之后算出每条边所对应的圆心角，而t就是他们的公有的因数了（gcd）

#include <iostream>
#include <cmath>
#include <cstdio>
using namespace std;
#define  pi acos(-1.0)
#define eps 1e-2
double dis(double x1, double y1, double x2, double y2)
{
    double d=(x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
    return sqrt(d);
}
double getcos(double a, double r)
{
    double cos1=(r*r+r*r-a*a)/(2*r*r);
    return cos1;
}
double gcd(double a, double b)
{
    if(fabs(b)<eps)return a;
    return gcd(b,fmod(a,b));//double类型的gcd
}
int main()
{
    double x1,y1,x2,y2,x3,y3;
    cin>>x1>>y1>>x2>>y2>>x3>>y3;
    double d1,d2,d3;
    d1=dis(x1,y1,x2,y2);
    d2=dis(x1,y1,x3,y3);
    d3=dis(x2,y2,x3,y3);
    double p=(d1+d2+d3)/2.0;
    double s=sqrt(p*(p-d1)*(p-d2)*(p-d3));
    double R=d1*d2*d3/s/4.0;
    double angle1=acos(getcos(d1,R));
    double angle2=acos(getcos(d2,R));
    double angle3=2*pi-angle1-angle2;//为了避免误差，第三条边的
    double angle=gcd(angle3,gcd(angle1,angle2));
    printf("%.6f",pi*R*R*sin(angle)/angle);
    return 0;
}
/*
double mod(double a,double b)//小数取模
{
  return a-(int)(a/b)*b;
}
*/