login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A261495 Number of necklaces with n white beads and n^2 black beads. 2
1, 1, 3, 19, 245, 4751, 124936, 4141383, 166237161, 7847250409, 426342182761, 26219808548110, 1801378010581175, 136784412621194274, 11378390032696241010, 1029218687419565103111, 100592759623604055645649, 10565465772302876757883823, 1186893721789951847976898669 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
LINKS
F. Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc. [Cached copy, with permission, pdf format only]
Eric Weisstein's World of Mathematics, Necklace
FORMULA
a(n) = 1/((n+1)*n) * Sum_{d|n} C((n+1)*n/d,n/d) * A000010(d) for n>0, a(0) = 1.
a(n) ~ exp(n+1/2) * n^(n-5/2) / sqrt(2*Pi). - Vaclav Kotesovec, Aug 22 2015
MAPLE
with(numtheory):
a:= n-> `if`(n=0, 1, add(binomial((n^2+n)/d, n/d)
*phi(d), d=divisors(n))/(n^2+n)):
seq(a(n), n=0..20);
MATHEMATICA
a[n_] := If[n==0, 1, DivisorSum[n, Binomial[(n^2+n)/#, n/#]*EulerPhi[#]&]/ (n^2 + n)];
Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Mar 25 2017, translated from Maple *)
PROG
(PARI) a(n) = if(n<1, 1, sumdiv(n, d, binomial((n^2 + n)/d, n/d) * eulerphi(d)) / (n^2 + n));
for(n=0, 20, print1(a(n), ", ")) \\ Indranil Ghosh, Mar 25 2017
CROSSREFS
Main diagonal of A261494.
Sequence in context: A340225 A118023 A054590 * A069344 A305562 A356672
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Aug 21 2015
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 05:39 EDT 2024. Contains 371235 sequences. (Running on oeis4.)