numb 0.138.0 → 0.152.0

data/spec/numb/integer_p_spec.rb

@@ -13,15 +13,15 @@ describe Float, "#integer?" do
     end
   end
 
-  [23.5, 4509.0000001, 2222981011.020].each do |n|
+  [23.5, 4509.00001, 2222981011.020].each do |n|
     it "should return false for float #{n}" do
       n.should_not be_integer
     end
   end
 
-  it "ignores digits past the tenth" do
-    3.00000000001.should be_integer
+  it "ignores digits past the seventh" do
+    3.00000001.should be_integer
     123.0000000000032.should be_integer
-    123.00000000032.should_not be_integer
+    123.000032.should_not be_integer
   end
 end

data/spec/numb/knuth_spec.rb

@@ -0,0 +1,29 @@
+# coding: utf-8
+# A007448
+KNUTH = [1,3,3,4,7,7,7,9,9,10,13,13,13,15,15,19,19,19,19,
+         21,21,22,27,27,27,27,27,28,31,31,31,39,39,39,39,
+         39,39,39,39,40,43,43,43,45,45,46,55,55,55,55,55,
+         55,55,55,55,57,57,58,63,63,63,63,63,64,67,67,67,
+         79,79,79,79]
+
+describe Integer, "#knuth?" do
+  KNUTH.uniq.each do |n|
+    it "returns true for Knuth number #{n}" do
+      n.should be_knuth
+    end
+  end
+
+  KNUTH.uniq.to_seq.invert.sample(10).each do |n|
+    it "returns false for non-Knuth number #{n}" do
+      n.should_not be_knuth
+    end
+  end
+end
+
+describe Integer, "#knuth" do
+  KNUTH.each_with_index do |k, n|
+    it "returns #{k} as the #{n}#{n.ordinal} Knuth number" do
+      n.knuth.should == k
+    end
+  end
+end

data/spec/numb/lah_spec.rb

@@ -0,0 +1,61 @@
+# coding: utf-8
+describe Integer, "#lah" do
+  @seq = {
+    # A000142
+    1 => [1,1,2,6,24,120,720,5040,40320,362880,3628800,
+          39916800,479001600,6227020800,87178291200,
+          1307674368000,20922789888000,355687428096000,
+          6402373705728000,121645100408832000,
+          2432902008176640000][1..-1],
+    # A001286
+    2 => [1,6,36,240,1800,15120,141120,1451520,16329600,
+          199584000,2634508800,37362124800,566658892800,
+          9153720576000,156920924160000,2845499424768000,
+          54420176498688000,1094805903679488000,
+          23112569077678080000],
+    # A001754
+    3 => [0,0,1,12,120,1200,12600,141120,1693440,21772800,
+          299376000,4390848000,68497228800,1133317785600,
+          19833061248000,366148823040000,7113748561920000,
+          145120470663168000,3101950060425216000,
+          69337707233034240000,1617879835437465600000].reject(&:zero?),
+    # A001755 
+    4 => [1,20,300,4200,58800,846720,12700800,199584000,
+          3293136000,57081024000,1038874636800,
+          19833061248000,396661224960000,8299373322240000,
+          181400588328960000,4135933413900288000,
+          98228418580131840000],
+    # A001777
+    5 => [1,30,630,11760,211680,3810240,69854400,
+          1317254400,25686460800,519437318400,
+          10908183686400,237996734976000,5394592659456000,
+          126980411830272000,3101950060425216000,
+          78582734864105472000],
+    # A001778
+    6 => [1,42,1176,28224,635040,13970880,307359360,
+          6849722880,155831195520,3636061228800,
+          87265469491200,2157837063782400,55024845126451200,
+          1447576694865100800,39291367432052736000,
+          1100158288097476608000],
+    # A111597
+    7 => [1,56,2016,60480,1663200,43908480,1141620480,
+          29682132480,779155977600,20777492736000,
+          565147802419200,15721384321843200,
+          448059453172531200,13097122477350912000,
+          392913674320527360000,12101741169072242688000],
+    }
+
+  @seq.each do |k, values|
+    (1...k).each do |n|
+      it "returns 0 for Lah(#{n}, #{k})" do
+        n.lah(k).should be_zero 
+      end
+    end
+
+    values.to_enum.with_index(k).each do |lah, n|
+      it "returns #{lah} for Lah(#{n}, #{k})" do
+        n.lah(k).should == lah
+      end
+    end
+  end
+end

