login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A073267 Number of compositions (ordered partitions) of n into exactly two powers of 2. 18

%I #29 Jul 28 2016 21:54:12

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

%T 2,0,2,0,0,0,2,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,2,2,0,

%U 2,0,0,0,2,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0

%N Number of compositions (ordered partitions) of n into exactly two powers of 2.

%C Starting with 1 = self-convolution of A036987, the characteristic function of the powers of 2. [_Gary W. Adamson_, Feb 23 2010]

%H Reinhard Zumkeller, <a href="/A073267/b073267.txt">Table of n, a(n) for n = 0..10000</a>

%H Sen-Peng Eu, Shu-Chung Liu and Yeong-Nan Yeh, <a href="http://dx.doi.org/10.1016/j.ejc.2007.06.019">Catalan and Motzkin numbers modulo 4 and 8</a>, Eur. J. Combinat. 29 (2008) 1449-1466.

%F G.f.: (Sum_{k>=0} x^(2^k) )^2. - _Vladeta Jovovic_, Mar 28 2005

%F a(n+1) = A000108(n) mod 4, n>=1 [Theorem 2.3 of Eu et al.]. - _R. J. Mathar_, Feb 27 2008

%F a(n) = sum (A209229(k)*A036987(n-k): k = 0..n), convolution of characteristic functions of 2^n and 2^n-1. [_Reinhard Zumkeller_, Mar 07 2012]

%F a(n+2) = A000168(n) mod 4. - _John M. Campbell_, Jul 07 2016

%e For 2 there is only composition {1+1}, for 3 there is {1+2, 2+1}, for 4 {2+2}, for 5 {1+4, 4+1}, for 6 {2+4,4+2}, for 7 none, thus a(2)=1, a(3)=2, a(4)=1, a(5)=2, a(6)=2 and a(7)=0.

%p f:= proc(n) local d;

%p d:= convert(convert(n,base,2),`+`);

%p if d=2 then 2 elif d=1 then 1 else 0 fi

%p end proc:

%p 0, 0, seq(f(n),n=2..100); # _Robert Israel_, Jul 07 2016

%t Table[Count[Map[{#, n - #} &, Range[0, n]], k_ /; Times @@ Boole@ Map[IntegerQ@ Log2@ # &, k] == 1], {n, 0, 88}] (* _Michael De Vlieger_, Jul 08 2016 *)

%o (Haskell)

%o a073267 n = sum $ zipWith (*) a209229_list $ reverse $ take n a036987_list

%o -- _Reinhard Zumkeller_, Mar 07 2012

%o (PARI)

%o N=166; x='x+O('x^N);

%o v=Vec( 'a0 + sum(k=0,ceil(log(N)/log(2)), x^(2^k) )^2 );

%o v[1] -= 'a0; v

%o /* _Joerg Arndt_, Oct 21 2012 */

%Y The second row of the table A073265. The essentially same sequence 1, 1, 2, 1, 2, 2, 0, 1, ... occurs for first time in A073202 as row 105 (the fix count sequence of A073290). The positions of 1's for n > 1 is given by the characteristic function of A000079, i.e. A036987 with offset 1 instead of 0 and the positions of 2's is given by A018900. Cf. also A023359.

%Y Cf. A036987. [_Gary W. Adamson_, Feb 23 2010]

%K nonn

%O 0,4

%A _Antti Karttunen_, Jun 25 2002

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 24 18:17 EDT 2024. Contains 371962 sequences. (Running on oeis4.)