OFFSET
1,1
COMMENTS
Proper subset of A024406.
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..400
John D. Cook, Sparse binary Pythagorean triples (2025).
H. S. Uhler, A Colossal Primitive Pythagorean Triangle, The American Mathematical Monthly, Vol. 57, No. 5 (May, 1950), pp. 331-332.
Wikipedia, Pythagorean triple.
Index entries for linear recurrences with constant coefficients, signature (340,-22848,348160,-1048576).
FORMULA
a(n) = (2^(4n) + 2^(2n+1)) * (2^(4n) - 2^(4n-2) - 2^(2n) - 1) / 2.
G.f.: 12*(7 - 280*x - 24832*x^2 + 163840*x^3)/((1 - 4*x)*(1 - 16*x)*(1 - 64*x)*(1 - 256*x)). - Stefano Spezia, Feb 13 2025
MATHEMATICA
LinearRecurrence[{340, -22848, 348160, -1048576}, {84, 25200, 6350784, 1614708480}, 20] (* Paolo Xausa, Feb 26 2025 *)
PROG
(PARI) a(n) = (2^(4*n) + 2^(2*n+1)) * (2^(4*n) - 2^(4*n-2) - 2^(2*n) - 1) / 2
(Magma) [(2^(4*n) + 2^(2*n+1)) * (2^(4*n) - 2^(4*n-2) - 2^(2*n) - 1) / 2: n in [1..20]];
(Python)
def A381009(n): return (m:=1<<(n<<1)-1)*(m-1)*(m+1)*(3*m+1)<<1 # Chai Wah Wu, Feb 13 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Robert C. Lyons, Feb 12 2025
STATUS
approved
