numb 0.68.0 → 0.72.0

data/VERSION

@@ -1 +1 @@
-0.68.0
+0.72.0

data/lib/numb.rb

@@ -1,10 +1,11 @@
 # coding: utf-8
 
-libs = %w{abundancy abundant achilles almost_perfect amicable aspiring 
-          automorphic balanced_prime biquadratic breeder brown carmichael carol
-          centered_n_gonal centered_triangular congruum composite coprime
-          core cototient cube d decagonal deficient dodecagonal dihedral_prime
-          dudeney economical emrip equidigital extravagant factorial factorion 
+libs = %w{abundancy abundant achilles almost_perfect almost_prime amicable
+          apocalyptic aspiring augmented_amicable automorphic balanced_prime 
+          binomial biquadratic breeder brown carmichael carol centered_n_gonal 
+          centered_triangular congruum composite coprime core cototient 
+          cube d decagonal deficient dodecagonal dihedral_prime dudeney 
+          economical emrip equidigital extravagant factorial factorion 
           fermat_pseudoprime fibonacci friendly frugal happy harshad 
           heptagonal hexagonal highly_composite highly_abundant hilbert 
           hyperperfect idoneal impolite integer_p interprime 

data/lib/numb/almost_prime.rb

@@ -0,0 +1,6 @@
+# coding: utf-8
+class Integer
+  def almost_prime?(k)
+    Ω == k
+  end
+end

data/lib/numb/apocalyptic.rb

@@ -0,0 +1,5 @@
+class Integer
+  def apocalyptic?
+    (2**self).to_s.include?('666')
+  end
+end

data/lib/numb/augmented_amicable.rb

@@ -0,0 +1,7 @@
+# coding: utf-8
+class Integer
+  def augmented_amicable?(n)
+    m = self
+    [m.σ, n.σ].all?{|sigma| sigma == m + n - 1}
+  end
+end

data/lib/numb/binomial.rb

@@ -0,0 +1,15 @@
+class Integer
+  def binomial?(exp=4)
+    x = self
+    return true if (0..2).include? x
+    (2..exp).each do |n|
+      (1...x).each do |a|
+        an = a**n
+        sign, *terms = an > x ? [:-, an, x] : [:+, x, an]
+        b = (terms.reduce(:-))**(1.0/n.to_f)
+        return true if b.integer? and x == an.send(sign, b**n)
+      end
+    end
+    false
+  end
+end

data/spec/almost_prime_spec.rb

@@ -0,0 +1,73 @@
+describe Integer, "#almost_prime?" do
+
+  @seq = {
+    # 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],
+    # A001358      
+    2 => [4,6,9,10,14,15,21,22,25,26,33,34,35,38,39,46,49,
+          51,55,57,58,62,65,69,74,77,82,85,86,87,91,93,94,
+          95,106,111,115,118,119,121,122,123,129,133,134,
+          141,142,143,145,146,155,158,159,161,166,169,177,
+          178,183,185,187],
+    # A014612 
+    3 => [8,12,18,20,27,28,30,42,44,45,50,52,63,66,68,70,
+          75,76,78,92,98,99,102,105,110,114,116,117,124,125,
+          130,138,147,148,153,154,164,165,170,171,172,174,
+          175,182,186,188,190,195,207,212,222,230,231,236,
+          238,242,244],
+    # A014613      
+    4 => [16,24,36,40,54,56,60,81,84,88,90,100,104,126,132,
+          135,136,140,150,152,156,184,189,196,198,204,210,
+          220,225,228,232,234,248,250,260,276,294,296,297,
+          306,308,315,328,330,340,342,344,348,350,351,364,
+          372,375,376],
+    # A014614      
+    5 => [32,48,72,80,108,112,120,162,168,176,180,200,208,
+          243,252,264,270,272,280,300,304,312,368,378,392,
+          396,405,408,420,440,450,456,464,468,496,500,520,
+          552,567,588,592,594,612,616,630,656,660,675,680,
+          684,688,696],
+    # A046306      
+    6 => [64,96,144,160,216,224,240,324,336,352,360,400,
+          416,486,504,528,540,544,560,600,608,624,729,736,
+          756,784,792,810,816,840,880,900,912,928,936,992,
+          1000,1040,1104,1134,1176,1184,1188,1215,1224,1232,
+          1260,1312,1320],
+    # A046308
+    7 => [128,192,288,320,432,448,480,648,672,704,720,800,
+          832,972,1008,1056,1080,1088,1120,1200,1216,1248,
+          1458,1472,1512,1568,1584,1620,1632,1680,1760,1800,
+          1824,1856,1872,1984,2000,2080,2187,2208,2268,2352,
+          2368,2376],
+    # A046310       
+    8 => [256,384,576,640,864,896,960,1296,1344,1408,1440,
+          1600,1664,1944,2016,2112,2160,2176,2240,2400,2432,
+          2496,2916,2944,3024,3136,3168,3240,3264,3360,3520,
+          3600,3648,3712,3744,3968,4000,4160,4374,4416,4536,
+          4704,4736],
+    # A046312 
+    9 => [512,768,1152,1280,1728,1792,1920,2592,2688,2816,
+          2880,3200,3328,3888,4032,4224,4320,4352,4480,4800,
+          4864,4992,5832,5888,6048,6272,6336,6480,6528,6720,
+          7040,7200,7296,7424,7488,7936,8000,8320,8748,8832,
+          9072,9408],
+  }
+
+  @seq.each do |k, members|
+    members.sample(10).each do |n|
+      it "should return true for #{k}-almost-prime #{n}" do
+        n.almost_prime?(k).should be_true
+      end
+    end
+
+    members.to_seq.invert.sample(10).each do |n|
+      it "should return false for non-#{k}-almost-prime #{n}" do
+        n.almost_prime?(k).should be_false
+      end
+    end
+  end
+end

data/spec/apocalyptic_spec.rb

@@ -0,0 +1,19 @@
+describe Integer, "#apocalyptic?" do
+  # A007356 
+  @seq = [157,192,218,220,222,224,226,243,245,247,251,278,
+          285,286,287,312,355,361,366,382,384,390,394,411,
+          434,443,478,497,499,506,508,528,529,539,540,541,
+          564,578,580,582,583,610].to_seq
+
+  @seq.each do |n|
+    it "should return true for apocalyptic number #{n}" do
+      n.should be_apocalyptic
+    end
+  end
+
+  @seq.invert.sample(100).each do |n|
+    it "should return false for non-apocalyptic number #{n}" do
+      n.should_not be_apocalyptic
+    end
+  end
+end

data/spec/augmented_amicable_spec.rb

@@ -0,0 +1,38 @@
+describe Integer, "#augmented_amicable?" do
+
+         # A007992, A015630
+  @seq = [6160,12220,23500,68908,249424,425500,434784,
+          649990,660825,1017856,1077336,1238380,1252216,
+          1568260,1754536,2166136,2362360,2482536,2537220,
+          2876445,3957525,4177524,4287825,5224660,5559510,
+          5641552].zip(
+         [11697,16005,28917,76245,339825,570405,871585,
+          697851,678376,1340865,2067625,1823925,1483785,
+          1899261,2479065,2580105,4895241,4740505,5736445,
+          3171556,4791916,6516237,4416976,7524525,9868075,
+          7589745])
+
+  @seq.each do |m, n|
+    it "should return true for augmented amicable pair (#{m}, #{n})" do
+      m.augmented_amicable?(n).should be_true
+    end
+
+    it "should return true for augmented amicable pair (#{n}, #{m})" do
+      n.augmented_amicable?(m).should be_true
+    end
+  end
+
+  [[234,234], [78982,72101], [13,7], [100019,276281], [1,0]].each do |m, n|
+    it "should return false for non-augmented amicable pair (#{m}, #{n})" do
+      m.augmented_amicable?(n).should be_false
+    end
+
+    it "should return false for non-augmented amicable pair (#{n}, #{m})" do
+      n.augmented_amicable?(m).should be_false
+    end
+  end
+
+  it "should return false for amicable but non-augmented pair (220, 284)" do
+    220.augmented_amicable?(284).should be_false
+  end
+end

data/spec/binomial_spec.rb

@@ -0,0 +1,20 @@
+describe Integer, "#binomial?" do
+  # A079299
+  @seq = {
+    false => [6,14,22,30,38,42,46,62,66,70,78,86,94,102,110,
+              114,118,134,138,142,150,154,158,166,174,182,186,
+              190,198,206,210,214,222,230,238,246,254,258,262,
+              266,270,278,282,286,294,302,310,318,322,326,330,
+              334,350,354,358].to_seq
+  }
+  @seq[true] = (0..5).to_a + @seq[false].invert.to_a
+
+
+  @seq.each do |key, members|
+    members.sample(10).each do |n|
+      it "should return #{key} for #{'non-' unless key}binomial number #{n}" do
+        n.binomial?.should == key
+      end
+    end
+  end
+end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification 
 name: numb
 version: !ruby/object:Gem::Version 
-  version: 0.68.0
+  version: 0.72.0
 platform: ruby
 authors: 
 - Run Paint Run Run
@@ -9,7 +9,7 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2010-01-15 00:00:00 +00:00
+date: 2010-01-17 00:00:00 +00:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -43,10 +43,14 @@ files:
 - lib/numb/abundant.rb
 - lib/numb/achilles.rb
 - lib/numb/almost_perfect.rb
+- lib/numb/almost_prime.rb
 - lib/numb/amicable.rb
+- lib/numb/apocalyptic.rb
 - lib/numb/aspiring.rb
+- lib/numb/augmented_amicable.rb
 - lib/numb/automorphic.rb
 - lib/numb/balanced_prime.rb
+- lib/numb/binomial.rb
 - lib/numb/biquadratic.rb
 - lib/numb/breeder.rb
 - lib/numb/brown.rb
@@ -163,10 +167,14 @@ files:
 - spec/abundant_spec.rb
 - spec/achilles_spec.rb
 - spec/almost_perfect_spec.rb
+- spec/almost_prime_spec.rb
 - spec/amicable_spec.rb
+- spec/apocalyptic_spec.rb
 - spec/aspiring_spec.rb
+- spec/augmented_amicable_spec.rb
 - spec/automorphic_spec.rb
 - spec/balanced_prime_spec.rb
+- spec/binomial_spec.rb
 - spec/biquadratic_spec.rb
 - spec/breeder_spec.rb
 - spec/brown_spec.rb
@@ -354,6 +362,7 @@ test_files:
 - spec/kynea_spec.rb
 - spec/centered_n_gonal_spec.rb
 - spec/highly_abundant_spec.rb
+- spec/augmented_amicable_spec.rb
 - spec/rough_spec.rb
 - spec/emrip_spec.rb
 - spec/nth_prime_spec.rb
@@ -410,6 +419,7 @@ test_files:
 - spec/interprime_spec.rb
 - spec/number_of_distinct_prime_factors_spec.rb
 - spec/balanced_prime_spec.rb
+- spec/binomial_spec.rb
 - spec/digital_sum_spec.rb
 - spec/brown_spec.rb
 - spec/factorion_spec.rb
@@ -417,6 +427,7 @@ test_files:
 - spec/extravagant_spec.rb
 - spec/self_descriptive_spec.rb
 - spec/decagonal_spec.rb
+- spec/apocalyptic_spec.rb
 - spec/kaprekar_spec.rb
 - spec/hilbert_spec.rb
 - spec/keith_spec.rb
@@ -424,6 +435,7 @@ test_files:
 - spec/myriagonal_spec.rb
 - spec/achilles_spec.rb
 - spec/coprime_spec.rb
+- spec/almost_prime_spec.rb
 - spec/pronic_spec.rb
 - spec/politeness_spec.rb
 - spec/noncototient_spec.rb