Cod sursă (job #199227)

Utilizator avatar calinalexandru Alexandru Calin calinalexandru IP ascuns
Problemă Înfășurătoare convexă Compilator cpp | 1,39 kb
Rundă Arhiva de probleme Status evaluat
Dată 13 feb. 2016 17:23:35 Scor 0
#include <fstream>
#include <iostream>
#include <algorithm>
#include <iomanip>
using namespace std;
ifstream in("infasuratoare.in");
ofstream out("infasuratoare.out");
struct punct
{
    double x,y;
} pct[120005],minp,rasp[120005];
bool det(punct a,punct b,punct c)
{
    if((a.x*b.y + b.x*c.y + c.x*a.y - a.y*b.x - b.y*c.x - c.y*a.x)>0) return true;
    else return false;
}
bool comp(punct a,punct b)
{
    if(det(minp,a,b)) return false;
    else return true;
}
int main()
{
    int n,i,poz,poz2,j,k=2;
    in>>n;
    for(i=1; i<=n; i++)
        in>>pct[i].x>>pct[i].y;
    minp.y=pct[1].y;
    poz=1;
    for(i=2; i<=n; i++)
        if(pct[i].y<minp.y)
        {
            minp.y=pct[i].y;
            poz=i;
        }
    minp.x=pct[poz].x;
    for(i=2; i<=n; i++)
        if(pct[i].y==minp.y && pct[i].x<minp.x)
        {
            minp.x=pct[i].x;
            poz=i;
        }
    swap(pct[poz].y,pct[n].y);
    swap(pct[poz].x,pct[n].x);
    n--;
    sort(pct+1,pct+n+1,comp);
    rasp[1]=minp;
    rasp[2]=pct[1];
    for(i=2; i<=n; i++)
    {
        while(det(rasp[k-1],rasp[k],pct[i]))
            k--;
        rasp[++k].x=pct[i].x;
        rasp[k].y=pct[i].y;
    }
    out<<k<<'\n';
    out<<fixed<<setprecision(6)<<rasp[1].x<<" "<<rasp[1].y<<'\n';
    for(i=k;i>1;i--)
        out<<fixed<<setprecision(6)<<rasp[i].x<<" "<<rasp[i].y<<'\n';
}