cp's OEIS Frontend

A227580 Number of lattice paths from {n}^3 to {0}^3 using steps that decrement one component such that for each point (p_1,p_2,p_3) we have p_1<=p_2<=p_3.

%I A227580 #12 Dec 20 2020 07:48:07
%S A227580 1,1,14,290,7680,238636,8285506,312077474,12509563082,526701471002,
%T A227580 23076216957520,1044813920439200,48630132961189400,
%U A227580 2317337976558074760,112689430179458971738,5577655817793682738378,280392321290875174774106,14290804691034216155457274
%N A227580 Number of lattice paths from {n}^3 to {0}^3 using steps that decrement one component such that for each point (p_1,p_2,p_3) we have p_1<=p_2<=p_3.
%H A227580 Alois P. Heinz, <a href="/A227580/b227580.txt">Table of n, a(n) for n = 0..200</a>
%F A227580 a(n) ~ 2^(6*n+10)/(sqrt(3)*Pi*(5*n)^4). - _Vaclav Kotesovec_, Jul 18 2013
%e A227580 a(2) = 14: [(2,2,2),(0,2,2),(0,0,2),(0,0,0)], [(2,2,2),(0,2,2),(0,0,2),(0,0,1),(0,0,0)], [(2,2,2),(0,2,2),(0,1,2),(0,0,2),(0,0,0)], [(2,2,2),(0,2,2),(0,1,2),(0,0,2),(0,0,1),(0,0,0)], [(2,2,2),(0,2,2),(0,1,2),(0,1,1),(0,0,1),(0,0,0)], [(2,2,2),(1,2,2),(0,2,2),(0,0,2),(0,0,0)], [(2,2,2),(1,2,2),(0,2,2),(0,0,2),(0,0,1),(0,0,0)], [(2,2,2),(1,2,2),(0,2,2),(0,1,2),(0,0,2),(0,0,0)], [(2,2,2),(1,2,2),(0,2,2),(0,1,2),(0,0,2),(0,0,1),(0,0,0)], [(2,2,2),(1,2,2),(0,2,2),(0,1,2),(0,1,1),(0,0,1),(0,0,0)], [(2,2,2),(1,2,2),(1,1,2),(0,1,2),(0,0,2),(0,0,0)], [(2,2,2),(1,2,2),(1,1,2),(0,1,2),(0,0,2),(0,0,1),(0,0,0)], [(2,2,2),(1,2,2),(1,1,2),(0,1,2),(0,1,1),(0,0,1),(0,0,0)], [(2,2,2),(1,2,2),(1,1,2),(1,1,1),(0,1,1),(0,0,1),(0,0,0)].
%p A227580 a:= proc(n) option remember; `if`(n<3, [1, 1, 14][n+1],
%p A227580       ((n+1)*(665*n^3-1433*n^2+980*n-204) *a(n-1)
%p A227580        -(n-2)*(1615*n^3-3218*n^2+1521*n-342) *a(n-2)
%p A227580        +192*(5*n-1)*(n-3)*(n-2)^2 *a(n-3)) /
%p A227580        (2*(n+2)*(5*n-6)*(n+1)^2))
%p A227580     end:
%p A227580 seq(a(n), n=0..30);
%t A227580 b[l_] := b[l] = If[l[[-1]] == 0, 1, Sum[Sum[b[ReplacePart[l, i -> j]], {j, If[i == 1, 0, l[[i-1]]], l[[i]] - 1}], {i, 1, Length[l]}]];
%t A227580 a[n_] := b[Array[n&, 3]];
%t A227580 a /@ Range[0, 30] (* _Jean-François Alcover_, Dec 20 2020, after _Alois P. Heinz_ in A227578 *)
%Y A227580 Column k=3 of A227578.
%K A227580 nonn
%O A227580 0,3
%A A227580 _Alois P. Heinz_, Jul 16 2013