nyoj 78 圈水池

http://acm.nyist.net/JudgeOnline/problem.php?pid=78

分析：简单的凸包问题，很适合入门练习。

思路也相对简单，就是graham扫描算法。

这是我的程序，初次写凸包代码，还很粗糙，见谅见谅。

#include<iostream>
#include<algorithm>

using namespace std;

struct point{
double x,y;
};

bool cmp(const point &l,const point &r)
{
return l.x<r.x||(l.x==r.x&&l.y<r.y);
}

bool mult( point sp ,point ep ,point op )
{
return (sp.x-op.x)*(ep.y-op.y)>=(ep.x-op.x)*(sp.y-op.y);
}

/*bool operator < (const point &l,const point &r)
{
return l.y<r.y||(l.y==r.y&&l.x<r.x);
}*/

int graham(point pnt[],int n,point res[])
{
int i,len,top=1;
sort(pnt,pnt+n,cmp);
if(n==0) return 0;res[0]=pnt[0];
if(n==1) return 1; res[1]=pnt[1];
for(i=2;i<n;i++)
{
while(top&&mult(pnt[i],res[top],res[top-1]))
top--;
res[++top]=pnt[i];
}
len=top;
res[++top]=pnt[n-2];
for(i=n-3;i>=0;i--)
{
while(top!=len&&mult(pnt[i],res[top],res[top-1]))
top--;
res[++top]=pnt[i];
}
return top;
}

int main()
{
int K,m,i,j,tmp;
point pnt[103],res[103];
cin>>K;
while(K--)
{
cin>>m;
for(i=0;i<m;i++)
{
cin>>pnt[i].x>>pnt[i].y;
}
tmp=graham(pnt,m,res);
sort(res,res+tmp,cmp);
for(i=0;i<tmp;i++)
cout<<res[i].x<<" "<<res[i].y<<endl;
}
return 0;
}