data/spec/numb/motzkin_spec.rb

@@ -0,0 +1,30 @@
+# coding: utf-8
+
+# A001006
+MOTZKIN = [1,1,2,4,9,21,51,127,323,835,2188,5798,15511,
+           41835,113634,310572,853467,2356779,6536382,
+           18199284,50852019,142547559,400763223,1129760415,
+           3192727797,9043402501,25669818476,73007772802,
+           208023278209,593742784829].first(15)
+
+describe Integer, "#motzkin" do
+  MOTZKIN.each_with_index.each do |m, n|
+    it "returns #{m} for Motzkin number #{n}" do
+      n.motzkin.should == m
+    end
+  end
+end
+
+describe Integer, "#motzkin?" do
+  MOTZKIN.sample(10).each do |n|
+    it "returns true for Motzkin number #{n}" do
+      n.should be_motzkin
+    end
+  end
+
+  MOTZKIN.to_seq.invert.sample(10).each do |n|
+    it "returns false for non-Motzkin number #{n}" do
+      n.should_not be_motzkin
+    end
+  end
+end

data/spec/numb/multiamicable_spec.rb

@@ -0,0 +1,26 @@
+# coding: utf-8
+describe Integer, "#multiamicable?" do
+  @seq = [
+    [1, 6, 76455288, 183102192],
+    [1, 7, 52920, 152280],
+    [1, 7, 16225560, 40580280],
+    [1, 7, 90863136, 227249568],
+    [1, 7, 16225560, 40580280],
+    [1, 7, 70821324288, 177124806144],
+    [1, 7, 199615613902848, 499240550375424],
+  ].first(5)
+
+  @seq.each do |a, b, m, n|
+    it "returns true for multi-amicable numbers #{n} and #{m} (α = #{a}, β = #{b})" do
+      n.should be_multiamicable(m, a, b)
+    end
+  end
+
+  @seq.each do |a, b, m, n|
+    a += 3
+    n -= 100
+    it "returns false for non-multi-amicable numbers #{n} and #{m} (α = #{a}, β = #{b})" do
+      n.should_not be_multiamicable(m, a, b)
+    end
+  end
+end

data/spec/numb/m/303/251nage_spec.rb

@@ -0,0 +1,31 @@
+# coding: utf-8
+# A000179
+MENAGE = [1,0,0,1,2,13,80,579,4738,43387,439792,4890741,
+          59216642,775596313,10927434464,164806435783,
+          2649391469058,45226435601207,817056406224416,
+          15574618910994665,312400218671253762,
+          6577618644576902053]
+
+describe Integer, "#ménage" do
+  MENAGE.first(20).each_with_index do |m, n|
+    it "returns #{m} for ménage number #{n}" do
+      n.ménage.should == m
+    end
+  end
+end
+
+describe Integer, "#ménage" do
+  @seq = MENAGE.uniq.first(15)
+  
+  @seq.each do |n|
+    it "returns true for ménage number #{n}" do
+      n.should be_ménage
+    end
+  end
+
+  @seq.to_seq.invert.sample(10).each do |n|
+    it "returns false for non-ménage number #{n}" do
+      n.should_not be_ménage
+    end
+  end
+end

data/spec/numb/nexus_spec.rb

