OFFSET
0,6
FORMULA
A(0,k) = 0, A(1,k) = 1; A(n,k) = A002203(k) * A(n-1,k) - (-1)^k * A(n-2,k) for n > 1.
A(n,k) = Pell(k*n)/Pell(k) for k > 0.
EXAMPLE
Square array begins:
0, 0, 0, 0, 0, 0, ...
1, 1, 1, 1, 1, 1, ...
2, 2, 6, 14, 34, 82, ...
3, 5, 35, 197, 1155, 6725, ...
4, 12, 204, 2772, 39236, 551532, ...
5, 29, 1189, 39005, 1332869, 45232349, ...
MATHEMATICA
A[n_, k_] := Fibonacci[k*n, 2]/Fibonacci[k, 2]; A[n_, 0] := n; Table[A[k, n - k], {n, 0, 10}, {k, 0, n}] // Flatten (* Amiram Eldar, May 08 2025 *)
PROG
(PARI) pell(n) = ([2, 1; 1, 0]^n)[2, 1];
a(n, k) = if(k==0, n, pell(k*n)/pell(k));
KEYWORD
AUTHOR
Seiichi Manyama, May 07 2025
STATUS
approved
