|
|
A033125
|
|
Base-7 digits are, in order, the first n terms of the periodic sequence with initial period 1,0,1.
|
|
1
|
|
|
1, 7, 50, 351, 2457, 17200, 120401, 842807, 5899650, 41297551, 289082857, 2023580000, 14165060001, 99155420007, 694087940050, 4858615580351, 34010309062457, 238072163437200, 1666505144060401, 11665536008422807, 81658752058959650, 571611264412717551
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
a(n) = 7*a(n-1) + a(n-3) - 7*a(n-4).
G.f.: x*(1 + x^2)/((1 - x)*(1 - 7*x)*(1 + x + x^2)). - Colin Barker, Dec 24 2015
E.g.f.: exp(-x/2)*(exp(3*x/2)*(25*exp(6*x) - 19) - 6*cos(sqrt(3)*x/2) + 8*sqrt(3)*sin(sqrt(3)*x/2))/171.
a(n) ~ 25*7^n/171. (End)
|
|
MATHEMATICA
|
Table[FromDigits[PadRight[{}, n, {1, 0, 1}], 7], {n, 30}] (* or *) LinearRecurrence[ {7, 0, 1, -7}, {1, 7, 50, 351}, 30] (* Harvey P. Dale, Feb 04 2019 *)
|
|
PROG
|
(PARI) Vec(x*(1+x^2)/((1-x)*(1-7*x)*(1+x+x^2)) + O(x^30)) \\ Colin Barker, Dec 24 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|