OFFSET
1,2
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16385
PROG
(PARI)
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; \\ This function from Charles R Greathouse IV, Aug 17 2011
A048250(n) = if(n<1, 0, sumdiv(n, d, if(core(d)==d, d)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 10 2017
STATUS
approved