CollectionProduct.jl

# 2,3,4,5 -alkulukujokkojen tulot alkuluvuista testi
nnumbers = Int64[]  # sisältää kongruenssin toteuttamat joukot
for i in 2:5
    joukko = collect(combinations(testi,i))
    for j in 1 : length(joukko)
        tulo = prod(joukko[j])
        if powermod(2,tulo-2,tulo) == 1
            push!(nnumbers,tulo)
            println(tulo,"   ",factor(tulo))
        end
    end
end