|
%I #16 Oct 01 2021 10:28:15
%S 10,1,11,2,13,4,17,8,25,7,32,5,37,10,47,11,58,13,71,8,79,16,95,14,109,
%T 10,119,11,130,4,134,8,142,7,149,14,163,10,173,11,184,13,197,17,214,7,
%U 221,5,226,10,236,11,247,13,260,8,268,16,284,14,298,19,317,11,328,13
%N Start with 10; write down the sum of its digits; add last two terms; repeat.
%H Michael S. Branicky, <a href="/A129888/b129888.txt">Table of n, a(n) for n = 1..10000</a>
%F a(2n+1)=A007618(n+2). For 1<=n<=10: a(2n)=A065075(n+1). - _R. J. Mathar_, Jun 14 2007
%t nxt[n_]:=Module[{t=Total[n]},{t,Total[IntegerDigits[t]]}]; Flatten[ NestList[ nxt,{10,1},35]] (* _Harvey P. Dale_, Mar 29 2011 *)
%o (Python)
%o def next2(n): sd = sum(map(int, str(n))); return [sd, n+sd]
%o def aupton(terms):
%o alst = [10]
%o while len(alst) < terms: alst.extend(next2(alst[-1]))
%o return alst[:terms]
%o print(aupton(66)) # _Michael S. Branicky_, Oct 01 2021
%Y Cf. A007618, A065075.
%K nonn,base
%O 1,1
%A _N. J. A. Sloane_, May 26 2007
%E More terms from _R. J. Mathar_, Jun 14 2007
|
