numb 0.99.0 → 0.109.0

data/Rakefile

@@ -35,4 +35,7 @@ rescue LoadError
   end
 end
 
+task :irb do 
+  exec("irb -r./lib/numb")
+end
 task :default => :spec

data/VERSION

@@ -1 +1 @@
-0.99.0
+0.109.0

data/lib/numb/amenable.rb

@@ -0,0 +1,5 @@
+class Integer
+  def amenable?
+    self != 4 and modulo(4) < 2
+  end
+end

data/lib/numb/betrothed.rb

@@ -0,0 +1,8 @@
+# coding: utf-8
+class Integer
+  def betrothed?(m)
+    σ == m.σ and σ == self.succ + m
+  end
+  alias :quasi_amicable? :betrothed?
+  alias :reduced_amicable? :betrothed?
+end

data/lib/numb/blum.rb

@@ -0,0 +1,6 @@
+class Integer
+  def blum?
+    return false unless prime_factors.size == 2
+    prime_factors.uniq.size == 2 and prime_factors.all?{|p| p.modulo(4) == 3}
+  end
+end

data/lib/numb/euclid.rb

@@ -0,0 +1,5 @@
+class Integer
+  def euclid?
+    (self - 1).primorial?
+  end
+end

data/lib/numb/giuga.rb

@@ -0,0 +1,7 @@
+class Integer
+  def giuga?
+    composite? and prime_factors.uniq.all? do |p|
+      ((self / p) - 1).divides?(p)
+    end
+  end
+end

data/lib/numb/iban.rb

@@ -0,0 +1,6 @@
+class Integer
+  def iban?
+    not words.include?('i')
+  end
+  alias :blind? :iban?
+end

data/lib/numb/lychrel.rb

@@ -0,0 +1,15 @@
+class Integer
+  def palindrome?
+    to_s.reverse.to_i == self
+  end
+
+  def lychrel?
+    n = self
+    # This limit is as arbitrary as it looks
+    100.times do
+      return false if n.palindrome?
+      n += n.to_s.reverse.to_i
+    end
+    true
+  end
+end

data/lib/numb/oban.rb

@@ -0,0 +1,5 @@
+class Integer
+  def oban?
+    not words.include?(?o)
+  end
+end

data/lib/numb/primorial.rb

@@ -17,10 +17,6 @@ class Integer
   def primorial
     return nil if self < 1
     return 1 if self == 1
-    if prime?
-      self * (self - 1).primorial
-    else
-      (self - 1).primorial
-    end
+    (prime? ? self : prev_prime).downto(2).select(&:prime?).reduce(:*)
   end
 end

data/lib/numb/strictly_non_palindromic.rb

@@ -0,0 +1,10 @@
+class String
+  def palindrome?; self == reverse; end
+end
+
+class Integer
+  def strictly_non_palindromic?
+    return true if (0..4).include?(self) or self == 6
+    prime? and (2..(Math.sqrt(self).floor)).none?{|base| base(base).palindrome?}
+  end
+end

data/lib/numb/uban.rb

@@ -0,0 +1,5 @@
+class Integer
+  def uban?
+    not words.include?(?u)
+  end
+end

data/spec/numb/amenable_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#amenable?" do
+  # A100832
+  @seq = [1,5,8,9,12,13,16,17,20,21,24,25,28,29,32,33,36,
+          37,40,41,44,45,48,49,52,53,56,57,60,61,64,65,68,
+          69,72,73,76,77,80,81,84,85,88,89,92,93,96,97,100,
+          101,104,105,108,109,112,113,116,117,120,121,124,
+          125,128,129,132]
+
+  @seq.each do |n|
+    it "returns true for amenable number #{n}" do
+      n.should be_amenable
+    end
+  end
+
+  @seq.to_seq.invert.each do |n|
+    it "returns false for non-amenable number #{n}" do
+      n.should_not be_amenable
+    end
+  end
+end

data/spec/numb/betrothed_spec.rb

