login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

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