已知三个点 p0,p1,p2 的叉积函数 cross:
}
叉积的一个重要性质,判断两向量互相之间的顺逆时针关系。
若 P×Q>0,则 P 在 Q 的顺时针方向;
若 P×Q<0,则 P 在 Q 的逆时针方向;
double cross(point p0,point p1,point p2)
{
return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y) ;}
叉积的一个重要性质,判断两向量互相之间的顺逆时针关系。
若 P×Q>0,则 P 在 Q 的顺时针方向;
若 P×Q<0,则 P 在 Q 的逆时针方向;
若 P×Q=0,则 P 和 Q 共线,但可能同向也可能反向;
注:Distance 函数名字可以换掉,但是不能换成 distance 会编译出错,原因是 distance 是 stl 中 计算距离的函数.按道理会自动选用我们自己定义的,但是不明白为什么报错,总之是会报错了;
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <algorithm>
using namespace std;
#define maxn 50+3
struct point
{
int x;
int y;
int no;
};
point P[maxn+3],res[maxn];
int pos; //全局变量,记录当前的 基准点
int cross(const point &p0, const point &p1, const point &p2) //计算叉积
{
return (p1.x-p0.x)*(p2.y-p0.y) - (p2.x-p0.x)*(p1.y-p0.y);
}
int Distance(const point &a, const point &b) //两向量间的距离
{
return (b.x-a.x)*(b.x-a.x) + (b.y-a.y)*(b.y-a.y);
}
bool cmp(const point & p1,const point & p2) //sort 的比较函数. 极角排序函数
{
int ans = cross(P[pos],p1,p2); //计算叉积
if(ans > 0) //叉积>0 则以 p[pos] 位置为基准 p1 在 p2 的顺时针方向
return true;
else if(ans == 0 && Distance(P[pos],p1)<Distance(P[pos],p2)) //叉积=0 三点共线, 距离短的在前面
return true;
return false;
}
void swap(point &a, point &b) //交换函数
{
point temp = a;
a = b;
b = temp;
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("in","r",stdin);
#endif
int M,N;
cin>>M;
while(M--)
{
cin>>N;
for(int i = 1; i <= N; i++)
{
cin>>P[i].no>>P[i].x>>P[i].y;
if(P[i].y < P[1].y)
swap(P[i],P[1]);
}
pos = 1;
sort(P+pos+1,P+N+1,cmp);
res[pos++] = P[pos];
for(int i = 2; i <= N; i++)
{
sort(P+pos+1,P+N+1,cmp);
res[pos++] = P[pos];
}
cout<<pos-1;
for(int i = 1; i <= N; i++)
{
cout<<" "<<res[i].no;
}
cout<<endl;
}
}