Tag Archives: ruby

Project Euler: Problem 9 in Ruby

I spent way too long trying “clever” solutions that didn’t work or took twice as long as my eventual brute force solution. Screw it.

Problem #9

A Pythagorean triplet is a set of three natural numbers, a < b < c, for which, a^(2) + b^(2) = c^(2) For example, 3^(2) + 4^(2) = 9 + 16 = 25 = 5^(2). There exists exactly one Pythagorean triplet for which a + b + c = 1000. Find the product abc.

target = 1000
min    = Math.sqrt(target).to_i
max    = target / 2

(min..max).each do |a|
  ((a + 1)..max).each do |b|
    c   = Math.sqrt(a**2 + b**2)
    sum = a + b + c
    if sum > target
      break
    elsif sum == target && (a**2 + b**2 == c**2)
      puts a * b * c
      break
    end 
  end
end

Olathe’s is pretty:

 n = 1000
a = (1..n / 2).to_a.find { |a|
  (n * (n / 2 - a) % (n - a)).zero?
}
b = n * (n / 2 - a) / (n - a)
puts "Product is #{a * b * (n - a - b)}."

Project Euler: Problem 8 in Ruby

Problem 8 is to Problem 7 as Problem 6 is to Problem. My original plan involved storing the current sum of 5 numbers and then simply dividing by $big_n[i – 1] and multiplying by $big_n[i + 5], but the zeros made it a bit messy so I thought I’d just give the straight-forward an approach a go to see what I was up against. And it ran correctly, and more importantly quickly so I let it be.

Problem #8

Find the greatest product of five consecutive digits in the 1000-digit number.

73167176531330624919225119674426574742355349194934
96983520312774506326239578318016984801869478851843
85861560789112949495459501737958331952853208805511
12540698747158523863050715693290963295227443043557
66896648950445244523161731856403098711121722383113
62229893423380308135336276614282806444486645238749
30358907296290491560440772390713810515859307960866
70172427121883998797908792274921901699720888093776
65727333001053367881220235421809751254540594752243
52584907711670556013604839586446706324415722155397
53697817977846174064955149290862569321978468622482
83972241375657056057490261407972968652414535100474
82166370484403199890008895243450658541227588666881
16427171479924442928230863465674813919123162824586
17866458359124566529476545682848912883142607690042
24219022671055626321111109370544217506941658960408
07198403850962455444362981230987879927244284909188
84580156166097919133875499200524063689912560717606
05886116467109405077541002256983155200055935729725
71636269561882670428252483600823257530420752963450

Note: Global variables ahead, proceed with caution.

$big_n ="731671765313306249192251196744265747423553
491949349698352031277450632623957831801698480186947
885184385861560789112949495459501737958331952853208
805511125406987471585238630507156932909632952274430
435576689664895044524452316173185640309871112172238
311362229893423380308135336276614282806444486645238
749303589072962904915604407723907138105158593079608
667017242712188399879790879227492190169972088809377
665727333001053367881220235421809751254540594752243
525849077116705560136048395864467063244157221553975
369781797784617406495514929086256932197846862248283
972241375657056057490261407972968652414535100474821
663704844031998900088952434506585412275886668811642
717147992444292823086346567481391912316282458617866
458359124566529476545682848912883142607690042242190
226710556263211111093705442175069416589604080719840
385096245544436298123098787992724428490918884580156
166097919133875499200524063689912560717606058861164
671094050775410022569831552000559357297257163626956
1882670428252483600823257530420752963450".split('')

$stack = []
$chars = 5
max    = 9

def find_product start
  total = 1
  $chars.times do |i|
    total *= $big_n[start + i].to_i
  end
  total
end

$chars.upto($big_n.length - $chars) do |i|
  current = find_product i
  max = [current,max].max
end

puts max

Blessed Olathe had a Ruby-er Solution:

class String
  def max_digit_product length
    # Get each string of 'length' or more nonzero digits.
    strings = scan(/((?:W*[1-9A-Z]){#{length},})/i).map { |string|
      string.first.scan(/([1-9A-Z])/i).map { |digit|
        digit.first.to_i 36
      }
    }
 
    # Find the maximum product of 'length' digits.
    max = 0
    strings.each { |string|
      product = string[0, length].inject(1) { |prod, digit|
        prod * digit
      }
      max = [max, product].max
      (0...(string.length - length)).each { |index|
        product = product / string[index] * string[index + length]
        max = [max, product].max
      }
    }
    max
  end
end

Breakout Knock Off in Ruby with Gosu

I Just finished up a Breakout knock off using Gosu and some basic extension classes I’ve been tinkering around with. There’s only “level” and I’m certain the paddle / ball bounce routine is totally incorrect but I’ve had enough fun with it for tonight. The samples were lovingly provided by http://freesoundfiles.tintagel.net/Audio/, but those swanky graphics were all me.

Download the Code!

Project Euler: Problem 7 in Ruby

It took a couple hours of a couple days for me to come up with the solution to Problem #7, heck of a time. Looking over it now, the solution seems laughably obvious to me given the problems leading up to it. I’m going to optimistically take this as a sign that I’m learning…or something.

When I first began the problem I spent a lot of time trying to wheel and sieve my data set while totally ignoring my redundant method of testing primality!. Now I’m actually utilizing the data I’d already computed to help me speed up my algorithm at the piddly expense of storing a couple thousand integers.

NOTE: I realize that it isn’t actually necessary to store all 10,001 numbers, but it’s my program and I’ll do what I want!

Problem #7

By listing the first six prime numbers: 2, 3, 5, 7, 11, and 13, we can see that the 6^(th) prime is 13.

What is the 10001^(st) prime number?

target = 10001
stack  = [2,3]
n      = 1

def is_prime stack, n
  max_check = Math.sqrt(n).floor
  stack.each do |i|
    if i > max_check
      return true
    end
    if n % i == 0
      return false
    end
  end
  true
end

while true do
  
  if is_prime stack, n += 4
    stack << n
  end

  if is_prime stack, n += 2
    stack << n
  end

  if stack.length >= target
    break
  end
end

puts stack[target - 1]

I’d show you Olathe’s, but he cheated again.

Project Euler: Problem 6 in Ruby

I thought this one was way too easy compared the the last one. There’s some mathematical trickery you can do to speed things up, but the straight-forward solution is running in .009s so I’m going to call this one “done”.

Problem #6

The sum of the squares of the first ten natural numbers is,
1^(2) + 2^(2) + … + 10^(2) = 385

The square of the sum of the first ten natural numbers is,
(1 + 2 + … + 10)^(2) = 55^(2) = 3025

Hence the difference between the sum of the squares of the first ten natural numbers and the square of the sum is 3025 − 385 = 2640.

Find the difference between the sum of the squares of the first one hundred natural numbers and the square of the sum.

max = 100

def sum_of_squares max
  total = 0
  max.times do |i|
    i += 1
    total += i ** 2
  end
  total
end

def square_of_sums max
  total = 0
  max.times do |i|
    i += 1
    total += i
  end
  total ** 2
end

puts square_of_sums(max) - sum_of_squares(max)

I think kode had a neat solution in the forums:

sum, variable = 0, 0
100.downto(1) { |x| sum+=x; variable += (x*x) }
sum_square = sum*sum
puts sum_square - variable

Project Euler: Problem 5 in Ruby

I was stuck on this one for way too long. I was able to figure out a brute force approach quickly enough, but it was taking forever. Finally after hours (told you I took way too long!) of scratching around in my notebook before I had my eureka moment. I basically just determined that the least common multiple between 2 contiguous numbers must be a factor of my final result, and from there I just stepped my way up the chain. It’s running in about 24ms on my machine!

I also determined that I didn’t need to check any numbers under 1/2 of the given number 20 since they were already factors of the higher numbers, but that doesn’t actually effect the Big O and the timing difference is negligible.

PS: Yes, I used the exact library that I accused Olathe and Mat of “cheating” for on Problem #4. Shush!

Problem #5

2520 is the smallest number that can be divided by each of the numbers from 1 to 10 without any remainder.

What is the smallest number that is evenly divisible by all of the numbers from 1 to 20?

 require 'mathn'

max = 20
min = (max / 2).floor + 1
n   = min

(min..max).each do |i|
  n = n.lcm i
end

puts n 

I was pretty proud of this solution and given all the incarnations this program’s been through I thought it came out really nice…until (you guessed it) I saw Olathes!

 require 'rational'
num = (1..20).inject(1) { |result, n| result.lcm n }
puts "Smallest evenly divisible number is #{ num }."

Project Euler: Problem 4 in Ruby

This one was my favorite, probably because it’s been the least mathy so far. I was also tickled pink by the idea of numbers being palindromes. I saw some other solutions on-line but I was really happy with how mine turned out.

Problem #4

A palindromic number reads the same both ways. The largest palindrome made from the product of two 2-digit numbers is 9009 = 91 × 99.

Find the largest palindrome made from the product of two 3-digit numbers.

My Solution:

def get_highest_palindrome high, low
  highest = 0
  high.downto low do |i|
    high.downto low do |j|
      sum = i * j
      if sum <= highest
        break
      end
      if is_palindrome(sum.to_s)
        highest = [highest, sum].max
      end
    end
  end
  highest
end

def is_palindrome str
  return str == str.reverse
end

puts get_highest_palindrome 999, 100

In this particular (and probably only) case I liked my solution better than Olathe's. I found that by counting down from 999 and breaking after any number less than my current max, I was able to save a lot of iterations.

max = 0
100.upto(999) { |a|
  a.upto(999) { |b|
    prod = a * b
    max = [max, prod].max if prod.to_s == prod.to_s.reverse
  }
}
puts "Maximum palindrome is #{ max }."

However, just before I smugged off to bed I saw a post by Begoner in Python. His explanation:

The palindrome can be written as:

abccba

Which then simpifies to:

100000a + 10000b + 1000c + 100c + 10b + a

And then:

100001a + 10010b + 1100c

Factoring out 11, you get:

11(9091a + 910b + 100c)

So the palindrome must be divisible by 11. Seeing as 11 is prime, at least one of the numbers must be divisible by 11. So brute force in Python, only with less numbers to be checked:

def c():
	max = maxI = maxJ = 0
	i = 999
	j = 990
	while (i > 100):
		j = 990
		while (j > 100):
			product = i * j
			if (product > max):
				productString = str(product)
				if (productString == productString[::-1]):
					max = product
					maxI = i
					maxJ = j
			j -= 11
		i -= 1
	return max, maxI, maxJ

I got PWNED!

Project Euler: Problem 3 in Ruby

It’s been a number of years since my last math class, so I had to bust out the pen and paper. I figured that for given number n you can stop checking for factors after n / 2, but I didn’t realize till much later (and by “realize”, I mean read up on prime numbers) that you can actually stop checking for prime factors at the square root. This saves considerable time when you’re dealing with numbers in the quadrabazillions. It’s the difference between < 2 seconds and I-did-the-dishes-and-it's-still-running. On a side note, I don't consider reading up on the subject "cheating". I consider looking up algorithms and other solutions "cheating". Problem #3

The prime factors of 13195 are 5, 7, 13 and 29.

What is the largest prime factor of the number 600851475143 ?

target = 600_851_475_143

def find_highest_prime_factor n
  (Math.sqrt(n).ceil).downto 2 do |i|
    if n % i == 0 && find_highest_prime_factor(i) == 1
      return i
    end
  end
  1
end

puts find_highest_prime_factor(target)

As expected, Olathe delivers. Runs roughly 3 times as fast and is easier on the eyes…even though he loaded an external library (cheater!)

 require 'mathn'
num, factor = 317_584_931_803, 0
primes = Prime.new
while num > 1
  factor = primes.next
  num /= factor while (num % factor).zero?
end
 
puts "Largest factor is #{ factor }."

However, if you’re going to cheat…cheat big! Posted by mat:

 require 'mathn'
317584931803.prime_division.last[0]

**UPDATE**
The two not-mine solutions use a different number (317584931803. instead of 600851475143) so the solution is different, but the algorithms are correct.

Project Euler: Problem 2 in Ruby

My solution for Project Euler problem #2 came out a little fugly mainly because of the global variables, but it works. The “secret” to getting the program to run in a reasonable amount of time is to store the values of the fibonacci sequence as you compute them instead of processing the same data over again.

Problem #2

Each new term in the Fibonacci sequence is generated by adding the previous two terms. By starting with 1 and 2, the first 10 terms will be:

1, 2, 3, 5, 8, 13, 21, 34, 55, 89, …

Find the sum of all the even-valued terms in the sequence which do not exceed four million.

My Solution

max    = 4_000_000
total  = 0
$stack = []

def fib n
  if n == 0
    return 1
  end
  if n == 1
    return 2
  end
  return $stack[n - 1] + $stack[n - 2]
end

max.times do |i|
  $stack[i] = fib i
  if $stack[i] > max
    break
  end
  if $stack[i] & 1 == 0
    total += $stack[i]
  end
end

puts total

As is quickly becoming the case, Olathe’s solution was much more elegant than mine:

x, y, sum = 1, 1, 0
while sum < 1_000_000
  sum += (x + y)
  x, y = x + 2*y, 2*x + 3*y
end
 
puts "Sum is #{ sum }."

Project Euler: Problem 1 in Ruby

I’ve been working through some of the Project Euler problems in Ruby and I thought I’d post my solutions. I refuse to “cheat”, but it’s really neat to see all the different ways other people solve the problems afterwards.

Problem #1

If we list all the natural numbers below 10 that are multiples of 3 or 5, we get 3, 5, 6 and 9. The sum of these multiples is 23.

Find the sum of all the multiples of 3 or 5 below 1000.

total = 0

1000.times do |i|
  if i % 3 == 0 || i % 5 == 0
    total += i
  end
end

puts total

I googled around a bit after I finished and particularly liked this solution, which seemed more ruby-ish to me:

answer = (0..999).select { |a| a%3 ==0 || a%5==0 }
puts answer.inject { |sum, n| sum+n }

And, from the forums I thought Olathe’s solution was really neat, although mathily over my head:

class Integer
  def sum_mod modulus
    n = self.div modulus
    modulus * n * (n + 1) / 2
  end
end
 
num = 999
sum = num.sum_mod(3) + num.sum_mod(5) - num.sum_mod(15)
 
puts "Sum is #{ sum }."