@@ -0,0 +1,42 @@
+# coding: utf-8
+describe Integer, "#nexus" do
+  @seq = {
+    0 => [1] * 60,
+    # A005408
+    1 => [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],
+    # A003215
+    2 => [1,7,19,37,61,91,127,169,217,271,331,397,469,547,
+          631,721,817,919,1027,1141,1261,1387,1519,1657,
+          1801,1951,2107,2269,2437,2611,2791,2977,3169,3367,
+          3571,3781,3997,4219,4447,4681,4921,5167,5419,5677,
+          5941,6211,6487],
+    # A005917
+    3 => [1,15,65,175,369,671,1105,1695,2465,3439,4641,
+          6095,7825,9855,12209,14911,17985,21455,25345,
+          29679,34481,39775,45585,51935,58849,66351,74465,
+          83215,92625,102719,113521,125055,137345,150415,
+          164289,178991],
+    # A022521
+    4 => [1,31,211,781,2101,4651,9031,15961,26281,40951,
+          61051,87781,122461,166531,221551,289201,371281,
+          469711,586531,723901,884101,1069531,1282711,
+          1526281,1803001,2115751,2467531],
+    # A022522
+    5 => [1,63,665,3367,11529,31031,70993,144495,269297,
+          468559,771561,1214423,1840825,2702727,3861089,
+          5386591,7360353,9874655,13033657,16954119,
+          21766121,27613783,34655985,43067087]
+  }
+
+  @seq.each do |d, values|
+    values.each_with_index do |nexus, n|
+      it "returns #{nexus} for #{n}" do
+        n.nexus(d).should == nexus
+      end
+    end
+  end
+end

data/spec/numb/nsw_spec.rb

@@ -0,0 +1,32 @@
+# coding: utf-8
+
+# A002315
+NSW = [1,7,41,239,1393,8119,47321,275807,1607521,
+       9369319,54608393,318281039,1855077841,10812186007,
+       63018038201,367296043199,2140758220993,
+       12477253282759,72722761475561,423859315570607,
+       2470433131948081]
+
+describe Integer, "#nsw" do
+  NSW.each_with_index do |nsw, n|
+    it "returns #{nsw} as the #{n}#{n.ordinal} NSW number" do
+      n.nsw.should == nsw
+    end
+  end
+end
+
+describe Integer, "#nsw?" do
+  @seq = NSW
+
+  @seq.each do |n|
+    it "returns true for NSW number #{n}" do
+      n.should be_nsw
+    end
+  end
+
+  @seq.to_seq.invert.sample(100).each do |n|
+    it "returns false for non-NSW number #{n}" do
+      n.should_not be_nsw
+    end
+  end
+end

data/spec/numb/pandigital_spec.rb

@@ -0,0 +1,29 @@
+# coding: utf-8
+describe Integer, "#pandigital?" do
+  # A050278  
+  @seq = [1023456789,1023456798,1023456879,1023456897,
+          1023456978,1023456987,1023457689,1023457698,
+          1023457869,1023457896,1023457968,1023457986,
+          1023458679,1023458697,1023458769,1023458796,
+          1023458967,1023458976] 
+  @seq.each do |n|
+    it "returns true for pandigital number #{n}" do
+      n.should be_pandigital
+    end
+  end
+  
+  # A050289
+  [123456789,123456798,123456879,123456897,123456978,123456987,123457689,
+   123457698,123457869,123457896,123457968,123457986,123458679,123458697,
+   123458769,123458796,123458967,123458976].each do |n|
+    it "returns false for zeroless pandigital number #{n}" do
+      n.should_not be_pandigital
+    end
+  end
+
+  [19, 123456780, 987653211, 474839302029273675, 100000001, 34, 334455].each do |n|
+    it "returns false for non-pandigital number #{n}" do
+      n.should_not be_pandigital
+    end
+  end
+end

data/spec/numb/persistent_spec.rb

@@ -0,0 +1,44 @@
+# coding: utf-8
+describe Integer, "#persistent?" do
+
+  @seq = {
+    # A051264
+    1 => [1023456798,1023456897,1023456978,1023456987,
+          1023457698,1023457896,1023457968,1023457986,
+          1023458697,1023458796,1023458967,1023458976,
+          1023459678,1023459687,1023459768,1023459786,
+          1023459867,1023459876],
+    # A051018      
+    2 => [1023456789,1023456879,1023457689,1023457869,
+          1023458679,1023458769,1023465789,1023465879,
+          1023467589,1023467859,1023468579,1023468759,
+          1023475689,1023475869,1023476589,1023476859,
+          1023478569,1023478659],
+    # A051019      
+    3 => [1052674893,1052687493,1052746893,1052748693,
+          1052867493,1052874693,1053267489,1053268749,
+          1053274869,1053286749,1053287469,1065273489,
+          1065287349,1067285493,1067328549,1068547293,
+          1068547329,1068549273],
+    # A051020      
+    4 => [1253094867,1267085493,1268547309,1269085473,
+          1273085469,1308547269,1308549267,1326854907,
+          1327068549,1328746905,1450687329,1450732869,
+          1450867293,1450928673,1452687309,1452690873,
+          1452730869,1453087269],
+  }
+
+  @seq.each do |n,ks|
+    ks.each do |k|
+      it "returns true for #{n}-persistent number #{k}" do
+        k.should be_persistent(n)
+      end
+    end
+
+    ks.to_seq.invert.sample(20).each do |k|
+      it "returns false for non-#{n}-persistent number #{k}" do
+        k.should_not be_persistent(n)
+      end
+    end
+  end
+end

