Yet another belphegor's prime script

belphegor.jl

#!/usr/bin/env julia

using DataStructures
using Primes
using Statistics

belphegor(x::Integer) = parse(BigInt,(z -> "1$(z)666$(z)1")("0"^x))

buf = CircularBuffer{Float64}(100)

for i in 0:10^4
    n = belphegor(i)
    t = @elapsed test = isprime(n)
    push!(buf, t)
    if test
        print(" "^32*"\r")
        println(i)
    else
        print("$i (10^$(round(log10(mean(buf)); digits = 2))s)"*" "^8*"\r")
    end
end

exit(0)