poj2189 凸包

这道题真的让我感受颇深啊！一直WA!我就一直以为是在算法那块写错了，于是我就在主要算法的那个几个函数一直在改，然后就一直wa，然后我就坚信是我的算法思想解决有问题，于是乎，我分析了别人的代码，然后把别人的部分不断的改我自己的，然后还是没有对，只到我改到了最后的一个函数，一个我从来没有检查也坚信不会错的，但是就是马虎出错的函数，让我从下午wa到晚上啊！！！！！以后检查程序的时候一定要全方位的细节的，别的就不多说了，经历了之后才发现认真是多么重要啊！

看到题的第一遍感觉是暴力，但是又一想如果暴力可以过的话就不会AC率不高了。。。于是就采用凸包的方法！

所谓凸包就是找到一个点集合中最外部的点，使得其余的点要么在这些点连线的里面或者是在线上！

对于凸包的算法可以参见：http://www.2cto.com/kf/201301/184235.html

首先对于坐标的排序，可以采用极坐标也可以直接求角度，然后建立栈不断的进行更新就可以，我WA的原因就是建立栈后对凸包外边的点求距离时没有用建立的栈的数据而是直接用了原来数组的，马虎啊马虎！

下面贴两个代码，一个是我本来就写的wa了好多次又改过来的，另一个是建于别人代码改过来的！

#include<iostream>
#include<cmath>
#include<cstdio>
#include<stack>
using namespace std;
#define exp 0.0000000001
typedef struct point
{
	int x,y;
	double degree;
}rr;
point a[50005];
int n;
int b[50005],Max;
int cmp(const void *a,const void *b)
{
	point *c=(point *)a;
	point *d=(point *)b;
	if(fabs(c->degree-d->degree)>exp)
	return c->degree>d->degree?1:-1;
	else//相等的情况
		return c->x-d->x;//升序排列
}
void Sort()
{
	int i,j,x=INT_MAX,y=INT_MAX;
	for(i=0; i<n; i++)
	{
		if(a[i].y<y)
		{
			y=a[i].y;j=i;x=a[i].x;
		}
		else if(a[i].y==y && a[i].x<x)
		{
			y=a[i].y;j=i;x=a[i].x;
		}
	}//找到最小的
	a[j].x=a[0].x;a[j].y=a[0].y;
	a[0].x=x;a[0].y=y;//j和0的位置进行互换的
	for(i=1; i<n; i++)
	{
        a[i].degree=acos(((double)(a[i].x-x))/(sqrt(pow((double)(a[i].x-x),2)+pow((double)(a[i].y-y),2))));
	}
	qsort(&a[1],n-1,sizeof(point),cmp);//对其进行升序排列的
//	for(i=0; i<n; i++)
//		printf("%d %d  %lf\n",a[i].x,a[i].y,a[i].degree);
}
double judge(point s,point t,point i)
{
	return (t.x-s.x)*(i.y-s.y)-(t.y-s.y)*(i.x-s.x);
}
void fun()
{
	int i,s,t;
	stack<int>q;
	q.push(0);
	q.push(1);//先进去两个的
	for(i=2; i<n; i++)
	{
		t=q.top();
		q.pop();
		if(!q.empty())
		{
		while(!q.empty())
		{
			s=q.top();
			if(judge(a[s],a[t],a[i])<-exp)//将t删除
			{
               q.pop();
			   t=s;
			   if(q.empty())//为空
			   {
                  q.push(t);//将t放入
				  break;
			   }
			   s=q.top();
			}
			else
			{
				q.push(t);//将原来的还放入的
				break;
			}
		}
		}
		else
		{
			q.push(t);
		}
		q.push(i);
	}
	i=0;
	while(!q.empty())
	{
		b[i++]=q.top();
		q.pop();
	}
	n=i;//返回的是啊
    return ;
}
void len()
{
	int i,j;
	for(i=0; i<n-1; i++)
	{
		for(j=i+1; j<n; j++)
		{
			int s=(a[b[j]].x-a[b[i]].x)*(a[b[j]].x-a[b[i]].x)+(a[b[j]].y-a[b[i]].y)*(a[b[j]].y-a[b[i]].y);
			if(s>Max)
				Max=s;
		}
	}
}
int main()
{
	int i;
	while(scanf("%d",&n)!=EOF)
	{
		for(i=0; i<n; i++) 
			scanf("%d%d",&a[i].x,&a[i].y); 
		Sort();
        fun();
		Max=0;
		len();
		printf("%d\n",Max);
	}
	return 0;
}

