cp's OEIS Frontend

A133457 Irregular triangle read by rows: row n gives exponents in expression for n as a sum of powers of 2.

%I A133457 #39 Dec 22 2021 02:20:12
%S A133457 0,1,0,1,2,0,2,1,2,0,1,2,3,0,3,1,3,0,1,3,2,3,0,2,3,1,2,3,0,1,2,3,4,0,
%T A133457 4,1,4,0,1,4,2,4,0,2,4,1,2,4,0,1,2,4,3,4,0,3,4,1,3,4,0,1,3,4,2,3,4,0,
%U A133457 2,3,4,1,2,3,4,0,1,2,3,4,5,0,5,1,5,0,1,5,2,5,0,2,5,1,2,5,0,1,2,5,3,5,0,3,5
%N A133457 Irregular triangle read by rows: row n gives exponents in expression for n as a sum of powers of 2.
%C A133457 This sequence contains every increasing finite sequence. For example, the finite sequence {0,2,3,5} arises from n = 45.
%C A133457 Essentially A030308(n,k)*k, then entries removed where A030308(n,k)=0. - _R. J. Mathar_, Nov 30 2007
%C A133457 In the corresponding irregular triangle {a(n)+1}, the m-th row gives all positive integer roots m_i of polynomial {m,k}. - see link [Shevelev]; see also A264613. - _Vladimir Shevelev_, Dec 13 2015
%H A133457 Reinhard Zumkeller, <a href="/A133457/b133457.txt">Rows n = 1..1024 of triangle, flattened</a>
%H A133457 Vladimir Shevelev, <a href="http://www.emis.de/journals/INTEGERS/papers/m1/m1.Abstract.html">The number of permutations with prescribed up-down structure as a function of two variables</a>, INTEGERS, 12 (2012), #A1. (See Section 3, Theorem 21 and Section 8, Theorem 50)
%F A133457 a(n) = A048793(n) - 1.
%e A133457 1 = 2^0.
%e A133457 2 = 2^1.
%e A133457 3 = 2^0 + 2^1.
%e A133457 4 = 2^2.
%e A133457 5 = 2^0 + 2^2.
%e A133457 etc. and reading the exponents gives the rows of the triangle.
%p A133457 A133457 := proc(n) local a,bdigs,i ; a := [] ; bdigs := convert(n,base,2) ; for i from 1 to nops(bdigs) do if op(i,bdigs) <> 0 then a := [op(a),i-1] ; fi ; od: a ; end: seq(op(A133457(n)),n=1..80) ; # _R. J. Mathar_, Nov 30 2007
%t A133457 Array[Join @@ Position[#, 1] - 1 &@ Reverse@ IntegerDigits[#, 2] &, 41] // Flatten (* _Michael De Vlieger_, Oct 08 2017 *)
%o A133457 (Haskell)
%o A133457 a133457 n k = a133457_tabf !! (n-1) !! n
%o A133457 a133457_row n = a133457_tabf !! (n-1)
%o A133457 a133457_tabf = map (fst . unzip . filter ((> 0) . snd) . zip [0..]) $
%o A133457                    tail a030308_tabf
%o A133457 -- _Reinhard Zumkeller_, Oct 28 2013, Feb 06 2013
%Y A133457 Cf. A003071, A030308, A048793, A067255, A264613.
%Y A133457 Cf. A073642 (row sums), A272011 (rows reversed).
%K A133457 base,tabf,easy,nonn,look
%O A133457 1,5
%A A133457 _Masahiko Shin_, Nov 27 2007
%E A133457 More terms from _R. J. Mathar_, Nov 30 2007