numb 0.152.0 → 0.170.0

data/spec/numb/jacobsthal_lucas_spec.rb

@@ -1,20 +1,28 @@
-describe Integer, "#jacobsthal_lucas?" do
-  # A014551
-  @seq = [2,1,5,7,17,31,65,127,257,511,1025,2047,4097,8191,
-          16385,32767,65537,131071,262145,524287,1048577,
-          2097151,4194305,8388607,16777217,33554431,
-          67108865,134217727,268435457,536870911,1073741825,
-          2147483647].to_seq
+# A014551
+JACOBSTHAL_LUCAS = [2,1,5,7,17,31,65,127,257,511,1025,2047,4097,8191,
+                    16385,32767,65537,131071,262145,524287,1048577,
+                    2097151,4194305,8388607,16777217,33554431,
+                    67108865,134217727,268435457,536870911,1073741825,
+                    2147483647]
 
-  @seq.each do |n|
-    it "should return true for Jacobsthal-Lucas number #{n}" do
+describe Integer, "#jacobsthal_lucas?" do
+  JACOBSTHAL_LUCAS.each do |n|
+    it "returns true for Jacobsthal-Lucas number #{n}" do
       n.should be_jacobsthal_lucas
     end
   end
 
-  @seq.invert.sample(10).each do |n|
-    it "should return false for non-Jacobsthal-Lucas number #{n}" do
+  JACOBSTHAL_LUCAS.to_seq.invert.sample(10).each do |n|
+    it "returns false for non-Jacobsthal-Lucas number #{n}" do
       n.should_not be_jacobsthal_lucas
     end
   end
 end
+
+describe Integer, "#jacobsthal_lucas" do
+  JACOBSTHAL_LUCAS.each_with_index do |l, n|
+    it "returns #{l} for the #{n}#{n.ordinal} Jacobsthal-Lucas number" do
+      n.jacobsthal_lucas.should == l
+    end
+  end
+end

data/spec/numb/lucas2_spec.rb

@@ -0,0 +1,30 @@
+# coding: utf-8
+describe Integer, "#lucas2" do
+  @seq = {
+    # A000204
+    [1, -1] => [2]+[1,3,4,7,11,18,29,47,76,123,199,322,521,843,1364,
+                2207,3571,5778,9349,15127,24476,39603,64079,
+                103682,167761,271443,439204,710647,1149851,
+                1860498,3010349,4870847,7881196,12752043],
+    # A002203
+    [2, -1] => [2,2,6,14,34,82,198,478,1154,2786,6726,16238,
+                39202,94642,228486,551614,1331714,3215042,7761798,
+                18738638,45239074,109216786,263672646,636562078,
+                1536796802,3710155682,8957108166,21624372014,
+                52205852194],
+    # A014551            
+    [1, -2] => [2,1,5,7,17,31,65,127,257,511,1025,2047,4097,8191,
+                16385,32767,65537,131071,262145,524287,1048577,
+                2097151,4194305,8388607,16777217,33554431,
+                67108865,134217727,268435457,536870911,1073741825,
+                2147483647]
+  }
+
+  @seq.each do |(p, q), seq|
+    seq.each_with_index do |l, n|
+      it "returns #{l} for the #{n}#{n.ordinal} Lucas number of the 2nd kind (#{p}, #{q})" do
+        n.lucas2(p, q).should == l
+      end
+    end
+  end
+end

data/spec/numb/lucas_spec.rb

@@ -1,20 +1,28 @@
-describe Integer, "#lucas?" do
-  # A000032
-  @seq = [2,1,3,4,7,11,18,29,47,76,123,199,322,521,843,
-          1364,2207,3571,5778,9349,15127,24476,39603,64079,
-          103682,167761,271443,439204,710647,1149851,
-          1860498,3010349,4870847,7881196,12752043,20633239,
-          33385282].to_seq
+# A000032
+LUCAS = [2,1,3,4,7,11,18,29,47,76,123,199,322,521,843,
+         1364,2207,3571,5778,9349,15127,24476,39603,64079,
+         103682,167761,271443,439204,710647,1149851,
+         1860498,3010349,4870847,7881196,12752043,20633239,
+         33385282]
 
