话说上回讲到海东集团推选老总的事情,最终的结果是XHD以微弱优势当选,从此以后,“徐队”的称呼逐渐被“徐总”所取代,海东集团(HDU)也算是名副其实了。
创业是需要地盘的,HDU向钱江肉丝高新技术开发区申请一块用地,很快得到了批复,据说这是因为他们公司研发的“海东牌”老鼠药科技含量很高,预期将占全球一半以上的市场。政府划拨的这块用地是一个多边形,为了描述它,我们用逆时针方向的顶点序列来表示,我们很想了解这块地的基本情况,现在请你编程判断HDU的用地是凸多边形还是凹多边形呢?
Input
输入包含多组测试数据,每组数据占2行,首先一行是一个整数n,表示多边形顶点的个数,然后一行是2×n个整数,表示逆时针顺序的n个顶点的坐标(xi,yi),n为0的时候结束输入。
Output
对于每个测试实例,如果地块的形状为凸多边形,请输出“convex”,否则输出”concave”,每个实例的输出占一行。
Sample Input
4 0 0 1 0 1 1 0 1 0
Sample Output
convex 海东集团终于顺利成立了!后面的路,他们会顺顺利利吗? 欲知后事如何,且听下回分解——
就是判断每个点形成的图形是不是凸包
#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<cstdio>
#include<math.h>
#include<stack>
#define inf 0x3f3f3f
using namespace std;
struct node
{
double x;
double y;
node friend operator -(node a,node b)
{
return {a.x-b.x,a.y-b.y};
}
}ac[105],s[105];
double dis(node a,node b)
{
node c=a-b;
return sqrt(c.x*c.x+c.y*c.y);
}
double X(node a,node b)
{
return a.x*b.y-a.y*b.x;
}
int cmp(node a,node b)
{
double x=X(a-ac[1],b-ac[1]);
if(x>0)return 1;
if(x==0&&dis(a,ac[1])<dis(b,ac[1])) return 1;
return 0;
}
double mulit(node a,node b,node c)
{
return X(b-a,c-a);
}
int main()
{
int n;
while(scanf("%d",&n),n)
{
for(int i=1;i<=n;i++)
{
scanf("%lf%lf",&ac[i].x,&ac[i].y);
}
if(n==1) {printf("0.00\n");continue;}
if(n==2) {printf("%.2f\n",dis(ac[1],ac[2]));continue;}
int k=1;
for(int i=2;i<=n;i++)
{
if(ac[i].y<ac[k].y||(ac[i].y==ac[k].y&&ac[i].x<ac[k].x))
k=i;
}
swap(ac[1],ac[k]);
sort(ac+2,ac+n+1,cmp);
s[1]=ac[1];
s[2]=ac[2];
int tot=2;
for(int i=2;i<=n;i++)
{
while(tot>=2&&mulit(s[tot-1],s[tot],ac[i])<=0) tot--;
s[++tot]=ac[i];
}
//s[tot+1]=s[1];
if(tot==n)
cout<<"convex"<<endl;
else
cout<<"concave"<<endl;
// printf("%.2f\n",sum);
}
return 0;
}