A063114 a(n) = n + product of the nonzero digits of n.
2, 4, 6, 8, 10, 12, 14, 16, 18, 11, 12, 14, 16, 18, 20, 22, 24, 26, 28, 22, 23, 26, 29, 32, 35, 38, 41, 44, 47, 33, 34, 38, 42, 46, 50, 54, 58, 62, 66, 44, 45, 50, 55, 60, 65, 70, 75, 80, 85, 55, 56, 62, 68, 74, 80, 86, 92, 98, 104, 66, 67, 74, 81, 88, 95, 102, 109, 116
Offset: 1
Examples
a(59) = 59 + 5*9 = 104. a(66) = 66 + 6*6 = 102.
Links
- Harry J. Smith, Table of n, a(n) for n = 1..1000
- P. A. Loomis, An Interesting Family of Iterated Sequences
- P. A. Loomis, An Introduction to Digit Product Sequences, J. Rec. Math., 32 (2003-2004), 147-151.
- P. A. Loomis, An Introduction to Digit Product Sequences, J. Rec. Math., 32 (2003-2004), 147-151. [Annotated archived copy]
- Index entries for Colombian or self numbers and related sequences
Programs
-
ARIBAS
var stk: stack; end; for n := 1 to 80 do s := itoa(n); for j := 0 to length(s) -1 do k := atoi(s[j..j]); if k > 0 then stack_push(stk,k); end; end; write(n + product(stack2array(stk))," "); end;
-
Haskell
a063114 n = n + a051801 n -- Reinhard Zumkeller, Jan 15 2012
-
Mathematica
Table[i+Times@@(IntegerDigits[i]/. 0->1), {i, 70}] -
PARI
a(n) = n + vecprod(select(x->(x!=0), digits(n))) \\ Harry J. Smith, Aug 19 2009
Formula
a(n) = n + A051801(n). - Reinhard Zumkeller, Jan 15 2012
Extensions
More terms from Robert G. Wilson v and Klaus Brockhaus, Aug 09 2001