@@ -0,0 +1,58 @@
+describe Integer, "#betrothed?" do
+  # A005276, augmented with http://djm.cc/reduced.fmtlist
+  @seq = [[48, 75],
+          [140, 195],
+          [1050, 1925],
+          [1575, 1648],
+          [2024, 2295],
+          [5775, 6128],
+          [8892, 16587],
+          [9504, 20735],
+          [62744, 75495],
+          [186615, 206504],
+          [196664, 219975],
+          [199760, 309135],
+          [266000, 507759],
+          [312620, 549219],
+          [526575, 544784],
+          [573560, 817479],  
+          [587460, 1057595],
+          [1000824, 1902215],
+          [1081184, 1331967],
+          [1139144, 1159095],
+          [1140020, 1763019],
+          [1173704, 1341495],
+          [1208504, 1348935],
+          [1233056, 1524831],
+          [1236536, 1459143],
+          [1279950, 2576945],
+          [1921185, 2226014],
+          [2036420, 2681019],
+          [2102750, 2142945],
+          [2140215, 2421704],
+  ]
+
+  @seq.each do |n, m|
+    it "returns true for betrothed pair #{n}, #{m}" do
+      n.should be_betrothed(m)
+      m.should be_betrothed(n)
+    end
+  end
+
+  @seq.each do |n, m|
+    n -= 1
+    it "returns false for non-betrothed pair #{n}, #{m}" do
+      n.should_not be_betrothed(m)
+      m.should_not be_betrothed(n)
+    end
+  end
+
+  10.times do
+    n, m = 2.times.map{ Random.new.rand(@seq.last.max) }
+    redo if @seq.include?([n, m])
+    it "returns false for non-betrothed pair #{n}, #{m}" do
+      n.should_not be_betrothed(m)
+      m.should_not be_betrothed(n)
+    end
+  end
+end

data/spec/numb/blum_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#blum?" do
+  # A016105
+  @seq = [21,33,57,69,77,93,129,133,141,161,177,201,209,
+          213,217,237,249,253,301,309,321,329,341,381,393,
+          413,417,437,453,469,473,489,497,501,517,537,553,
+          573,581,589,597,633,649,669,681,713,717,721,737,
+          749,753,781,789]
+
+  @seq.each do |n|
+    it "returns true for Blum integer #{n}" do
+      n.should be_blum
+    end
+  end
+
+  @seq.to_seq.invert.each do |n|
+    it "returns false for non-Blum integer #{n}" do
+      n.should_not be_blum
+    end
+  end
+end

data/spec/numb/euclid_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#euclid?" do
+  # A006862 
+  @seq = [2,3,7,31,211,2311,30031,510511,9699691,223092871,
+          6469693231,200560490131,7420738134811,
+          304250263527211,13082761331670031,
+          614889782588491411,32589158477190044731,
+          1922760350154212639071]
+
+  @seq.each do |n|
+    it "returns true for Euclid number #{n}" do
+      n.should be_euclid
+    end
+  end
+
+  @seq.to_seq.invert.first(100).each do |n|
+    it "returns false for Euclid number #{n}" do
+      n.should_not be_euclid
+    end
+  end
+end

data/spec/numb/giuga_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#giuga?" do
+  # A007850
+  @seq = [30,858,1722,66198,2214408306,24423128562,
+          432749205173838,14737133470010574,
+          550843391309130318,244197000982499715087866346,
+          554079914617070801288578559178,
+          1910667181420507984555759916338506].first(6)
+
+  @seq.each do |n|
+    it "returns true for Giuga number #{n}" do
+      n.should be_giuga
+    end
+  end
+
+  @seq.to_seq.invert.sample(10).each do |n|
+    it "returns false for non-Giuga number #{n}" do
+      n.should_not be_giuga
+    end
+  end
+end

data/spec/numb/iban_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#iban?" do
+  # A089589
+  @seq = [0,1,2,3,4,7,10,11,12,14,17,20,21,22,23,24,27,40,
+          41,42,43,44,47,70,71,72,73,74,77,100,101,102,103,
+          104,107,110,111,112,114,117,120,121,122,123,124,
+          127,140,141,142,143,144,147,170,171,172,173,174,
+          177,200,201]
+
+  @seq.each do |n|
+    it "returns true for iban number #{n}" do
+      n.should be_iban
+    end
+  end
+
+  @seq.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-iban number #{n}" do
+      n.should_not be_iban
+    end
+  end
+end

data/spec/numb/lychrel_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#lychrel?" do
+  # A023108 
+  @seq = [196,295,394,493,592,689,691,788,790,879,887,978,
+          986,1495,1497,1585,1587,1675,1677,1765,1767,1855,
+          1857,1945,1947,1997,2494,2496,2584,2586,2674,2676,
+          2764,2766,2854,2856,2944,2946,2996,3493,3495,3583,
+          3585,3673,3675]
+
+  @seq.each do |n|
+    it "returns true for Lychrel number #{n}" do
+      n.should be_lychrel
+    end
+  end
+
+  @seq.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-Lychrel number #{n}" do
+      n.should_not be_lychrel
+    end
+  end
+end

