import Pkg; Pkg.add("Primes")

using Primes  # 用于处理素数和因式分解的软件包

# 该函数检查数字n是否几乎是阶数k的素数
# 第k阶的几乎素数是具有恰好k个素数因子的数（考虑到多重性）

isalmostprime(n::Integer, k::Integer) = sum(values(factor(n))) == k

# 该函数返回前N个几乎k素数
function almostprimes(N::Integer, k::Integer) 
    # 创建用于存储类型为k的大小为N的结果的向量
    P = Vector{typeof(k)}(undef,N)
    # 初始化已找到数字的计数器和当前已检查数字
    i = 0; n = 2
    # 循环继续，直到我们找到N个数字。
    while i < N
        # 如果当前数在k的量级上几乎是素数
        if isalmostprime(n, k) 
            # 我们将其保存到results数组并递增计数器。
            P[i += 1] = n 
        end
        # 继续下一个数字
        n += 1
    end
    # 我们返回找到的几乎素数的向量
    return P
end

# 对于k从1到5的每个值，我们生成并输出前10个几乎素数
for k in 1:5
    # 我们生成阶数k的前10个几乎素数
    primes_k = almostprimes(10, k)
    # 将数字转换为字符串并用逗号连接
    result_str = join(primes_k, ", ")
    # 我们以格式化的形式输出结果
    println("$k-几乎素数：resul result_str。..")
end