A096482 - OEIS

%I #5 Oct 15 2013 22:32:26

%S 3,67,31,401,241,211,773,2221,1913,7649,3229,1669,2477,10009,5749,

%T 33647,9973,14107,60821,130729,16141,15683,113233,86629,95651,74959,

%U 35617,388403,214993,557093,248909,637003,296843,448451,186481,1145899

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

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

%e a[n]+2n=A000040[1+A000040(A096480[n])]

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

%o (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

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

%K nonn

%O 1,1

%A _Labos Elemer_, Jun 23 2004

%E a(31) - a(36) from _Klaus Brockhaus_, Jun 27 2004