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 21000. Lo que hago es esto:
#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;
}
21000 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 21000 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.
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:
mpz_init(num)
Oye una cosa¿ la condición del while en pseudocódigo cual es?
Muchas gracias por responder.
Lo he inicializado, solo que dándole otro valor distinto de 0. Así:
mpz_set_ui(num, 2); //Inicializo num a 2
Además, debe estar bien, porque esa variable la utilizo para hacer 21000 y dicho resultado me da bien.
Saludos.
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í :
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)
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:
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...
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 21000¿no?
Sí, sería sumando los dígitos de 21000, pero ahora me he quedado con la curiosidad de por qué sale mal lo que he hecho... Seguiré probando, a ver si saco algo...
Estoy mirando el código y ten mucho cuidado con las funciones que utilizas , mira la función esta lo que retorna:
mpz_cdiv_ui(resultado, 10)
No retorna el número dividido , sino el resto de la división...
Ya está listo dame un momento y te publico lo que está mal...
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?
#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!
¡Muchas gracias! Creía que era esa función la de dividir...
La suma no da... Así que he puesto esto:
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_tdiv_q_ui(resultado, 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.
getchar();
}
Para ver qué pasaba, y los primeros números los suma bien, pero a partir de 565 sale un número muy grande, números negativos, etc. (no sé por qué salen negativos, cuando "suma" es unsigned long int...)...
Pues es que unsigned long int , es de 4 bytes de longitud 4*8= 32 bits , o sea que solo puede contener hasta 232=4,294,967,295 , pero eso no tiene nada que ver porque esas cifras no suman tal cantidad de numeros, así que no sé , cuanto debe de dar la suma entonces?
El problema es de una página muy conocida de retos de programación, por lo que buscando en Internet acabo de ver que da 1366. No sé si será el resultado correcto...
Eso me dá a mí , a ti ¿cuanto te dá? ¿Que página es?
Qué raro... A mí me da 1620828886. La página es Project Euler.
A ver vamos a hacer una cosa , vamos a guardarla en la variable suma pero está ahora vá a ser de tipo mpz_t ahora te paso el codigo a ver que tal , pero me parece rarísimo.Toma aquí tienes , ahora lo que he hecho es que la variable suma es de tipo mpz_t
#include <stdio.h>
#include <gmp.h>
int main()
{
mpz_t num, resultado, mod,suma; //Variables
mpz_inits(resultado,mod,suma,NULL);//Resultado suma 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
mpz_add_ui(suma,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("Suma: %Zd\n\n", suma); //Además de que el bucle es infinito, he comprobado que la suma la hace mal.
//getchar();
}
gmp_printf("Suma final: %Zd\n\n", suma);
mpz_clears(resultado,num,mod,suma,NULL); //Liberamos memoria
getchar();
return 0;
}
Vamos a ver una cosa compila este código y dime lo que te sale también...
#include <stdio.h>
int main()
{
printf("El tamano de un unsigned long int en tu maquina es de: %d bytes",sizeof(unsigned long int));
getchar();
return 0;
}
Tanto en tu código como en uno que he hecho yo (que ha resultado ser igual que el tuyo), me da 43123298370... :(
Me sale 4 bytes.
Pues ni idea la verdad , igual es que mis librerias como las compile para mi ordenador no sirven para el tuyo...
Yo creo que sí sirven, porque sino me daría error al compilar, o me daría mal 21000... Veré que puedo hacer.
Gracias por haberme ayudado!
Saludos.
No sé yo eh , sobre que Windows estás compilando x64 o x32 ? Yo compilé mis librerias desde x64 a ver si no sirven , compila las tuyas.
Desde Windows 7 x64
Yo que tú compilaba las librerías , con el tutorial ese que posteé.
¡Mira lo que me ha pasado ahora! :-(
Para probar si me hacía la división bien, he puesto:
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_tdiv_q_ui(resultado, resultado, 10); //Dividimos resultado entre 10
gmp_printf("Resultado: %Zd\n\n", resultado);
getchar();
}
(el getchar() para pausar y poder ver que todo va bien)
Y con eso me ha ido haciendo esto:
Resultado: 107150860718626732094842504906000181056140481170553360744375038837035105112493612
24931983788156958581275946729175531468251871452856923140435984577574698574803934567774824230
98542107460506237114187795418215304647498358194126739876755916554394607706291457119647768654
2167660429831652624386837205668069376
Resultado: 107150860718626732094842504906000181056140481170553360744375038837035105112493612
24931983788156958581275946729175531468251871452856923140435984577574698574803934567774824230
98542107460506237114187795418215304647498358194126739876755916554394607706291457119647768654
216766042983165262438683720566806937
Resultado: 107150860718626732094842504906000181056140481170553360744375038837035105112493612
24931983788156958581275946729175531468251871452856923140435984577574698574803934567774824230
98542107460506237114187795418215304647498358194126739876755916554394607706291457119647768654
21676604298316526243868372056680693
Resultado: GNU MP: Cannot reallocate memory (old_size=291 new_size=300)
This application has requested the Runtime to terminate it in an unusual way.
Please contact the application's support team for more information.
No creo que sea por lo de la librería (ya que las operaciones me las hace bien, o hasta cierto punto), pero igualmente voy a intentar compilarla, por si acaso...
EDITO: Ya he compilado la librería por mí mismo y sigue saliendo mal...
Creo que lo mejor será dejarlo, porque si en un ordenador sale y en otro no... :-(
Pues es muy raro , mira intenta inicializar la variable resultado así , mpz_init2(resultado,10000000);
Tampoco va... Sigo probando cosas y leyendo un poco el manual, pero dudo que venga algo allí, porque si el programa está bien...
Voy a pasarte el exe , si no funciona , está claro que no son las librerías y que tu SO es el problema supongo.
Aquí lo tienes: http://www.mediafire.com/?fpt5yyatz64yz47
Siento si te he molestado ya mucho, pero... No me deja descargarlo :(
Que no molestas hombre :) , será el antivirus , prueba con este:
http://www.mediafire.com/?mbw2v6ywhndpstv
Me dice:
"El programa no puede iniciarse porque falta libgcc_s_dw2-1.dll en el equipo. Intente reinstalar el programa para corregir este problema." :huh:
No sabes como odio los problemas con las .dll's :xD prueba con este que lleva la dll:
http://www.mediafire.com/?kld0s4vbqd22e6e
Muchas gracias. No da el resultado correcto... Por lo que no debería ser del código, y me he llevado toda la tarde mirándolo xD
Lo único que se me ocurre es que intentes compilar una versión anterior de la librería. A ver si va a ser un bug...
Si fuera un bug a ti también te tendría que salir mal, supongo...
¿Y si depende del ordenador o algo de eso? Es que no te puedo orientar más , hay muy poca información al respecto y yo no sé mucho ... De todas maneras prueba con esta que está precompilada http://www.cs.nyu.edu/exact/core/gmp/gmp-static-mingw-4.1.tar.gz , tienes que hacer lo mismo que con la que instalaste , pero quita la que te pasé yo antes!
EDITADO: Es posible que si intentas compilar tu código modificado por mi con esa librería , obtengas errores refiriendose a la funcion mpz_inits y mpz_clears supongo que en antiguas versiones no estaba implementada.
¡Muchísimas gracias, Avesudra! Me has ayudado muchísimo! Debe ser un bug como bien has dicho ya, porque ahora sí me funciona!!
Reitero, muchísimas gracias!!
Saludos.
¡Muchas de nadas! Debe de ser un bug que dependa del tipo de hardware del ordenador(cpu , cantidad de ram...) a lo mejor, es que no tengo ni idea la verdad, porque si mis librerías estuviesen mal compiladas no me saldría a mí tampoco... ¡Pues bueno problema solucionado! ¡A programar tranquilamente!
¡Un saludo crack!
Eso sí, no me gustaría que le pasara a otro, porque yo mismo he vivido este sufrimiento que se tiene cuando ves que una cosa está bien y no sale lo que debe salir.
¿Ves bien que intente contactar a los creadores de la librería para informar del bug?
Eso te iba a decir , no estaría mal , pero hay una especie de documentación para reportarlo mira http://gmplib.org/manual/Reporting-Bugs.html , una cosa sería mejor que comentasemos esto ya por privado porque el tema principal del hilo se ha zanjado.