numb 0.184.0 → 0.185.0

data/VERSION

@@ -1 +1 @@
-0.184.0
+0.185.0

data/lib/numb/divisors.rb

@@ -472,7 +472,7 @@ class Integer
   #     56.xgcd(72) #=> [8, 4, -3]
   #
   # @param [Integer] b the number `self` is divided by
-  # @returns [Array<Integer>] GCD, `x`, and `y`
+  # @return [Array<Integer>] GCD, `x`, and `y`
   def xgcd b
     a = self
     x, y, u, v = 0, 1, 1, 0

data/lib/numb/kronecker.rb

@@ -0,0 +1,113 @@
+# -*- coding: utf-8 -*-
+class Integer
+  # Returns the Kronecker symbol for `self`/`b`
+  #
+  # The Kronecker symbol is an extension of the Jacobi symbol to all
+  # integers.
+  #
+  #    2.kronecker 3 #=> -1
+  #    -8.kronecker 3 #=> 1
+  #    6.kronecker 7 #=> -1
+  #
+  # @param [Integer] b The "denominator"
+  # @return [Integer] -1, 0, 1 
+  def kronecker b
+    # The following algorithm is excerpted from Henri Cohen's A Course
+    # in Computational Algebraic Number Theory, 3rd. ed. It appears on
+    # page 48, and is entitled "Algorithm 1.4.12 (Kronecker-Binary).
+
+    # Given a, b ∈ ℤ, this algorithm computes the Kronecker symbol
+    # (a/b) (hence the Legendre symbol when b is an odd prime).
+    raise ArgumentError unless b.is_a?(Integer)
+
+    a = self
+    # The following Array is a lookup table for computing (-1)^(a^2 -
+    # 1)/8 (Cohen, p. 29)
+    tab2 = [0, 1, 0, -1, 0, -1, 0, 1]
+
+    # 1. [Test b = 0] If b = 0, then output 0 if |a| ≠ 1, 1 if |a| = 1
+    # and terminate the algorithm
+    return a.abs == 1 ? 1 : 0 if b.zero?
+
+    # 2. [Remove 2's from b] If a and b are both even, output 0 and
+    # terminate the algorithm. Otherwise, set v ← 0 and while b is
+    # even set v ← v + 1 and b ← b/2. Then if v is even set k ← 1,
+    # otherwise set k ← (-1)^(a^2 - 1)/8 (by table lookup, not by
+    # computing (a^2 - 1)/8). Finally, if b < 0 set b ← -b, and if in
+    # addition a < 0 set k ← -k.
+    return 0 if a.even? and b.even?
+    v, k = 0, nil
+    while b.even?
+      v += 1
+      b = b.quo(2).to_i
+    end
+    k = v.even? ? 1 : tab2[a&7]
+    if b < 0
+      b = -b
+      k =-k if a < 0
+    end
+
+    # 3. [Reduce size once] (Here b is odd and b > 0.) Set a ← a mod b
+    raise unless b.odd? and b > 0
+    a %= b
+
+    loop do
+      # 4. [Finished?] If a = 0, output 0 if b > 1, k if b = 1, and
+      # terminate the algorithm
+      return (b > 1) ? 0 : k if a.zero?
+
+      # 5. [Remove powers of 2] Set v ← 0 and, while a is even, set v ←
+      # v + 1 and a ← a/2. If v is odd, set k ← (-1)^((b^2 - 1)/8) * k
+      v = 0
+      while a.even?
+        v += 1
+        a = a.quo(2).to_i
+      end
+      k = ((-1)**((b**2).pred.quo 8) * k).to_i if v.odd?
+      
+      # 6. [Subtract and apply reciprocity] (Here a and b are odd.) Set
+      # r ← b - a. If r > 0, then set k ← (-1)^((a-1)(b-1)/4) * k (using
+      # if statements), b ← a and a ← r, else set a ← -r. Go to step 4.
+      raise unless a.odd? and b.odd?
+      r = b - a
+      if r > 0
+        # The optimisation below is from Cohen, p. 29
+        k = -k unless (a & b & 2).zero?
+        b = a
+        a = r
+      else
+        a = -r
+      end   
+    end 
+  end
+
+  # Returns the Legendre symbol for `self`/`p`
+  #
+  # The Legendre symbol is a multiplicative function with values 1,
+  # −1, 0: its value on a (nonzero) quadratic residue mod `p` is 1 and
+  # on a quadratic non-residue is −1. It is defined only when `p` is
+  # an odd prime.
+  # 
+  #    12345.legendre 331 #=> -1
+  #    0.legendre 83 #=> 0
+  #    55.legendre 10 #=> nil
+  #    6.legendre 7 #=> -1
+  #
+  # @param [Integer] p The "denominator"
+  # @return [Integer, nil] -1, 0, 1, or `nil` 
+  def legendre p
+    kronecker(p) if p.odd? and p > 2 and p.prime?
+  end
+
+  # Returns the Jacobi symbol for `a`/`n`
+  #
+  # For any integer `a` and any positive odd integer `n` the Jacobi
+  # symbol is defined as the product of the Legendre symbols
+  # corresponding to the prime factors of `n`.
+  # 
+  # @param [Integer] p The "denominator"
+  # @return [Integer, nil] -1, 0, 1, or `nil` 
+  def jacobi n
+    kronecker(n) if n.odd? and n > 0
+  end
+end

