|
%I #15 Nov 28 2020 06:56:29
%S 1,2,1,4,2,1,7,14,7,16,8,4,2,1,5,10,5,12,6,3,3,5,9,15,23,33,45,59,75,
%T 93,93,96,48,24,12,6,3,24,12,6,3,7,15,27,43,63,87,115,147,183,183,188,
%U 94,47,67,92,46,23,63,108,54,27,39,57,81,111,147,189,237,291,291,298,149
%N a(n) = a(n-1) + A007954(n-1) if a(n-1) is odd, otherwise a(n) = a(n-1) / 2.
%H Harvey P. Dale, <a href="/A134586/b134586.txt">Table of n, a(n) for n = 1..1000</a>
%t nxt[{n_,a_}]:={n+1,If[EvenQ[a],a/2,a+Times@@IntegerDigits[n+1]]}; Transpose[ NestList[nxt,{0,1},80]][[2]] (* _Harvey P. Dale_, Sep 11 2015 *)
%o (PARI) lista(nn) = {my(va = vector(nn)); va[1] = 1; for (n=2, nn, va[n] = if (va[n-1] % 2, va[n-1] + vecprod(digits(n-1)), va[n-1]/2);); va;} \\ _Michel Marcus_, Nov 28 2020
%Y Cf. A007954.
%K base,nonn
%O 1,2
%A _Ctibor O. Zizka_, Jan 24 2008
%E Name corrected by _Michel Marcus_, Nov 28 2020
|
