Cod sursa(job #2098579)

Utilizator BlueLightningMihai Alexandru BlueLightning Data 3 ianuarie 2018 01:55:22
Problema Infasuratoare convexa Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 2.38 kb
///Acoperire convexa. Graham's Scan
#include <iostream>
#include <iomanip>
#include <fstream>
#include <algorithm>
#define N 120005
using namespace std;

ifstream f("infasuratoare.in");
ofstream g("infasuratoare.out");

struct punct
{
    double x,y;
};

punct stiva[N],v[N];///v = vectorul cu punctele date; stiva = stiva in care se va afla acoperirea convexa
int n,vf;///n = numarul de puncte; vf = numar puncte de pe acoperirea convexa

int test(punct a,punct b,punct c)///test orientare
{
    double o=(b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
    if(o==0)///punctele sunt coliniare
        return 0;
    else
        if(o>0)///viraj la stanga
            return 1;
        else///viraj la dreapta
            return 2;
}

int dist(punct a,punct b)///distanta la patrat dintre doua puncte
{
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}

bool cmp(punct a,punct b)///sortare in functie de unghiul si distanta polara
{
    int o=test(v[1],a,b);
    if(o==0)
        if(dist(v[1],b)<=dist(v[1],a))
            return false;
        else
            return true;
    else
        if(o==1)
            return true;
        else
            return false;
}

int main()
{
    int i;
    ///citire
    f>>n;
    for(i=1;i<=n;i++)
        f>>v[i].x>>v[i].y;
    ///se alege cel mai din stanga-jos punct; va fi primul punct al acoperirii convexe
    int poz=1;
    for(i=2;i<=n;i++)
        if((v[i].x<v[poz].x) or (v[i].x==v[poz].x and v[i].y<v[poz].y))
            poz=i;
    swap(v[1],v[poz]);
    ///sortarea punctelor dupa unghiul polar/distanta polara
    sort(v+2,v+n+1,cmp);
    ///primele doua varfuri vor face parte din acoperirea convexa
    stiva[1]=v[1];
    stiva[2]=v[2];
    vf=2;
    for(i=3;i<=n;i++)///verificam pentru restul varfurilor daca fac parte din acoperirea convexa
    {
        while(vf>=2 and test(stiva[vf-1],stiva[vf],v[i])!=1)///cat timp sunt cel putin doua varfuri in acoperirea convexa
            ///se scot varfuri cat timp varful curent nu formeaza cu ultimele doua varfuri un viraj la stanga
            vf--;
        stiva[++vf]=v[i];///se adauga varful curent la stiva ce va contine acoperirea convexa
    }
    ///afisare
    if(vf<3)
        g<<"Nu e posibil.";
    else
    {
        g<<vf<<"\n";
        for(i=1;i<=vf;i++)
            g<<fixed<<setprecision(6)<<stiva[i].x<<" "<<stiva[i].y<<"\n";
    }
}