Language: C++

#include <cassert>
#include <iostream>
#include <set>
#include <vector>
#include <cmath>
#include <iomanip>
using namespace std;

// Definiciones iniciales.
typedef long double Double;
typedef long long Long;
typedef pair<int, int> pii;

const Double INF = 1ll << 60;
const Double ERROR = 1e-9;
const Double M_2PI = 2.0 * M_PI;

// Tolerancia a flotantes
#define EQ(a, b) (fabs((a) - (b)) < ERROR)

// Punto en 2D.
struct Punto {
	Double x, y;
	Punto() : x(), y() {}
	Punto(Double X, Double Y) : x(X), y(Y) {}
	// Izquierda a derecha, abajo a arriba.
	bool operator<(const Punto& q) const {
		if (!EQ(x, q.x)) return x < q.x;
		return EQ(y, q.y)? false: y < q.y;
	}
	bool operator==(const Punto& q) const {
		return EQ(x, q.x) && EQ(y, q.y);
	}

	//Sumar los vectores actual y q.
	Punto operator +(const Punto& p)const { return Punto(x+p.x, y+p.y); }
	//Restar los vectores (o trasladarlos), tomando p como el origen.
	Punto operator -(const Punto& p)const { return Punto(x-p.x, y-p.y); }
	//Obtener el vector opuesto.
	Punto operator -() const { return Punto(-x, -y); }
	//Obtener el vector por un escalar.
	Punto operator*( Double c ) const { return Punto( x*c, y*c ); }
	//Obtener el vector dividido por un escalar.
	Punto operator/( Double c ) const { return Punto( x/c, y/c ); }
	//Magnitud del vector.
	Double Abs() const { return hypot(x, y); }
	//Magnitud elevada al cuadrado del vector.
	Double Norm() const{ return x*x + y*y; }
	//Distancia entre el punto actual y q.
	Double Dist( const Punto& q ) const { return (*this - q ).Abs(); }
	//Distancia elevada al cuadrado entre el punto actual y q.
	Double Dist2(const Punto& q ) const { return ( *this - q).Norm(); }
	//Producto punto entre el punto actual y q.
	Double Dot( const Punto& q ) const { return x*q.x + y*q.y; }
	//Producto cruz entre el punto actual y q.
	Double Cruz( const Punto& q ) const { return x*q.y - y*q.x; }

	//Rotar 90 grados en CCW.
	Punto rot90() const { return Punto(-y, x); }

	//Rotar t radianes en cw respecto al origen.
	Punto RotarCW(const Double& t) const {
		return Punto( x*cos(t) + y*sin(t),
						y*cos(t) - x*sin(t));
	}
	//Rotar t radianes en ccw respecto al origen.
	Punto RotarCCW(const Double& t) const {
		return Punto( x*cos(t) - y*sin(t), 
						y*cos(t) + x*sin(t));
	}
};

//Conversion de grados a radianes.
Double GradARad(Double grd) {
	return (grd * M_PI) / 180.0;
}

// Conversion de radianes a grados.
Double RadAGrad(Double rad) {
	return (rad * 180.0) / M_PI;
}

//Angulo mas chico formado por los puntos aob en radianes.
Double Angulo(const Punto& a, const Punto& o, const Punto& b){
	Punto u(a - o), v(b - o);
	return acos( u.Dot(v) / (u.Abs() * v.Abs()));
}

// Test de mano derecha: CCW = 1, CW = -1, Colineal = 0.
int ManoDerecha(const Punto& o, const Punto& p, const Punto& q) {
	Double ccw = (p-o).Cruz(q-o);
	return EQ(ccw, 0)? 0: (ccw < 0)? -1: 1;
}

