cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A089398 a(n) = n-th column sum of binary digits of k*2^(k-1), where summation is over k>=1, without carrying between columns.

This page as a plain text file.
%I A089398 #16 Mar 29 2015 13:56:43
%S A089398 1,0,2,1,1,1,3,2,2,0,3,2,2,2,4,3,3,1,2,2,2,2,4,3,3,1,4,3,3,3,5,4,4,2,
%T A089398 3,1,2,2,4,3,3,1,4,3,3,3,5,4,4,2,3,3,3,3,5,4,4,2,5,4,4,4,6,5,5,3,4,2,
%U A089398 1,2,4,3,3,1,4,3,3,3,5,4,4,2,3,3,3,3,5,4,4,2,5,4,4,4,6,5,5,3,4,2,3,3,5,4,4
%N A089398 a(n) = n-th column sum of binary digits of k*2^(k-1), where summation is over k>=1, without carrying between columns.
%C A089398 sum(k=1,n, a(k)*2^(k-1)) = 2^A089399(n)+1 for n>2, with a(1)=a(2)=1.
%C A089398 Row sums of triangular arrays in A103588 and in A103589. - _Philippe Deléham_, Apr 04 2005
%C A089398 a(k) = 0 for k = 2, 10, 2058, 2058 + 2^2059, ..., that is, for k = A034797(n) - 1, n>=2. - _Philippe Deléham_, Nov 16 2007
%F A089398 a(2^n)=n-1 (for n>0), a(2^n-1)=n (for n>0), a(2^n+1)=n-1 (for n>1), a(2^n-k)=n-A089400(k) (for n>k>0), a(2^n+k)=n-A089401(k) (for n>k>0), where sequences have limits: A089400={0, 2, 2, 2, 1, 4, 2, 2, 1, 3, 3, ...} and A089401={1, 1, 3, 2, 4, 5, 6, 5, 7, 8, 11, 9, ...},
%e A089398 Binary expansions of k*2^(k-1), with bits in ascending order by powers of 2, are:
%e A089398 1
%e A089398 001
%e A089398 0011
%e A089398 000001
%e A089398 0000101
%e A089398 00000011
%e A089398 000000111
%e A089398 00000000001
%e A089398 000000001001
%e A089398 0000000000101
%e A089398 00000000001101
%e A089398 000000000000011
%e A089398 0000000000001011
%e A089398 .................
%e A089398 Giving column sums:
%e A089398 10211132203222433...
%t A089398 f[n_] := Block[{lg = Floor[Log[2, n]] + 1}, Sum[ Join[ Reverse[ IntegerDigits[n - i + 1, 2]], {0}][[i]], {i, lg}]]; Table[ f[n], {n, 105}] (* _Robert G. Wilson v_, Mar 26 2005 *)
%o A089398  (PARI) /* Prints initial 1000 terms: */
%o A089398 {A=vector(1000);for(n=1,#A,Bn=binary(n*2^(n-1));for(k=1,min(#Bn,#A),A[k]=A[k]+Bn[#Bn-k+1]) );print(A)}
%Y A089398 Cf. A089399, A089400, A089401.
%K A089398 base,nonn
%O A089398 1,3
%A A089398 _Paul D. Hanna_, Oct 30 2003