OFFSET
1,1
FORMULA
a(n) = A131581(n-1)^2.
MATHEMATICA
Table[Prime[1 + PrimePi[Sqrt[10^n]]]^2, {n, 0, 25}] (* T. D. Noe, Aug 19 2011 *)
PROG
(Python)
from math import isqrt
from sympy import nextprime
def a(n): return nextprime(isqrt(10**(n-1)))**2
print([a(n) for n in range(1, 21)]) # Michael S. Branicky, Apr 25 2021, edited Jan 11 2026
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Jaroslav Krizek, Aug 24 2010
STATUS
approved
