Petando la pila. Problemas y retos usando recursividad.

[kutcher]

Iniciado por leosansan

¡¡¡ Duda, duda ¡¡¡. ¿Y cúal es la duda?. Sencillo el valor que se obtiene supera los límites del tipo "int" y como los resultados obtenidos con el código son correctos me lleva a plantearme dos posibilidades:

** El programa trunca el valor, pero al hacerlo tanto en la cifra del número como en la del factor donde aparece dicho dígito, se compensan los errores.

** El programa al encontrarse un valor tan alto pasa a tomarlo como float, manteniendo todas sus cifras, y tampoco se produce error.

Excelente observación leosansan y ahora paso a explicar porque funciona,
como pueden notar nos encontramos ante un comportamiento indefinido en el lenguaje C en este caso con respecto al operador de desplazamiento a la izquierda:

"El resultado es indefinido si el operando de la derecha es negativo, o mayor o igual que el número de bits en el tipo de la expresión izquierda."

Como sabemos al desplazar:

Código C++:

1 << 54
 

Se espera como resultado el cero, debido a que los bits que salen por la izquierda se pierden, los que entran por la derecha se rellenan con ceros pero en la expresión que utilizo eso no sucede ocurre una especie de rotación aquí dejo una pagina que explica porque pasa eso en los procesadores Intel arquitectura x86:

http://choorucode.com/2004/12/20/c-shift-operator-mayhem/

Ocurre el siguiente desplazamiento con el ejemplo de leosansan el numero 39:

Código C++:

1 << (39 % 10) * 6 = 1 << 9 * 6 = 1 << 54 = 4194304
 

Código C++:

00000000 01000000 00000000 00000000
 

Y en vista de todo esto expongo otra opción pero utilizando otro método para identificar los dígitos, como van a ver hago uso del producto de los números primos:

Código C++:

#define MAX(a,b) ((a)>(b)?(a):(b))
#define MIN(a,b) ((a)<(b)?(a):(b))
 
int dec[] = { 1, 10, 100, 1000, 10000, 100000 };
int primos[10] = {2, 3, 5, 7, 11, 13, 17, 23, 29, 31};
 
int fact(int num, int t)
{
    if (!num) return t;
    t *= primos[num % 10];
    return fact(num/10, t);
}
int comprobar(int min, int max, int num, int b, int t)
{
    if(min >= max) return 0;
 
    b = num / min;
    if (min * b == num && ((min % 10) || (b%10))
        && t == fact(min,1) * fact(b,1) )
        return 1;
 
    return comprobar(min + 1, max, num, b, t);
}
int esVampiro(int num)
{
    int min, max, b = 0, nd = log10(num)+1, t = fact(num, 1);
 
    if (nd % 2) return 0;
    nd /= 2;
    min = MAX(dec[nd - 1], (num + dec[nd] - 2)/(dec[nd] - 1));
    max = MIN(num/min, sqrt(num));
 
    return comprobar(min, max, num, b, t);
}
 

Un saludo