The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A023533 a(n) = 1 if n is of the form m(m+1)(m+2)/6, and 0 otherwise. 45

%I

%S 1,1,0,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,

%T 0,1,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,

%U 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

%N a(n) = 1 if n is of the form m(m+1)(m+2)/6, and 0 otherwise.

%C a(n) is the characteristic function of tetrahedral numbers. - _Mikael Aaltonen_, Mar 28 2015

%H Danny Rorabaugh, <a href="/A023533/b023533.txt">Table of n, a(n) for n = 0..10000</a>

%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>

%F a(A000292(n))=1; a(A145397(n))=0; a(n)=1-A014306(n). - _Reinhard Zumkeller_, Oct 14 2008

%F For n > 0, a(n) = floor(t(n) + 1/(3 * t(n)) - 1) - floor(t(n-1) + 1/(3 * t(n-1)) - 1), where t(n) = ( sqrt(243*n^2-1)/3^(3/2) + 3*n )^(1/3). - _Mikael Aaltonen_, Mar 28 2015

%t With[{ms=Table[m(m+1)(m+2)/6,{m,0,20}]},Table[If[MemberQ[ms,n],1,0], {n,0,100}]] (* _Harvey P. Dale_, Jul 25 2011 *)

%o (Sage) # Generates an array with at least N terms.

%o def A023533_list(N):

%o A = []

%o for m in range(ceil((6*N)^(1/3))):

%o A.extend(*(binomial(m+2, 3) - len(A)) + )

%o return A

%o print(A023533_list(40))

%o # _Danny Rorabaugh_, Mar 16 2015

%o (PARI) lista(nn) = {v = vector(nn); for (n=0, nn, i = 1+n*(n+1)*(n+2)/6; if (i > nn, break); v[i] = 1;); v;} \\ _Michel Marcus_, Mar 16 2015

%Y Cf. A000292, A014306, A145397.

%K nonn,easy

%O 0,1

%A _Clark Kimberling_, Jun 14 1998

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified April 8 14:52 EDT 2020. Contains 333314 sequences. (Running on oeis4.)