OFFSET
1,2
COMMENTS
When interspersed with 0's, exponents in expansion of A065481 as a product zeta(n)^(-a(n)).
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
Juggling Information Service, Site Swap FAQs
M. Macauley, Braids and Juggling patterns, Thesis (Harvey Mudd College) 2003, eq. (A1)
G. Niklasch, Some number theoretical constants: 1000-digit values [Cached copy]
FORMULA
a(n) ~ 3^n/n. - Vaclav Kotesovec, Mar 05 2016
Inverse Euler transform of A133494. - Alois P. Heinz, Jun 23 2018
EXAMPLE
We have one period 1 (2), two period 2 (31/13 and 40/04) and six period three 2-ball siteswaps (312, 330, 411, 420, 501, 600) (The average of the digits is always 2).
MAPLE
[seq(DistSS(p, 2), p=1..60)];
A065178 := proc(n)
add( mobius(n/d)*(3^d-2^d), d=numtheory[divisors](n)) /n ;
end proc:
seq(A065178(n), n=1..30) ; # R. J. Mathar, Aug 05 2015
MATHEMATICA
a[n_] := DivisorSum[n, MoebiusMu[n/#] * (3^#-2^#)&] / n; Array[a, 30] (* Jean-François Alcover, Mar 05 2016, after R. J. Mathar *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 19 2001
STATUS
approved