使用不同整数kind值时Fortran程序出现异常行为

问题根源 #

你的代码存在两个核心问题：

1. 累加变量s未初始化：Fortran中局部变量默认不会自动初始化为0，程序启动时s会取内存中的随机垃圾值，后续累加操作会基于这个错误值计算，导致结果完全偏离预期。这和整数kind无关——哪怕用kind=8，只是偶然随机值刚好让结果看起来正常，本质逻辑依然错误。
2. kind值写法依赖编译器：直接用kind=8/16/32是编译器相关的非标准写法，不同编译器对这些数字的定义可能不同（比如有些编译器中kind=8对应64位，但kind=16不一定支持128位）。正确的做法是用selected_int_kind()函数获取对应精度的整数kind。

修正后的代码 #

program multiples_of_3_and_5
    implicit none
    
    ! 使用selected_int_kind获取指定精度的整数类型，38对应128位整数（可表示10^38以内的数）
    integer, parameter :: int128 = selected_int_kind(38)
    
    integer(int128) :: n
    ! 初始化累加变量s为0
    integer(int128) :: s = 0
    
    do n = 1, 999
        if (mod(n, 3) == 0 .or. mod(n, 5) == 0) then
            s = s + n
        end if
    end do
    
    print *,"The result is: ", s
    
end program multiples_of_3_and_5

关键修改说明 #

• 用selected_int_kind(38)定义128位整数的kind参数int128，保证代码在不同编译器下的兼容性。
• 显式初始化s = 0，避免垃圾值干扰累加计算。
• 简化循环条件为1, 999，和1000-1等价但更清晰。

正确结果 #

运行修正后的代码，输出应为：The result is: 233168，这和数学计算的结果一致（1到999中能被3或5整除的数之和：(3×333×334/2)+(5×199×200/2)-(15×66×67/2) = 166833 + 99500 - 33165 = 233168）。

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