OFFSET
1,4
COMMENTS
After zero, sequence contains only terms of A048103.
Note that there are fixed points, like a(15) = 15.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..1728
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
PROG
(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A129251(n) = { my(f = factor(n)); sum(k=1, #f~, (f[k, 2]>=f[k, 1])); };
A276086(n) = { my(i=0, m=1, pr=1, nextpr); while((n>0), i=i+1; nextpr = prime(i)*pr; if((n%nextpr), m*=(prime(i)^((n%nextpr)/pr)); n-=(n%nextpr)); pr=nextpr); m; };
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 07 2019
STATUS
approved