L3-021 神坛 (30分)

在古老的迈瑞城，巍然屹立着 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

 

样例解释

输出的数值等于图中红色或紫色框线的三角形的面积。

 

思路：面积需要用到叉积来计算，大体分为两个步骤： ①极角排序 ②先判断象限，再比较角度 ③最后叉积算最小面积

注意数据范围要开到long long，不然会卡样例

 

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#define ll long long
using namespace std;
const int N=5050;
int n;
struct Node{
	ll x;
	ll y;
	int x_x;//象限 
}a[N],b[N];
int judge(Node x) {	
	if(x.x > 0 && x.y > 0)  return 1;	
	if(x.x > 0 && x.y < 0)  return 2;	
	if(x.x < 0 && x.y < 0)  return 3;	
	if(x.x < 0 && x.y > 0)  return 4;
}
bool cmp(Node a1,Node b1) {
	if(a1.x_x != b1.x_x) return a1.x_x < b1.x_x;
	return a1.x*b1.y - a1.y*b1.x < 0; 
}
int main() {
	double ans=-1;
	int cnt;
	cin >> n;
	for(int i=1;i<=n;i++)	cin >> a[i].x>>a[i].y;
	
	for(int i=1;i<=n;i++) {
		cnt=1;
		for(int j=1;j<=n;j++) {
			if(i==j)
				continue;
			b[cnt].x=a[j].x-a[i].x;
			b[cnt].y=a[j].y-a[i].y;
			b[cnt].x_x = judge(b[cnt]);
			cnt++;
		}
		sort(b+1,b+n,cmp); //极角排序
		for(int j=1;j<n-1;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("%.3f\n",ans);
	return 0;
}