data/spec/numb/oban_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#oban?" do
+  # A008521
+  @seq = [3,5,6,7,8,9,10,11,12,13,15,16,17,18,19,20,23,25,
+          26,27,28,29,30,33,35,36,37,38,39,50,53,55,56,57,
+          58,59,60,63,65,66,67,68,69,70,73,75,76,77,78,79,
+          80,83,85,86,87,88,89,90,93,95,96,97,98,99,300,303,
+          305,306,307]
+
+  @seq.each do |n|
+    it "returns true for oban number #{n}" do
+      n.should be_oban
+    end
+  end
+
+  @seq.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-oban number #{n}" do
+      n.should_not be_oban
+    end
+  end
+end

data/spec/numb/primorial_spec.rb

@@ -22,7 +22,7 @@ describe Integer, "#primorial?" do
     end
   end
 
-  ((1..PRIMORIAL.last).first(20) - PRIMORIAL).uniq.each do |n|
+  PRIMORIAL.to_seq.invert.first(20).each do |n|
     it "should return false for non-primorial number #{n}" do
       n.should_not be_primorial
     end

data/spec/numb/strictly_non_palindromic_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#strictly_non_palindromic?" do
+  # A016038
+  @seq = [0,1,2,3,4,6,11,19,47,53,79,103,137,139,149,163,
+          167,179,223,263,269,283,293,311,317,347,359,367,
+          389,439,491,563,569,593,607,659,739,827,853,877,
+          977,983,997,1019,1049,1061,1187,1213,1237,1367,
+          1433,1439,1447,1459]
+
+  @seq.each do |n|
+    it "returns true for strictly-non-palindromic number #{n}" do
+      n.should be_strictly_non_palindromic
+    end
+  end
+
+  @seq.to_seq.invert.each do |n|
+    it "returns false for non-strictly-non-palindromic number #{n}" do
+      n.should_not be_strictly_non_palindromic
+    end
+  end
+end

data/spec/numb/uban_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#uban?" do
+  # A089590 
+  @seq = [1,2,3,5,6,7,8,9,10,11,12,13,15,16,17,18,19,20,21,
+          22,23,25,26,27,28,29,30,31,32,33,35,36,37,38,39,
+          40,41,42,43,45,46,47,48,49,50,51,52,53,55,56,57,
+          58,59,60,61,62,63,65,66,67,68,69,70,71,72,73,75,
+          76,77,78,79]
+
+  @seq.each do |n|
+    it "returns true for uban number #{n}" do
+      n.should be_uban
+    end
+  end
+
+  @seq.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-uban number #{n}" do
+      n.should_not be_uban
+    end
+  end
+end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification 
 name: numb
 version: !ruby/object:Gem::Version 
-  version: 0.99.0
+  version: 0.109.0
 platform: ruby
 authors: 
 - Run Paint Run Run
@@ -46,6 +46,7 @@ files:
 - lib/numb/aliquot_sequence.rb
 - lib/numb/almost_perfect.rb
 - lib/numb/almost_prime.rb
+- lib/numb/amenable.rb
 - lib/numb/amicable.rb
 - lib/numb/apocalyptic.rb
 - lib/numb/aspiring.rb
@@ -53,8 +54,10 @@ files:
 - lib/numb/automorphic.rb
 - lib/numb/balanced_prime.rb
 - lib/numb/base.rb
+- lib/numb/betrothed.rb
 - lib/numb/binomial.rb
 - lib/numb/biquadratic.rb
+- lib/numb/blum.rb
 - lib/numb/breeder.rb
 - lib/numb/brown.rb
 - lib/numb/carmichael.rb
@@ -78,6 +81,7 @@ files:
 - lib/numb/economical.rb
 - lib/numb/emrip.rb
 - lib/numb/equidigital.rb
+- lib/numb/euclid.rb
 - lib/numb/evil.rb
 - lib/numb/extravagant.rb
 - lib/numb/factorial.rb
@@ -86,6 +90,7 @@ files:
 - lib/numb/fibonacci.rb
 - lib/numb/friendly.rb
 - lib/numb/frugal.rb
+- lib/numb/giuga.rb
 - lib/numb/goldbach.rb
 - lib/numb/hamming.rb
 - lib/numb/happy.rb
