A172102 Primes in A118482.

3, 11, 29, 59, 101, 239, 619, 809, 4253, 5323, 5923, 6551, 29131, 37277, 48341, 54413, 58711, 60937, 70537, 101063, 110533, 214993, 224603, 417203, 445069, 466537, 473867, 511391, 519089, 534629, 633449, 686269, 713771, 741913, 770767, 1000537

Offset: 1

Author

Jonathan Vos Post, Jan 25 2010

Keywords

• easy
• nonn

Crossrefs

Cf. A118482.

Programs

Maple

isA001358 := proc(n) return ( numtheory[bigomega](n) = 2 ); end proc:
isA109611 := proc(n) isprime(n) and ( isprime(n+2) or isA001358(n+2) ); end proc:
A109611 := proc(n) option remember; local a; if n = 1 then 2; else a := nextprime( procname(n-1) ) ; while not isA109611(a) do a := nextprime(a) ; end do ; return a; end if; end proc:
A118482 := proc(n) option remember ; 1+add( A109611(j),j=1..n) ; end proc:
isA172102 := proc(n) if isprime(n) then for j from 1 do if A118482(j) > n then return false; elif A118482(j) = n then return true; end if; end do ; else false ; end if; end proc:
for n from 1 to 10000000 do if isA172102(n) then printf("%d,\n",n) ; end if; end do ;
# R. J. Mathar, Feb 07 2010

Extensions

Extended by R. J. Mathar, Feb 07 2010