data/spec/numb/prime_power_spec.rb

@@ -0,0 +1,66 @@
+# coding: utf-8
+describe Integer, "#prime_power?" do
+  @seq = {
+    # A000961
+    0 => [1,2,3,4,5,7,8,9,11,13,16,17,19,23,25,27,29,31,32,
+          37,41,43,47,49,53,59,61,64,67,71,73,79,81,83,89,
+          97,101,103,107,109,113,121,125,127,128,131,137,
+          139,149,151,157,163,167,169,173,179,181,191,193,
+          197,199,211,223,227],
+    # A000040
+    1 => [2,3,5,7,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,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],
+    # A001248      
+    2 => [4,9,25,49,121,169,289,361,529,841,961,1369,1681,
+          1849,2209,2809,3481,3721,4489,5041,5329,6241,6889,
+          7921,9409,10201,10609,11449,11881,12769,16129,
+          17161,18769,19321,22201],
+    # A030078
+    3 => [8,27,125,343,1331,2197,4913,6859,12167,24389,
+          29791,50653,68921,79507,103823,148877,205379,
+          226981,300763,357911,389017,493039,571787,704969,
+          912673,1030301,1092727,1225043,1295029,1442897],
+    # A030514          
+    4 => [16,81,625,2401,14641,28561,83521,130321,279841,
+          707281,923521,1874161,2825761,3418801,4879681,
+          7890481,12117361,13845841,20151121,25411681,
+          28398241,38950081,47458321,62742241,88529281,
+          104060401],
+    5 => [32,243,3125,16807,161051,371293,1419857,2476099,
+          6436343,20511149,28629151,69343957,115856201,
+          147008443,229345007,418195493,714924299,844596301,
+          1350125107,1804229351,2073071593,3077056399,
+          3939040643],
+  }
+
+  @seq[0].each do |n|
+    it "returns true for prime-power #{n}" do
+      n.should be_prime_power
+    end
+  end
+
+  @seq[0].to_seq.invert.sample(10).each do |n|
+    it "returns false for non-prime-power #{n}" do
+      n.should_not be_prime_power
+    end
+  end
+
+  @seq.each do |k, values|
+    next if k.zero?
+
+    values.each do |n|
+      it "returns true for #{k}#{k.ordinal}-prime-power #{n}" do
+        n.should be_prime_power(k)
+      end
+    end
+
+    values.to_seq.invert.sample(10).each do |n|
+      it "returns false for non-#{k}#{k.ordinal}-prime-power #{n}" do
+        n.should_not be_prime_power(k)
+      end
+    end
+  end
+end

data/spec/numb/rough_spec.rb

@@ -1,48 +1,48 @@
+# coding: utf-8
 describe Integer, "#rough?" do
-
-  ROUGH = {
+  @seq = {
     # 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],
+    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],
+    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       
+    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       
+    # 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      
+    # 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
+    # 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,
@@ -51,19 +51,19 @@ describe Integer, "#rough?" do
     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]
+           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
+  @seq.each do |k, values|
+    values.each do |n|
+      it "returns true for #{k}-rough number #{n}" do
+        n.should be_rough(k)
+      end
+    end
+
+    values.to_seq.invert.sample(5).each do |n|
+      it "returns false for non-#{k}-rough number #{n}" do
+        n.should_not be_rough(k)
       end
     end
   end

data/spec/numb/square_spec.rb

@@ -1,24 +1,28 @@
 describe Integer, "#square?" do
