OFFSET
1,4
COMMENTS
For p prime, a(p^k) = p^(k-1). For p and q distinct primes, a(pq) = min(p,q).
a(n) >= n/sopf(n), where sopf is A008472; when the right hand side is an integer, this is an equality.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..5000
EXAMPLE
For n = 12, the sum 2+2+2+3+3 uses each prime factor at most 3 times, so a(12) = 3.
PROG
(PARI) a(n)=if(n<=1, return(0)); my(fm=factor(n), p=prod(i=1, #fm~, 1+x^fm[i, 1] + O(x*x^n))); for(k=0, n, if(polcoef(p^k, n)!=0, return(k)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Franklin T. Adams-Watters, Aug 29 2009
STATUS
approved
