OFFSET
1,2
COMMENTS
A046523(a(n))=A046523(n); A046523(n)<=a(n)<=n; A001221(a(n))=A001221(n), A001222(a(n))=A001222(n); A020639(a(n))=2, A006530(a(n))=A000040(A001221(n))<=A006530(n); A000005(a(n))=A000005(n);
LINKS
Daniel Forgues, Table of n, a(n) for n=1..100000
FORMULA
In prime factorization of n, replace least prime by 2, next least by 3, etc.
EXAMPLE
a(105875) = a(5*5*5*7*11*11) = 2*2*2*3*5*5 = 600.
MATHEMATICA
Table[ e = Last /@ FactorInteger[n]; Product[Prime[i]^e[[i]], {i, Length[e]}], {n, 88}] (* Ray Chandler, Sep 23 2005 *)
PROG
(Haskell)
a071364 = product . zipWith (^) a000040_list . a124010_row
-- Reinhard Zumkeller, Feb 19 2012
(PARI) a(n) = f = factor(n); for (i=1, #f~, f[i, 1] = prime(i)); factorback(f); \\ Michel Marcus, Jun 13 2014
(Python)
from math import prod
from sympy import prime, factorint
def A071364(n): return prod(prime(i+1)**p[1] for i, p in enumerate(sorted(factorint(n).items()))) # Chai Wah Wu, Sep 16 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, May 21 2002
EXTENSIONS
Extended by Ray Chandler, Sep 23 2005
STATUS
approved