-  @seq.each do |n|
+describe Integer, "#lucas?" do
+  LUCAS.each do |n|
     it "returns true for Lucas number #{n}" do
       n.should be_lucas
     end
   end
 
-  @seq.invert.sample(10).each do |n|
+  LUCAS.to_seq.invert.sample(10).each do |n|
     it "returns false for non-Lucas number #{n}" do
       n.should_not be_lucas
     end
   end
 end
+
+describe Integer, "#lucas" do
+  LUCAS.each_with_index do |l, n|
+    it "returns #{l} for the #{n}#{n.ordinal} Lucas number" do
+      n.lucas.should == l
+    end
+  end
+end

data/spec/numb/minimal_spec.rb

@@ -5,7 +5,7 @@ describe Integer, "#minimal?" do
           2880,3072,3600,4096,5040,5184,6300,6480,6720,7560,
           9216,10080].to_seq
 
-  @seq.first(20).sample(10).each do |n|
+  @seq.first(25).each do |n|
     it "should return true for minimal number #{n}" do
       n.should be_minimal
     end

data/spec/numb/near_square_spec.rb

@@ -0,0 +1,54 @@
+# coding: utf-8
+describe Integer, "#near_square?" do
+  @seq = {
+    # A028875
+    -5 => [4,11,20,31,44,59,76,95,116,139,164,191,220,251,
+           284,319,356,395,436,479,524,571,620,671,724,779,
+           836,895,956,1019,1084,1151,1220,1291,1364,1439,
+           1516,1595,1676,1759,1844,1931,2020,2111,2204,2299,
+           2396],
+    # A028347       
+    -4 => [0,5,12,21,32,45,60,77,96,117,140,165,192,221,252,
+          285,320,357,396,437,480,525,572,621,672,725,780,
+          837,896,957,1020,1085,1152,1221,1292,1365,1440,
+          1517,1596,1677,1760,1845],
+    # A028872
+    -3 => [1,6,13,22,33,46,61,78,97,118,141,166,193,222,253,
+          286,321,358,397,438,481,526,573,622,673,726,781,
+          838,897,958,1021,1086,1153,1222,1293,1366,1441,
+          1518,1597,1678,1761,1846,1933,2022,2113,2206,
+          2301],
+    # A117950      
+    3 => [3,4,7,12,19,28,39,52,67,84,103,124,147,172,199,
+          228,259,292,327,364,403,444,487,532,579,628,679,
+          732,787,844,903,964,1027,1092,1159,1228,1299,1372,
+          1447,1524,1603,1684,1767,1852,1939,2028,2119,2212,
+          2307,2404,2503],
+    # A087475
+    4 => [4,5,8,13,20,29,40,53,68,85,104,125,148,173,200,
+          229,260,293,328,365,404,445,488,533,580,629,680,
+          733,788,845,904,965,1028,1093,1160,1229,1300,1373,
+          1448,1525,1604,1685,1768,1853,1940,2029,2120,2213,
+          2308,2405,2504],
+    # A117951       
+    5 => [5,6,9,14,21,30,41,54,69,86,105,126,149,174,201,
+          230,261,294,329,366,405,446,489,534,581,630,681,
+          734,789,846,905,966,1029,1094,1161,1230,1301,1374,
+          1449,1526,1605,1686,1769,1854,1941,2030,2121,2214,
+          2309,2406,2505]
+  }
+
+  @seq.each do |k, values|
+    values.each do |n|
+      it "returns true for #{n}^2 #{'+' if k > 0}#{k}" do
+        n.should be_near_square(k)
+      end
+    end
+
+    values.to_seq.invert.sample(100).each do |n|
+      it "returns false for #{n}^2 #{'+' if k > 0}#{k}" do
+        n.should_not be_near_square(k)
+      end
+    end
+  end
+end

data/spec/numb/octahedral_spec.rb