// Par de puntos mas cercanos en un conjunto de puntos P.
pair<Punto, Punto> ParMasCercano(vector<Punto> P) {
	// Si ya esta ordenado, no usar sort.
	sort(P.begin(), P.end());

	set<Punto> rect;
	pair<Punto, Punto> par;
	int prev = 0; Double delta = 1e9;
	for (int i = 0; i < P.size(); ++i) {
		while (P[i].x - P[prev].x > delta)
			rect.erase(Punto(P[prev].y, P[prev].x)), prev++;

		set<Punto>::iterator it = rect.lower_bound(
			Punto(P[i].y - delta, P[0].x));

		for (; it != rect.end() && it->x <= P[i].y + delta; ++it) {
			Double dist = hypot(P[i].x - it->y, P[i].y - it->x);
			if (dist < delta) delta = dist, par = make_pair(
				Punto(it->y, it->x), P[i]);
		}
		rect.insert(Punto(P[i].y, P[i].x));
	}
	return par; // Alternativamente puede devolver delta.
}
// Linea en 2D.
struct Linea {
	Punto p, q;
	Double a, b, c;
	Linea() : p(), q(), a(), b(), c() {}
	Linea(Long A, Long B, Long C)
		: p(), q(), a(A), b(B), c(C) {
		if (EQ(a, 0)) {
			c /= -b; b = -1;
			p = Punto(0, c);
			q = Punto(1, c);
		} else if (EQ(b, 0)) {
			c /= -a; a = -1;
			p = Punto(c, 0);
			q = Punto(c, 1);
		} else {
			p = Punto(-c/a, 0);
			q = Punto(-(b+c)/a, 0);
		} if (q < p) swap(p, q);
	}
	Linea(const Punto& P, const Punto& Q)
		: p(P), q(Q), a(), b(), c() {
		// Asegura p como punto menor.
		if (q < p) swap(p, q);
		a = q.y - p.y;
		b = p.x - q.x;
		if (!a) c = p.y, b = -1;
		else if (!b) c = p.x, a = -1;
		else c = -a*p.x - b*p.y;
	}
	// PELIGRO! Ordena por ecuacion de recta.
	bool operator<(const Linea& cmp) const {
		if (!EQ(a, cmp.a)) return a < cmp.a;
		if (!EQ(b, cmp.b)) return b < cmp.b;
		return EQ(c, cmp.c)? false: c < cmp.c;
	}
};

// Saber si un punto p esta en la recta r.
bool PuntoEnRecta(const Punto& p, const Linea& r){
	return !ManoDerecha( r.p, r.q, p );
}
// Saber si un punto p esta en el segmento s.
bool PuntoEnSegmento(const Punto& p, const Linea& s) {
	return PuntoEnRecta(p, s) && !(p < s.p || s.q < p);
}

// Saber si dos lineas l y m son paralelas.
bool LineasParalelas(const Linea& l, const Linea& m) {
	return EQ( 0, ( l.q - l.p ).Cruz( m.q - m.p ) );
}

// Saber si dos lineas l y m son iguales.
bool LineasIguales(const Linea& l, const Linea& m) {
	return LineasParalelas(l, m)
		&& EQ( 0, ( l.p - l.q ).Cruz( l.p - m.q ) )
		&& EQ( 0, ( m.q - m.p ).Cruz( m.q - l.p ) );
}
// Saber si dos lineas l y m son perpendiculares.
bool LineasPerpendiculares(const Linea& l, const Linea& m) {
	return EQ( 0, ( l.q - l.p ).Dot( m.q - m.p ) );
}

// Obtener una linea paralela a l que pase por p.
Linea ParalelaEnPunto(const Linea& l, const Punto& p) {
	return Linea(p, p + ( l.q - l.p ));
}

// Obtener una linea perpendicular a l que pase por p.
Linea PerpendicularEnPunto(const Linea& l, const Punto& p) {
	return Linea(p, p + ( l.q - l.p ).rot90());
}

// Saber si dos rectas r y s se intersectan.
// No intersectan = 0, Interseccion en un punto = 1,
// Interseccion paralela en infinitos puntos = -1.
int InterseccionRectas(const Linea& r, const Linea& s) {
	if (LineasIguales(r, s)) return -1;
	return LineasParalelas(r, s)? 0: 1;
}

