X-Git-Url: https://git.llucax.com/z.facultad/75.29/dale.git/blobdiff_plain/3939466e952e19024a4856cfa7c6ca2b07bc297a..ed4673952ba5ec554d3a045cfeced013644919e7:/src/number.h?ds=inline

diff --git a/src/number.h b/src/number.h
index a079e8f..27d19a8 100644
--- a/src/number.h
+++ b/src/number.h
@@ -50,21 +50,21 @@ struct number
 	// Constructores (después de construído, el chunk siempre tiene al
 	// menos un elemento).
 	// Constructor default (1 'átomo con valor 0)
-	number(): chunk(1, 0) {}
+	number(): chunk(1, 0), sign(positive) {}
 
 	// Constructor a partir de buffer (de 'átomos') y tamaño
 	// Copia cada elemento del buffer como un 'átomo' del chunk
 	// (el átomo menos significativo es el chunk[0] == buf[0])
-	number(native_type* buf, size_type len, sign_type sign = positive):
-		chunk(buf, buf + len), sign(sign)
+	number(native_type* buf, size_type len, sign_type s = positive):
+		chunk(buf, buf + len), sign(s)
 	{
 		fix_empty();
 	}
 
 	// Constructor a partir de un 'átomo' (lo asigna como único elemento
 	// del chunk). Copia una vez N en el vector.
-	number(native_type n, sign_type sign = positive):
-		chunk(1, n), sign(sign) {}
+	number(native_type n, sign_type s = positive):
+		chunk(1, n), sign(s) {}
 
 	number(const std::string& str);
 
@@ -80,7 +80,7 @@ struct number
 	number& operator<<= (const size_type n);
 	number& operator-= (const number& n);
 	bool	operator< (const number& n);
-	bool operator==(const number& n) const;
+	bool    operator==(const number& n) const;
 
 	// Devuelve referencia a 'átomo' i del chunk (no debería ser necesario
 	// si la multiplicación es un método de este objeto).
@@ -132,19 +132,18 @@ struct number
 	// borrow)
 	void borrow(size_type i)
 	{
-		if (chunk.size() >= i)
+		// para poder pedir prestado debo tener uno a la izquierda
+		assert (chunk.size() >= i);
+		
+		if (chunk[i] == 0)
 		{
-			if (chunk[i] == 0)
-			{
-				borrow(i+1); // Overflow, pido prestado
-				chunk[i] = ~((N)0); //quedo con el valor máximo
-			}
-			else
-			{
-				--chunk[i]; //tengo para dar, pero pierdo uno yo
-			}
+			borrow(i+1); // Overflow, pido prestado
+			chunk[i] = ~((N)0); //quedo con el valor máximo
+		}
+		else
+		{
+			--chunk[i]; //tengo para dar, pero pierdo uno yo
 		}
-		//else ERROR, están haciendo a-b con a>b
 	}
 	// Verifica si es un número par
 	bool es_impar() const
@@ -202,6 +201,23 @@ number< N, E >::number(const std::string& origen)
 template < typename N, typename E >
 number< N, E >& number< N, E >::operator+= (const number< N, E >& n)
 {
+	// Si tienen distinto signo, restamos...
+	if (sign != n.sign)
+	{
+		if (sign == positive) // n es negativo
+		{
+			number< N, E > tmp = n;
+			tmp.sign = positive;
+			*this -= tmp;
+		}
+		else // n es positivo, yo negativo
+		{
+			sign = positive;
+			*this = n - *this;
+		}
+		return *this;
+	}
+
 	native_type c = 0;
 	size_type ini = 0;
 	size_type fin = std::min(chunk.size(), n.chunk.size());
@@ -297,7 +313,6 @@ number< N, E >& number< N, E >::operator-= (const number< N, E >& n)
 		minuend.sign = positive;
 	}
 
-	native_type c = 0;
 	size_type ini = 0;
 	size_type fin = std::min(minuend.chunk.size(), subtrahend.chunk.size());
 	size_type i; //problema de VC++, da error de redefinición
@@ -310,11 +325,17 @@ number< N, E >& number< N, E >::operator-= (const number< N, E >& n)
 		// si no alcanza para restar pido prestado
 		if ((minuend.chunk[i] < subtrahend.chunk[i]))
 		{
-			minuend.borrow(i);
+			// no puedo pedir si soy el más significativo ...
+			assert (i != fin);
+
+			// le pido uno al que me sigue
+			minuend.borrow(i+1);
 		}
 		
-		// resto el chunk i-ésimo
-		minuend.chunk[i] -= subtrahend.chunk[i];
+		// es como hacer 24-5: el 4 pide prestado al 2 (borrow(i+1)) y después
+		// se hace 4 + (9-5) + 1
+
+		minuend.chunk[i] += (~((N)0) - subtrahend.chunk[i]) + 1; 
 	}
 
 	//retorno el minuendo ya restado
@@ -381,17 +402,17 @@ number< N, E > operator<< (const number< N, E >& n, typename number< N, E >::siz
 	return tmp;
 }
 
-template < typename N, typename E >
-std::ostream& operator<< (std::ostream& os, const number< N, E >& n)
+template < typename NN, typename EE >
+std::ostream& operator<< (std::ostream& os, const number< NN, EE >& n)
 {
 	// FIXME sacar una salida bonita en ASCII =)
 	if (n.sign == positive)
 		os << "+ ";
 	else
 		os << "- ";
-	for (typename number< N, E >::const_reverse_iterator i = n.chunk.rbegin();
+	for (typename number< NN, EE >::const_reverse_iterator i = n.chunk.rbegin();
 			i != n.chunk.rend(); ++i)
-		os << std::setfill('0') << std::setw(sizeof(N) * 2) << std::hex
+		os << std::setfill('0') << std::setw(sizeof(NN) * 2) << std::hex
 			<< *i << " ";
 	return os;
 }
@@ -540,8 +561,7 @@ number < N, E > naif(const number< N, E > &u, const number< N, E > &v)
 
 	sign_type sign;
 
-	if ( (u.sign == positive && v.sign == positive) ||
-			(u.sign == negative && v.sign == negative) ) {
+	if (u.sign == v.sign) {
 		sign = positive;
 	} else {
 		sign = negative;
@@ -621,8 +641,7 @@ number < N, E > karatsuba(const number< N, E > &u, const number< N, E > &v)
 
 	sign_type sign;
 
-	if ( (u.sign == positive && v.sign == positive) ||
-			(u.sign == negative && v.sign == negative) ) {
+	if (u.sign == v.sign) {
 		sign = positive;
 	} else {
 		sign = negative;
@@ -631,7 +650,7 @@ number < N, E > karatsuba(const number< N, E > &u, const number< N, E > &v)
 	if (chunk_size == 1) {
 		E tmp;
 		tmp = static_cast< E >(u.chunk[0]) * static_cast< E >(v.chunk[0]);
-		num_type tnum = num_type(static_cast< N* >(&tmp), 2, sign);
+		num_type tnum = num_type(reinterpret_cast< N* >(&tmp), 2, sign);
 		return tnum;
 	}