Cod sursa(job #1002387)

Utilizator CosminRusuCosmin Rusu CosminRusu Data 27 septembrie 2013 17:44:00
Problema Infasuratoare convexa Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 2 kb
#include <fstream>
#include <vector>
#include <bitset>
#include <queue>
#include <algorithm>
#include <utility>
#include <cstring>
#include <string>
#include <stack>
#include <deque>
#include <iomanip>
#include <set>
#include <map>
#include <cassert>
#include <ctime>
#include <list>
#include <iomanip>

#define x first
#define y second

using namespace std;

string file = "infasuratoare";

ifstream cin( (file + ".in").c_str() );
ofstream cout( (file + ".out").c_str() );

const int MAXN = 120005;
const int oo = 0x3f3f3f3f;

typedef pair<double, double> pdd;
const inline int min(const int &a, const int &b) { if( a > b ) return b;   return a; }
const inline int max(const int &a, const int &b) { if( a < b ) return b;   return a; }
const inline void Get_min(int &a, const int b)    { if( a > b ) a = b; }
const inline void Get_max(int &a, const int b)    { if( a < b ) a = b; }


pdd P[MAXN];
int N, p, k, st[MAXN];
double xm = oo, ym = oo, X, Y;

struct ClassComp {
    inline bool operator () (const pdd &a, const pdd &b) const {
        return (a.x > X) * (b.y - Y) < (b.x - X) * (a.y - Y);
    }
};

inline bool Match(pdd x, pdd y, pdd z) {
    double sol=x.x*y.y+y.x*z.y+z.x*x.y-x.y*y.x-y.y*z.x-z.y*x.x;
    return sol > 0;
}

int main() {
    cin >> N;
    for(int i = 1 ; i <= N ; ++ i) {
        cin >> P[i].x >> P[i].y;
        if(P[i].x < xm || ( P[i].x == xm && P[i].y < ym) ) {
            xm = P[i].x;
            ym = P[i].y;
            p = i;
        }
    }
    X = P[p].x;
    Y = P[p].y;
    swap(P[p], P[N]);
    sort(P + 1, P + N, ClassComp());
    st[ ++ k ] = N;
    for(int i = 1 ; i <= N ; ++ i) {
        if(p == i)
            continue;
        while( k >= 2 && Match(P[st[k-1]], P[st[k]], P[i]))
            -- k;
        st[++k] = i;
    }
    cout << k << fixed << setprecision(6) << '\n';
    for(int i = k ; i ; -- i)
        cout << P[st[i]].x << ' ' << P[st[i]].y << '\n';
    cin.close();
    cout.close();
    return 0;
}