OFFSET
1,3
COMMENTS
Note that A130158 is a list of the positive divisors of the terms of this sequence.
EXAMPLE
MAPLE
A130157 := proc(nmax) local a, a130158, n, anext, i ; a := [1] ; a130158 := [] ; while nops(a) < nmax do n := nops(a)+1 ; a130158 := [op(a130158), op(numtheory[divisors](op(-1, a)))] ; anext :=0 ; for i from 1 to nops(a) do if gcd(op(i, a), op(n-1, a130158)) = 1 then anext := anext+1 ; fi ; od ; a := [op(a), anext] ; od ; RETURN(a) ; end: A130157(80) ; # R. J. Mathar, Jun 12 2007
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, May 13 2007
EXTENSIONS
More terms from R. J. Mathar, Jun 12 2007
Edited by Charles R Greathouse IV, Apr 27 2010
Edited by N. J. A. Sloane, May 09 2010
STATUS
approved