@@ -0,0 +1,29 @@
+# coding: utf-8
+# A005900
+OCTAHEDRAL = [0,1,6,19,44,85,146,231,344,489,670,891,1156,1469,
+              1834,2255,2736,3281,3894,4579,5340,6181,7106,8119,
+              9224,10425,11726,13131,14644,16269,18010,19871,
+              21856,23969,26214,28595,31116,33781,36594,39559,
+              42680]
+
+describe Integer, "#octahedral" do
+  OCTAHEDRAL.each_with_index do |o,n|
+    it "returns #{o} as the #{n}#{n.ordinal} octahedral number" do
+      n.octahedral.should == o
+    end
+  end
+end
+
+describe Integer, "#octahedral?" do
+  OCTAHEDRAL.each do |n|
+    it "returns true for octahedral number #{n}" do
+      n.should be_octahedral
+    end
+  end
+
+  OCTAHEDRAL.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-octahedral number #{n}" do
+      n.should_not be_octahedral
+    end
+  end
+end

data/spec/numb/palindrome_spec.rb

@@ -0,0 +1,43 @@
+# coding: utf-8
+describe Integer, "#palindrome?" do
+  @seq = {
+    # A057148
+    2 => [0,1,0b11,0b101,0b111,0b1001,0b1111,0b10001,0b10101,0b11011,0b11111,
+          0b100001,0b101101,0b110011,0b111111,0b1000001,0b1001001,
+          0b1010101,0b1011101,0b1100011,0b1101011,0b1110111,0b1111111,
+          0b10000001,0b10011001,0b10100101,0b10111101,0b11000011,
+          0b11011011],
+    # A014190      
+    3 => [0,1,2,4,8,10,13,16,20,23,26,28,40,52,56,68,80,82,
+          91,100,112,121,130,142,151,160,164,173,182,194,
+          203,212,224,233,242,244,280,316,328,364,400,412,
+          448,484,488,524,560,572,608,644,656,692,728,730,
+          757,784,820,847,874,910],
+    # A014192 
+    4 => [0,1,2,3,5,10,15,17,21,25,29,34,38,42,46,51,55,59,
+          63,65,85,105,125,130,150,170,190,195,215,235,255,
+          257,273,289,305,325,341,357,373,393,409,425,441,
+          461,477,493,509,514,530,546,562,582,598,614,630,
+          650,666,682,698,718,734],
+    # A002113
+    10 => [0,1,2,3,4,5,6,7,8,9,11,22,33,44,55,66,77,88,99,
+           101,111,121,131,141,151,161,171,181,191,202,212,
+           222,232,242,252,262,272,282,292,303,313,323,333,
+           343,353,363,373,383,393,404,414,424,434,444,454,
+           464,474,484,494,505,515]
+  }
+
+  @seq.each do |base, ns|
+    ns.each do |n|
+      it "returns true for base-#{base} palindrome #{n}" do
+        n.should be_palindrome(base)
+      end
+    end
+
+    ns.to_seq.invert.sample(10).each do |n|
+      it "returns false for non-base-#{base} palindrome #{n}" do
+        n.should_not be_palindrome(base)
+      end
+    end
+  end
+end

data/spec/numb/pell_lucas_spec.rb

@@ -0,0 +1,28 @@
+# A002203
+PELL_LUCAS = [2,2,6,14,34,82,198,478,1154,2786,6726,16238,
+              39202,94642,228486,551614,1331714,3215042,7761798,
+              18738638,45239074,109216786,263672646,636562078,
+              1536796802,3710155682,8957108166,21624372014,
+              52205852194]
+
+describe Integer, "#pell_lucas?" do
+  PELL_LUCAS.each do |n|
+    it "returns true for Pell-Lucas number #{n}" do
+      n.should be_pell_lucas
+    end
+  end
+
+  PELL_LUCAS.to_seq.invert.sample(10).each do |n|
+    it "returns false for non-Pell-Lucas number #{n}" do
+      n.should_not be_pell_lucas
+    end
+  end
+end
+
+describe Integer, "#pell_lucas" do
+  PELL_LUCAS.each_with_index do |l, n|
+    it "returns #{l} for the #{n}#{n.ordinal} Pell-Lucas number" do
+      n.pell_lucas.should == l
+    end
+  end
+end

