Cod sursa(job #693085)

Utilizator Oancea.CatalinOancea Catalin Oancea.Catalin Data 27 februarie 2012 09:19:19
Problema Infasuratoare convexa Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.69 kb
#include <algorithm>
#include <stdio.h>
#define INF 0x3f3f3f
#define NMAX 120001
using namespace std;
int n;
struct puncte{
    double x,y,pante;
};
puncte a[NMAX];
double pxmin,pymin;
int pozmin;
int stiva[NMAX],nr;
double calcularepanta(int xx)
{
    if (a[0].x==a[xx].x)
        return INF;
    return ((a[xx].y-a[0].y)/(a[xx].x-a[0].x));
}
bool cmp(puncte x1,puncte x2)
{
    return (x1.pante<x2.pante);
}
void citire()
{
    freopen("infasuratoare.in","r",stdin);
    scanf("%d",&n);
    scanf("%lg%lg",&a[1].x,&a[1].y);
    pxmin=a[1].x;
    pymin=a[1].y;pozmin=1;
    for (int i=2;i<=n;++i)
    {
        scanf("%lg%lg",&a[i].x,&a[i].y);
        if (a[i].x<pxmin)
            pxmin=a[i].x,pymin=a[i].y,pozmin=i;
        else if (a[i].x==pxmin && a[i].y<pymin)
            pxmin=a[i].x,pymin=a[i].y,pozmin=i;
    }
    a[0].x=pxmin;
    a[0].y=pymin;
    a[pozmin].x=a[n].x;
    a[pozmin].y=a[n].y;
    --n;
    for (int i=1;i<=n;++i)
        a[i].pante=calcularepanta(i);
    sort(a+1,a+n+1,cmp);
}
double determinant(int x1,int y1,int z1)
{
    double s1=0;
    s1+=(a[x1].x*a[y1].y);s1+=(a[y1].x*a[z1].y);s1+=(a[x1].y*a[z1].x);
    s1-=(a[y1].y*a[z1].x);s1-=(a[x1].x*a[z1].y);s1-=(a[y1].x*a[x1].y);
    return s1;

}
void rezolvare()
{
    stiva[0]=0;
    stiva[1]=1;
    nr=1;
    for (int i=2;i<=n;i++)
    {
        while (nr>0 && determinant(stiva[nr-1],stiva[nr],i)<0)
                nr--;
        stiva[++nr]=i;
    }
    freopen("infasuratoare.out","w",stdout);
    printf("%d\n",nr+1);
    for (int i=0;i<=nr;i++)
        printf("%.6lf %.6lf\n",a[stiva[i]].x,a[stiva[i]].y);
}
int main()
{
    citire();
    rezolvare();
    return 0;
}