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
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