|
|
A060221
|
|
Number of orbits of length n under the full 18-shift (whose periodic points are counted by A001027).
|
|
1
|
|
|
18, 153, 1938, 26163, 377910, 5667681, 87460002, 1377481950, 22039920504, 357046533675, 5842582734474, 96402612275775, 1601766528128550, 26772383354990049, 449776041098370870, 7589970692848393200, 128583032925805678350, 2185911559727674682148, 37275544492386193492506
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Number of Lyndon words (aperiodic necklaces) with n beads of 18 colors. - Andrew Howroyd, Dec 10 2017
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (1/n)* Sum_{d|n} mu(d)*A001027(n/d).
G.f.: Sum_{k>=1} mu(k)*log(1/(1 - 18*x^k))/k. - Ilya Gutkovskiy, May 20 2019
|
|
EXAMPLE
|
a(2)=153 since there are 324 points of period 2 in the full 18-shift and 18 fixed points, so there must be (324-18)/2 = 153 orbits of length 2.
|
|
PROG
|
(PARI) a001027(n) = 18^n;
a(n) = (1/n)*sumdiv(n, d, moebius(d)*a001027(n/d)); \\ Michel Marcus, Sep 11 2017
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|