程序设计天梯赛L3-21 (极角排序)

题目
题意: 在二维空间种给定n个点，选出任意三个点组成一个三角形，使得三角形的面积最小。(若不能组成三角形，输出0) n <= 5000
思路: n^3暴力(doge)
  看网上的说法是可以以某个点为A，花O(n)的时间求出与其他点的向量，之后按照斜率从小到大排序，求出所有相邻两个向量对应边的最小值，以该点为A的三角形的最小面积即为其中的最小值。但是这样不严谨，因为有可能这个点不能为A，所以要将所有点作为A枚举一遍。

以左下角为A就不够用，还得枚举到左上角的点为A才能找到真正的最小值。

时间复杂度: O(nnlogn)
代码:

#include<bits/stdc++.h>
using namespace std;
const int N = 5002;
#define int long long
int n,m,k,T;
struct node{
	int x,y;
	bool operator<(const node&rhs)const
	{
		return (y*rhs.x) < (x*rhs.y);
	}
}a[N],va[N];
typedef pair<int,int> PII;
map<PII,int> mp;
double ans = 1e18;
void solve()
{
	cin>>n;
	bool flag = false;
	for(int i=0;i<n;++i)
	{
		cin>>a[i].x>>a[i].y;
	    if(mp[{a[i].x,a[i].y}]) flag = true;
		mp[{a[i].x,a[i].y}]++;
	}
	for(int i=0;i<n&&!flag;++i)
	{
		int tot = 0;
		for(int j = 0; j < n; j++)
		{
			if(i==j)continue;
			va[tot].x = a[j].x-a[i].x;
			va[tot].y = a[j].y-a[i].y;
			tot++;
		}
		sort(va,va+tot);
		for(int j=0;j<tot;++j)
		{
			int x1 = va[j].x;
			int y1 = va[j].y;
			int x2 = va[(j+1)%tot].x;
			int y2 = va[(j+1)%tot].y;
			ans = min(ans,0.5 * fabs(x1*y2-x2*y1)  );
		}
	}
	if(flag) ans = 0;
	printf("%.3lf",ans);
}
signed main(void)
{
	solve();
	return 0;
}