login
A136346
Octagonal numbers which are the sums of exactly two positive octagonal numbers.
2
560, 736, 1541, 3201, 5461, 6816, 7400, 9976, 11041, 11408, 13333, 14981, 15408, 15841, 19521, 21000, 21505, 25761, 28616, 30401, 41536, 45141, 50440, 51221, 52008, 54405, 56856, 61920, 63656, 65416, 69008, 75525, 76480, 81345, 82336, 85345, 87381, 89441
OFFSET
1,1
COMMENTS
For sums of two positive octagonal numbers, see A136345. This is to octagonal numbers A000567 as A089982 is to triangular numbers A000217, as A009000 is to squares A000290, as A136117 is to pentagonal numbers A000326, as A133215 is to hexagonal numbers A000384, and as A117104 is to heptagonal numbers A000566. If Oc(a) + Oc(b) = Oc(c) then a(3a-2) + b(3b+2) = c(3c+2), so solving the quadratic equations for c we have (when an integer): c = (2 + sqrt(4 + 36a^2 + 36b^2 - 24a - 24b))/6.
LINKS
Eric Weisstein's World of Mathematics, Octagonal Number
FORMULA
A000567 INTERSECTION {A000567(i) + A000567(j), i, j > 0}. {i*(3*i-2)} INTERSECTION {i*(3*i-2) + j(3*j-2), i > 0}.
EXAMPLE
Where Oc(n) = A000567(n) = n-th octagonal number:
a(1) = 560 = Oc(14) = 280 + 280 = Oc(10) + Oc(10).
a(2) = 736 = Oc(16) = 560 + 176 = Oc(14) + Oc(8).
a(3) = 1541 = Oc(23) = 1408 + 133 = Oc(22) + Oc(7).
a(4) = 3201 = Oc(33) = 2465 + 736 = Oc(29) + Oc(16).
a(5) = 5461 = Oc(43) = 2821 + 2640 = Oc(31) + Oc(30).
MATHEMATICA
Module[{nn=300, ono}, ono=PolygonalNumber[8, Range[nn]]; Union[Select[ Total/@ Tuples[ono, 2], MemberQ[ono, #]&]]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Oct 26 2019 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jonathan Vos Post, Dec 25 2007
EXTENSIONS
Corrected and edited by B. D. Swan (bdswan(AT)gmail.com), Dec 20 2008
STATUS
approved