From d9e13dbc20cd9df48d849fac8e36923f8af5c5b1 Mon Sep 17 00:00:00 2001
From: Rafael Fourquet <fourquet.rafael@gmail.com>
Date: Sat, 24 Jun 2017 16:53:31 +0200
Subject: [PATCH] fix prodfactor(factor(1)) (was an error)

---
 src/Primes.jl    | 2 +-
 test/runtests.jl | 3 ++-
 2 files changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/Primes.jl b/src/Primes.jl
index 8676625..47e6451 100644
--- a/src/Primes.jl
+++ b/src/Primes.jl
@@ -362,7 +362,7 @@ julia> prodfactors(factor(100))
 """
 function prodfactors end
 
-prodfactors(factors::Associative) = prod(p^i for (p, i) in factors)
+prodfactors{K}(factors::Associative{K}) = isempty(factors) ? one(K) : prod(p^i for (p, i) in factors)
 prodfactors(factors::Union{AbstractArray, Set, IntSet}) = prod(factors)
 
 """
diff --git a/test/runtests.jl b/test/runtests.jl
index e73287c..1c7e4ec 100644
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -341,7 +341,8 @@ end
 end
 
 for T in (Int, UInt, BigInt)
-    for n in rand(T(1):T(100000), 10)
+    for n in [T(1); rand(T(2):T(100000), 10)]
+        # for n=T(1), must not error out (#51)
         for C = (Factorization, Vector, Dict)
             @test prodfactors(factor(C, n)) == n
         end