data/spec/numb/jacobi_spec.rb

@@ -0,0 +1,65 @@
+# -*- coding: utf-8 -*-
+describe Integer, "#jacobi" do
+
+  @seq = {
+    # A165591
+    59701 => [0,1,-1,1,1,1,-1,-1,-1,1,-1,1,1,-1,1,1,1,-1,-1,-1,
+              1,-1,-1,1,-1,1,1,1,-1,-1,-1,-1,-1,1,1,-1,1,-1,1,
+              -1,-1,1,1,1,1,1,-1,-1,1,1,-1,-1,-1,-1,-1,1,1,-1,1,
+              -1,1,-1,1,-1,1,-1,-1,1,-1,1,1,-1,-1,-1,1,1,-1,-1,
+              1,-1,1],
+    # A165596
+    59881 => [0,1,1,1,1,1,1,-1,1,1,1,-1,1,1,-1,1,1,-1,1,-1,1,
+              -1,-1,1,1,1,1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1,1,
+              -1,1,-1,1,1,1,1,1,1,-1,1,1,1,-1,-1,-1,1,-1,1,-1,1,
+              -1,1,1,-1,-1,-1,1,-1,-1,1,-1,-1,1,-1,1,1,-1,1,1,1,
+              1,-1,-1],
+  }
+
+  @seq.each do |a, ns|
+    ns.each_with_index do |j, n|
+      it "returns #{j} for #{a}.jacobi(#{n})" do
+        n.jacobi(a).should == j
+      end
+    end
+  end
+
+  # A077009
+  [1,-1,1,-1,0,-1,1,1,1,-1,0,-1,0,0,1,-1,-1,-1,1,0,
+   1,-1,0,-1,0,-1,1,0,0,-1,1,0,-1,-1,1,-1,1,0,1,-1,0,
+   -1,1,1,1,-1,0,1,1,0,1,-1,0,-1,1,0,1,1,0,1,0,-1,0,
+   -1,0,-1,1,0,1,-1,1,-1,-1,0,1,-1,0,0,1,1,1,-1,0,-1,
+   0,0,1,0,-1,-1,1,0,1,1,0,-1,1,0,1,-1].each_with_index do |j, n|
+    n = (2 * n).succ
+    a = n.φ
+    it "returns #{j} for #{a}.jacobi(#{n})" do
+      a.jacobi(n).should == j
+    end
+  end
+
+  # A077010
+  [1,1,1,1,1,1,1,0,1,1,-1,1,1,1,1,1,0,1,1,-1,1,1,0,
+   1,1,0,1,1,-1,1,1,-1,-1,1,0,1,1,1,1,1,1,1,1,0,1,0,
+   -1,0,1,0,1,1,0,1,1,-1,1,1,0,1,1,0,1,1,-1,1,1,0,1,
+   1,0,1,0,0,1,1,0,1,1,0,-1,1,0,1,1,-1,1,-1,0,1,1,-1,
+   1,1,1,1,1,0,1,1].each_with_index do |j, n|
+    n = (2 * n).succ
+    a = n.σ
+    it "returns #{j} for #{a}.jacobi(#{n})" do
+      a.jacobi(n).should == j
+    end
+  end
+  
+  [2, 10, -4, 0].each do |n|
+    10.times do |a|
+      it "returns nil for #{a}.jacobi(#{n})" do
+        a.jacobi(n).should be_nil
+      end
+    end
+  end
+end
+  
+
+
+
+