// Saber si una recta r y un segmento s se intersectan.
// No intersectan = 0, Interseccion en un punto = 1,
// Interseccion paralela en infinitos puntos = -1.
int IntersecRectaSegmen(const Linea& r, const Linea& s) {
	if (LineasIguales(r, s)) return -1;
	if (LineasParalelas(r, s)) return 0;
	int t1 = ManoDerecha(r.p, r.q, s.p);
	int t2 = ManoDerecha(r.p, r.q, s.q);
	return (t1 != t2)? 1: 0;
}

// Saber si dos segmentos s y t se intersectan.
// No intersectan = 0, Interseccion en un punto = 1,
// Interseccion paralela en infinitos puntos = -1.
int InterseccionSegmentos(const Linea& s, const Linea& t) {
	int t1 = ManoDerecha(s.p, s.q, t.p);
	int t2 = ManoDerecha(s.p, s.q, t.q);

	if (t1 == t2) return t1? 0:
		(PuntoEnSegmento(s.p, t) ||
		 PuntoEnSegmento(s.q, t) ||
		 PuntoEnSegmento(t.p, s) ||
		 PuntoEnSegmento(t.q, s))? -1: 0;

	return (ManoDerecha(t.p, t.q, s.p) !=
			ManoDerecha(t.p, t.q, s.q))? 1: 0;
}

// Obtener punto de interseccion entre lineas l y m.
Punto PuntoInterseccion(const Linea& l, const Linea& m) {
	Punto a = l.p, b = l.q, c = m.p, d = m.q;
	b = b-a; d = c-d; c =c-a;
	assert( !LineasParalelas( l, m ) );
	return a + b*(c.Cruz( d )/b.Cruz( d ));
}
// Obtener proyeccion del vector v en la recta r.
Punto ProyeccionEnRecta(const Punto& v, const Linea& r) {
	return r.p + ( r.q - r.p ) * 
			(( v - r.p ).Dot( r.q - r.p )/( r.q  - r.p ).Norm());
}

// Distancia entre un punto p y una recta r.
Double DistanciaPuntoRecta(const Punto& p, const Linea& r) {
	return p.Dist(ProyeccionEnRecta(p, r));
}

// Distancia entre un punto p y un segmento s.
Double DistanciaPuntoSegmento(const Punto& p, const Linea& s) {
	Punto proy = ProyeccionEnRecta(p, s);
	if (proy < s.p) return p.Dist(s.p);
	if (s.q < proy) return p.Dist(s.q);
	return p.Dist(proy);
}

// Distancia entre dos lineas l y m.
Double DistanciaRectaRecta(const Linea& l, const Linea& m) {
	return LineasParalelas(l, m)? DistanciaPuntoRecta(l.p, m): 0;
}

// Distancia entre dos segmentos s y r.
Double DistanciaSegmenSegmen(const Linea& s, const Linea& r) {
	if (InterseccionSegmentos(s, r)) return 0;
	return min(min(DistanciaPuntoSegmento(s.p, r),
					 DistanciaPuntoSegmento(s.q, r)),
				 min(DistanciaPuntoSegmento(r.p, s),
					 DistanciaPuntoSegmento(r.q, s)));
}

// Un poligono es una serie de
// vertices conectados por aristas.
// P = p1 -> p2 -> p3 -> ... -> pn -> p1.
typedef vector<Punto> Poligono;

// Saber si un punto esta en el perimetro de un poligono.
bool PuntoEnPerimetro(const Punto& p, const Poligono& P) {
	for (int i = 1; i < P.size(); ++i) {
		if( PuntoEnSegmento( p, Linea( P[i-1], P[i] ) ) ) return true;
	}
	return false;
}

// Prueba de punto en poligono convexo.
// En el perimetro = -1, Fuera = 0, Dentro = 1.
int PuntoEnConvexo(const Punto& p, const Poligono& P) {
	if (PuntoEnPerimetro(p, P)) return -1;
	int dir = ManoDerecha(P[0], P[1], p);
	for (int i = 2; i < P.size(); ++i)
		if (ManoDerecha(P[i - 1], P[i], p)
			!= dir) return 0; // Fuera.
	return 1; // Dentro.
}