#include<iostream>
#include<cmath>
#include<cstdio>
#include<stack>
using namespace std;
#define exp 0.0000000001
typedef struct point
{
	int x,y;
}rr;
point a[50005];
int n;
int b[50005],Max;
point *p;
int cross(point p1,point p2,point p3,point p4)
{
	return (p2.x-p1.x)*(p4.y-p3.y)-(p2.y-p1.y)*(p4.x-p3.x);
}
int distance(point p1,point p2,point p3,point p4)//a大了要排在后面的
{
	return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)>(p3.x-p4.x)*(p3.x-p4.x)+(p3.y-p4.y)*(p3.y-p4.y)?-1:1;
}
int cmp(const void *a,const void *b)
{
	point *c=(point *)a;
	point *d=(point *)b;
    int temp=cross(*p,*c,*p,*d);
	if(temp<0)
		return 1;
	else if(temp==0)
		return distance(*p,*c,*p,*d); 
	else
	return -1; 
}
void Sort()
{
	int i,j,x=INT_MAX,y=INT_MAX;
	for(i=0; i<n; i++)
	{
		if(a[i].y<y)
		{
			y=a[i].y;j=i;x=a[i].x;
		}
		else if(a[i].y==y && a[i].x<x)
		{
			y=a[i].y;j=i;x=a[i].x;
		}
	}//找到最小的
	a[j].x=a[0].x;a[j].y=a[0].y;
	a[0].x=x;a[0].y=y;//j和0的位置进行互换的
    p=&a[0];//这个是原点坐标的
//	p.x=x;p.y=y;
	qsort(&a[1],n-1,sizeof(point),cmp);//对其进行升序排列的
//	for(i=0; i<n; i++)
//		printf("%d %d  \n",a[i].x,a[i].y);
}
double judge(point s,point t,point i)
{
	return (t.x-s.x)*(i.y-s.y)-(t.y-s.y)*(i.x-s.x);
}
void fun()
{
	int i,s,t;
	stack<int>q;
	q.push(0);
	q.push(1);//先进去两个的
	for(i=2;; i++)
	{
		i=i%n;
		if(i==1)
			break;
		if(i!=0 && a[i].x==a[i-1].x && a[i].y==a[i-1].y)//去除重点
			continue;
		t=q.top();
		q.pop();
		if(!q.empty())
		{
		while(!q.empty())
		{
			s=q.top();
			if(judge(a[s],a[t],a[i])<-exp)//将t删除
			{
               q.pop();
			   t=s;
			   if(q.empty())//为空
			   {
                  q.push(t);//将t放入
				  break;
			   }
			   s=q.top();
			}
			else
			{
				q.push(t);//将原来的还放入的
				break;
			}
		}
		}
		else
		{
			q.push(t);
		}
		q.push(i);
	}
	i=0;
	while(!q.empty())
	{
		b[i++]=q.top();
		q.pop();
	}
	n=i;//返回的是啊
    return ;
}
/*void fun()
{
	int i,pc;
	b[1]=0;
	b[2]=1;
	pc=2;
	b[0]=2;
	for(i=2; i<=n;)
	{
		if(cross(a[b[pc-1]],a[b[pc]],a[b[pc]],a[i])<0)
		{
			pc--;
			b[0]--;
		}
		else
		{
			b[++pc]=i++;
			b[0]++;
		}
	}
	n=b[0];
}*/
void len()
{
	int i,j;
	for(i=0; i<n-1; i++)
	{
		for(j=i+1; j<n; j++)
		{
			int s=(a[b[j]].x-a[b[i]].x)*(a[b[j]].x-a[b[i]].x)+(a[b[j]].y-a[b[i]].y)*(a[b[j]].y-a[b[i]].y);
			if(s>Max)
				Max=s;
		}
	}
}
int main()
{
	int i;
	while(scanf("%d",&n)!=EOF)
	{
		for(i=0; i<n; i++) 
			scanf("%d%d",&a[i].x,&a[i].y); 
		Sort();
        fun();
		Max=0;
		len();
		printf("%d\n",Max);
	}
	return 0;
}