@@ -96,6 +101,7 @@ files:
 - lib/numb/highly_composite.rb
 - lib/numb/hilbert.rb
 - lib/numb/hyperperfect.rb
+- lib/numb/iban.rb
 - lib/numb/idoneal.rb
 - lib/numb/impolite.rb
 - lib/numb/integer_p.rb
@@ -110,6 +116,7 @@ files:
 - lib/numb/leyland.rb
 - lib/numb/lucas.rb
 - lib/numb/lucas_carmichael.rb
+- lib/numb/lychrel.rb
 - lib/numb/mersenne_prime.rb
 - lib/numb/minimal.rb
 - lib/numb/mms_pair.rb
@@ -123,6 +130,7 @@ files:
 - lib/numb/nonhypotenuse.rb
 - lib/numb/nth_prime.rb
 - lib/numb/number_of_divisors.rb
+- lib/numb/oban.rb
 - lib/numb/octagonal.rb
 - lib/numb/odious.rb
 - lib/numb/ordinal.rb
@@ -165,6 +173,7 @@ files:
 - lib/numb/square_free.rb
 - lib/numb/square_triangular.rb
 - lib/numb/squared_triangular.rb
+- lib/numb/strictly_non_palindromic.rb
 - "lib/numb/st\xC3\xB8rmer.rb"
 - lib/numb/sublime.rb
 - lib/numb/sum_of_squares.rb
@@ -176,6 +185,7 @@ files:
 - lib/numb/triangular.rb
 - lib/numb/trimorphic.rb
 - lib/numb/twin_prime.rb
+- lib/numb/uban.rb
 - lib/numb/undulating.rb
 - lib/numb/unhappy.rb
 - lib/numb/unitary_amicable.rb
@@ -197,6 +207,7 @@ files:
 - spec/numb/achilles_spec.rb
 - spec/numb/almost_perfect_spec.rb
 - spec/numb/almost_prime_spec.rb
+- spec/numb/amenable_spec.rb
 - spec/numb/amicable_spec.rb
 - spec/numb/apocalyptic_spec.rb
 - spec/numb/aspiring_spec.rb
@@ -204,8 +215,10 @@ files:
 - spec/numb/automorphic_spec.rb
 - spec/numb/balanced_prime_spec.rb
 - spec/numb/base_spec.rb
+- spec/numb/betrothed_spec.rb
 - spec/numb/binomial_spec.rb
 - spec/numb/biquadratic_spec.rb
+- spec/numb/blum_spec.rb
 - spec/numb/breeder_spec.rb
 - spec/numb/brown_spec.rb
 - spec/numb/carmichael_spec.rb
@@ -230,6 +243,7 @@ files:
 - spec/numb/economical_spec.rb
 - spec/numb/emrip_spec.rb
 - spec/numb/equidigital_spec.rb
+- spec/numb/euclid_spec.rb
 - spec/numb/evil_spec.rb
 - spec/numb/extravagant_spec.rb
 - spec/numb/factorial_spec.rb
@@ -238,6 +252,7 @@ files:
 - spec/numb/fibonacci_spec.rb
 - spec/numb/friendly_spec.rb
 - spec/numb/frugal_spec.rb
+- spec/numb/giuga_spec.rb
 - spec/numb/goldbach_spec.rb
 - spec/numb/hamming_spec.rb
 - spec/numb/happy_spec.rb
@@ -248,6 +263,7 @@ files:
 - spec/numb/highly_composite_spec.rb
 - spec/numb/hilbert_spec.rb
 - spec/numb/hyperperfect_spec.rb
+- spec/numb/iban_spec.rb
 - spec/numb/idoneal_spec.rb
 - spec/numb/integer_p_spec.rb
 - spec/numb/interprime_spec.rb
@@ -261,6 +277,7 @@ files:
 - spec/numb/leyland_spec.rb
 - spec/numb/lucas_carmichael_spec.rb
 - spec/numb/lucas_spec.rb
+- spec/numb/lychrel_spec.rb
 - spec/numb/maris_mcgwire_sosa_pair_spec.rb
 - spec/numb/mersenne_prime_spec.rb
 - spec/numb/minimal_spec.rb
@@ -273,6 +290,7 @@ files:
 - spec/numb/nth_prime_spec.rb
 - spec/numb/number_of_distinct_prime_factors_spec.rb
 - spec/numb/number_of_prime_factors_spec.rb
