OFFSET
1,1
COMMENTS
All terms are even.
EXAMPLE
4 = 2^2 is a member because is an Ulam number and 3, 5 are a couple of twin primes.
72 = 2^3 * 3^2 is a member because is an Ulam number and 71, 73 are a couple of twin primes.
MAPLE
N:= 10^5: # for terms <= N
U:= [1, 2]: V:= Vector(N): V[3]:= 1: R:= NULL: count:= 0:
for i from 3 do
for k from U[-1]+1 to N do
if V[k] = 1 then
J:= select(`<=`, U +~ k, N);
V[J]:= V[J] +~ 1;
U:= [op(U), k];
if isprime(k+1) and isprime(k-1) then R:= R, k; count:= count+1; fi;
break
fi
od;
if k > N then break fi;
od:
R; # Robert Israel, Jan 24 2025
MATHEMATICA
seq[numUlams_] := Module[{ulams = {1, 2}}, Do[AppendTo[ulams, n = Last[ulams]; While[n++; Length[DeleteCases[Intersection[ulams, n - ulams], n/2, 1, 1]] != 2]; n], {numUlams}]; Select[ulams, And @@ PrimeQ[# + {-1, 1}] &]]; seq[2500] (* Amiram Eldar, Jan 24 2025, after Jean-François Alcover at A002858 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Massimo Kofler, Jan 24 2025
STATUS
approved
