A033922 Base-2 digital convolution sequence.
1, 1, 1, 2, 1, 2, 2, 3, 2, 3, 3, 4, 3, 4, 4, 5, 1, 2, 2, 3, 2, 3, 3, 4, 3, 4, 4, 5, 4, 5, 5, 6, 2, 3, 3, 4, 3, 4, 4, 5, 4, 5, 5, 6, 5, 6, 6, 7, 3, 4, 4, 5, 4, 5, 5, 6, 5, 6, 6, 7, 6, 7, 7, 8, 2, 3, 3, 4, 3, 4, 4, 5, 4, 5, 5, 6, 5, 6, 6, 7, 3, 4, 4, 5, 4, 5
Offset: 0
Examples
For example, 6 = 2^2 + 2^1, so a(6) = a(2) + a(1) = 2.
Links
- Alois P. Heinz, Table of n, a(n) for n = 0..10000
- Benoit Cloitre, Plot of a(n) for n=1 up to 100000
Crossrefs
Programs
-
Maple
a:= proc(n) option remember; local c, m, t; if n=0 then 1 else m:= n; c:=0; for t from 0 while m<>0 do c:= c+ `if`(irem(m, 2, 'm')=1, a(t), 0) od; c fi end: seq(a(n), n=0..120); # Alois P. Heinz, Jul 13 2011
-
PARI
al(n)=local(v,k,e);v=vector(n+1);v[1]=1;for(m=1,n,k=m;e=0;while(k>0,if(k%2,v[m+1]+=v[e+1]);e++;k\=2));v /* Benoit Cloitre, Jan 10 2013 */
-
PARI
/* to compute quickly 2^m terms of the sequence */ m=10;v=[0,1];for(n=2,m,v=concat(v,vector(2^n/2,i,v[i]+v[n])));a(n)=if(n<2,1,v[n]) /* Benoit Cloitre, Jan 16 2013 */
Formula
From Flávio V. Fernandes, Jul 31 2025: (Start)
a(n) = a(2^n).
Extensions
Edited by Franklin T. Adams-Watters, Jul 13 2011
Comments