OFFSET
1,2
PROG
(PARI) A = matrix(50, 50); for (i = 1, 50, A[1, i] = 2^(i - 1)); for (j = 2, 50, count = 0; p = prime(j); for (d = 1, j, for (x = 1, d, y = d + 1 - x; count++; if (count <= 50, A[j, count] = p*A[x, y]))); d = j; while (count < 50, d++; for (x = 1, j, y = d + 1 - x; count++; if (count <= 50, A[j, count] = p*A[x, y])))); for (d = 1, 50, s = 0; for (y = 1, d, x = d + 1 - y; s += A[x, y]); print1(s, ", ")); \\ David Wasserman, Dec 13 2004
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, May 25 2003
EXTENSIONS
Corrected and extended by David Wasserman, Dec 13 2004
STATUS
approved