data/spec/numb/kronecker_spec.rb

@@ -0,0 +1,76 @@
+describe Integer, "#kronecker" do
+  @seq = {
+    # A109017
+    -6 => [0,1,0,0,0,1,0,1,0,0,0,1,0,-1,0,0,0,-1,0,-1,0,0,0,
+           -1,0,1,0,0,0,1,0,1,0,0,0,1,0,-1,0,0,0,-1,0,-1,0,0,
+           0,-1,0,1,0,0,0,1,0,1,0,0,0,1,0,-1,0,0,0,-1,0,-1,0,
+           0,0,-1,0,1,0,0,0,1,0,1,0,0,0,1,0,-1,0,0,0,-1,0,-1,
+           0,0,0,-1,0,1,0,0,0,1,0,1,0].map.with_index,
+    # -5..-2 are extracted from <http://mathworld.wolfram.com/KroneckerSymbol.html>
+    -5 => [1, -1, 1, 1, 0, -1, 1, -1, 1, 0, -1, 1, -1, -1, 0, 1, -1,
+           -1, -1, 0].map.with_index(1),
+    -4 => [1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0,
+           -1, 0].map.with_index(1),
+    -3 => [1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0, 1, -1, 0,
+           1, -1].map.with_index(1),
+    -2 => [1, 0, 1, 0, -1, 0, -1, 0, 1, 0, 1, 0, -1, 0, -1, 0, 1, 0,
+           1, 0].map.with_index(1),
+    # A034947
+    -1 => [1,1,-1,1,1,-1,-1,1,1,1,-1,-1,1,-1,-1,1,1,1,-1,1,
+           1,-1,-1,-1,1,1,-1,-1,1,-1,-1,1,1,1,-1,1,1,-1,-1,1,
+           1,1,-1,-1,1,-1,-1,-1,1,1,-1,1,1,-1,-1,-1,1,1,-1,
+           -1,1,-1,-1,1,1,1,-1,1,1,-1,-1,1,1,1,-1,-1,1,-1,-1,
+           1,1].map.with_index(1),
+    # 0..1 are extracted from <http://mathworld.wolfram.com/KroneckerSymbol.html>
+    0 => [1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+          0].map.with_index(1),
+    1 => [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1].map.with_index(1),
+    # A091337
+    2 => [1,0,-1,0,-1,0,1,0,1,0,-1,0,-1,0,1,0,1,0,-1,0,-1,
+          0,1,0,1,0,-1,0,-1,0,1,0,1,0,-1,0,-1,0,1,0,1,0,-1,
+          0,-1,0,1,0,1,0,-1,0,-1,0,1,0,1,0,-1,0,-1,0,1,0,1,
+          0,-1,0,-1,0,1,0,1,0,-1,0,-1,0,1,0,1,0,-1,0,-1,0,1,
+          0,1,0,-1,0,-1,0,1,0,1,0,-1,0,-1,0,1,0,1].map.with_index(1),
+    # A091338
+    3 => [1,-1,0,1,-1,0,-1,-1,0,1,1,0,1,1,0,1,-1,0,-1,-1,0,
+          -1,1,0,1,-1,0,-1,-1,0,-1,-1,0,1,1,0,1,1,0,1,-1,0,
+          -1,1,0,-1,1,0,1,-1,0,1,-1,0,-1,1,0,1,1,0,1,1,0,1,
+          -1,0,-1,-1,0,-1,1,0,1,-1,0,-1,-1,0,-1,-1,0,1,1,0,
+          1,1,0,-1,-1,0,-1,1,0,-1,1,0,1,-1,0,1,-1,0].map.with_index(1),
+    # A080891
+    5 => [0,1,-1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,
+          0,1,-1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0,
+          1,-1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0,1,
+          -1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0,1,
+          -1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0,1,-1,-1,1,0].map.with_index,
+    # The below is extracted from <http://mathworld.wolfram.com/KroneckerSymbol.html>
+    6 => [1, 0, 0, 0, 1, 0, -1, 0, 0, 0, -1, 0, -1, 0, 0, 0, -1, 0, 1,
+          0].map.with_index(1)
+  }
+
+  @seq.each do |a, bs|
+    bs.each do |k, b|
+      it "returns #{k} for #{a}.kronecker(#{b})" do
+        a.kronecker(b).should == k
+      end
+    end
+  end
+
+  # A071935, A071936 
+  [1,-1,1,1,1,-1,1,1,1,-1,-1,1,1,-1,1,1,1,-1,1,1,1,
+  -1,-1,1,1,-1,-1,1,1,-1,1,1,1,-1,1,1,1,-1,1,1,1,-1,
+  -1,1,1,-1,-1,1,1,-1,1,1,1,-1,-1,1,1,-1,-1,1,1,-1,
+  1,1,1,-1,1,1,1,-1,1,1,1,-1,-1,1,1,-1,1,1,1,-1,1,1,
+  1,-1,-1,1,1,-1,-1,1,1,-1,-1,1,1,-1,1,1,1,-1,1,1,1,
+  -1,-1].map.with_index(1){|k, n| [n, n.succ, k]} + \
+  [1,-1,1,1,1,1,1,1,1,-1,1,1,1,1,1,1,1,-1,1,1,1,1,1,
+  1,1,-1,1,1,1,1,1,1,1,-1,1,1,1,1,1,1,1,-1,1,1,1,1,
+  1,1,1,-1,1,1,1,1,1,1,1,-1,1,1,1,1,1,1,1,-1,1,1,1,
+  1,1,1,1,-1,1,1,1,1,1,1,1,-1,1,1,1,1,1,1,1,-1,1,1,
+  1,1,1,1,1,-1,1,1,1,1,1,1,1,-1,1,1,1,1,1,1,1,-1,1,
+  1,1,1,1,1].map.with_index(1){|k, n| [n.succ, n, k]}.each do |a, b, k| 
+    it "returns #{k} for #{a}.kronecker(#{b})" do
+      a.kronecker(b).should == k 
+    end 
+  end 
+end

