codeforces-1C Ancient Berland Ci

题意：给定一个正多边形三点坐标，求能组成的正多边形的最小面积。
思路：
1、根据三点坐标计算三边长度，然后计算该三点所组成三角形的外接圆半径
2、根据三边与半径求出圆心角
3、找到三个圆心角最大公约数
4、求正多边形的面积
正多边形一共有2PI / angle个小三角形
每个小三角形的面积为R2sin(angle) / 2

#include <stdio.h>
#include <math.h>
#define pi acos(-1.0)  
#define eps 1e-4

typedef struct Point
{
    double x,y;
}Point;

double dist(Point p1, Point p2){
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double Angle(double a,double b,double c){
    return acos((a*a+b*b-c*c)/(2*b*a));
}
double gcd(double a,double b){
    if(b+eps>0&&b-eps<0) return a;
    if(a+eps>0&&a-eps<0) return b;
    return gcd(b,fmod(a,b));
}

int main()
{
    /* code */

    Point p1,p2,p3;

    while(~scanf("%lf%lf%lf%lf%lf%lf",&p1.x,&p1.y,&p2.x,&p2.y,&p3.x,&p3.y)){

            double a,b,c,p,s,r;
            double angle1,angle2,angle3,angle;
            a=dist(p1,p2);
            b=dist(p1,p3);
            c=dist(p2,p3);
            //半周长
            p=(a+b+c)/2;
            //三角形面积
            s=sqrt(p*(p-a)*(p-b)*(p-c));
            //外接圆半径
            r=a*b*c/(4*s);
            angle1=Angle(r,r,a);
            angle2=Angle(r,r,b);
            angle3=2*pi - angle1 - angle2;
            angle=gcd(angle1,gcd(angle2,angle3));
            double ans;
            ans=0.5*r*r*sin(angle)*(2*pi/angle);
            printf("%.6lf\n",ans);
    }
    return 0;
}