%I #23 Jan 03 2022 20:39:15
%S 1,1,2,12,107,1178,14805,203885,3002973,46573347,752521980,
%T 12571607865,215925120675,3796546970232,68106673339365,
%U 1243210765414512,23041656826384341
%N a(n) is the n-th prime 3-dimensional Catalan number.
%C "Prime" here is used in the sense of "primitive" or "irreducible".
%H Manuel Wettstein, <a href="http://arxiv.org/abs/1602.07235">Trapezoidal Diagrams, Upward Triangulations, and Prime Catalan Numbers</a>, arXiv:1602.07235 [cs.CG], 2016.
%F Lemma 15 of Wettstein (2016) gives a formula in terms of the 3-dimensional Catalan numbers (A005789).
%p A005789 := proc(n)
%p 2*(3*n)!/(n+2)!/(n+1)!/n! ;
%p end proc:
%p maxn := 30 :
%p Cx := add(A005789(i)*x^i,i=0..maxn) ;
%p d := 3:
%p for i from 0 to maxn do
%p coeftayl(1/Cx^(d*i-1),x=0,i) ;
%p %/(1-d*i) ;
%p printf("%d,",%) ;
%p end do: # _R. J. Mathar_, Feb 27 2018
%t A005789[n_] := 2*(3*n)!/(n+2)!/(n+1)!/n!; Maxn = 30; Cx = Sum[A005789[i]* x^i, {i, 0, Maxn}]; d = 3; Reap[For[i = 0, i <= Maxn, i++, sc = SeriesCoefficient[1/Cx^(d*i-1), {x, 0, i}]; Sow[sc/(1-d*i)]]][[2, 1]] (* _Jean-François Alcover_, Mar 24 2018, after _R. J. Mathar_ *)
%Y Primitive terms from A000108, A005789.
%K nonn
%O 0,3
%A _N. J. A. Sloane_, Feb 24 2016
%E 7 more terms from _R. J. Mathar_, Feb 27 2018
|