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!)
A052394 Number of unlabeled 4-ary cacti having n polygons. 4

%I #22 Jan 09 2024 10:54:23

%S 1,1,4,10,44,197,1228,7692,52828,373636,2735952,20506258,156922676,

%T 1221179926,9642496488,77092885016,623121750844,5085013101164,

%U 41850600967984,347060754685888,2897800158952304,24344668688424333,205667187527660076,1746375819789491996,14898241072028602276

%N Number of unlabeled 4-ary cacti having n polygons.

%H Andrew Howroyd, <a href="/A052394/b052394.txt">Table of n, a(n) for n = 0..200</a>

%H Miklos Bona, Michel Bousquet, Gilbert Labelle, and Pierre Leroux, <a href="https://doi.org/10.1006/aama.1999.0665">Enumeration of m-ary cacti</a>, Advances in Applied Mathematics, 24 (2000), 22-56.

%H <a href="/index/Ca#cacti">Index entries for sequences related to cacti</a>

%F a(n) = (1/n)*(Sum_{d|n} phi(n/d)*binomial(4*d, d)) - 3*binomial(4*n, n)/(3*n+1) for n > 0. - _Andrew Howroyd_, May 02 2018

%F a(n) ~ 2^(8*n + 1/2) / (sqrt(Pi) * n^(5/2) * 3^(3*n + 3/2)). - _Vaclav Kotesovec_, Jul 17 2017

%t a[n_] := If[n == 0, 1, (Binomial[4*n, n]/(3*n + 1) + DivisorSum[n, Binomial[4*#, #]*EulerPhi[n/#]*Boole[# < n] & ])/n]; Table[a[n], {n, 0, 24}] (* _Jean-François Alcover_, Jul 17 2017 *)

%o (PARI) a(n) = if(n==0, 1, sumdiv(n, d, eulerphi(n/d)*binomial(4*d, d))/n - 3*binomial(4*n, n)/(3*n+1)) \\ _Andrew Howroyd_, May 02 2018

%Y Column k=4 of A303912.

%Y Cf. A054362, A052395.

%K nonn

%O 0,3

%A _Simon Plouffe_

%E More terms from _Jean-François Alcover_, Jul 17 2017

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 20:05 EDT 2024. Contains 371254 sequences. (Running on oeis4.)