整数精度[重复]

Question

我正在使用 Fortran 解决 Project Euler 问题；问题是创建一个斐波那契数列并找到所有小于 400 万的偶数的总和。这是我写的

implicit none
integer*4::a(1:4000000),sum
integer*4::i,maxc
maxc = 3999999

a(1) = 1
a(2) = 2
do i = 3,maxc,1
   a(i) = a(i-1) + a(i-2)
end do
sum = 0
do i = 1,maxc
   if (mod(a(i),2)==0) then
      sum = sum + a(i)
   end if
end do
print*,sum


end

输出为-1833689714 知道出了什么问题吗？

Answer 1

由于您选择的整数类型的大小，您可以表示的数字是有限的。 在您的代码中，它是2147483647（带有gfortran，由print *,huge(sum) 获得）。 可以证明，在您的实现中，i=59 超出了此限制。 然后，您会得到一个整数溢出，并且该值变为负数。

简单地使用浮点表示求和，即

real :: sum

成功了。