+- spec/numb/oban_spec.rb
 - spec/numb/octagonal_spec.rb
 - spec/numb/odious_spec.rb
 - spec/numb/ordinal_spec.rb
@@ -316,6 +334,7 @@ files:
 - spec/numb/square_spec.rb
 - spec/numb/square_triangular_spec.rb
 - spec/numb/squared_triangular_spec.rb
+- spec/numb/strictly_non_palindromic_spec.rb
 - "spec/numb/st\xC3\xB8rmer_spec.rb"
 - spec/numb/sublime_spec.rb
 - spec/numb/sum_of_squares_spec.rb
@@ -327,6 +346,7 @@ files:
 - spec/numb/triangular_spec.rb
 - spec/numb/trimorphic_spec.rb
 - spec/numb/twin_prime_spec.rb
+- spec/numb/uban_spec.rb
 - spec/numb/undulating_spec.rb
 - spec/numb/unhappy_spec.rb
 - spec/numb/unitary_amicable_spec.rb
@@ -381,6 +401,8 @@ test_files:
 - spec/numb/unitary_divisor_spec.rb
 - spec/numb/primitive_pseudoperfect_spec.rb
 - spec/numb/polydivisible_spec.rb
+- spec/numb/strictly_non_palindromic_spec.rb
+- spec/numb/iban_spec.rb
 - spec/numb/equidigital_spec.rb
 - spec/numb/primorial_spec.rb
 - spec/numb/superabundant_spec.rb
@@ -406,6 +428,7 @@ test_files:
 - spec/numb/carmichael_spec.rb
 - spec/numb/woodall_spec.rb
 - spec/numb/abundancy_spec.rb
+- spec/numb/giuga_spec.rb
 - spec/numb/polite_spec.rb
 - spec/numb/unitary_amicable_spec.rb
 - spec/numb/squared_triangular_spec.rb
@@ -416,6 +439,7 @@ test_files:
 - spec/numb/jacobsthal_lucas_spec.rb
 - spec/numb/aban_spec.rb
 - spec/numb/hyperperfect_spec.rb
+- spec/numb/uban_spec.rb
 - spec/numb/zeisel_spec.rb
 - spec/numb/smooth_spec.rb
 - spec/numb/d_spec.rb
@@ -436,6 +460,7 @@ test_files:
 - spec/numb/safe_prime_spec.rb
 - spec/numb/sublime_spec.rb
 - spec/numb/refactorable_spec.rb
+- spec/numb/betrothed_spec.rb
 - spec/numb/smith_spec.rb
 - spec/numb/unusual_spec.rb
 - spec/numb/octagonal_spec.rb
@@ -466,10 +491,12 @@ test_files:
 - spec/numb/frugal_spec.rb
 - spec/numb/divides_spec.rb
 - spec/numb/odious_spec.rb
+- spec/numb/euclid_spec.rb
 - spec/numb/semiprime_spec.rb
 - spec/numb/deficient_spec.rb
 - spec/numb/almost_perfect_spec.rb
 - spec/numb/amicable_spec.rb
+- spec/numb/amenable_spec.rb
 - spec/numb/happy_spec.rb
 - spec/numb/twin_prime_spec.rb
 - spec/numb/hexagonal_spec.rb
@@ -479,6 +506,7 @@ test_files:
 - spec/numb/integer_p_spec.rb
 - spec/numb/semi_perfect_spec.rb
 - spec/numb/pentagonal_spec.rb
+- spec/numb/blum_spec.rb
 - spec/numb/smarandache_wellin_spec.rb
 - spec/numb/maris_mcgwire_sosa_pair_spec.rb
 - spec/numb/core_spec.rb
@@ -504,6 +532,7 @@ test_files:
 - spec/numb/factorion_spec.rb
 - spec/numb/undulating_spec.rb
 - spec/numb/extravagant_spec.rb
+- spec/numb/lychrel_spec.rb
 - spec/numb/self_descriptive_spec.rb
 - spec/numb/decagonal_spec.rb
 - spec/numb/apocalyptic_spec.rb
@@ -513,6 +542,7 @@ test_files:
 - spec/numb/keith_spec.rb
 - spec/numb/narcissistic_spec.rb
 - spec/numb/myriagonal_spec.rb
+- spec/numb/oban_spec.rb
 - spec/numb/achilles_spec.rb
 - spec/numb/coprime_spec.rb
 - spec/numb/almost_prime_spec.rb