data/spec/numb/pentatope_spec.rb

@@ -0,0 +1,15 @@
+# coding: utf-8
+describe Integer, "#pentatope" do
+  # A000332 
+  @seq = [0,0,0,0,1,5,15,35,70,126,210,330,495,715,1001,
+          1365,1820,2380,3060,3876,4845,5985,7315,8855,
+          10626,12650,14950,17550,20475,23751,27405,31465,
+          35960,40920,46376,52360,58905,66045,73815,82251,
+          91390,101270,111930,123410][3..-1]
+
+  @seq.each_with_index do |pe, n|
+    it "returns #{pe} for the #{n}#{n.ordinal} pentatope number" do
+      n.pentatope.should == pe
+    end
+  end
+end

data/spec/numb/primitive_abundant_spec.rb

@@ -0,0 +1,21 @@
+# coding: utf-8
+describe Integer, "#primitive_abundant?" do
+  # A071395
+  @seq = [20,70,88,104,272,304,368,464,550,572,650,748,836,
+          945,1184,1312,1376,1430,1504,1575,1696,1870,1888,
+          1952,2002,2090,2205,2210,2470,2530,2584,2990,3128,
+          3190,3230,3410,3465,3496,3770,3944,4030,4070,4095,
+          4216,4288]
+
+  @seq.each do |n|
+    it "returns true for primitive abundant number #{n}" do
+      n.should be_primitive_abundant
+    end
+  end
+
+  @seq.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-primitive abundant number #{n}" do
+      n.should_not be_primitive_abundant
+    end
+  end
+end

data/spec/numb/pyramidal_spec.rb

@@ -0,0 +1,42 @@
+# coding: utf-8
+describe Integer, "#pyramidal" do
+  @seq = {
+    # A002411
+    5 => [0,1,6,18,40,75,126,196,288,405,550,726,936,1183,
+          1470,1800,2176,2601,3078,3610,4200,4851,5566,6348,
+          7200,8125,9126,10206,11368,12615,13950,15376,
+          16896,18513,20230,22050,23976,26011,28158,30420,
+          32800,35301][1..-1],
+    # A002413
+    7 => [1,8,26,60,115,196,308,456,645,880,1166,1508,1911,
+          2380,2920,3536,4233,5016,5890,6860,7931,9108,
+          10396,11800,13325,14976,16758,18676,20735,22940,
+          25296,27808,30481,33320,36330,39516,42883,46436,
+          50180,54120],
+    # A002414      
+    8 => [1,9,30,70,135,231,364,540,765,1045,1386,1794,
+          2275,2835,3480,4216,5049,5985,7030,8190,9471,
+          10879,12420,14100,15925,17901,20034,22330,24795,
+          27435,30256,33264,36465,39865,43470,47286,51319,
+          55575,60060,64780],
+    # A007586      
+    11 => [0,1,12,42,100,195,336,532,792,1125,1540,2046,
+           2652,3367,4200,5160,6256,7497,8892,10450,12180,
+           14091,16192,18492,21000,23725,26676,29862,33292,
+           36975,40920,45136,49632,54417,59500][1..-1],
+    # A130566       
+    47 => [1,48,186,460,915,1596,2548,3816,5445,7480,9966,
+           12948,16471,20580,25320,30736,36873,43776,51490,
+           60060,69531,79948,91356,103800,117325,131976,
+           147798,164836,183135,202740,223696,246048,269841,
+           295120,321930,350316,380323]
+  }
+
+  @seq.each do |k, ns|
+    ns.to_enum.with_index(1).each do |py, n|
+      it "returns #{py} for #{k}-pyramidal number #{n}" do
+        n.pyramidal(k).should == py
+      end
+    end
+  end
+end

data/spec/numb/q_spec.rb

