cp's OEIS Frontend

A047750 If n mod 2 = 0 then m := n/2 and a(n) = (3*m)!*(5*m+1)/((m+1)!*(2*m+1)!); otherwise m := (n-1)/2, a(n) = 6*(3*m+2)!/(m!*(2*m+3)!).

%I A047750 #43 Mar 21 2024 08:32:02
%S A047750 1,2,3,6,11,24,48,110,231,546,1183,2856,6324,15504,34884,86526,197087,
%T A047750 493350,1134705,2861430,6633315,16829280,39268320,100134216,234930276,
%U A047750 601661144,1418201268,3645533040,8627761528,22249511328
%N A047750 If n mod 2 = 0 then m := n/2 and a(n) = (3*m)!*(5*m+1)/((m+1)!*(2*m+1)!); otherwise m := (n-1)/2, a(n) = 6*(3*m+2)!/(m!*(2*m+3)!).
%H A047750 L. W. Beineke and R. E. Pippert, <a href="http://dx.doi.org/10.4153/CJM-1974-006-x">Enumerating dissectable polyhedra by their automorphism groups</a>, Canad. J. Math., 26 (1974), 50-67.
%F A047750 From _Gary W. Adamson_, Jul 14 2011: (Start)
%F A047750 a(n) = sum of top row terms in M^n, M = the infinite square production matrix:
%F A047750   1, 1, 0, 0, 0, 0, ...
%F A047750   0, 0, 1, 0, 0, 0, ...
%F A047750   1, 1, 0, 1, 0, 0, ...
%F A047750   0, 0, 1, 0, 1, 0, ...
%F A047750   1, 1, 0, 1, 0, 1, ...
%F A047750   ... (End)
%F A047750 8*n*(n+2)*a(n) + 4*(7*n^2 - 7*n - 17)*a(n-1) + 6*(-9*n^2 + 9*n - 17)*a(n-2) - 21*(3*n-5)*(3*n-7)*a(n-3) = 0. - _R. J. Mathar_, Jul 10 2013
%F A047750 From _Robert A. Russell_, Mar 20 2024: (Start)
%F A047750 a(n) = V(n) in the Beineke and Pippert link.
%F A047750 G.f.: 2*(G(z^2) - 1)/z + 2*G(z^2)^2 - G(z^2), where G(z) = 1 + z*G(z)^3 is the g.f. for A001764. (End)
%p A047750 series(RootOf(x*A^3-2*A^2+3*A-1, A)^2, x=0, 30);  # _Mark van Hoeij_, May 16 2013
%t A047750 a[0] = 1; a[1] = 2; a[n_] := a[n] = 3(2n+3)(3n-4)(3n-2)a[n-2]/(4n(n+2)(2n+1)) + (3(18n+16)a[n-1])/(4n(n+2)(2n+1)); Table[a[n], {n, 0, 29}] (* _Jean-François Alcover_, Dec 02 2016 *)
%t A047750 Table[If[OddQ[n],6Binomial[(3n+1)/2,n+1]/(n+2),(5n+2)Binomial[3n/2,n/2] / ((n+1)(n+2))],{n,0,30}] (* _Robert A. Russell_, Feb 16 2024 *)
%o A047750 (PARI) a047750(n)={if(n%2,my(m=(n-1)/2);6*(3*m+2)!/(m!*(2*m+3)!),my(m=n/2);(3*m)!*(5*m+1)/((m+1)!*(2*m+1)!))};
%o A047750 for(k=0,29,print1(a047750(k),", ")) \\ _Hugo Pfoertner_, Mar 07 2020
%Y A047750 Cf. A001764, A047749, A047760, A047773.
%K A047750 nonn
%O A047750 0,2
%A A047750 _N. J. A. Sloane_