data/spec/numb/legendre_spec.rb

@@ -0,0 +1,142 @@
+describe Integer, "#legendre" do
+
+  @seq = {
+    # A011582
+    11 => [0,1,-1,1,1,1,-1,-1,-1,1,-1,0,1,-1,1,1,1,-1,-1,-1,
+           1,-1,0,1,-1,1,1,1,-1,-1,-1,1,-1,0,1,-1,1,1,1,-1,
+           -1,-1,1,-1,0,1,-1,1,1,1,-1,-1,-1,1,-1,0,1,-1,1,1,
+           1,-1,-1,-1,1,-1,0,1,-1,1,1,1,-1,-1,-1,1,-1,0,1,-1,
+           1],
+    # A011583
+    13 => [0,1,-1,1,1,-1,-1,-1,-1,1,1,-1,1,0,1,-1,1,1,-1,-1,
+           -1,-1,1,1,-1,1,0,1,-1,1,1,-1,-1,-1,-1,1,1,-1,1,0,
+           1,-1,1,1,-1,-1,-1,-1,1,1,-1,1,0,1,-1,1,1,-1,-1,-1,
+           -1,1,1,-1,1,0,1,-1,1,1,-1,-1,-1,-1,1,1,-1,1,0,1,
+           -1],
+    # A011584
+    17 => [0,1,1,-1,1,-1,-1,-1,1,1,-1,-1,-1,1,-1,1,1,0,1,1,
+           -1,1,-1,-1,-1,1,1,-1,-1,-1,1,-1,1,1,0,1,1,-1,1,-1,
+           -1,-1,1,1,-1,-1,-1,1,-1,1,1,0,1,1,-1,1,-1,-1,-1,1,
+           1,-1,-1,-1,1,-1,1,1,0,1,1,-1,1,-1,-1,-1,1,1,-1,-1,
+           -1],
+    # A011585
+    19 => [0,1,-1,-1,1,1,1,1,-1,1,-1,1,-1,-1,-1,-1,1,1,-1,0,
+           1,-1,-1,1,1,1,1,-1,1,-1,1,-1,-1,-1,-1,1,1,-1,0,1,
+           -1,-1,1,1,1,1,-1,1,-1,1,-1,-1,-1,-1,1,1,-1,0,1,-1,
+           -1,1,1,1,1,-1,1,-1,1,-1,-1,-1,-1,1,1,-1,0,1,-1,-1,
+           1],
+
+    # A011586
+    23 => [0,1,1,1,1,-1,1,-1,1,1,-1,-1,1,1,-1,-1,1,-1,1,-1,
+           -1,-1,-1,0,1,1,1,1,-1,1,-1,1,1,-1,-1,1,1,-1,-1,1,
+           -1,1,-1,-1,-1,-1,0,1,1,1,1,-1,1,-1,1,1,-1,-1,1,1,
+           -1,-1,1,-1,1,-1,-1,-1,-1,0,1,1,1,1,-1,1,-1,1,1,-1,
+           -1],
+    # A011587
+    29 => [0,1,-1,-1,1,1,1,1,-1,1,-1,-1,-1,1,-1,-1,1,-1,-1,
+           -1,1,-1,1,1,1,1,-1,-1,1,0,1,-1,-1,1,1,1,1,-1,1,-1,
+           -1,-1,1,-1,-1,1,-1,-1,-1,1,-1,1,1,1,1,-1,-1,1,0,1,
+           -1,-1,1,1,1,1,-1,1,-1,-1,-1,1,-1,-1,1,-1,-1,-1,1,
+           -1,1],
+    # A011588
+    31 => [0,1,1,-1,1,1,-1,1,1,1,1,-1,-1,-1,1,-1,1,-1,1,1,1,
+           -1,-1,-1,-1,1,-1,-1,1,-1,-1,0,1,1,-1,1,1,-1,1,1,1,
+           1,-1,-1,-1,1,-1,1,-1,1,1,1,-1,-1,-1,-1,1,-1,-1,1,
+           -1,-1,0,1,1,-1,1,1,-1,1,1,1,1,-1,-1,-1,1,-1,1,-1,
+           1],
+    # A011589
+    37 => [0,1,-1,1,1,-1,-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,
+           -1,-1,1,-1,-1,-1,1,1,1,1,-1,1,-1,-1,1,1,-1,1,0,1,
+           -1,1,1,-1,-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,-1,-1,
+           1,-1,-1,-1,1,1,1,1,-1,1,-1,-1,1,1,-1,1,0,1,-1,1,1,
+           -1,-1],
+    # A011590
+    41 => [0,1,1,-1,1,1,-1,-1,1,1,1,-1,-1,-1,-1,-1,1,-1,1,
+           -1,1,1,-1,1,-1,1,-1,-1,-1,-1,-1,1,1,1,-1,-1,1,1,
+           -1,1,1,0,1,1,-1,1,1,-1,-1,1,1,1,-1,-1,-1,-1,-1,1,
+           -1,1,-1,1,1,-1,1,-1,1,-1,-1,-1,-1,-1,1,1,1,-1,-1,
+           1,1,-1,1],
+    # A011591
+    43 => [0,1,-1,-1,1,-1,1,-1,-1,1,1,1,-1,1,1,1,1,1,-1,-1,
+           -1,1,-1,1,1,1,-1,-1,-1,-1,-1,1,-1,-1,-1,1,1,-1,1,
+           -1,1,1,-1,0,1,-1,-1,1,-1,1,-1,-1,1,1,1,-1,1,1,1,1,
+           1,-1,-1,-1,1,-1,1,1,1,-1,-1,-1,-1,-1,1,-1,-1,-1,1,
+           1,-1],
+    # A011626
+    227 => [0,1,-1,1,1,-1,-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,
+            -1,1,-1,1,-1,1,1,1,1,1,1,-1,-1,1,1,-1,1,-1,-1,-1,
+            1,-1,-1,1,1,-1,-1,1,1,1,-1,-1,-1,1,-1,-1,-1,1,-1,
+            1,-1,-1,1,1,1,1,-1,-1,-1,1,1,1,-1,1,1,1,1,1,1,1,
+            -1],
+    # A011627
+    229 => [0,1,-1,1,1,1,-1,-1,-1,1,-1,1,1,-1,1,1,1,1,-1,1,1,
+            -1,-1,-1,-1,1,1,1,-1,-1,-1,-1,-1,1,-1,-1,1,1,-1,
+            -1,-1,-1,1,1,1,1,1,-1,1,1,-1,1,-1,1,-1,1,1,1,1,-1,
+            1,1,1,-1,1,-1,-1,-1,1,-1,1,1,-1,-1,-1,1,1,-1,1,-1,
+            1],
+    # A011628
+    233 => [0,1,1,-1,1,-1,-1,1,1,1,-1,-1,-1,1,1,1,1,-1,1,1,
+            -1,-1,-1,1,-1,1,1,-1,1,1,1,1,1,1,-1,-1,1,1,1,-1,
+            -1,-1,-1,-1,-1,-1,1,-1,-1,1,1,1,1,-1,-1,1,1,-1,1,
+            -1,1,-1,1,1,1,-1,1,-1,-1,-1,-1,1,1,-1,1,-1,1,-1,
+            -1,-1,-1],
+    # A165573
+    257 => [0,1,1,-1,1,-1,-1,-1,1,1,-1,1,-1,1,-1,1,1,1,1,-1,
+            -1,1,1,1,-1,1,1,-1,-1,1,1,1,1,-1,1,1,1,-1,-1,-1,
+            -1,-1,1,-1,1,-1,1,-1,-1,1,1,-1,1,-1,-1,-1,-1,1,1,
+            1,1,1,1,-1,1,-1,-1,1,1,-1,1,-1,1,1,-1,-1,-1,-1,-1,
+            1,-1,1],
+    # A165574
+    263 => [0,1,1,1,1,-1,1,-1,1,1,-1,1,1,1,-1,-1,1,1,1,-1,-1,
+            -1,1,1,1,1,1,1,-1,-1,-1,1,1,1,1,1,1,1,-1,1,-1,-1,
+            -1,1,1,-1,1,-1,1,1,1,1,1,-1,1,-1,-1,-1,-1,-1,-1,1,
+            1,-1,1,-1,1,-1,1,1,1,-1,1,-1,1,1,-1,-1,1,-1,-1,1,
+            -1,1],
+    # A165471
+    65537 => [0,1,1,-1,1,-1,-1,-1,1,1,-1,-1,-1,1,-1,1,1,1,1,1,
+              -1,1,-1,-1,-1,1,1,-1,-1,-1,1,-1,1,1,1,1,1,1,1,-1,
+              -1,-1,1,-1,-1,-1,-1,-1,-1,1,1,-1,1,1,-1,1,-1,-1,
+              -1,-1,1,-1,-1,-1,1,-1,1,-1,1,1,1,1,1,-1,1,-1,1,1,
+              -1,1,-1,1],
+    # A165476
+    131071 => [0,1,1,-1,1,1,-1,1,1,1,1,1,-1,-1,1,-1,1,1,1,-1,1,
+               -1,1,1,-1,1,-1,-1,1,1,-1,1,1,-1,1,1,1,-1,-1,1,1,
+               -1,-1,1,1,1,1,1,-1,1,1,-1,-1,-1,-1,1,1,1,1,1,-1,
+               -1,1,1,1,-1,-1,-1,1,-1,1,-1,1,1,-1,-1,-1,1,1,-1,1,
+               1,-1,1,-1,1],
+    # A165481
+    28657 => [0,1,1,1,1,-1,1,-1,1,1,-1,-1,1,-1,-1,-1,1,-1,1,1,
+              -1,-1,-1,-1,1,1,-1,1,-1,1,-1,-1,1,-1,-1,1,1,-1,1,
+              -1,-1,1,-1,-1,-1,-1,-1,1,1,1,1,-1,-1,1,1,1,-1,1,1,
+              -1,-1,1,-1,-1,1,1,-1,-1,-1,-1,1,-1,1,1,-1,1,1,1,
+              -1,-1,-1,1],
+     # A165586
+    514229 => [0,1,-1,-1,1,1,1,1,-1,1,-1,1,-1,1,-1,-1,1,1,-1,-1,
+               1,-1,-1,1,1,1,-1,-1,1,1,1,1,-1,-1,-1,1,1,1,1,-1,
+               -1,-1,1,1,1,1,-1,1,-1,1,-1,-1,1,-1,1,1,-1,1,-1,-1,
+               -1,1,-1,1,1,1,1,1,1,-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,
+               1,1,1,1],
+   # A165581
+    524287 => [0,1,1,-1,1,-1,-1,-1,1,1,-1,-1,-1,1,-1,1,1,-1,1,
+               -1,-1,1,-1,-1,-1,1,1,-1,-1,1,1,1,1,1,-1,1,1,1,-1,
+               -1,-1,1,1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,-1,1,-1,1,1,
+               1,1,-1,1,-1,1,-1,1,1,-1,1,1,1,1,1,1,-1,-1,1,-1,1,
+               -1,1,1,-1],
+  }
+
+  @seq.each do |p, as|
+    as.each_with_index do |l, a|
+      it "returns #{l} for #{a}.legendre(#{p})" do
+        a.legendre(p).should == l
+      end
+    end
+  end
+
+  [1, 2, 4].each do |b|
+    100.times do |a|
+      it "returns nil for #{a}.legendre(#{b})" do
+        a.legendre(b).should be_nil
+      end
+    end
+  end
+end

