Cod sursa(job #1003666)

Utilizator Adrian1997Radulescu Adrian Adrian1997 Data 1 octombrie 2013 11:07:40
Problema Infasuratoare convexa Scor 50
Compilator cpp Status done
Runda Arhiva educationala Marime 1.19 kb
#include <fstream>
#include <stdio.h>
#include <limits.h>
#include <algorithm>
using namespace std;
ifstream f("infasuratoare.in");
FILE *g=fopen("infasuratoare.out","w");
int n,S[120010];

struct puncte{
    double x,y,tg;
};
puncte v[120010],t,ct;

double det(puncte x,puncte y,puncte z){
    return (x.x*y.y+y.x*z.y+z.x*x.y-z.x*y.y-y.x*x.y-x.x*z.y);
}

int cmp(puncte a,puncte b){
    return (det(v[1],a,b)>0);
}

int main(void){
    register int i,j,p=0,k;

    f>>n;
    v[0].x=INT_MAX,v[0].y=INT_MAX;
    for(i=1;i<=n;i++){
        f>>v[i].x>>v[i].y;
        if(v[i].x<v[p].x ||(v[i].x==v[p].x && v[i].y<v[p].y))
            p=i;
    }

    ct=v[p];
    v[p]=v[1],v[1]=ct;
    v[1].x=0,v[1].y=0;
    for(i=2;i<=n;i++){
        v[i].x-=ct.x;
        v[i].y-=ct.y;
        v[i].tg=(v[i].x==0?INT_MAX:(v[i].y/v[i].x));
    }

    sort(v+2,v+n+1,cmp);

    S[1]=1,S[2]=2,k=2;
    for(i=3;i<=n;i++){
        while(k>=2 && det(v[S[k-1]],v[S[k]],v[i])<=0)
            k--;
        S[++k]=i;
    }

    fprintf(g,"%d\n",k);
    for(i=1;i<=k;i++){
        fprintf(g,"%.6f %.6f\n",v[S[i]].x+ct.x,v[S[i]].y+ct.y);
    }
    f.close();
    fclose(g);
    return 0;
}