A134586 a(n) = a(n-1) + A007954(n-1) if a(n-1) is odd, otherwise a(n) = a(n-1) / 2.

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, 93, 93, 96, 48, 24, 12, 6, 3, 24, 12, 6, 3, 7, 15, 27, 43, 63, 87, 115, 147, 183, 183, 188, 94, 47, 67, 92, 46, 23, 63, 108, 54, 27, 39, 57, 81, 111, 147, 189, 237, 291, 291, 298, 149

Offset

1,2

Links

Harvey P. Dale, Table of n, a(n) for n = 1..1000

Mathematica

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 *)

Prog

(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

Crossrefs

Cf. A007954.

Sequence in context: A210034 A074586 A277812 * A135287 A359803 A290650

Adjacent sequences: A134583 A134584 A134585 * A134587 A134588 A134589

Keyword

base,nonn

Author

Ctibor O. Zizka, Jan 24 2008

Extensions

Name corrected by Michel Marcus, Nov 28 2020

Status

approved