numb 0.21.0 → 0.63.0

data/spec/repunit_spec.rb

@@ -0,0 +1,25 @@
+describe Integer, "#repunit?" do
+  REPUNITS = {
+    10 => [1, 11, 111, 1111, 11111, 111111, 1111111, 11111111, 111111111],
+    9 =>  [1, 10, 91, 820, 7381, 66430, 597871, 5380840, 48427561, 435848050],
+    8 =>  [1, 9, 73, 585, 4681, 37449, 299593, 2396745, 19173961, 153391689],
+    7 =>  [1, 8, 57, 400, 2801, 19608, 137257, 960800, 6725601, 47079208],
+    6 =>  [1, 7, 43, 259, 1555, 9331, 55987, 335923, 2015539, 12093235],
+    5 =>  [1, 6, 31, 156, 781, 3906, 19531, 97656, 488281, 2441406],
+    4 =>  [1, 5, 21, 85, 341, 1365, 5461, 21845, 87381, 349525],
+    3 =>  [1, 4, 13, 40, 121, 364, 1093, 3280, 9841, 29524],
+    2 =>  [1, 3, 7, 15, 31, 63, 127, 255, 511, 1023],
+  }
+  REPUNITS.each do |base, members|
+    members.each do |n|
+      it "should return true for #{n} when base = #{base}" do
+        n.repunit?(base).should be_true
+      end
+
+      m = n.odd? ? n - 1 : n + 2
+      it "should return false for #{m} when base = #{base}" do
+        m.repunit?(base).should be_false
+      end
+    end
+  end
+end

data/spec/rough_spec.rb

@@ -0,0 +1,70 @@
+describe Integer, "#rough?" do
+
+  ROUGH = {
+    # A000027
+    2 =>  [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,
+           20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,
+           36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,
+           52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,
+           68,69,70,71,72,73,74,75,76,77],
+    # A005408
+    3 =>  [1,3,5,7,9,11,13,15,17,19,21,23,25,27,29,31,33,35,
+           37,39,41,43,45,47,49,51,53,55,57,59,61,63,65,67,
+           69,71,73,75,77,79,81,83,85,87,89,91,93,95,97,99,
+           101,103,105,107,109,111,113,115,117,119,121,123,
+           125,127,129,131],
+    # A007310
+    5 =>  [1,5,7,11,13,17,19,23,25,29,31,35,37,41,43,47,49,
+           53,55,59,61,65,67,71,73,77,79,83,85,89,91,95,97,
+           101,103,107,109,113,115,119,121,125,127,131,133,
+           137,139,143,145,149],
+    # A007775       
+    7 =>  [1,7,11,13,17,19,23,29,31,37,41,43,47,49,53,59,61,
+           67,71,73,77,79,83,89,91,97,101,103,107,109,113,
+           119,121,127,131,133,137,139,143,149,151,157,161,
+           163,167,169,173,179,181,187,191,193,197,199,203,
+           209],
+    # A008364       
+    11 => [1,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,
+           71,73,79,83,89,97,101,103,107,109,113,121,127,131,
+           137,139,143,149,151,157,163,167,169,173,179,181,
+           187,191,193,197,199,209,211,221,223,227,229,233,
+           239,241,247],
+    # A008365       
+    13 => [1,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,
+           73,79,83,89,97,101,103,107,109,113,127,131,137,
+           139,149,151,157,163,167,169,173,179,181,191,193,
+           197,199,211,221,223,227,229,233,239,241,247,251,
+           257,263,269],
+    # A008366      
+    17 => [1,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,
+           79,83,89,97,101,103,107,109,113,127,131,137,139,
+           149,151,157,163,167,173,179,181,191,193,197,199,
+           211,223,227,229,233,239,241,251,257,263,269,271,
+           277,281,283],
+    # A166061
+    19 => [1,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,
+           83,89,97,101,103,107,109,113,127,131,137,139,149,
+           151,157,163,167,173,179,181,191,193,197,199,211,
+           223,227,229,233,239,241,251,257,263,269,271],
+    # A166063
+    23 => [1,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,
+           89,97,101,103,107,109,113,127,131,137,139,149,151,
+           157,163,167,173,179,181,191,193,197,199,211,223,
+           227,229,233,239,241,251,257,263,269,271]
+  }
+
+  ROUGH.each do |base, members|
+    (members.first..members.last).each do |n|
+      if members.include?(n)
+        it "should return true for rough number #{n} when k=#{base}" do
+          n.rough?(base).should be_true
+        end
+      else
+        it "should return false for non-rough number #{n} when k=#{base}" do
+          n.rough?(base).should be_false
+        end
+      end
+    end
+  end
+end

