login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A054365 Number of unlabeled 5-gonal cacti having n polygons. 5
1, 1, 1, 3, 17, 102, 811, 6626, 58385, 532251, 5011934, 48344880, 475982471, 4766639628, 48434621610, 498363430232, 5184274255789, 54451326151253, 576810990484823, 6156943228387305, 66170786572330174, 715564777086617766 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

COMMENTS

Also, the number of noncrossing partitions up to rotation composed of n blocks of size 5. - Andrew Howroyd, May 04 2018

LINKS

Andrew Howroyd, Table of n, a(n) for n = 0..200

Miklos Bona, Michel Bousquet, Gilbert Labelle and Pierre Leroux, Enumeration of m-ary cacti, Advances in Applied Mathematics, 24 (2000), 22-56 (pdf, dvi).

Luke Lippstreu, Jorge Mago, Marcus Spradlin, Anastasia Volovich, Weak Separation, Positivity and Extremal Yangian Invariants, arXiv:1906.11034 [hep-th], 2019.

Index entries for sequences related to cacti

FORMULA

a(n) = ((Sum_{d|n} phi(n/d)*binomial(5*d, d)) + (Sum_{d|gcd(n-1, 5)} phi(d)*binomial(5*n/d, (n-1)/d)))/(5*n) - binomial(5*n, n)/(4*n+1) for n > 0. - Andrew Howroyd, May 04 2018

MAPLE

with(combinat): with(numtheory): m := 5: for p from 2 to 28 do s1 := 0: s2 := 0: for d from 1 to p do if p mod d = 0 then s1 := s1+phi(p/d)*binomial(m*d, d) fi: od: for d from 1 to p-1 do if gcd(m, p-1) mod d = 0 then s2 := s2+phi(d)*binomial((p*m)/d, (p-1)/d) fi: od: printf(`%d, `, (s1+s2)/(m*p)-binomial(m*p, p)/(p*(m-1)+1)) od: # Zerinvary Lajos, Dec 01 2006

MATHEMATICA

a[0] = 1;

a[n_] := (DivisorSum[n, EulerPhi[n/#] Binomial[5#, #]&] + DivisorSum[GCD[n - 1, 5], EulerPhi[#] Binomial[5n/#, (n-1)/#]&])/(5n) - Binomial[5n, n]/ (4n+1);

Table[a[n], {n, 0, 21}] (* Jean-Fran├žois Alcover, Jul 01 2018, after Andrew Howroyd *)

PROG

(PARI) a(n) = {if(n==0, 1, (sumdiv(n, d, eulerphi(n/d)*binomial(5*d, d)) + sumdiv(gcd(n-1, 5), d, eulerphi(d)*binomial(5*n/d, (n-1)/d)))/(5*n) - binomial(5*n, n)/(4*n+1))} \\ Andrew Howroyd, May 04 2018

CROSSREFS

Column k=5 of A303694.

Cf. A054363, A054364, A303871.

Sequence in context: A074565 A339565 A241768 * A116886 A163064 A020069

Adjacent sequences:  A054362 A054363 A054364 * A054366 A054367 A054368

KEYWORD

nonn

AUTHOR

Simon Plouffe

EXTENSIONS

More terms from Zerinvary Lajos, Dec 01 2006

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 27 11:07 EST 2021. Contains 341649 sequences. (Running on oeis4.)