metadata

@@ -4,9 +4,9 @@ version: !ruby/object:Gem::Version
   prerelease: false
   segments: 
   - 0
-  - 184
+  - 185
   - 0
-  version: 0.184.0
+  version: 0.185.0
 platform: ruby
 authors: 
 - Run Paint Run Run
@@ -14,7 +14,7 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2010-10-17 00:00:00 +01:00
+date: 2010-10-18 00:00:00 +01:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -110,6 +110,7 @@ files:
 - lib/numb/keith.rb
 - lib/numb/knodel.rb
 - lib/numb/knuth.rb
+- lib/numb/kronecker.rb
 - lib/numb/kynea.rb
 - lib/numb/lah.rb
 - lib/numb/leonardo.rb
@@ -286,14 +287,17 @@ files:
 - spec/numb/interprime_spec.rb
 - spec/numb/inv_mod_spec.rb
 - spec/numb/isqrt_spec.rb
+- spec/numb/jacobi_spec.rb
 - spec/numb/jacobsthal_lucas_spec.rb
 - spec/numb/k_perfect_spec.rb
 - spec/numb/kaprekar_spec.rb
 - spec/numb/keith_spec.rb
 - spec/numb/knodel_spec.rb
 - spec/numb/knuth_spec.rb
+- spec/numb/kronecker_spec.rb
 - spec/numb/kynea_spec.rb
 - spec/numb/lah_spec.rb
