XMU 1401 寻找最大内角

思路：枚举 利用余弦定理和反三角函数求出最大角

易错点： 浮点型比大小： a > b  &&  fabs( a - b ) > 1e-6

描述

请定义Point结构体，Point结构体有两个int域代表点的坐标。

现在给定N个点的坐标，问在这些点中任意选取3个点作为三角形，有可能得到的最大内角是多少。

注意：你提交的代码将被插入到以下程序框架中一起编译，所以请不要提交全部代码。允许添加自定义函数。

//固定头部开始
#include <stdio.h>
#include <math.h>
//固定头部结束

//你的代码开始
double GetMaxAngle(struct Point pts[], int n)
{

}
//你的代码结束

//固定尾部开始
int main()
{
    int i, n;
    struct Point pts[50];

    scanf("%d", &n);
    for (i = 0; i < n; i++)
        scanf("%d%d", &pts[i].x, &pts[i].y);

    printf("%.6f", GetMaxAngle(pts, n));

    return 0;
}
//固定尾部结束

输入

第一行输入一个整数N，代表点的个数。（N \le 50N≤50）

接下来N行每行两个整数X Y，代表每个点的坐标，（0 \le X,Y \le 10000≤X,Y≤1000），保证没有重叠点以及三点共线。

输出

输出一个实数，表示最大内角的度数，弧度制，保留6位小数。

输入样例 1 

4
0 0
5 5
3 5
4 6

输出样例 1

2.896614

提示

余弦定理：c^2 = a^2 + b^2 - 2ab \cos Cc2=a2+b2−2abcosC

其中a、b、c是三角形的三条边，边c是角C的对边。

//你的代码开始
struct Point
    {
    	int x;
    	int y;
	};
 
	
	double dis(int a, int b, int c, int d)
	{
	    return sqrt((a-c)*(a-c)*1.0 + (b-d)*(b-d)*1.0);	
	}
	
	double angle(double a, double b, double c)
	{
		return acos((a*a + b*b -c*c)/(2*a*b));
	}
	
	double max(double a, double b, double c)
	{
		a = (a>b && fabs(a-b)>1e-6) ? a : b;
		a = (a>c && fabs(a-c)>1e-6) ? a : c;
		return a;
	}
 
double GetMaxAngle(struct Point pts[], int n)
{
    
	
	double maxa = 0;
	
	for(int i=0; i<n; i++)
	  for(int j=i+1; j<n; j++)
	    for(int k=j+1; k<n; k++)
	    {
	    	double a = dis(pts[i].x, pts[i].y, pts[j].x, pts[j].y);
	    	double b = dis(pts[j].x, pts[j].y, pts[k].x, pts[k].y);
	    	double c = dis(pts[i].x, pts[i].y, pts[k].x, pts[k].y);
	    	double amax = max(angle(a,b,c),angle(a,c,b),angle(b,c,a));
	    	if( amax - maxa > 0 && fabs(amax - maxa)>1e-6) maxa = amax;
		}
	
	return maxa;
	
}
//你的代码结束