OFFSET
1,1
FORMULA
MAPLE
A075319 := proc(nmax) local r, k, a, pairs ; a := [4] ; pairs := [1, 3] ; k := 2 ; r := 5 ; while nops(a) < nmax do while r in pairs do r := r+2 ; od ; if r+2*k in pairs then printf("inconsistency", k) ; fi ; a := [op(a), 2*(r+k)] ; pairs := [op(pairs), r, r+2*k] ; k := k+1 ; od ; RETURN(a) ; end: a := A075319(200) : for n from 1 to nops(a) do printf("%d, ", op(n, a)) ; od ; # R. J. Mathar, Nov 12 2006
PROG
(Python)
from math import isqrt
def A075319(n): return (n+isqrt(5*n**2)&-2)+n-1<<1 # Chai Wah Wu, Aug 16 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Sep 14 2002
EXTENSIONS
More terms from R. J. Mathar, Nov 12 2006
STATUS
approved