POJ 1696 Space Ant

题目大意求一旋转向内的连线，可以想象成蜗牛的壳

不得不说这个我弄的狠纠结，一开是用的Graham，带上模板，一WA到死，怒了，代码全删了，自己看过包裹法，又认真看了遍叉积，点积，自己写了个只用叉积的算法。

算法思路：

step 1: 第一次取y最小，y相等取x最小的点。p0

step 2: 每次取第一个未标记的点pk，尝试拓展 p0-pk ，若发现存在pi 使得 p0-pi 在p0-pk 的左边（他们的叉积  < 0），则 k = i ；

step 3: 重复执行step 2，直到取得所有点

参考代码：

struct POINT {
    int x,y;
    int idx;
    bool mark;
    POINT(int i,int a,int b) {
        idx=i;x=a;y=b;mark=false;
    }
    POINT(){}
}p[MAXN],res[MAXN];

int n,idx;

inline double multiply(POINT o,POINT s,POINT e) {
    return ((s.x - o.x) * (e.y - o.y) - (s.y - o.y) * (e.x - o.x));
}


void ground() {
    int cnt = 0;
    for(int i = 1;i<n;i++) {
        if(p[i].y < p[0].y ||
         ((p[i].y == p[0].y) &&
          (p[i].x < p[0].x)))
          swap(p[0],p[i]);
    }
    res[cnt++] = p[0];
    p[0].mark = true;
    while(cnt < n) {
        int k = 0;
        for(int i = 0;i<n;i++) {
            if(!p[i].mark) {
                k = i;
                break;
            }
        }
        for(int i = k + 1;i<n;i++) {
            if(!p[i].mark) {
                if(multiply(res[cnt-1],p[k],p[i]) < 0)
                    k = i;
            }
        }
        res[cnt++] = p[k];
        p[k].mark = true;
    }
}

void init() {
    SCF(n);
    int id,a,b;
    F(i,n) {
        SCFT(id,a,b);
        p[i] = POINT(id,a,b);
    }
    ground();
    PCF(n);
    F(i,n) printf(" %d",res[i].idx);
    puts("");
}


int main() {
    FOPEN;
    int cse;SCF(cse);
    while(cse --) {
        init();
    }
}