A217623 a(n) = prime(prime(n^2)).

3, 17, 83, 241, 509, 877, 1433, 2063, 2897, 3911, 4943, 6353, 8011, 9661, 11909, 13693, 16141, 18787, 21727, 24781, 28307, 32261, 35801, 40093, 44621, 49139, 54251, 59417, 64853, 70621, 77047, 83617, 90203, 97039, 103991, 112097, 120223, 128683, 136813, 145903

Offset: 1

Author

Vincenzo Librandi, Oct 13 2012

Keywords

• nonn
• easy

Links

• Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A000040, A011757.

Subsequence of A006450.

Programs

Magma

[NthPrime(NthPrime(n^2)): n in [1..50]];

Maple

a:= n-> (ithprime@@2)(n^2):
seq(a(n), n=1..40);  # Alois P. Heinz, Mar 17 2021

Mathematica

Table[Prime[Prime[n^2]], {n, 100}]

PARI

a(n) = prime(prime(n^2)); \\ Michel Marcus, Mar 17 2021

Python

from sympy import prime
def a(n): return prime(prime(n**2))
print([a(n) for n in range(1, 41)]) # Michael S. Branicky, Mar 17 2021