OFFSET
0,2
COMMENTS
Equivalently, output sequences with primitive period n from a simple cycling shift register.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..400
E. N. Gilbert and J. Riordan, Symmetry types of periodic sequences, Illinois J. Math., 5 (1961), 657-665.
FORMULA
a(n) = Sum mu(d)*6^(n/d); d|n.
a(0) = 1, a(n) = n * A032164(n).
G.f.: 1 + 6 * Sum_{k>=1} mu(k) * x^k / (1 - 6*x^k). - Ilya Gutkovskiy, Apr 14 2021
MAPLE
with(numtheory):
a:= n-> `if`(n=0, 1, add(mobius(d)*6^(n/d), d=divisors(n))):
seq(a(n), n=0..30); # Alois P. Heinz, Oct 21 2012
MATHEMATICA
a[0] = 1; a[n_] := Sum[MoebiusMu[d]*6^(n/d), {d, Divisors[n]}]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Mar 11 2014 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Apr 20 2000
STATUS
approved