OFFSET
1,3
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
Antti Karttunen, First 105368 terms (from 0 to 9699690)
FORMULA
For all n >= 1, A328307(a(n)) = n.
PROG
(PARI)
A003415checked(n) = if(n<=1, 0, my(f=factor(n), s=0); for(i=1, #f~, if(f[i, 2]>=f[i, 1], return(0), s += f[i, 2]/f[i, 1])); (n*s));
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; };
isA099308(n) = { while(n>1, n = A003415checked(n)); (n); };
isA328116(n) = isA099308(A276086(n));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 08 2019
STATUS
approved