OFFSET
1,1
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
FORMULA
T(n,k) = 2*n + (n+k-1)*(2*n+2*k-1), k>=1, n>=1.
EXAMPLE
Northwest corner:
5....12...23...38...57
14...25...40...59...82
27...42...61...84...111
44...63...86...113..144
MATHEMATICA
f[n_, k_]:=2n+(n+k-1)(2n+2k-1);
TableForm[Table[f[n, k], {n, 1, 10}, {k, 1, 15}]]
Table[f[n-k+1, k], {n, 14}, {k, n, 1, -1}]//Flatten
PROG
(Python)
from math import comb, isqrt
def A185871(n):
a = (m:=isqrt(k:=n<<1))+(k>m*(m+1))
x = n-comb(a, 2)
y = a-x+1
return y*((y+(c:=x<<1)<<1)-3)+x*(c-1)+1 # Chai Wah Wu, Jun 18 2025
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Clark Kimberling, Feb 05 2011
STATUS
approved
