|
%I #34 Apr 14 2021 22:25:07
%S 1,3,6,24,72,240,696,2184,6480,19656,58800,177144,530640,1594320,
%T 4780776,14348640,43040160,129140160,387400104,1162261464,3486725280,
%U 10460350992,31380882456,94143178824,282428998560,847288609200,2541864234000,7625597465304
%N Number of ternary sequences with primitive period n.
%C Equivalently, output sequences with primitive period n from a simple cycling shift register.
%H Seiichi Manyama, <a href="/A054718/b054718.txt">Table of n, a(n) for n = 0..2095</a> (terms 0..650 from Alois P. Heinz)
%H E. N. Gilbert and J. Riordan, <a href="http://projecteuclid.org/euclid.ijm/1255631587">Symmetry types of periodic sequences</a>, Illinois J. Math., 5 (1961), 657-665.
%F a(n) = Sum_{d|n} mu(d)*3^(n/d).
%F a(0) = 1, a(n) = n * A027376(n).
%F a(n) = 3 * A034741(n).
%F G.f.: 1 + 3 * Sum_{k>=1} mu(k) * x^k / (1 - 3*x^k). - _Ilya Gutkovskiy_, Apr 14 2021
%p with(numtheory):
%p a:= n-> `if`(n=0, 1, add(mobius(d)*3^(n/d), d=divisors(n))):
%p seq(a(n), n=0..30); # _Alois P. Heinz_, Oct 21 2012
%t 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_ *)
%o (PARI) a(n) = if(n==0,1,sumdiv(n,d, moebius(d) * 3^(n/d) )); \\ _Joerg Arndt_, Apr 14 2013
%Y Column k=3 of A143324.
%Y Cf. A027375, A054719, A054720, A054721.
%K nonn
%O 0,2
%A _N. J. A. Sloane_, Apr 20 2000
|
