julia - 有没有办法阻止 Julia 使用 BigFloats 的科学(指数)表示法？

Question

好吧，从表面上看，我有一个简单的问题。我想将无理数的小数部分取到指定的位数并将其视为整数。例如，如果我的无理数是 2.657829……而我想要五位数字，那么我正在寻找 65782(尽管我实际上是在处理“大”数字)。

这可以使用字符串轻松实现，例如，如果我想要根 3 到 50 位数字的小数部分:

function main_1(n::Int, m::Int)::BigInt
    setprecision(Int(trunc((m + 3) * log2(10))))
    sr = string(sqrt(big(n)))
    ff = findfirst(sr, '.')
    dp = parse(BigInt, sr[ff + 1:ff + m])
    return dp
end

@time main_1(3, 50)

输出为 73205080756887729352744634150587236694280525381038 .

但是，当我只处理数字时，我讨厌使用字符串!我想要做的是从 BigFloat 开始，减去整数部分，将结果乘以适当的 10 倍，将结果四舍五入为零，然后将其转换为 BigInt。问题是 Julia 使用科学/指数表示法，所以我似乎无法仅使用数字来实现我想要的。以下(部分)代码显示了问题:

function main_2(n::Int, m::Int)::BigFloat
    setprecision(Int(trunc((m + 3) * log2(10))))
    sr = sqrt(big(n))
    tr = trunc(sr)
    dp = (sr - tr) * big(10.0) ^ 50
    return dp
end

@time main_2(3, 50)

这种情况下的输出是 7.32050807568877293527446341505872366942805253810380625e+49 (在舍入阶段会删除一些额外的数字)。

所以我的问题是，有没有办法在不诉诸字符串的情况下实现我的目标？

Answer 1

实现此目的的一种方法是在不使用字符串的情况下将结果及其整数部分转换为 BigInt，然后再进行减法(并将函数类型从 BigFloat 更改为 BigInt):

function main_2(n::Int, m::Int)::BigInt
    setprecision(Int(trunc((m + 3) * log2(10))))

    # Calc the sqrt
    result = sqrt(big(n))

    # Convert the whole number to BigInt to the specified precision
    sr = convert(BigInt, trunc(result*big(10)^m))

    # Convert the integer part to BigInt
    tr = convert(BigInt, trunc(result)*big(10)^m)

    dp = sr - tr
    return dp
end

将上述实现与 main_1 进行比较功能，有一点改进:

julia> @time main_1(3, 50)
  0.000042 seconds (36 allocations: 5.254 KiB)
73205080756887729352744634150587236694280525381038

julia> @time main_2(3, 50)
  0.000028 seconds (51 allocations: 1.617 KiB)
73205080756887729352744634150587236694280525381038

编辑:

其他方式(如@Bill 所评论)只是截断结果(以摆脱 InexactError() )并将函数类型更改为 BigInt:

function main_2(n::Int, m::Int)::BigInt
    setprecision(Int(trunc((m + 3) * log2(10))))
    sr = sqrt(big(n))
    tr = trunc(sr)
    dp = (sr - tr) * big(10.0) ^ 50
    return trunc(dp)
end

测试后:

julia> @time main_2(3,50)
  0.000026 seconds (28 allocations: 1.016 KiB)
73205080756887729352744634150587236694280525381038