POJ 1696 极角排序

已知三个点 p0,p1,p2 的叉积函数 cross:

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;
    }
}