OFFSET
0,2
COMMENTS
Also the sum of the first n^2+1 terms of A004016.
LINKS
Peter Kagey, Table of n, a(n) for n = 0..1000
Code Golf Stack Exchange, Triangular Lattice Points close to the Origin
FORMULA
a(n) = A053416(2*n).
PROG
(PARI) a(n)=1+2*vecsum(Vec(qfrep([2, 1; 1, 2], n^2, 1))) /* Georg Fischer, Apr 07 2022 (from codegolf link) */
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Kagey, Jun 28 2019
STATUS
approved