Cod sursa(job #656046)

Utilizator Marius_mFMI-M2 Marius Melemciuc Marius_m Data 3 ianuarie 2012 20:54:53
Problema Infasuratoare convexa Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.62 kb
#include <stdio.h>
#include <algorithm>
#define Nmax 120001
#define inf 1000000002

using namespace std;

struct punct
{
    double x, y, p;
} a[Nmax],s[Nmax],z;

int h, n = 2, poz;

bool cmp(punct a,punct b)
{
    if (a.p>b.p)
        return false;
    return true;
}

void citire()
{
	int i;
    z.x =z.y=inf;
    scanf ("%d ",&h);
    for (i=0;i<h;++i)
    {
        scanf ("%lf %lf ", &a[i].x, &a[i].y);
        if (a[i].x < z.x || a[i].x == z.x && a[i].y < z.y)
        {
            z.x = a[i].x;
            z.y = a[i].y;
            poz = i;
        }
    }
    a[poz]=a[h-1];
    --h;
}

void afisare()
{
	int i;
    printf ("%d\n",n);
    for (i=0;i<n;++i)
        printf ("%lf %lf\n", s[i].x, s[i].y);
}

void panta()
{
	int i;
    for (i=0;i<h;++i)
        if (a[i].x==z.x)
            a[i].p=inf;
        else
            a[i].p=(a[i].y-z.y)/(a[i].x-z.x);
}

bool determinant(punct a, punct b, punct c)
{
    double det = a.x * b.y + a.y * c.x + b.x * c.y - a.y * b.x - b.y * c.x - c.y * a.x;
    if (det >= 0)
        return true;
    return false;
}

void infasuratoare()
{
	int i;
    s[0] = z;
    s[1] = a[0];
    for (i=1;i<h;++i)
        if (determinant(s[n-2],s[n-1],a[i]))
            s[n++]=a[i];
        else
        {
            while (!determinant(s[n-2],s[n-1],a[i]))
                s[--n] = a[i];
            ++n;
        }
}

int main()
{
    freopen ("infasuratoare.in", "r", stdin);
    freopen ("infasuratoare.out", "w", stdout);
    citire();
    panta();
    sort(a, a+h, cmp);
    infasuratoare();
    afisare();
    return 0;
}