A172102 - OEIS

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

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,1

LINKS

Table of n, a(n) for n=1..36.

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

CROSSREFS

Cf. A118482.

Sequence in context: A069350 A021005 A053845 * A242807 A188475 A072610

Adjacent sequences: A172099 A172100 A172101 * A172103 A172104 A172105

KEYWORD

easy,nonn

AUTHOR

Jonathan Vos Post, Jan 25 2010

EXTENSIONS

Extended by R. J. Mathar, Feb 07 2010

STATUS

approved