吃巧克力 计算几何

吃巧克力

Time Limit: 1000MS		Memory Limit: 1000K

Description

CYF 喜欢吃巧克力，有一天，CYF收到了一份镶嵌在平面直角坐标系上的巧克力，这个巧克力的形状是一个正多边形。日子一天天过去，CYF吃完了这份巧克力，但是却忘记了自己吃的巧克力有几条边，只记得正多边形巧克力上三个点的坐标。

CYF知道送自己巧克力的人非常小气，所以正多边形巧克力的边数不会超过100。现在的CYF非常想知道自己吃的正多边形巧克力到底有多少条边，你能帮他算一下吗？

Input

每个样例包含一组输入。

一个输入有三行每行是一个点的坐标

Xi Yi

坐标之间用空格隔开

Output

输出CYF吃的巧克力最少有多少个角

Sample Input

0.000000 0.000000

1.000000 1.000000

0.000000 1.000000

Sample Output

4

题解：算出来三角形的外心  然后for i

算三个角的方法是   算出来两个小于π的角   然后第三个角就等于2π减去两个角

这个方法是防止夹角大于π

因为acos返回的是[0,π]之间的角

然后for i 求出i边形的内角度 判断这个角度的某个倍数与三个角的差是否小于eps  如果是  输出i

ps：本题精度有点坑  eps 1e-8过不了  1e-6就能过

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define pi acos(-1.0)
#define eps 1e-6
struct Point
{
    double x;
    double y;
};
struct Traingle
{
    struct Point p[3];
};
double Dis(struct Point p, struct Point q)
{
    double dx=p.x-q.x;
    double dy=p.y-q.y;
    return sqrt(dx*dx+dy*dy);
}
double solve(Point a,Point b,Point cen){//求两个向量的夹角
	double xx=a.x-cen.x,yy=a.y-cen.y,xxx=b.x-cen.x,yyy=b.y-cen.y;
	return acos((xx*xxx+yy*yyy)/Dis(a,cen)/Dis(b,cen));
}
int jud(double ci,int t){//判断t边形下  内角的某个倍数是否等于ci
	double sd=ci*t/2.0/pi;
	if(fabs(sd-(int)(sd+0.5))<eps)return 1;
	return 0;
}
int main(){
	Traingle s;
	scanf("%lf%lf%lf%lf%lf%lf",&s.p[0].x,&s.p[0].y,&s.p[1].x,&s.p[1].y,&s.p[2].x,&s.p[2].y);
	double x1=s.p[0].x,y1=s.p[0].y,x2=s.p[1].x,y2=s.p[1].y,x3=s.p[2].x,y3=s.p[2].y;
	double disd=0,ac1,ac2,ac3;
	//求外接圆圆心  
    double t1=x1*x1+y1*y1;  
    double t2=x2*x2+y2*y2;  
    double t3=x3*x3+y3*y3;  
    double temp=x1*y2+x2*y3+x3*y1-x1*y3-x2*y1-x3*y2;  
    double x=(t2*y3+t1*y2+t3*y1-t2*y1-t3*y2-t1*y3)/temp/2;  
    double y=(t3*x2+t2*x1+t1*x3-t1*x2-t2*x3-t3*x1)/temp/2;  
    Point center;
    center.x=x;
    center.y=y;
	disd=Dis(s.p[0],s.p[1]);
	disd=max(disd,Dis(s.p[1],s.p[2]));
	disd=max(disd,Dis(s.p[0],s.p[2]));
	if(fabs(disd-Dis(s.p[0],s.p[1]))<eps){
		ac1=solve(s.p[0],s.p[2],center);
		ac2=solve(s.p[1],s.p[2],center);
		ac3=2*pi-ac1-ac2;
	}
	else if(fabs(disd-Dis(s.p[1],s.p[2]))<eps){
		ac1=solve(s.p[0],s.p[2],center);
		ac2=solve(s.p[1],s.p[0],center);
		ac3=2*pi-ac1-ac2;
	}
	else{
		ac1=solve(s.p[1],s.p[0],center);
		ac2=solve(s.p[1],s.p[2],center);
		ac3=2*pi-ac1-ac2;
	}
	for(int i=3;i<=100;i++){
		if(jud(ac1,i)&&jud(ac2,i)&&jud(ac3,i)){
			printf("%d\n",i);
			return 0;
		}
	}
	return 0;
}