-  SQUARES = [0,1,4,9,16,25,36,49,64,81,100,121,144,169,196,
-             225,256,289,324,361,400,441,484,529,576,625,676,
-             729,784,841,900,961,1024,1089,1156,1225,1296,1369,
-             1444,1521,1600,1681,1764,1849]
+  @seq = [0,1,4,9,16,25,36,49,64,81,100,121,144,169,196,
+          225,256,289,324,361,400,441,484,529,576,625,676,
+          729,784,841,900,961,1024,1089,1156,1225,1296,1369,
+          1444,1521,1600,1681,1764,1849]
   
-  it "returns true if the number is square" do
-    SQUARES.each do |number|
-      number.should be_square
+  @seq.each do |n|
+    it "returns true for square number #{n}" do
+      n.should be_square
     end
   end
 
-  it "returns false if the number isn't square" do
-    ((1..1850).to_a - SQUARES).each do |number|
-      number.should_not be_square
+  @seq.to_seq.invert.each do |n|
+    it "returns false for non-square number #{n}" do
+      n.should_not be_square
     end
   end
 
   it "returns false for negative numbers" do
-    [-1, -378, -9].each do |number|
-      number.should_not be_square
+    [-1, -378, -9].each do |n|
+      n.should_not be_square
     end
   end
+
+  it "returns nil when the square root is infinite" do
+    9999999880000000359859986946294623781838875946012712893365986548632235701967954558670969253503997913415008896487205684793069342787691170258832794954572186327945677037037154860793120217960493002317298159645342431034751955276791509386257991014312642060444444445753597283950617283950617283950617283950617283950617283950617283950617283950617284297678218281.square?.should be_nil
+  end
 end

data/spec/numb/sum_of_infinitary_divisors_spec.rb

@@ -0,0 +1,15 @@
+# coding: utf-8
+describe Integer, "#σ∞" do
+  # A049417 
+  @seq = [1,3,4,5,6,12,8,15,10,18,12,20,14,24,24,17,18,30,
+          20,30,32,36,24,60,26,42,40,40,30,72,32,51,48,54,
+          48,50,38,60,56,90,42,96,44,60,60,72,48,68,50,78,
+          72,70,54,120,72,120,80,90,60,120,62,96,80,85,84,
+          144,68,90]
+
+  @seq.to_enum.with_index(1).each do |sum, n|
+    it "returns #{sum} for #{n}" do
+      n.σ∞.should == sum
+    end
+  end
+end

data/spec/numb/superabundant_spec.rb

@@ -1,18 +1,18 @@
 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]
+  @seq = [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].first(17)
 
-  SUPERABUNDANT.first(20).each do |n|
+  @seq.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|
+  @seq.to_seq.invert.sample(20).each do |n|
     it "should return false for non-superabundant number #{n}" do
       n.should_not be_superabundant
     end

data/spec/numb/vampire_spec.rb

@@ -1,19 +1,19 @@
 describe Integer, "#vampire?" do
-  # http://www.research.att.com/~njas/sequences/A014575
-  VAMPIRE = [1260,1395,1435,1530,1827,2187,6880,102510,104260,
-             105210,105264,105750,108135,110758,115672,116725,
-             117067,118440,120600,123354,124483,125248,125433,
-             125460,125500,126027,126846,129640]
+  #A014575
+  @seq = [1260,1395,1435,1530,1827,2187,6880,102510,104260,
+          105210,105264,105750,108135,110758,115672,116725,
+          117067,118440,120600,123354,124483,125248,125433,
+          125460,125500,126027,126846,129640]
   
-  it "returns true for vampire numbers" do
-    VAMPIRE.each do |number|
-      number.should be_vampire 
+  @seq.each do |n|
+    it "returns true for vampire number #{n}" do
+      n.should be_vampire 
     end
   end
   
-  it "returns false for non-vampire numbers" do
-    [1259,22,1396].each do |number|
-      number.should_not be_vampire 
+  @seq.to_seq.invert.sample(10).each do |n|
+    it "returns false for non-vampire number #{n}" do
+      n.should_not be_vampire 
     end
   end
 end

data/spec/spec_helper.rb

@@ -25,7 +25,7 @@ class Seq
   end
 
   def invert
-    return [] if @include.size < 2
+    return [] if @include.size < 2 or @include.size == self.end.succ - self.begin
     Seq.new(self.begin..self.end).tap{|s| s.exclude = @include }
   end