当前位置: 代码迷 >> 综合 >> XMU 1401 寻找最大内角
  详细解决方案

XMU 1401 寻找最大内角

热度:62   发布时间:2023-12-13 23:35:03.0

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

易错点: 浮点型比大小: 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;}
//你的代码结束