+- spec/numb/legendre_spec.rb
 - spec/numb/leonardo_spec.rb
 - spec/numb/leyland_spec.rb
 - spec/numb/liouville_spec.rb
@@ -477,6 +481,7 @@ test_files:
 - spec/numb/woodall_spec.rb
 - spec/numb/unitary_sociable_spec.rb
 - spec/numb/refactorable_spec.rb
+- spec/numb/kronecker_spec.rb
 - spec/numb/maris_mcgwire_sosa_pair_spec.rb
 - spec/numb/zeisel_spec.rb
 - spec/numb/q_spec.rb
@@ -589,6 +594,7 @@ test_files:
 - spec/numb/amicable_spec.rb
 - spec/numb/bell_spec.rb
 - spec/numb/super_catalan_spec.rb
+- spec/numb/legendre_spec.rb
 - spec/numb/motzkin_spec.rb
 - spec/numb/mobius_spec.rb
 - spec/numb/chen_prime_spec.rb
@@ -628,6 +634,7 @@ test_files:
 - spec/numb/nsw_spec.rb
 - spec/numb/knuth_spec.rb
 - spec/numb/first_with_n_divisors_spec.rb
+- spec/numb/jacobi_spec.rb
 - spec/numb/isqrt_spec.rb
 - spec/numb/inrt_spec.rb
 - spec/numb/nexus_spec.rb