Fortran 95如何计算大数值 解决Silverfrost 95整数溢出报错问题

报错原因说明 #

1. 整数溢出报错是因为你要计算的5^5432854数值量级远超过Fortran原生整数类型的存储上限：
   • 64位有符号整数最大值仅为9.2e18，而5^27就已经超过该阈值，5^5432854的十进制位数接近380万位，不可能用原生整数类型存储。
2. Kind参数报错是因为Silverfrost FTN95的整数Kind取值规则和GCC等其他编译器不同，它支持的整数Kind仅为1、2、3、4、7，分别对应1字节到8字节的整数类型，不存在kind=8的配置，其中kind=7对应8字节（64位）整数，但如上文所述，就算用64位整数也远远不够存储目标数值。

对应解决方案 #

分两种场景选择实现方式：

场景1：需要输出完整的精确整数结果 #

Silverfrost FTN95本身没有内置任意精度整数支持，需要手动实现大整数运算逻辑：

• 用数组存储大整数的每一位或分段数值
• 配合快速幂算法实现大数乘法和进位处理，降低运算耗时
  参考代码示例：

! 示例：数组存储逐位计算5的N次方，数组长度预留400万位足够容纳结果
integer, parameter :: MAX_BIT = 4000000
integer :: num_arr(MAX_BIT), curr_len, exponent, base, i

! 初始化数值为1（对应5^0）
num_arr(1) = 1
curr_len = 1
base = 5
exponent = 5432854

do i = 1, exponent
    call mult_by(base, num_arr, curr_len)
end do

! 倒序输出完整结果
do i = curr_len, 1, -1
    write(*, '(I1)', advance='no') num_arr(i)
end do
write(*, *)

contains
! 大整数乘单精度整数的子routine
subroutine mult_by(factor, arr, len)
    integer, intent(in) :: factor
    integer, intent(inout) :: arr(*), len
    integer :: j, carry, prod
    carry = 0
    do j = 1, len
        prod = arr(j) * factor + carry
        arr(j) = mod(prod, 10)
        carry = prod / 10
    end do
    do while(carry > 0)
        len = len + 1
        arr(len) = mod(carry, 10)
        carry = carry / 10
    end do
end subroutine mult_by

如果需要进一步提升运算速度，可以把逐位存储改成按万进制、十万进制分段存储，同时把顺序乘改成快速幂实现，可大幅降低运算次数。

场景2：仅需要近似值，不需要完整精确整数 #

直接用浮点数计算对数值即可得到科学计数法表示的近似结果，Silverfrost FTN95中用real(kind=2)声明双精度浮点数即可满足精度需求：

real(kind=2) :: log_val, coeff
integer :: exp_val
! 计算常用对数值
log_val = 5432854.0d0 * log10(5.0d0)
! 拆分系数和指数部分
exp_val = int(aint(log_val))
coeff = 10.0d0 ** (log_val - real(exp_val, kind=2))
write(*, *) "近似结果为：", coeff, "E", exp_val

计算可得结果约为3.027E3797302。

内容的提问来源于stack exchange，提问作者Mohammad Sh