data/spec/smooth_spec.rb

@@ -0,0 +1,72 @@
+describe Integer, "#smooth?" do
+  SMOOTH = {
+   # A000079
+   2 =>  [1,2,4,8,16,32,64,128,256,512,1024,2048,4096,8192,
+          16384,32768,65536,131072,262144,524288,1048576,
+          2097152,4194304,8388608,16777216,33554432,
+          67108864,134217728,268435456,536870912,1073741824,
+          2147483648,4294967296,8589934592],
+   # A003586 
+   3 =>  [1,2,3,4,6,8,9,12,16,18,24,27,32,36,48,54,64,72,
+          81,96,108,128,144,162,192,216,243,256,288,324,384,
+          432,486,512,576,648,729,768,864,972,1024,1152,
+          1296,1458,1536,1728,1944,2048,2187,2304,2592,2916,
+          3072,3456,3888],
+   # A051037       
+   5 =>  [1,2,3,4,5,6,8,9,10,12,15,16,18,20,24,25,27,30,32,
+          36,40,45,48,50,54,60,64,72,75,80,81,90,96,100,108,
+          120,125,128,135,144,150,160,162,180,192,200,216,
+          225,240,243,250,256,270,288,300,320,324,360,375,
+          384,400,405],
+   # A002473
+   7 =>  [1,2,3,4,5,6,7,8,9,10,12,14,15,16,18,20,21,24,25,
+          27,28,30,32,35,36,40,42,45,48,49,50,54,56,60,63,
+          64,70,72,75,80,81,84,90,96,98,100,105,108,112,120,
+          125,126,128,135,140,144,147,150,160,162,168,175,
+          180,189,192],
+   # A051038
+   11 => [1,2,3,4,5,6,7,8,9,10,11,12,14,15,16,18,20,21,22,
+          24,25,27,28,30,32,33,35,36,40,42,44,45,48,49,50,
+          54,55,56,60,63,64,66,70,72,75,77,80,81,84,88,90,
+          96,98,99,100,105,108,110,112,120,121,125,126,128,
+          132,135,140],
+   # A080197
+   13 => [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,18,20,21,
+          22,24,25,26,27,28,30,32,33,35,36,39,40,42,44,45,
+          48,49,50,52,54,55,56,60,63,64,65,66,70,72,75,77,
+          78,80,81,84,88,90,91,96,98,99,100,104,105,108,110,
+          112,117,120],
+   # A080681       
+   17 => [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,20,
+          21,22,24,25,26,27,28,30,32,33,34,35,36,39,40,42,
+          44,45,48,49,50,51,52,54,55,56,60,63,64,65,66,68,
+          70,72,75,77,78,80,81,84,85,88,90,91,96,98,99,100,
+          102,104,105],
+   # A080682
+   19 => [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,
+          20,21,22,24,25,26,27,28,30,32,33,34,35,36,38,39,
+          40,42,44,45,48,49,50,51,52,54,55,56,57,60,63,64,
+          65,66,68,70,72,75,76,77,78,80,81,84,85,88,90,91,
+          95,96,98,99,100],
+   # A080683
+   23 => [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,
+          20,21,22,23,24,25,26,27,28,30,32,33,34,35,36,38,
+          39,40,42,44,45,46,48,49,50,51,52,54,55,56,57,60,
+          63,64,65,66,68,69,70,72,75,76,77,78,80,81,84,85,
+          88,90,91,92,95]
+  }
+
+  SMOOTH.each do |base, members|
+    20.times.map{ Random.new.rand(members.first..members.last) }.each do |n|
+      if members.include?(n)
+        it "should return true for #{base}-smooth number #{n}" do
+          n.smooth?(base).should be_true
+        end
+      else
+        it "should return false for non-#{base}-rough number #{n}" do
+          n.smooth?(base).should be_false
+        end
+      end
+    end
+  end
+end

data/spec/spec_helper.rb

@@ -1,9 +1,59 @@
+# coding: utf-8
 require 'spec'
 
 $LOAD_PATH.unshift(File.dirname(__FILE__))
 $LOAD_PATH.unshift(File.join(File.dirname(__FILE__), '..', 'lib'))
 require 'numb'
 
