Evenpseudoprimes.jl

#Even pseudoprimes, integers n such that 2^n = 2 ( mod n)
using Primes
using DelimitedFiles
function evenpseudo(nmin::Int64,nmax::Int64)
    if nmin % 2 != 0
        nmin = nmin +1  # nmin must be even
    s:Int64 = 0
    evenpseudos = Int64[]
    for n::Int64 in nmin : 2 : nmax
        if powermod(2,n,n) == 2
            s += 1
            println(s," ",n," ",factor(n))
            push!(evenpseudos,n)
        end
    end
    #open("evenpseudos.dat", "w") do io writedlm(io,evenpseudos) end
    return evenpseudos
end