OFFSET
1,1
COMMENTS
If for some k >= 0 we have q=p+2*(k+1)*(2*k+1) and r=p+2*(k+1)*(2*k+3) then p is in the sequence. This includes A022004 in the case k=0.
It appears that all terms are of the above form.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 17 is in the sequence because with p = 17, q = 19, r = 23, q*r = 437 == 17 (mod 19+23=42).
MAPLE
q:= 2: r:= 3: count:= 0: R:= NULL:
while count < 100 do
p:= q: q:= r: r:= nextprime(r);
if (q*r) mod (q+r) = p then count:= count+1; R:= R, p fi;
od:
R;
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Feb 09 2021
STATUS
approved