OFFSET
1,1
EXAMPLE
sqrt(16) = 4. replace 16 with prime(4) = 7.
MATHEMATICA
snp[n_]:=Module[{s=Sqrt[n]}, If[IntegerQ[s], Prime[s], n]]; Array[snp, 80] (* Harvey P. Dale, Jun 21 2013 *)
PROG
(PARI) replsqpx(n) = { for(x=1, n, if(issquare(x), y=prime(floor(sqrt(x))), y=x); print1(floor(y)", ")) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Cino Hilliard, Aug 23 2004
STATUS
approved