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

1, 2, 12, 106, 238, 612, 1020, 1732, 2088, 2902, 4336, 5278, 7102, 8500, 9282, 10566, 12640, 15268, 15678, 19488, 21142, 22062, 24766, 27430, 30108, 33892, 37056, 38706, 40786, 41980, 43626, 52360, 54982, 57900, 60496, 66808, 68398, 72670

Offset

0,2

Links

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

Formula

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

Mathematica

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

Prog

(Magma)
p:=NthPrime;
A141133:= func< n | n eq 0 select 1 else p(p(p(p(n))-1)-1) -1 >;
[A141133(n): n in [0..50]]; // G. C. Greubel, Aug 05 2024
(SageMath)
p=nth_prime
def A141133(n): return 1 if n==0 else p(p(p(p(n))-1)-1)-1
[A141133(n) for n in range(51)] # G. C. Greubel, Aug 05 2024

Crossrefs

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

Sequence in context: A125031 A258230 A085867 * A217801 A036077 A275765

Adjacent sequences: A141130 A141131 A141132 * A141134 A141135 A141136

Keyword

nonn

Author

Juri-Stepan Gerasimov, Jul 31 2008

Extensions

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

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

Status

approved