A141136 a(n) = prime(prime(A028815(n) - 1) - 1) - 1.

1, 2, 6, 28, 52, 112, 172, 270, 316, 420, 592, 700, 910, 1060, 1150, 1290, 1510, 1782, 1830, 2212, 2376, 2472, 2740, 2998, 3256, 3630, 3930, 4078, 4270, 4390, 4546, 5350, 5590, 5866, 6100, 6658, 6802, 7186, 7602, 7828, 8218, 8520, 8712, 9310, 9438, 9732

Offset: 0

Author

Juri-Stepan Gerasimov, Jul 31 2008

Keywords

• nonn

Links

• G. C. Greubel, Table of n, a(n) for n = 0..10000

Crossrefs

Cf. A000040, A028815, A141130, A141132, A141133, A141138.

Programs

Magma

p:=NthPrime;
A141136:= func< n | n eq 0 select 1 else p(p(p(n))-1) -1 >;
[A141136(n): n in [0..50]]; // G. C. Greubel, Aug 05 2024

Mathematica

A141136[n_]:= With[{p=Prime}, If[n==0, 1, p[p[p[n]]-1]-1]];
Table[A141136[n], {n,0,60}] (* G. C. Greubel, Aug 05 2024 *)

SageMath

p=nth_prime
def A141136(n): return 1 if n==0 else p(p(p(n))-1)-1
[A141136(n) for n in range(51)] # G. C. Greubel, Aug 05 2024

Formula

a(n) = A000040(A000040(A000040(n)) - 1) - 1, with a(0) = 1. - G. C. Greubel, Aug 05 2024

Extensions

Corrected and extended by D. S. McNeil, Mar 21 2009

Offset changed by G. C. Greubel, Aug 05 2024