Serabe Reloaded » RubyEuler http://www.serabe.com Conecto ergo sum. Non conecto ergo urgueo. Sun, 09 Oct 2011 21:30:53 +0000 en hourly 1 http://wordpress.org/?v=3.2.1 Problem 3 http://www.serabe.com/2008/10/12/problem-3/ http://www.serabe.com/2008/10/12/problem-3/#comments Sun, 12 Oct 2008 14:13:42 +0000 Serabe http://www.serabe.com/?p=259 Wording: (Original) The prime factors of 13195 are 5, 7, 13 and 29.
What is the largest prime factor of the number 600851475143 ?

Solution:
First the code:

PLAIN TEXT
RUBY:
  1. include Math
  2. def prime_factors(num, factor=2)
  3.   return []  if num <= 1
  4.   next_pf = (factor..(sqrt(num).ceil)).find(lambda {num}){ |x| num%x == 0 }
  5.   return [next_pf] + prime_factors(num/next_pf, next_pf)
  6. end
  7.  
  8. puts prime_factors(600851475143).max

The code is a direct port from the one in PyEuler. It is based mainly in the idea that the first factor of a number is always prime.

Line 3 breaks recursion and line 4 works out the next prime factor (actually, it finds the smallest factor). For doing this, it uses find method. If no element match the criteria, then it returns the number (that is what the lambda block is for).

Finally, line 5 merges the results and line 8 print the result.

This time just one code is showed.

article clipper Problem 3
 

share save 171 16 Problem 3

]]> http://www.serabe.com/2008/10/12/problem-3/feed/ 0 Problema 2 http://www.serabe.com/2008/03/23/problema-2/ http://www.serabe.com/2008/03/23/problema-2/#comments Sun, 23 Mar 2008 21:37:14 +0000 Serabe http://www.serabe.com/2008/03/23/problema-2/ Enunciado:
Suma todos los números pares de la sucesión de Fibonacci menores que un cuatro millón.

Primera solución:

La primera solución es muy simple. Hay un método que devuelve un array con todos los números de la sucesión de Fibonacci menores que un máximo max tomando como inicio un array de dos elementos arr.

PLAIN TEXT
RUBY:
  1. def fib1(max,arr=[1,1])
  2.   while(arr.last <max) do
  3.     arr <<(arr.last+arr[arr.size-2])
  4.   end
  5.   arr.pop
  6.   arr
  7. end

Después simplemente se eliminan los impares y se suman los que quedan.

PLAIN TEXT
RUBY:
  1. def euler2a(max, ini=[1,1])
  2.   fib1(max,ini).delete_if{|x| x%2 == 1}.inject{|memo,obj| memo+obj}
  3. end
  4.  
  5. puts euler2a(4e6)

Segunda solución:

Esta solución es un poco mejor. Se va a hacer un nuevo método que calcule todos los números de la sucesión de Fibonacci menores que un máximo max pero que sólo almacene aquellos que pasen una condición que se le pasa como bloque:

PLAIN TEXT
RUBY:
  1. def fib2(max,a=1,b=1)
  2.   arr = []
  3.   arr <<a if yield(a)
  4.   while(b<max) do
  5.     arr <<b if yield(b)
  6.     a,b=b,a+b
  7.   end
  8.   arr
  9. end

Después sólo hace falta pasar el filtro adecuado, y sumar los resultados:

PLAIN TEXT
RUBY:
  1. def euler2b(max,a=1,b=1)
  2.   fib2(max,a,b){|x| x%2==0}.inject{|memo,obj| memo+obj}
  3. end
  4.  
  5. puts euler2b(4e6)

Tercera solución:

Esta tercera solución es una modificación directa de la segunda. En vez de almacenar los datos, se suman directamente.

PLAIN TEXT
RUBY:
  1. def fib3(max,a=1,b=1)
  2.   res = 0
  3.   res += a if yield(a)
  4.   while(b<max) do
  5.     res += b if yield(b)
  6.     a,b=b,a+b
  7.   end
  8.   res
  9. end
  10.  
  11. def euler2c(max,a=1,b=1)
  12.   fib3(max,a,b){|x| x%2==0}
  13. end
  14.  
  15. puts euler2c(4e6)

Y eso es todo.

Actualización: Hay diferencias entre la página de PyEuler y la del Proyecto Euler. El problema ha sido actualizado para corresponderse con esta última.

article clipper Problema 2
 

share save 171 16 Problema 2

]]> http://www.serabe.com/2008/03/23/problema-2/feed/ 0 Problema 1 http://www.serabe.com/2008/02/27/problema-1/ http://www.serabe.com/2008/02/27/problema-1/#comments Wed, 27 Feb 2008 12:20:56 +0000 Serabe http://www.serabe.com/2008/02/27/problema-1/ Empiezo aquí una serie de artículos de periodicidad variable que lo único que pretende es ser algo similar (de una manera muy amplia) a PyEuler. La lista completa de problemas está en el Proyecto Euler y, como no podía ser de otra manera, empiezo por el primero.

Enunciado:
Halla la suma de todos los números menores que 1000 y múltiplos de 3 o de 5.

Solución 1:

PLAIN TEXT
RUBY:
  1. def euler1a(e, numbers)
  2.   (1..e).select{|x| numbers.any?{|y| (x%y == 0)}}.inject{|memo,o| memo+=o}
  3. end
  4.  
  5. puts euler1a(1000,[3,5])

La explicación es simple. Tenemos un rango desde 1 hasta el límite, filtramos los que son múltiplos de algún elemento del array y después se suman.

Solución 2:

PLAIN TEXT
RUBY:
  1. def euler1b(e,numbers)
  2.   (1..e).inject(0){|memo,o| (numbers.any?{|x| (o%x)==0}) ? memo+o : memo }
  3. end
  4.  
  5. puts euler1b(1000,[3,5])

Éste es similar al anterior, sólo que se suma según se recorre el rango.

Solución 3:
Para este necesitamos antes un par de funciones. La primera, halla el máximo común divisor de dos números. Para ello, básicamente usa el algoritmo de Euclides.

PLAIN TEXT
RUBY:
  1. def mcd(n1,n2)
  2.   if n1 <n2
  3.     mcd(n2,n1)
  4.   elsif (n1%n2) == 0
  5.     n2
  6.   else
  7.     mcd(n2,n1%n2)
  8.   end
  9. end

Por otra parte, está la función que calcula el mínimo común múltiplo. Para ello usa su relación con el m.c.d..

PLAIN TEXT
RUBY:
  1. def mcm(n1,n2)
  2.   n1*n2/mcd(n1,n2)
  3. end

Con esto, ya podemos pasar a la tercera solución, que es diferente de las otras dos en que sólo sirve si utilizamos dos números como filtro.

PLAIN TEXT
RUBY:
  1. def euler1c(limit, n1, n2)
  2.   (n1*((limit/n1)*((limit/n1)+1)/2)+n2*((limit/n2)*((limit/n2)+1)/2)-mcm(n1,n2)*((limit/mcm(n1,n2))*((limit/mcm(n1,n2))+1)/2))
  3. end
  4.  
  5. puts euler1c(1000,3,5)

La explicación es un poco más complicada, pues implica teoría básica de conjuntos y la conocida fórmula de la suma de 1 a n.

Ahora es vuestro turno. ¿Cómo se os ocurre hacerlo?

article clipper Problema 1
 

share save 171 16 Problema 1

]]> http://www.serabe.com/2008/02/27/problema-1/feed/ 0