¿Uso mal las funciones de la librería GMP?

[DickGumshoe]

Hola.

Gracias a Avesudra, ayer pude configurar la librería GMP en Codeblocks. Llevo toda la mañana intentando calcular la suma de todos los dígitos del resultado de hacer 2¹⁰⁰⁰. Lo que hago es esto:

Código (c) [Seleccionar] Expandir

#include <stdio.h>
#include <gmp.h>

unsigned long int suma=0, i=0;
//Lo pongo como variable global porque si lo meto en main() me dice "Suma dígitos.exe dejó de funcionar. Windows está buscando una solución"

int main()
{

    mpz_t num, resultado, mod; //Variables

    mpz_set_ui(num, 2); //Inicializo num a 2
    mpz_init(resultado); //resultado = 0
    mpz_init(mod); //mod=0

    mpz_pow_ui(resultado, num, 1000); // 2^1000
    gmp_printf("Resultado: %Zd\n\n", resultado); //Imprimo 2^1000, y el resultado es correcto

    while(mpz_cmp_ui(resultado, 0) > 0)
    {//BUCLE INFINITO
        suma += mpz_mod_ui(mod, resultado, 10); //suma = suma + el resto de dividir resultado entre 10
        mpz_set_ui(resultado, mpz_cdiv_ui(resultado, 10)); //Dividimos resultado entre 10
        printf("%d\n", suma); //Además de que el bucle es infinito, he comprobado que la suma la hace mal.

    }


    printf("Suma: %d", suma);

    mpz_clear(resultado); //Liberamos memoria
    mpz_clear(num);
    mpz_clear(mod);

    return 0;
}


2¹⁰⁰⁰ lo hace bien, pero se queda en el while(), no sé por qué...

Mi idea era inicializar suma a 0, hallar el resto de dividir 2¹⁰⁰⁰ entre 10, y sumárselo, y así hasta que "resultado = 0", puesto que si llega a 0 significa que ya hemos sumado todas sus cifras.

¿Alguien tiene una idea de qué hago mal?

Muchas gracias.

Saludos.

[avesudra]

Emmm... me parece que se te ha olvidado inicializar mpz_t num , no sé si ese es el único error , lo he mirado así por encima.

EDITO: Tienes que inicializar las variables antes de asignarles cualquier otro tipo de datos así que te falta:

Código [Seleccionar] Expandir

mpz_init(num)
Oye una cosa¿ la condición del while en pseudocódigo  cual es?

[DickGumshoe]

Muchas gracias por responder.

Lo he inicializado, solo que dándole otro valor distinto de 0. Así:

Código (c) [Seleccionar] Expandir

mpz_set_ui(num, 2); //Inicializo num a 2

Además, debe estar bien, porque esa variable la utilizo para hacer 2¹⁰⁰⁰ y dicho resultado me da bien.

Saludos.

[avesudra]

Ah si perdona es que no había visto las funciones de inicialización y asignación combinadas en la documentación , de todas maneras es así :

Código (c) [Seleccionar] Expandir

mpz_init_set_ui (num, 2);
En cuanto al bucle que condición es la que quieres poner, que resultado sea igual a 0?
entonces sería así while(mpz_cmp_ui(resultado,0)!=0)

[DickGumshoe]

Muchas gracias por corregir esa función, con la que había puesto yo tenía que crear la variable resultado también porque sino me daba error.

Sí, lo que quiero poner en el while() es que mientras resultado > 0, haga todo eso.

Saludos.

EDITO: También puse esa condición y no salía...

EDITO 2: Creo que el error no es la condición, sino lo de dentro, ya que he probado a poner:

Código (c) [Seleccionar] Expandir

while(mpz_cmp_ui(resultado, 0) !=0)
    {//BUCLE INFINITO
        suma += mpz_mod_ui(mod, resultado, 10); //suma = suma + el resto de dividir resultado entre 10
        mpz_set_ui(resultado, mpz_cdiv_ui(resultado, 10)); //Dividimos resultado entre 10
        gmp_printf("%Zd\n", resultado); //Además de que el bucle es infinito, he comprobado que la suma la hace mal.

    }

Y no da nunca resultado/10...

[avesudra]

Ya he visto que no me funciona ...Voy a ver si te lo hago desde 0 , entonces recapitulamos , tenemos que sumar todos los dígitos de 2¹⁰⁰⁰¿no?

[DickGumshoe]

Sí, sería sumando los dígitos de 2¹⁰⁰⁰, pero ahora me he quedado con la curiosidad de por qué sale mal lo que he hecho... Seguiré probando, a ver si saco algo...

[avesudra]

Estoy mirando el código y ten mucho cuidado con las funciones que utilizas , mira la función esta lo que retorna:

Código [Seleccionar] Expandir

mpz_cdiv_ui(resultado, 10)
No retorna el número dividido , sino el resto de la división...

[avesudra]

Ya está listo dame un momento y te publico lo que está mal...

[avesudra]

Pues mira te explico lo que tenías mal:

• Primero de todo la función mpz_set_ui(num,2); no inicializaba resultado , era mpz_init_set_ui(num,2)
• Segundo la función mpz_cdiv_ui(resultado, 10); no retornaba el número dividido si no el resto de la división
• Tercero el especificador de conversion para un tipo unsigned long int es %lu

Y creo que no se me olvida más nada , mmm ah si el codigo ! ¿Que valor daba la suma?

Código (c) [Seleccionar] Expandir

#include <stdio.h>
#include <gmp.h>

unsigned long int suma=0, i=0;
//Lo pongo como variable global porque si lo meto en main() me dice "Suma dígitos.exe dejó de funcionar. Windows está buscando una solución"

int main()
{

   mpz_t num, resultado, mod; //Variables
   mpz_inits(resultado,mod,NULL);//Resultado y mod inicializado a 0;
   mpz_init_set_ui(num, 2); //Inicializo num a 2
   mpz_pow_ui(resultado, num, 1000); // 2^1000
   gmp_printf("Resultado: %Zd\n\n", resultado); //Imprimo 2^1000, y el resultado es correcto

   while(mpz_cmp_ui(resultado, 0) > 0)
   {//BUCLE NO INFINITO
       suma += mpz_mod_ui(mod, resultado, 10); //suma = suma + el resto de dividir resultado entre 10
       //printf("SUMA :%d \n", suma); Probando la suma...
       mpz_tdiv_q_ui(resultado,resultado, 10);//Realizamos una division entre 10 y no redondeamos, lo que sobra lo truncamos.
       //mpz_set_ui(resultado, mpz_cdiv_ui(resultado, 10)); Esto está mal , la función mpz_cdiv_ui devuelve el resto de la division!!
       gmp_printf("Resultado dividiendo: %Zd\n\n", resultado);
       printf("%lu\n", suma); //Además de que el bucle es infinito, he comprobado que la suma la hace mal.

   }
   printf("Suma: %lu\n", suma);
   mpz_clears(resultado,num,mod,NULL); //Liberamos memoria
   getchar();
   return 0;
}

¡Un saludo!