This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
a(3)=286 because 259^2+286^2=148877=53^3.
a(3)=234 because 234^2+415^2=226981=61^3.
Table[Select[PowersRepresentations[61^n, 2, 2], CoprimeQ[#[[1]], #[[2]]] &][[1, 1]], {n, 22}] (* T. D. Noe, Nov 04 2013 *)
Table[Select[PowersRepresentations[61^n, 2, 2], CoprimeQ@@#&][[1, 1]], {n, 1, 40}] (* Vincenzo Librandi, Mar 02 2014 *)
a(3) = 296 because 296^2 + 549^2 = 389017 = 73^3.
f:=n -> min([abs@Re,abs@Im]((3+8*I)^n)): map(f, [$1..50]); # Robert Israel, Mar 31 2017
Table[Select[PowersRepresentations[73^n, 2, 2], CoprimeQ@@#&][[1, 1]], {n, 1, 40}] (* Vincenzo Librandi, Mar 02 2014 *)
Table[Min[Abs[Re[(3 + 8I)^n]], Abs[Im[(3 + 8I)^n]]], {n, 30}] (* Indranil Ghosh, Mar 31 2017, after formula by Robert Israel *)
from sympy import I, re, im print([min(abs(re((3 + 8*I)**n)), abs(im((3 + 8*I)**n))) for n in range(1, 31)]) # Indranil Ghosh, Mar 31 2017, after formula by Robert Israel
Comments