OFFSET
0,2
COMMENTS
Equivalently, output sequences with primitive period n from a simple cycling shift register.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..2095 (terms 0..650 from Alois P. Heinz)
E. N. Gilbert and J. Riordan, Symmetry types of periodic sequences, Illinois J. Math., 5 (1961), 657-665.
FORMULA
a(n) = Sum_{d|n} mu(d)*3^(n/d).
a(0) = 1, a(n) = n * A027376(n).
a(n) = 3 * A034741(n).
G.f.: 1 + 3 * Sum_{k>=1} mu(k) * x^k / (1 - 3*x^k). - Ilya Gutkovskiy, Apr 14 2021
MAPLE
with(numtheory):
a:= n-> `if`(n=0, 1, add(mobius(d)*3^(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]*3^(n/d), {d, Divisors[n]}]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Mar 11 2014, after Alois P. Heinz *)
PROG
(PARI) a(n) = if(n==0, 1, sumdiv(n, d, moebius(d) * 3^(n/d) )); \\ Joerg Arndt, Apr 14 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Apr 20 2000
STATUS
approved