// Punto en poligono concavo por ray casting.
// En el perimetro = -1, Fuera = 0, Dentro = 1.
int RayCasting(const Punto& p, const Poligono& P) {
	if (PuntoEnPerimetro(p, P)) return -1;
	Punto o = *min_element(P.begin(), P.end());
	Linea rayo(p, Punto(o.x - M_PI, o.y - M_E));
	int cruces = 0;
	for (int i = 1; i < P.size(); ++i)
		if (InterseccionSegmentos(rayo,
			Linea(P[i - 1], P[i]))) ++cruces;
	return cruces & 1;
}

// Punto en poligono concavo por angle summation.
// En el perimetro = -1, Fuera = 0, Dentro = 1.
int AngleSummation(const Punto& p, const Poligono& P) {
	if (PuntoEnPerimetro(p, P)) return -1;
	Double angulo = 0;
	for (int i = 1; i < P.size(); ++i)
		angulo += ManoDerecha(p, P[i - 1], P[i]) * Angulo(
			P[i - 1], p,P[i]);
	return (fabs(angulo) > M_PI)? 1: 0;
}

// Area de un poligono.
Double Area(const Poligono& P) {
	Double area = 0;
	for (int i = 1; i < P.size(); ++i)
		area += P[i - 1].Cruz(P[i]);
	return fabs(area) / 2.0;
}
// Perimetro de un poligono.
Double Perimetro(const Poligono& P) {
	Double perimetro = 0;
	for (int i = 1; i < P.size(); ++i)
		perimetro += P[i].Dist(P[i - 1]);
	return perimetro;
}

// Cerco convexo de un conjunto de puntos.
template<class It> Poligono CercoConvexo(It lo, It hi) {
	int k = 0;
	if (hi - lo <= 1) return Poligono(lo, hi);    
	Poligono H(2 * (int)(hi - lo));
	sort(lo, hi);

	for (It it = lo; it != hi; ++it) {
		while (k >= 2 && ManoDerecha(H[k - 2], H[k - 1], *it) < 1) k--;
		H[k++] = *it;
	}
	int t = k + 1;
	for (It it = hi - 2; it != lo - 1; --it) {
		while (k >= t && ManoDerecha(H[k - 2], H[k - 1], *it) < 1) k--;
		H[k++] = *it;
	}
	H.resize(k);
	return H;
}
// Centroide de un poligono.
Punto Centroide(const Poligono& P) {
	Double x = 0, y = 0, k = 0;
	for (int i = 1; i < P.size(); ++i) {
		Double cruz = (P[i-1]).Cruz(P[i]);
		x += cruz * (P[i - 1].x + P[i].x);
		y += cruz * (P[i - 1].y + P[i].y);
		k += 3*cruz;
	}
	return Punto(x/k, y/k);
}

// Saber si una recta corta un poligono.
bool RectaCortaPoligono(
	const Linea& r, const Poligono& P) {
	for (int i = 0, prim = 0; i < P.size(); ++i) {
		int lado = ManoDerecha(r.p, r.q, P[i]);
		if (!lado) continue; if (!prim) prim = lado;
		else if (lado != prim) return true;
	}
	return false;
}

// Obtiene los poligonos resultantes de
// cortar un poligono convexo con una recta.
vector<Poligono> CortarPoligono(
	const Poligono& P, const Linea& r) {
	if (!RectaCortaPoligono(r, P))
		return vector<Poligono>(1, P);
	int ind = 0;
	vector<Poligono> Ps(2);
	for (int i = 1; i < P.size(); ++i) {
		Linea s(P[i - 1], P[i]);
		if (IntersecRectaSegmen(r, s)) {
			Punto p = PuntoInterseccion(r, s);
			if (P[i - 1] == p) continue;
			Ps[ind].push_back(P[i - 1]);
			Ps[1 - ind].push_back(p);
			Ps[ind].push_back(p);
			ind = 1 - ind;
		}
		else Ps[ind].push_back(P[i - 1]);
	}
	Ps[0].push_back(Ps[0][0]);
	Ps[1].push_back(Ps[1][0]);
	return Ps;
}











