多倍長計算が面倒なので、既存のライブラリ FMLIB を利用することにします。MPFUN とどっちがいいのかよくわかりません。 数字列はファイルから読み取った方が楽だと思いましたが、まぁ整形してソース・プログラム中にべたうちすることにしました。 ソース・…

i番目の三角数は i(i+1)/2 素因数分解は i, i+1, 2 についてやって足し引きする。約数の数は素因数に1足して積。 ソース・プログラム program PEuler012 implicit none integer :: i, k1, k2, n integer, allocatable :: ipfac(:), ipfac0(:), ipfac1(:) i =…

めんどくさいｗ 美しくない。 ソース・プログラム この入力だと配列の列と行が反対になってる。 program PEuler011 implicit none integer, parameter :: n = 20 integer, parameter :: itab(n, n) = [ & ! transpose [08,02,22,97,38,15,00,40,00,75,04,05,…