OFFSET
1,1
MAPLE
A074137 := proc(nmax) local a, ipiv, apiv, incsum, nexta ; a := [1, 2, 3] ; ipiv := 2 ; while nops(a) < nmax do apiv := op(ipiv, a) ; while nops(a) < ipiv+apiv-1 do a := [op(a), op(-1, a)+1] ; od ; incsum := add(op(i, a), i=ipiv+1..nops(a)) ; nexta := apiv - ( incsum mod apiv ) ; while nexta <= op(-1, a) do nexta := nexta+apiv ; od ; a := [op(a), nexta] ; ipiv := ipiv+1 ; od ; a ; end ; a074137 := A074137(80) ; for n from 1 to op(-1, a074137) do if not n in a074137 then printf("%d, ", n) ; fi ; od ; # R. J. Mathar, Mar 30 2007
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Aug 28 2002
EXTENSIONS
More terms from R. J. Mathar, Mar 30 2007
STATUS
approved