%I #25 Jul 17 2020 18:11:14
%S 1,3,5,4,5,3,1,1,3,5,4,5,3,1,1,3,5,4,5,3,1,1,3,5,4,5,3,1,1,3,5,4,5,3,
%T 1,1,3,5,4,5,3,1,1,3,5,4,5,3,1,1,3,5,4,5,3,1,1,3,5,4,5,3,1,1,3,5,4,5,
%U 3,1,1,3,5,4,5,3,1,1,3,5,4,5,3,1,1,3,5,4,5,3,1
%N a(n) = A001517(n) mod 7.
%C Periodic with period length 7.
%H Colin Barker, <a href="/A272647/b272647.txt">Table of n, a(n) for n = 0..1000</a>
%H D. H. Lehmer, <a href="http://www.jstor.org/stable/1968107">Arithmetical periodicities of Bessel functions</a>, Annals of Mathematics, 33 (1932): 143-150. The sequence is on page 149.
%H <a href="/index/Rec#order_07">Index entries for linear recurrences with constant coefficients</a>, signature (0,0,0,0,0,0,1).
%F G.f.: (1 + 3*x + 5*x^2 + 4*x^3 + 5*x^4 + 3*x^5 + x^6) / ((1 - x)*(1 + x + x^2 + x^3 + x^4 + x^5 + x^6)). - _Colin Barker_, May 10 2016
%F a(n) = (3*m^6 - 54*m^5 + 365*m^4 - 1140*m^3 + 1582*m^2 - 636*m + 60)/60, where m = n mod 7. - _Luce ETIENNE_, Oct 18 2018
%p f:=proc(n) option remember; if n = 0 then 1 elif n=1 then 3 else f(n-2)+(4*n-2)*f(n-1); fi; end;
%p [seq(f(n) mod 7, n=0..120)];
%t PadRight[{},120,{1,3,5,4,5,3,1}] (* _Harvey P. Dale_, Jul 17 2020 *)
%o (PARI) Vec((1+3*x+5*x^2+4*x^3+5*x^4+3*x^5+x^6)/((1-x)*(1+x+x^2+x^3+x^4+x^5+x^6)) + O(x^50)) \\ _Colin Barker_, May 10 2016
%Y Cf. A001517, A002119, A272648.
%Y Cf. A010876.
%K nonn,easy
%O 0,2
%A _N. J. A. Sloane_, May 09 2016