login
A217623
a(n) = prime(prime(n^2)).
1
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,1
LINKS
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}]
PROG
(Magma) [NthPrime(NthPrime(n^2)): n in [1..50]];
(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
(PARI) a(n) = prime(prime(n^2)); \\ Michel Marcus, Mar 17 2021
CROSSREFS
Subsequence of A006450.
Sequence in context: A225342 A362390 A194596 * A083217 A225307 A318769
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Oct 13 2012
STATUS
approved