login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A256432
Characteristic function of octahedral numbers.
2
1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1
OFFSET
0
COMMENTS
Octahedral numbers are of the form (2*n^3 + n)/3.
FORMULA
For n > 0, a(n) = floor(t(n) - 1/(6 * t(n))) - floor(t(n-1) - 1/(6 * t(n-1))), where t(n) = ( sqrt(243*n^2+2)/(4*3^(3/2)) + 3*n/4 )^(1/3).
MATHEMATICA
Module[{upto=100, max, oct}, max=(Surd[6, 3](27upto+Sqrt[6+729upto^2])^ (2/3))/ (Surd[6, 3]^2 (27upto+Sqrt[6+729upto^2])^(1/3)); oct=Table[ (2n^3+n)/ 3, {n, 0, max}]; Table[If[MemberQ[oct, k], 1, 0], {k, 0, upto}]] (* Harvey P. Dale, Dec 01 2017 *)
CROSSREFS
Cf. A005900 (octahedral numbers).
Sequence in context: A266253 A113052 A369967 * A218171 A362130 A361018
KEYWORD
nonn
AUTHOR
Mikael Aaltonen, Mar 28 2015
STATUS
approved