//Dado un poligono convexo, devuelve todos los pares antipodales en orden
vector<pii> RotatingCalipers( Poligono &P ){
	int k = 1, i = P.size()-2, j = 0;
	while( fabs( ( P[j] - P[i] ).Cruz( P[k+1] - P[i] ) ) > 
		fabs( ( P[j] - P[i] ).Cruz( P[k] - P[i] )) ) k++;

	i = 0, j = k;

	vector< pii > res;
	while( i <= k && j < P.size()-1 ){
		res.push_back( pii(i,j) );
		while( j < P.size()-1 && fabs( ( P[i+1] - P[i] ).Cruz( P[j+1] - P[i] ) ) > 
		fabs( ( P[i+1] - P[i] ).Cruz( P[j] - P[i] ) ) ){
			j++;
			res.push_back( pii( i, j ));
		}
		i++;
	}
	return res;
}

//Para ancho de poligono
// vector< pii> tmp = RotatingCalipers( pol );
// tmp.insert(tmp.end(), tmp.begin(), tmp.end() );
// Double res = INF;
// for( int i = 0; i < tmp.size()-1; i++){
//  if( tmp[i].first == tmp[i+1].first && tmp[i].second+1 == tmp[i+1].second ){
//      res = min( res, DistanciaPuntoRecta( pol[tmp[i].first], Linea( pol[tmp[i].second], pol[tmp[i+1].second] ) ) );
//  }
// }

// for( int i = 0; i < tmp.size()-1; i++){
//  if( tmp[i].second == tmp[i+1].second && tmp[i].first+1 == tmp[i+1].first ){
//      res = min( res, DistanciaPuntoRecta( pol[tmp[i].second], Linea( pol[tmp[i].first], pol[tmp[i+1].first] ) ) );
//  }
// }    


// Circulo en 2D.
struct Circulo {
	Punto c; Double r; Circulo() : c(), r() {}
	Circulo(const Punto& C, Double R) : c(C), r(R) {}

	bool operator<(const Circulo& cmp) const {
		if (!(c == cmp.c)) return c < cmp.c;
		return EQ(r, cmp.r)? false: r < cmp.r;
	}
};

// Circunferencia de un circulo.
Double Circuferencia(const Circulo& c) {
	return M_2PI * c.r;
}

// Area de un circulo.
Double Area(const Circulo& c) {
	return M_PI * c.r * c.r;
}

// Saber si un punto esta dentro de un circulo.
// En circunferencia = -1, Fuera = 0, Dentro = 1.
int PuntoEnCirculo(const Punto& p, const Circulo& c) {
	Double dist = p.Dist(c.c);
	if (EQ(dist, c.r)) return -1;
	return (dist < c.r)? 1: 0;
}

// Distancia de un punto p a un circulo c
Double DistanciaPuntoCirculo(const Punto& p, const Circulo& c) {
	Double dist = p.Dist( c.c) - c.r;
	return (dist < 0)? 0: dist;
}

// Proyecta un punto fuera de un circulo en su circunferencia.
Punto ProyPuntoCircunferencia(const Punto& p, const Circulo& c) {
	Punto v = c.c - p;
	Double prop = DistanciaPuntoCirculo(p, c) / v.Abs();
	return p + v*prop;
}




// Obtiene dos puntos que, desde el punto p, forman
// lineas tangentes a la circunferencia del circulo c.
Linea ProyTangentes(const Punto& p, const Circulo& c) {
	Double a = acos(c.r / p.Dist( c.c));
	Punto p_ = ProyPuntoCircunferencia(p, c) - c.c;
	return Linea(c.c + p_.RotarCCW(a),
				 c.c + p_.RotarCW(a) );
}

// Saber si se intersecta un circulo c y una recta r.
// Tangente = -1, No se intersectan = 0, Cuerda = 1.
int IntersecCirculoRecta(const Circulo& c, const Linea& r) {
	Double dist = DistanciaPuntoRecta(c.c, r);
	if (EQ(dist, c.r)) return -1;
	return (dist < c.r)? 1: 0;
}

