POJ 1696 Space Ant (向量的叉积)

一只很特殊的蚂蚁，只能向坐转，并且不能经过已经走过的路。一张地图上有n个食物让蚂蚁去采集，求蚂蚁经过食物的顺序。

一开始想偷懒，直接甩graham凸包模版上去，加了一条判断和搜索点的重排。。果断WA

好吧果然还是不能偷懒啊！每确定一个点以它为起点找逆时针一个个碰到的点。。其实大概或许就是卷包裹求凸包了吧。。不太了解的说。。

判断点用的就是叉积了

//Memory: 260K//Time: 0MS#include <iostream>#include <string.h>#include <math.h>#include <algorithm>using namespace std;bool use[100];struct POINT{double x,y;int num;};POINT PointSet[100],ch[100];double dist(POINT p1,POINT p2)              { return( sqrt( (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y) ) ); }/******************************************************************************\r=multiply(sp,ep,op),得到(sp-op)和(ep-op)的叉积 r>0：ep在矢量opsp的逆时针方向； r=0：opspep三点共线； r<0：ep在矢量opsp的顺时针方向 \******************************************************************************/ double multiply(POINT sp,POINT ep,POINT op) { return((sp.x-op.x)*(ep.y-op.y)-(ep.x-op.x)*(sp.y-op.y)); }void getans(int n) { int i,k=0,sum=1,minp=50,p=0;bool flag;while(sum<=n){minp=100000;flag=true;for(i=1;i<=n;i++){if(!use[ PointSet[i].num ]){if(flag){minp=i;flag=false;}else{if(multiply(PointSet[minp],PointSet[i],PointSet[p])<0||// 极角更小  (multiply(PointSet[minp],PointSet[i],PointSet[p])==0) &&/* 极角相等，距离更短 */        dist(PointSet[p],PointSet[minp])>dist(PointSet[p],PointSet[i]))minp=i;}}}ch[sum++]=PointSet[minp];use[ PointSet[minp].num ]=true;p=minp;}}int main(){int n,cas,i;cin>>cas;while(cas--){cin>>n;memset(use,0,sizeof(use));double miny=1000000;for(i=1;i<=n;i++){cin>>PointSet[i].num>>PointSet[i].x>>PointSet[i].y;miny=min(miny,PointSet[i].y);}PointSet[0].x=0;PointSet[0].y=miny;getans(n);cout<<n;for(i=1;i<=n;i++)cout<<" "<<ch[i].num;cout<<endl;}return 0;}