OFFSET
1,1
MATHEMATICA
start = {2, 4}; primes = Table[Prime[n], {n, 1, 10000}];
f[x_, y_] := If[MemberQ[primes, x*y + 1], x*y + 1]
b[x_] :=
Block[{w = x},
Select[Union[
Flatten[AppendTo[w,
Table[f[w[[i]], w[[j]]], {i, 1, Length[w]}, {j, 1, i}]]]], # <
50000 &]];
t = FixedPoint[b, start] (* A192581 *)
CROSSREFS
KEYWORD
nonn,fini,full
AUTHOR
Clark Kimberling, Jul 04 2011
STATUS
approved