Cod sursa(job #280171)

Utilizator mircea_infoSuciu Mircea-Gabriel mircea_info Data 13 martie 2009 11:25:52
Problema Infasuratoare convexa Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.2 kb
#include <stdio.h>
#define nmax 120001

struct punct {double x,y;} a[nmax],pct;
int n,s[nmax],m;

double prod_vec(punct pct, punct a, punct b)
{
	return (a.x-pct.x)*(b.y-pct.y) - (b.x-pct.x)*(a.y-pct.y);
}

void swap(punct &x, punct &y)
{
	punct aux=x;
	x=y;
	y=aux;
}

void sort(punct arr[], int beg, int end)
{
	if(end>beg+1)
	{
		punct piv=arr[beg];
		int l=beg+1, r=end;
		while(l<r)
		{
			if(prod_vec(pct,arr[l],piv)>0)
				  l++;
			else
				swap(arr[l], arr[--r]);
		}
		swap(arr[--l], arr[beg]);
		sort(arr, beg, l);
		sort(arr, r, end);
	}
}

int main()
{
	freopen("infasuratoare.in","r",stdin);
	freopen("infasuratoare.out","w",stdout);
	scanf("%d",&n);
	scanf("%lf %lf",&a[1].x,&a[1].y);
	pct=a[1];
	int i,j=1;
	for(i=2;i<=n;i++)
	{
		scanf("%lf %lf",&a[i].x,&a[i].y);
		if(a[i].y<pct.y || (a[i].y==pct.y && a[i].x<pct.x))
		{
			pct=a[i];
			j=i;
		}
	}
	punct aux=a[1];
	a[1]=a[j];
	a[j]=aux;
	sort(a,2,n+1);
	m=2;
	s[1]=1;
	s[2]=2;
	for(i=3;i<=n;i++)
	{
		while(m>=2 && prod_vec(a[s[m]],a[i],a[s[m-1]])<=0)
			m--;
		s[++m]=i;
	}
	printf("%d\n",m);
	for(i=1;i<=m;i++)
		printf("%lf %lf\n",a[s[i]].x,a[s[i]].y);
	fclose(stdout);
	return 0;
}