OFFSET
1,2
COMMENTS
For all n > 1, a(n) is semiprime; n is semiprime iff a(n) = n; n is prime iff a(n) = n^2;
For all primes p and k > 0, a(p^k) = p^2.
LINKS
R. Zumkeller, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[n_] := Transpose[ FactorInteger[n]][[1]]; Table[ First[f[n]]*Last[f[n]], {n, 2, 100} ]
p[n_]:=Module[{fctrs=Transpose[FactorInteger[n]][[1]]}, First[fctrs] Last[fctrs]]; Array[p, 70] (* Harvey P. Dale, Jun 11 2011 *)
PROG
(PARI) { for (n = 1, 10000, if (n==1, a=1, f=factor(n)~; a=f[1, 1] * f[1, length(f)]); write("b066048.txt", n, " ", a) ) } \\ Harry J. Smith, Nov 08 2009
CROSSREFS
KEYWORD
nonn,nice,easy
AUTHOR
EXTENSIONS
More terms from Robert G. Wilson v, Jan 03 2002
STATUS
approved