在古老的迈瑞城,巍然屹立着 n 块神石。长老们商议,选取 3 块神石围成一个神坛。因为神坛的能量强度与它的面积成反比,因此神坛的面积越小越好。特殊地,如果有两块神石坐标相同,或者三块神石共线,神坛的面积为 0.000。
长老们发现这个问题没有那么简单,于是委托你编程解决这个难题。
输入格式:
输入在第一行给出一个正整数 n(3 ≤ n ≤ 5000)。随后 n 行,每行有两个整数,分别表示神石的横坐标、纵坐标(−10
9
≤ 横坐标、纵坐标 <10
9
)。
输出格式:
在一行中输出神坛的最小面积,四舍五入保留 3 位小数。
输入样例:
8
3 4
2 4
1 1
4 1
0 3
3 0
1 3
4 2
输出样例:
0.500
样例解释
输出的数值等于图中红色或紫色框线的三角形的面积。
#include <bits/stdc++.h>
using namespace std;
int n;
struct node
{
long long int x;
long long int y;
int w;
}a[5001],b[5001];
//先根据极角进行排序,这样可以把三维降成二维
bool cmp(node a,node b)
{
if(a.w!=b.w)
{
return a.w<b.w;
}
return a.x*b.y-b.x*a.y>0;
//根据右手螺旋定则,向量a叉乘向量b,右手四指指向向量a的方向,然后向着手心方向转向b的方向,拇指指向的就是叉积向量的方向
//这里我们为了让向量按照极角进行排序,极角是向量方向与ox方向的夹角,所以只有当向量a叉乘向量b小于0时才表明a的极角比b的小
}
int judge(node a)
{
if(a.x>0&&a.y>0)
{
return 1;
}
else if(a.x<0&&a.y>0)
{
return 2;
}
else if(a.x<0&&a.y<0)
{
return 3;
}
else if(a.x>0&&a.y<0)
{
return 4;
}
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
{
scanf("%lld %lld",&a[i].x,&a[i].y);
}
double ans=-1;
int cnt;
for(int i=1;i<=n;i++)
{
cnt=1;
for(int j=1;j<=n;j++)
{
if(i==j)
{
continue;
}
b[cnt].x=a[i].x-a[j].x;
b[cnt].y=a[i].y-a[j].y;
b[cnt].w=judge(b[cnt]);
cnt++;
}
sort(b+1,b+n,cmp);//注意b只有n-1个元素
/* for(int j=1;j<=n;j++)
{
cout<<b[j].x<<" "<<b[j].y<<endl;
}
cout<<endl;*/
for(int j=2;j<n;j++)
{
if(ans==-1||fabs(b[j-1].x*b[j].y-b[j-1].y*b[j].x)*0.5<ans)
{
ans=fabs(b[j-1].x*b[j].y-b[j-1].y*b[j].x)*0.5;
}
}
}
printf("%.3lf\n",ans);
return 0;
}