login
A026594
a(n) = T(2*n-1, n-2), where T is given by A026584.
17
1, 2, 13, 42, 225, 802, 4235, 15478, 82425, 304156, 1634435, 6064389, 32819839, 122244344, 665162897, 2484851486, 13577768505, 50841782786, 278745377821, 1045763359942, 5749240499515, 21603797860416, 119040956286133, 447922312642212, 2472886893122590, 9315646385012666, 51514464212546865, 194255376492836212
OFFSET
2,2
LINKS
FORMULA
a(n) = A026584(2*n-1, n-2).
MATHEMATICA
T[n_, k_]:= T[n, k]= If[k==0 || k==2*n, 1, If[k==1 || k==2*n-1, Floor[n/2], If[EvenQ[n+k], T[n-1, k-2] + T[n-1, k], T[n-1, k-2] + T[n-1, k-1] + T[n-1, k]]]]; (*T=A026584*)
Table[T[2*n-1, n-2], {n, 2, 40}] (* G. C. Greubel, Dec 13 2021 *)
PROG
(Sage)
@CachedFunction
def T(n, k): # T = A026584
if (k==0 or k==2*n): return 1
elif (k==1 or k==2*n-1): return (n//2)
else: return T(n-1, k-2) + T(n-1, k) if ((n+k)%2==0) else T(n-1, k-2) + T(n-1, k-1) + T(n-1, k)
[T(2*n-1, n-2) for n in (2..40)] # G. C. Greubel, Dec 13 2021
KEYWORD
nonn
EXTENSIONS
Terms a(19) onward added by G. C. Greubel, Dec 13 2021
STATUS
approved