A268226 Complement of A056991.

2, 3, 5, 6, 8, 11, 12, 14, 15, 17, 20, 21, 23, 24, 26, 29, 30, 32, 33, 35, 38, 39, 41, 42, 44, 47, 48, 50, 51, 53, 56, 57, 59, 60, 62, 65, 66, 68, 69, 71, 74, 75, 77, 78, 80, 83, 84, 86, 87, 89, 92, 93, 95, 96, 98, 101, 102, 104, 105, 107, 110, 111, 113, 114

Offset

1,1

Links

Colin Barker, Table of n, a(n) for n = 1..1000

H. I. Okagbue, M. O. Adamu, S. A. Iyase, A. A. Opanuga, Sequence of Integers Generated by Summing the Digits of their Squares, Indian Journal of Science and Technology, Vol 8(15), DOI: 10.17485/ijst/2015/v8i15/69912, July 2015.

Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,1,-1).

Formula

From Colin Barker, Feb 06 2016: (Start)

a(n) = a(n-1)+a(n-5)-a(n-6) for n>6.

G.f.: x*(2+x)*(1-x+x^2)*(1+x+x^2) / ((1-x)^2*(1+x+x^2+x^3+x^4)).

(End)

Mathematica

Select[Range@ 100, MemberQ[{0, 2, 3, 5, 6, 8}, #] &@ FixedPoint[Total@ IntegerDigits@ # &, #] &] (* Michael De Vlieger, Feb 01 2016 *)

Prog

(PARI) Vec(x*(2+x)*(1-x+x^2)*(1+x+x^2)/((1-x)^2*(1+x+x^2+x^3+x^4)) + O(x^100)) \\ Colin Barker, Feb 06 2016
(Magma) I:=[2, 3, 5, 6, 8, 11]; [n le 6 select I[n]  else Self(n-1) +Self(n-5) -Self(n-6): n in [1..70]]; // Vincenzo Librandi, Feb 06 2016

Crossrefs

Cf. A056991.

Sequence in context: A028794 A028966 A028767 * A138529 A266409 A035057

Adjacent sequences: A268223 A268224 A268225 * A268227 A268228 A268229

Keyword

nonn,base,easy

Author

N. J. A. Sloane, Jan 30 2016

Extensions

Warning: the Indian Journal of Science and Technology is run by the 'Science & Knowledge Publishing Corporation Limited', which is listed in Jeffrey Beall's list (http://scholarlyoa.com/publishers/) of predatory publishers. - Michel Marcus, Jan 31 2016

Status

approved