OFFSET
1,2
COMMENTS
The sequence is equal to A053837 up to the 488th term.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
a(345)=2 because 3+4+5=12 and 1*2=2.
MAPLE
P:=proc(n) local i, k, w; for i from 1 by 1 to n do w:=0; k:=i; while k>0 do w:=w+k-(trunc(k/10)*10); k:=trunc(k/10); od; k:=w; w:=1; while k>0 do w:=w*(k-(trunc(k/10)*10)); k:=trunc(k/10); od; print(w); od; end: P(500);
# alternative
f:= n -> convert(convert(convert(convert(n, base, 10), `+`), base, 10), `*`):
map(f, [$1..100]); # Robert Israel, Dec 09 2016
MATHEMATICA
sdpd[n_]:=Module[{s=Total[IntegerDigits[n]]}, Times@@IntegerDigits[s]]; Array[sdpd, 110] (* Harvey P. Dale, Dec 17 2013 *)
CROSSREFS
KEYWORD
easy,nonn,base
AUTHOR
Paolo P. Lava and Giorgio Balzarotti, May 03 2007
EXTENSIONS
Offset corrected by Robert Israel, Dec 09 2016
STATUS
approved