Pagini recente » Cod sursa (job #1219902) | Cod sursa (job #437459) | Cod sursa (job #1685011) | Cod sursa (job #1935811) | Cod sursa (job #2098579)
///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";
}
}