@@ -0,0 +1,31 @@
+# coding: utf-8
+# A005185
+Q = [1,1,2,3,3,4,5,5,6,6,6,8,8,8,10,9,10,11,11,12,12,
+     12,12,16,14,14,16,16,16,16,20,17,17,20,21,19,20,
+     22,21,22,23,23,24,24,24,24,24,32,24,25,30,28,26,
+     30,30,28,32,30,32,32,32,32,40,33,31,38,35,33,39,
+     40,37,38,40,39]
+
+describe Integer, "#q" do
+  Q.to_enum.with_index(1).each do |q, n|
+    it "returns #{q} for #{n}" do
+      n.q.should == q
+    end
+  end
+end
+
+describe Integer, "#q?" do
+  Q.each do |n|
+    it "returns true for Q-number #{n}" do
+      n.should be_q
+    end
+  end
+
+  # This is unsatisfactory because the abscence of a term in A005185 does not
+  # imply the sequence does not contain it...
+  Q.to_seq.invert.each do |n|
+    it "returns false for non-Q-number #{n}" do
+      n.should_not be_q
+    end
+  end
+end

data/spec/numb/schr/303/266der_spec.rb

@@ -0,0 +1,28 @@
+# coding: utf-8
+SCHRODER = [1,2,6,22,90,394,1806,8558,41586,206098,1037718,
+            5293446,27297738,142078746,745387038,3937603038,
+            20927156706,111818026018,600318853926,
+            3236724317174,17518619320890,95149655201962,
+            518431875418926]
+
+describe Integer, "#schröder" do
+  SCHRODER.each_with_index do |s, n|
+    it "returns #{s} as the #{n}#{n.ordinal} Schröder number" do
+      n.schröder.should == s
+    end
+  end
+end
+
+describe Integer, "#schröder?" do
+  SCHRODER.each do |n|
+    it "returns true for Schröder number #{n}" do
+      n.should be_schröder
+    end
+  end
+
+  SCHRODER.to_seq.invert.sample(100) do |n|
+    it "returns false for non-Schröder number #{n}" do
+      n.should_not be_schröder
+    end
+  end
+end

data/spec/numb/segmented_spec.rb

@@ -0,0 +1,29 @@
+# coding: utf-8
+# A002048 
+SEGMENTED = [1,2,4,5,8,10,14,15,16,21,22,25,26,28,33,34,35,36,
+             38,40,42,46,48,49,50,53,57,60,62,64,65,70,77,80,
+             81,83,85,86,90,91,92,100,104,107,108,116,119,124,
+             127,132,133,137,141,144,145,148,150,151,154,158,
+            159,163,165]
+
+describe Integer, "#segmented" do
+  SEGMENTED.to_enum.with_index(1).each do |s, n|
+    it "returns #{s} as the #{n}#{n.ordinal} segmented number" do
+      n.segmented.should == s
+    end
+  end
+end
+
+describe Integer, "#segmented?" do
+  SEGMENTED.each do |n|
+    it "returns true for segmented number #{n}" do
+      n.should be_segmented
+    end
+  end
+
+  SEGMENTED.to_seq.invert.sample(20).each do |n|
+    it "returns false for non-segmented number #{n}" do
+      n.should_not be_segmented
+    end
+  end
+end

data/spec/numb/square_part_spec.rb

@@ -0,0 +1,15 @@
+# coding: utf-8
+describe Integer, "#square_part" do
+  # A008833
+  @seq = [1,1,1,4,1,1,1,4,9,1,1,4,1,1,1,16,1,9,1,4,1,1,1,4,
+          25,1,9,4,1,1,1,16,1,1,1,36,1,1,1,4,1,1,1,4,9,1,1,
+          16,49,25,1,4,1,9,1,4,1,1,1,4,1,1,9,64,1,1,1,4,1,1,
+          1,36,1,1,25,4,1,1,1,16,81,1,1,4,1,1,1,4,1,9,1,4,1,
+          1,1,16,1]
+
+  @seq.to_enum.with_index(1).each do |sp, n|
+    it "returns #{sp} as the square part of #{n}" do
+      n.square_part.should == sp
+    end
+  end
+end

data/spec/numb/star_spec.rb

