A096482 a(n) = A000040[A000040[A096480(n)]] = A006450[A096480(n)] = A000040[A096481(n)].

3, 67, 31, 401, 241, 211, 773, 2221, 1913, 7649, 3229, 1669, 2477, 10009, 5749, 33647, 9973, 14107, 60821, 130729, 16141, 15683, 113233, 86629, 95651, 74959, 35617, 388403, 214993, 557093, 248909, 637003, 296843, 448451, 186481, 1145899

Offset

1,1

Links

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

Example

Both a[n] and a[n]+2n are primes while Pi[a(n)]=A096481(n),Pi[Pi[a(n)]]=A096480(n).
a[n]+2n=A000040[1+A000040(A096480[n])]

Mathematica

Prime[Prime[Min[Flatten[Position[Table[Prime[Prime[n]+1]- Prime[Prime[n]], {n, 1, 5000}], 2*j]]], {j, 1, 100}]]]

Prog

(PARI) {m=36; for(n=1, m, k=1; while((prime(prime(k)+1)-prime(prime(k)))!=2*n, k++); print1(prime(prime(k)), ", "))} - Klaus Brockhaus, Jun 27 2004

Crossrefs

Cf. A000040, A006450, A073124, A072677, A096477-A096481.

Sequence in context: A339599 A105443 A065425 * A139868 A110716 A142035

Adjacent sequences: A096479 A096480 A096481 * A096483 A096484 A096485

Keyword

nonn

Author

Labos Elemer, Jun 23 2004

Extensions

a(31) - a(36) from Klaus Brockhaus, Jun 27 2004

Status

approved