reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
1, 5, 25, 125, 65, 3125, 15625, 325, 390625, 1953125, 1625, 48828125, 4225, 1105, 6103515625, 30517578125, 40625, 21125, 3814697265625, 203125, 95367431640625, 476837158203125, 5525, 11920928955078125, 274625, 5078125, 1490116119384765625, 528125, 25390625, 186264514923095703125
approved
editing
proposed
approved
editing
proposed
Ray Chandler, <a href="/A006339/b006339.txt">Table of n, a(n) for n = 0..1438</a> (a(1439) exceeds 1000 digits).
approved
editing
_David W. Wilson (davidwwilson(AT)comcast.net)_
editing
approved
primeCon1Mod4 = Select[ Prime[ Range[ 50 ] ], (Mod[ #1, 4 ] == 1) & ]; (Times @@ (Take[ primeCon1Mod4, Length[ #1 ] ]^Reverse[ (#1 - 1)/2 ])) & /@ ((Join @@ (Table[ #1[ [ 1 ] ], {#1[ [ 2 ] ]} ] & /@ FactorInteger[ #1*2 + 1 ])) & /@ Range[ 0, 50 ]) (* Albert Mao, Dec 03 2006 *)
oneModFourPrimes[1] = 5;
oneModFourPrimes[n_] := oneModFourPrimes[n] = NestWhile[NextPrime, NextPrime[oneModFourPrimes[n - 1]], Mod[#, 4] != 1 & ];
factorizations[1, limit_] = {{}};
factorizations[n_, limit_] := factorizations[n, limit] = Join @@ Table[Prepend[#, d]& /@ factorizations[n/d, d], {d, Select[Rest[Divisors[n]], # <= limit & ]}];
leastHypotenuse[n_] := Min[(Times @@ (Array[oneModFourPrimes, Length[#]]^((# - 1)/2)) & ) /@ factorizations[2*n + 1, 2*n + 1]];
Array[leastHypotenuse, 30, 0]
(* Albert H. Mao, Jan 06 2012 *)
approved
editing