// Cuerda de interseccion entre un circulo y una recta.
vector<Punto> CuerdaInterseccion( const Circulo& c, const Linea& r ){
	vector<Punto> ret;
	Punto b = r.q - r.p, a = r.p - c.c;
	Double A = b.Dot(b);
	Double B = a.Dot(b);
	Double C = a.Dot(a) - c.r*c.r;
	Double D = B*B - A*C;
	if( D + ERROR < 0 ) return ret;
	ret.push_back( c.c + a + b * ((-B+sqrt(D+ERROR)) / A) );
	if( D > ERROR )
		ret.push_back( c.c + a + b * ((-B-sqrt(D)) / A) );
	return ret;
}

// Saber si un circulo c esta dentro de un circulo d.
bool CirculoEnCirculo(const Circulo& c, const Circulo& d) {
	return c.c.Dist(d.c) < d.r - c.r;
}

// Saber si el circulo c intersecta con el circulo d.
// Uno dentro del otro = -1, Disjuntos = 0, Intersectan = 1.
int IntersecCirculoCirculo(const Circulo& c, const Circulo& d) {
	Double dist = c.c.Dist(d.c);
	if (dist < fabs(c.r - d.r)) return -1;
	return (dist > c.r + d.r)? 0: 1;
}

//Dados 3 puntos, calcula el centro del circulo que los contiene.
Punto CalcularCentroCirculo( Punto a, Punto b, Punto c ){
	b = ( a + b ) * 0.5;
	c = ( a + c ) * 0.5;
	return PuntoInterseccion( Linea( b, b + ( a-b ).RotarCW(M_PI/2.0)), Linea( c, c+(a-c).RotarCW( M_PI/2.0 ) ) );
}






//Puntos de interseccion entre 2 circulos dados.
//Regresa un vector con los puntos.
vector<Punto> InterseccionCirculoCirculo( const Circulo& C, const Circulo& D ){
	vector<Punto> ret;
	Double d = ( C.c - D.c ).Abs();
	if( d > C.r + D.r || d + min( C.r, D.r ) < max( C.r, D.r ) ) return ret;
	Double x = ( d*d - D.r*D.r + C.r*C.r ) / (2*d);
	Double y = sqrt( C.r*C.r - x*x );
	Punto v = ( D.c - C.c )/d;
	ret.push_back( C.c+v*x + v.rot90()*y );
	if( y > 0 )
		ret.push_back( C.c+v*x - v.rot90()*y );
	return ret;
}

// Obtiene tangentes exteriores (las que NO se cruzan) entre dos circulos.
int TangenteExtCirculoCirculo(const Circulo& a, const Circulo& b, Linea &s, Linea &t) {
	// Circulos identicos. Tangentes infinitas (segun el problema).
	if ( EQ(a.r, b.r) && a.c == b.c) 
		return 0; // O etiqueta de caso especial.
	// Uno es circulo interior del otro. Comparten una tangente.
	// EL CALCULO PUEDE COPIARSE A TangenteInt SI SE REQUIERE.
	Punto u;
	bool unico = false;
	if (b.r < a.r && EQ(a.c.Dist(b.c) + b.r, a.r)) {
		u = a.c + (b.c-a.c) * (a.r / a.c.Dist(b.c));
		unico = true;
	}
	if (a.r < b.r && EQ(a.c.Dist(b.c) + a.r, b.r)) {
		u = b.c + (a.c-b.c)*(b.r / a.c.Dist(b.c));
		unico = true;
	}
	if (unico) {
		s = t = PerpendicularEnPunto(Linea(a.c, b.c), u); // Recta de tangencia; un punto es referencia.
		//s = t = Linea(u, u); // Punto de tangencia.
		return 1;
	}
	// Circulo en circulo. No hay tangentes.
	if (CirculoEnCirculo(a, b) || CirculoEnCirculo(b, a)) 
		return 0;

	// Calcular las 2 rectas tangentes.
	Linea proy;
	Punto v;
	if (EQ(a.r, b.r)) {
		proy = Linea(a.c, a.c);
		Linea perp = PerpendicularEnPunto(Linea(a.c, b.c), a.c);
		u = (perp.p-perp.q)* (b.r / perp.p.Dist(perp.q)); 
		v = -u;
	} else {
		Circulo c(a.c, fabs(a.r - b.r));
		proy = ProyTangentes(b.c, c);
		u = (proy.p-a.c)* (b.r / (a.r - b.r));
		v = (proy.q-a.c)* (b.r / (a.r - b.r));
	}
	s = Linea( proy.p+u, b.c+u);
	t = Linea( proy.q+v, b.c+v);
	return 2;
}


