|
%I #27 Sep 08 2022 08:46:15
%S 2,3,5,6,8,11,12,14,15,17,20,21,23,24,26,29,30,32,33,35,38,39,41,42,
%T 44,47,48,50,51,53,56,57,59,60,62,65,66,68,69,71,74,75,77,78,80,83,84,
%U 86,87,89,92,93,95,96,98,101,102,104,105,107,110,111,113,114
%N Complement of A056991.
%H Colin Barker, <a href="/A268226/b268226.txt">Table of n, a(n) for n = 1..1000</a>
%H H. I. Okagbue, M. O. Adamu, S. A. Iyase, A. A. Opanuga, <a href="http://www.indjst.org/index.php/indjst/article/view/69912">Sequence of Integers Generated by Summing the Digits of their Squares</a>, Indian Journal of Science and Technology, Vol 8(15), DOI: 10.17485/ijst/2015/v8i15/69912, July 2015.
%H <a href="/index/Rec#order_06">Index entries for linear recurrences with constant coefficients</a>, signature (1,0,0,0,1,-1).
%F From _Colin Barker_, Feb 06 2016: (Start)
%F a(n) = a(n-1)+a(n-5)-a(n-6) for n>6.
%F G.f.: x*(2+x)*(1-x+x^2)*(1+x+x^2) / ((1-x)^2*(1+x+x^2+x^3+x^4)).
%F (End)
%t Select[Range@ 100, MemberQ[{0, 2, 3, 5, 6, 8}, #] &@ FixedPoint[Total@ IntegerDigits@ # &, #] &] (* _Michael De Vlieger_, Feb 01 2016 *)
%o (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
%o (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
%Y Cf. A056991.
%K nonn,base,easy
%O 1,1
%A _N. J. A. Sloane_, Jan 30 2016
%E 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
|
