%I #16 Sep 01 2016 11:51:53
%S 2,14,101,708,4958,34706,242945,1700616,11904314,83330198,583311389,
%T 4083179724,28582258070,200075806490,1400530645433,9803714518032,
%U 68626001626226,480382011383582,3362674079685077,23538718557795540,164771029904568782,1153397209331981474
%N Base-7 digits are, in order, the first n terms of the periodic sequence with initial period 2,0,3,1.
%H Colin Barker, <a href="/A037726/b037726.txt">Table of n, a(n) for n = 1..1000</a>
%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (7,0,0,1,-7).
%F a(n) = 7*a(n-1) + a(n-4) - 7*a(n-5).
%F From _Colin Barker_, Dec 24 2015: (Start)
%F a(n) = 1/200*(-25*(-1)^n+(8+6*i)*(-i)^n+(8-6*i)*i^n+59*7^n-50) where i=sqrt(-1).
%F G.f.: x*(2+3*x^2+x^3) / ((1-x)*(1+x)*(1-7*x)*(1+x^2)).
%F (End)
%t Table[FromDigits[PadRight[{},n,{2,0,3,1}],7],{n,30}] (* or *) LinearRecurrence[{7,0,0,1,-7},{2,14,101,708,4958},30] (* _Harvey P. Dale_, Sep 01 2016 *)
%o (PARI) Vec(x*(2+3*x^2+x^3)/((1-x)*(1+x)*(1-7*x)*(1+x^2)) + O(x^30)) \\ _Colin Barker_, Dec 24 2015
%K nonn,base,easy
%O 1,1
%A _Clark Kimberling_
|