OFFSET
1,1
COMMENTS
We interpret the definition as saying that we find a(n+1) from a(n) by testing a(n)-+2*d (first the minus, then the plus, d=1,2,3,4...., d not used before) for primality and for not already being on the list. This creates a list of odd primes, but there is no proof (yet) that all odd primes are in the sequence. - R. J. Mathar, Jul 29 2007
MAPLE
A094747 := proc(nmax) local a, diffs, d; a := [3] ; diffs := {} ; while nops(a) < nmax do d := 2 ; while true do if isprime(op(-1, a)-d) and not op(-1, a)-d in a and not abs(d) in diffs then a := [op(a), op(-1, a)-d] ; diffs := diffs union {d} ; break ; elif isprime(op(-1, a)+d) and not op(-1, a)+d in a and not abs(d) in diffs then a := [op(a), op(-1, a)+d] ; diffs := diffs union {d} ; break ; else d := d+2 ; fi ; od ; od; RETURN(a) ; end: A094747(100) ; # R. J. Mathar, Jul 29 2007
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, May 24 2004
EXTENSIONS
More terms from R. J. Mathar, Jul 29 2007
STATUS
approved