@@ -0,0 +1,29 @@
+# coding: utf-8
+# A003154
+STAR = [1,13,37,73,121,181,253,337,433,541,661,793,937,
+        1093,1261,1441,1633,1837,2053,2281,2521,2773,3037,
+        3313,3601,3901,4213,4537,4873,5221,5581,5953,6337,
+        6733,7141,7561,7993,8437,8893,9361,9841,10333,
+        10837]
+
+describe Integer, "#star" do
+  STAR.to_enum.with_index(1).each do |s, n|
+    it "returns #{s} as the #{n}#{n.ordinal} star number" do
+      n.star.should == s
+    end
+  end
+end
+
+describe Integer, "#star?" do
+  STAR.each do |n|
+    it "returns true for star number #{n}" do
+      n.should be_star
+    end
+  end
+
+  STAR.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-star number #{n}" do
+      n.should_not be_star
+    end
+  end
+end

data/spec/numb/stella_octangula_spec.rb

@@ -0,0 +1,29 @@
+# coding: utf-8
+# A007588
+STELLA_OCTANGULA = [0,1,14,51,124,245,426,679,1016,1449,1990,2651,
+                    3444,4381,5474,6735,8176,9809,11646,13699,15980,
+                    18501,21274,24311,27624,31225,35126,39339,43876,
+                    48749,53970,59551,65504,71841,78574,85715,93276,
+                    101269]
+
+describe Integer, "#stella_octangula" do
+  STELLA_OCTANGULA.each_with_index do |o,n|
+    it "returns #{o} as the #{n}#{n.ordinal} stella octangula number" do
+      n.stella_octangula.should == o
+    end
+  end
+end
+
+describe Integer, "#stella_octangula?" do
+  STELLA_OCTANGULA.each do |n|
+    it "returns true for stella octangula number #{n}" do
+      n.should be_stella_octangula
+    end
+  end
+
+  STELLA_OCTANGULA.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-stella-octangula number #{n}" do
+      n.should_not be_stella_octangula
+    end
+  end
+end

data/spec/numb/stirling2_spec.rb

@@ -0,0 +1,23 @@
+# coding: utf-8
+describe Integer, "#stirling2" do
+  @seq = [
+          [1],
+          [1,   1],
+          [1,   3,    1],
+          [1,   7,    6,    1],
+          [1,  15,   25,   10,     1],
+          [1,  31,   90,   65,     15,     1],
+          [1,  63,  301,  350,    140,    21,    1],
+          [1, 127,  966, 1701,   1050,   266,   28,   1],
+          [1, 255, 3025, 7770,   6951,  2646,  462,  36,  1],
+          [1, 511, 9330, 34105, 42525, 22827, 5880, 750, 45, 1]
+  ]
+
+  @seq.to_enum.with_index(1).each do |row, n|
+    row.to_enum.with_index(1).each do |s, m|
+      it "returns #{s} as the #{n}#{n.ordinal} Stirling number of the 2nd kind (m = #{m})" do
+        n.stirling2(m).should == s
+      end
+    end
+  end
+end

data/spec/numb/stirling_spec.rb

@@ -0,0 +1,23 @@
+# coding: utf-8
+describe Integer, "#stirling" do
+  @seq = [
+          [    1                                                       ],
+          [   -1,         1                                            ],
+          [    2,        -3,      1                                    ],
+          [   -6,        11,     -6,      1                            ],
+          [   24,       -50,     35,    -10,     1                     ],
+          [ -120,       274,   -225,     85,   -15,     1              ],
+          [  720,     -1764,   1624,   -735,   175,   -21,   1         ],
+          [-5040,     13068, -13132,   6769, -1960,   322, -28,   1    ],
+          [40320,   -109584, 118124, -67284, 22449, -4536, 546, -36, 1 ],
+          [-362880, 1026576                                            ]
+  ]
+
+  @seq.to_enum.with_index(1).each do |row, n|
+    row.to_enum.with_index(1).each do |s, m|
+      it "returns #{s} as the #{n}#{n.ordinal} Stirling number of the 1st kind (m = #{m})" do
+        n.stirling(m).should == s
+      end
+    end
+  end
+end

data/spec/numb/super_catalan_spec.rb

