login
a(n) = wt(k) + f(j) if n = 6k+j, 0 <= j < 6, where wt() = A000120(), f() = A151899().
5

%I #5 Feb 16 2023 05:12:27

%S 0,0,1,1,1,2,1,1,2,2,2,3,1,1,2,2,2,3,2,2,3,3,3,4,1,1,2,2,2,3,2,2,3,3,

%T 3,4,2,2,3,3,3,4,3,3,4,4,4,5,1,1,2,2,2,3,2,2,3,3,3,4,2,2,3,3,3,4,3,3,

%U 4,4,4,5,2,2,3,3,3,4,3,3,4,4,4,5,3,3,4,4,4,5,4,4,5,5,5,6,1,1,2,2,2,3,2,2,3

%N a(n) = wt(k) + f(j) if n = 6k+j, 0 <= j < 6, where wt() = A000120(), f() = A151899().

%p f := proc(n) local j; j:=n mod 6; if (j<=1) then 0 elif (j<=4) then 1 else 2; fi; end;

%p wt := proc(n) local w,m,i; w := 0; m := n; while m > 0 do i := m mod 2; w := w+i; m := (m-i)/2; od; w; end;

%p A151902 := proc(n) local k,j; k:=floor(n/6); j:=n-6*k; wt(k)+f(j); end;

%t wt[n_] := DigitCount[n, 2, 1];

%t f[n_] := {0, 0, 1, 1, 1, 2}[[Mod[n, 6] + 1]];

%t a[n_] := wt[Floor[n/6]] + f[n - 6 Floor[n/6]];

%t Table[a[n], {n, 0, 104}] (* _Jean-François Alcover_, Feb 16 2023 *)

%Y Cf. A000120, A151899, A151904-A151907.

%K nonn

%O 0,6

%A _N. J. A. Sloane_, Jul 31 2009