// Obtiene tangentes interiores (las que SI se cruzan) entre dos circulos.
int TangenteIntCirculoCirculo(const Circulo& a, const Circulo& b, Linea &s, Linea &t) {
	// Circulos identicos. Tangentes infinitas (segun el problema)
	if (EQ(a.r, b.r) && a.c == b.c) {
		return 0; // O etiqueta de caso especial.
	}
	// Uno es circulo interior del otro. Comparten una tangente.
	// CALCULO HECHO EN TangenteInt. Copiar de este si se requiere.
	Punto u;
	// Circulos tangentes. Obtener recta tangente unica.
	if (EQ(a.c.Dist(b.c), a.r + b.r)) {
		u = b.c + (a.c-b.c) * (b.r / a.c.Dist(b.c));
		s = t = PerpendicularEnPunto(Linea(a.c, b.c), u); // Recta de tangencia.
		//s = t = Linea(u, u); // Punto de tangencia.
		return 1;
	}
	// Circulos se traslapan. no hay tangentes.
	if (!(a.r + b.r < a.c.Dist(b.c))) {
		return 0;
	}
	// Obtener 2 rectas tangentes.
	Linea proy;
	Punto v;
	Circulo c(a.c, a.r + b.r);
	proy = ProyTangentes(b.c, c);
	u = (proy.p - a.c) * (b.r / (a.r + b.r));
	v = (proy.q - a.c) * (b.r / (a.r + b.r));
	s = Linea( proy.p - u, b.c-u);
	t = Linea( proy.q - v, b.c-v);
	return 2;
}

int main(){
	int N;

	cin >> N;

	Poligono pol(N);
	for( int i = 0; i < N; i++ ){
		cin >> pol[i].x >> pol[i].y;
	}

	pol = CercoConvexo( pol.begin(), pol.end() );
	vector<pii> tmp = RotatingCalipers( pol );

	Double distancia = 0.0;
	for( int i = 0; i < tmp.size(); i++ ){
		distancia = max( distancia, pol[tmp[i].first].Dist( pol[tmp[i].second] ) );
	}
	cout << fixed << setprecision( 7 ) <<distancia << '\n';
	return 0;
}

Compiler output:
a.cc: In function 'std::pair<Punto, Punto> ParMasCercano(std::vector<Punto>)':
a.cc:98:9: error: 'sort' was not declared in this scope; did you mean 'sqrt'?
   98 |         sort(P.begin(), P.end());
      |         ^~~~
      |         sqrt
a.cc: In function 'int RayCasting(const Punto&, const Poligono&)':
a.cc:297:20: error: 'min_element' was not declared in this scope
  297 |         Punto o = *min_element(P.begin(), P.end());
      |                    ^~~~~~~~~~~
a.cc: In instantiation of 'Poligono CercoConvexo(It, It) [with It = __gnu_cxx::__normal_iterator<Punto*, std::vector<Punto> >; Poligono = std::vector<Punto>]':
a.cc:652:20:   required from here
  652 |         pol = CercoConvexo( pol.begin(), pol.end() );
      |               ~~~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~
a.cc:337:13: error: 'sort' was not declared in this scope; did you mean 'sqrt'?
  337 |         sort(lo, hi);
      |         ~~~~^~~~~~~~
      |         sqrt