|
|
A060219
|
|
Number of orbits of length n under the full 16-shift (whose periodic points are counted by A001025).
|
|
1
|
|
|
16, 120, 1360, 16320, 209712, 2795480, 38347920, 536862720, 7635496960, 109951057896, 1599289640400, 23456246655680, 346430740566960, 5146970983535160, 76861433640386288, 1152921504338411520
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Number of monic irreducible polynomials of degree n over GF(16). - Robert Israel, Jan 07 2015
Number of Lyndon words (aperiodic necklaces) with n beads of 16 colors. - Andrew Howroyd, Dec 10 2017
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (1/n)* Sum_{d|n} mu(d)*16^(n/d).
G.f.: Sum_{k>=1} mu(k)*log(1/(1 - 16*x^k))/k. - Ilya Gutkovskiy, May 19 2019
|
|
EXAMPLE
|
a(2)=120 since there are 256 points of period 2 in the full 16-shift and 16 fixed points, so there must be (256-16)/2 = 120 orbits of length 2.
|
|
MAPLE
|
f:= (n, p) -> add(numtheory:-mobius(d)*p^(n/d), d=numtheory:-divisors(n))/n:
|
|
PROG
|
(PARI) a(n) = sumdiv(n, d, moebius(d)*16^(n/d))/n; \\ Michel Marcus, Jan 07 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|