%I #11 Dec 01 2017 10:33:19
%S 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,
%T 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,
%U 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1
%N Characteristic function of octahedral numbers.
%C Octahedral numbers are of the form (2*n^3 + n)/3.
%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>
%F 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).
%t 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 *)
%Y Cf. A005900 (octahedral numbers).
%K nonn
%O 0
%A _Mikael Aaltonen_, Mar 28 2015