+
+class Seq
+  include Enumerable
+  attr_accessor :exclude, :include
+  def initialize(*args)
+    @include = args.first?
+    @exclude = []
+  end
+
+  def each
+    @include.each do |e|
+      next if @exclude.include?(e)
+      yield e
+    end
+  end
+
+  def invert
+    Seq.new(self.begin..self.end).tap{|s| s.exclude = @include }
+  end
+
+  def begin
+    @include.respond_to?(:begin) ? @include.begin : @include.min
+  end
+
+  def end
+    @include.respond_to?(:end) ? @include.end : @include.max
+  end
+
+  def sample(n=nil)
+    return @include.sample(n) if @include.respond_to?(:sample)
+    prng = Random.new
+    (n||1).times.map do 
+      r = prng.rand(@include) 
+      redo if @exclude.include?(r)
+      r
+    end.tap{|a| return a.first unless n}
+  end
+end
+
+class Array
+  def first?
+    size <= 1 ? first : self
+  end
+
+  def to_seq
+    Seq.new *self
+  end
+end
+
 Spec::Runner.configure do |config|
-  
+    
 end

data/spec/square_free_spec.rb

@@ -1,3 +1,5 @@
+# coding: utf-8
+
 describe Integer, "#square_free?" do
   # http://www.research.att.com/~njas/sequences/A005117
   SQUARE_FREE = [1,2,3,5,6,7,10,11,13,14,15,17,19,21,22,23,26,29,

data/spec/sublime_spec.rb

@@ -0,0 +1,18 @@
+describe Integer, "#sublime?" do
+  #  A081357
+  SUBLIME = [12,
+     6086555670238378989670371734243169622657830773351885970528324860512791691264]
+
+  SUBLIME.first(1).each do |n|
+    it "should return true for sublime number #{n}" do
+      n.should be_sublime
+    end
+  end
+
+  (1..100).each do |n|
+    next if SUBLIME.include?(n)
+    it "should return false for non-sublime number #{n}" do
+      n.should_not be_sublime
+    end
+  end
+end

data/spec/sum_of_squares.rb

@@ -0,0 +1,24 @@
+describe Integer, "#sum_of_squares" do
+
+  SOS = [1, # [0, 0, 0, 1]
+         2, # [0, 0, 1, 1]
+         3, # [0, 1, 1, 1]
+         7, # [1, 1, 1, 2]
+         9, # Square number
+         11,# Prime number
+         58, 9634, # Special cases
+         123456789, # Large number
+  ] + 1000.times.map { Random.new.rand(10..1000_000)} # Random numbers
+  SOS.each do |n|
+    squares = n.sum_of_squares
+    it "should return an Array of 4 Integers for #{n}" do
+      squares.should be_kind_of(Array)
+      squares.size.should == 4
+      squares.all?{|e| e.should be_kind_of(Integer)}
+    end
+
+    it "should return integers whose squares sum to give the receiver: #{n}" do
+      squares.map{|i| i**2}.reduce(:+).should == n
+    end
+  end
+end

data/spec/superabundant_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#superabundant?" do
+  # A004394
+  SUPERABUNDANT = [1,2,4,6,12,24,36,48,60,120,180,240,360,720,840,
+                  1260,1680,2520,5040,10080,15120,25200,27720,55440,
+                  110880,166320,277200,332640,554400,665280,720720,
+                  1441440,2162160,3603600,4324320,7207200,8648640,
+                  10810800]
+
+  SUPERABUNDANT.first(20).each do |n|
+    it "should return true for superabundant number #{n}" do
+      n.should be_superabundant
+    end
+  end
+
+  ((1..SUPERABUNDANT.last).first(20) - SUPERABUNDANT).each do |n|
+    it "should return false for non-superabundant number #{n}" do
+      n.should_not be_superabundant
+    end
+  end
+end

data/spec/superperfect_spec.rb

@@ -0,0 +1,17 @@
+describe Integer, "#superperfect?" do
+  # A019279
+  SUPERPERFECT = [2,4,16,64,4096,65536,262144,1073741824,
+                  1152921504606846976]
+
+  SUPERPERFECT[0..-2].each do |n|
+    it "returns true for superperfect number #{n}" do
+      n.should be_superperfect
+    end
+  end
+
+  [1, 3, 5, 78789,  4097].each do |n|
+    it "returns false for non-superperfect number #{n}" do
+      n.should_not be_superperfect
+    end
+  end
+end

data/spec/totient_spec.rb

@@ -0,0 +1,13 @@
+describe Integer, "#totient" do
+  # A000010
+  TOTIENT = [1,1,2,2,4,2,6,4,6,4,10,4,12,6,8,8,16,6,18,8,12,
+             10,22,8,20,12,18,12,28,8,30,16,20,16,24,12,36,18,
+             24,16,40,12,42,20,24,22,46,16,42,20,32,24,52,18,
+             40,24,36,28,58,16,60,30,36,32,48,20,66,32,44]
+
+  TOTIENT.each_with_index do |totient, i|
+    it "should return #{totient} for #{i.succ}" do
+      i.succ.totient.should == totient
+    end
+  end
+end

data/spec/unitary_divisor_spec.rb

@@ -0,0 +1,32 @@
+describe Integer, "#unitary_divisor?" do
+  # A034460 
+  UNITARY_DIVISORS = [0,1,1,1,1,6,1,1,1,8,1,8,1,10,9,1,1,12,1,10,11,14,
+                      1,12,1,16,1,12,1,42,1,1,15,20,13,14,1,22,17,14,1,
+                      54,1,16,15,26,1,20,1,28,21,18,1,30,17,16,23,32,1,
+                      60,1,34,17,1,19,78,1,22,27,74,1,18,1,40,29,24,19,
+                      90,1,22,1,44]
+
+  UNITARY_DIVISORS.each_with_index do |sum, i|
+    i = i.succ
+    it "should return true for the unitary divisors of #{i}" do
+      got = (i.proper_divisors.
+              select{|d| i.unitary_divisor?(d)}.
+              reduce(:+)) || 0
+      sum.should == got
+    end
+  end
+
+  it "should return true for 60 and 5" do
+    60.unitary_divisor?(5).should be_true
+  end
+
+  (10..20).each do |n|
+    it "returns true when the argument is 1" do
+      n.unitary_divisor?(1).should be_true
+    end
+  end
+
+  it "should return false for 60 and 6" do
+    60.unitary_divisor?(6).should be_false
+  end
+end

data/spec/unitary_perfect.rb

@@ -0,0 +1,14 @@
+describe Integer, "#unitary_perfect?" do
+  UNITARY_PERFECT = [6,60,90,87360,146361946186458562560000]
+
+  UNITARY_PERFECT.first(4).each do |n|
+    it "should return true for unitary perfect number #{n}" do
+      n.should be_unitary_perfect
+    end
+
+    m = n + 7
+    it "should return false for non-unitary-perfect number #{m}" do
+      m.should_not be_unitary_perfect
+    end
+  end
+end

data/spec/untouchable_spec.rb

@@ -0,0 +1,19 @@
+describe Integer, "#untouchable?" do
+  UNTOUCHABLE = [2,5,52,88,96,120,124,146,162,188,206,210,216,238,
+                 246,248,262,268,276,288,290,292,304,306,322,324,
+                 326,336,342,372,406,408,426,430,448,472,474,498,
+                 516,518,520,530,540,552,556,562,576,584,612,624,
+                 626,628,658]
+
+  UNTOUCHABLE.first(5).each do |n|
+    it "should return true for untouchable number #{n}" do
+      n.should be_untouchable
+    end
+  end
+    
+  ((1..UNTOUCHABLE.last).to_a - UNTOUCHABLE).first(5).each do |n|
+    it "should return false for non-untouchable number #{n}" do
+      n.should_not be_untouchable
+    end
+  end
+end

data/spec/zeisel_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#zeisel?" do
+  # A051015
+  ZEISEL = [105,1419,1729,1885,4505,5719,15387,24211,25085,
+            27559,31929,54205,59081,114985,207177,208681,
+            233569,287979,294409,336611,353977,448585,507579,
+            721907,982513,1012121,1073305,1242709,1485609,
+            2089257,2263811,2953711,3077705,3506371,3655861,
+            3973085,4648261,5069629,6173179,6253085,6985249,
+            7355239,7355671,7558219,8011459,8413179,8444431,
+            8712985,9271805,9773731,15411785,18175361,
+            18578113,19827641,20771801,23691481,26000605,
+            26758057]
+
+  # 721907 appears to be in this sequence by error...
+  (ZEISEL - [721907]).each do |n|
+    it "should return true for Zeisel number #{n}" do
+      n.should be_zeisel
+    end
+  end
+end