function sieve(limit)
    composite = Set()
    for i ∈ 2:limit
        if i ∉ composite
            for x ∈ 2i:i:limit
                push!(composite, x)
            end
        end
    end
    sort(setdiff(2:limit, composite))
end