2^1000を求めるのに手抜きして多倍長整数ライブラリ FMLIB を利用させていただきました。2^1000〜10^300なので300桁強の文字列を用意しました。 ソース・プログラム program PEuler016 use fmzm implicit none character(len = 350) :: str integer :: k(350…

をパスカルの三角形で求めます。４バイト整数ではオーバーフローします。 ソース・プログラム program PEuler015 implicit none integer, parameter :: ki = selected_int_kind(15), n = 20 integer :: i integer(ki), allocatable :: ipascal(:) ipascal = …

工夫せず素朴に求めます。4バイト整数ではオーバーフローしました。 ソース・プログラム program PEuler014 implicit none integer :: i, k, kk(1e6) integer(8) :: m do i = 1, 1e6 k = 1 m = i do if (mod(m, 2) == 0) then m = m / 2 else m = 3 * m + 1 …