@@ -0,0 +1,29 @@
+# coding: utf-8
+# A001003 
+SUPER_CATALAN = [1,1,3,11,45,197,903,4279,20793,103049,518859,
+                 2646723,13648869,71039373,372693519,1968801519,
+                 10463578353,55909013009,300159426963,
+                 1618362158587,8759309660445,47574827600981,
+                 259215937709463,1416461675464871]
+
+describe Integer, "#super_catalan" do
+  SUPER_CATALAN.to_enum.with_index(1).each do |s, n|
+    it "returns #{s} as the #{n}#{n.ordinal} super-Catalan number" do
+      n.super_catalan.should == s
+    end
+  end
+end
+
+describe Integer, "#super_catalan?" do
+  SUPER_CATALAN.each do |n|
+    it "returns true for super-Catalan number #{n}" do
+      n.should be_super_catalan
+    end
+  end
+
+  SUPER_CATALAN.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-super-Catalan number #{n}" do
+      n.should_not be_super_catalan
+    end
+  end
+end

data/spec/numb/super_d_spec.rb

@@ -0,0 +1,68 @@
+# coding: utf-8
+describe Integer, "#super_d?" do
+  @seq = {
+    # A032743
+    2 => [19,31,69,81,105,106,107,119,127,131,169,181,190,
+          219,231,247,269,281,310,318,319,331,332,333,334,
+          335,336,337,338,339,348,369,381,419,431,454,469,
+          481,511,519,531,558,569,581,601,619,631,669,679,
+          681,690,715],
+    # A014569       
+    3 => [261,462,471,481,558,753,1036,1046,1471,1645,1752,
+          1848,1923,1926,1968,2031,2231,2232,2363,2395,2471,
+          2591,2610,3058,3087,3148,3163,3172,3181,3471,3494,
+          3542,3851,3884,4143,4269,4314,4471,4527,4554,4620,
+          4710,4732],
+    # A032744      
+    4 => [1168,4972,7423,7752,8431,10267,11317,11487,11549,
+          11680,16588,16664,16837,18257,18597,19784,19933,
+          22217,22504,22819,22829,24078,24331,24514,25296,
+          25698,26685,26738,27812,27973,28988,32466,32467,
+          32735],
+    # A032745       
+    5 => [4602,5517,7539,12955,14555,20137,20379,26629,
+          32767,35689,35825,37706,46020,46715,51988,55170,
+          66344,73338,73974,75390,76157,86025,91497,105852,
+          114488,129550,132234,145550,146399,158651,160897,
+          171673,174782],
+    # A032746      
+    6 => [27257,272570,302693,323576,364509,502785,513675,
+          537771,676657,678146,731378,831122,836553,913797,
+          920456,921269,1045361,1144983,1169054,1283069,
+          1288697,1292673,1343642,1346117,1472078,1523993,
+          1640026],
+    # A032747      
+    7 => [140997,490996,1184321,1259609,1409970,1783166,
+          1886654,1977538,2457756,2714763,2750425,2980991,
+          3043607,3283057,3689639,4191601,4258476,4642725,
+          4909960,4973029,5242829,5349973,5444788,5523544,
+          5682065],
+    # A032748      
+    8 => [185423,641519,1551728,1854230,6415190,12043464,
+          12147605,15517280,16561735,18542300,26908132,
+          29242698,33491333,34982204,35866945,37584428,
+          44263715,45980752,54555936,56148739,60883944,
+          64151900],
+    # A032749      
+    9 => [17546133,32613656,93568867,107225764,109255734,
+          113315082,121251742,175461330,180917907,182557181,
+          190215626,227366552,267791979,297981691,298002027,
+          301321202,321541434,326136560,326687646,329178946,
+          351038353]
+  }
+
+  @seq.each do |d, values|
+    values.each do |n|
+      it "returns true for super-#{d} number #{n}" do
+        n.should be_super_d(d)
+      end
+    end
+
+    values.to_seq.invert.sample(100).each do |n|
+      it "returns false for non-super-#{d} number #{n}" do
+        n.should_not be_super_d(d)
+      end
+    end
+  end
+end
+