P2kongr_ver8.jl

Using Primes
using DelimitedFiles
function p2kongr(knumber1::Int64, knumber2::Int64)
    s::Int64 = 0
    # 3511^2 = 12327121
    nnumbers = Int128[]
    for k::Int64 in knumber1 :  knumber2
        residue::Int64 = (5+6*k) % 8
        if (residue == 1 || residue == 7)
            n::Int128 = 21658751597 + Int128(k)*43268194710
            if powermod(2,n-2,n) == 1
                s += 1
                println(s)
                push!(nnumbers,n)
            end
        end
    end
    #return nnumbers
    #open("kongr2.dat", "w") do io writedlm(io,nnumbers) end
end