OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
For n=21: 21->{3,7} 3+7=10, 10->{2,5} 2+5=7, 7->{7} 7; 3 terms found {21,10,7}, therefore a(21) = 3.
For n=2: 2->{2} 2, 1 term found {2}, therefore a(2) = 1.
For n=1: 1->{} 0, 2 term found {1,0}, therefore a(1) = 2.
MAPLE
f:= proc(n) option remember;
if isprime(n) then 1
else 1+procname(convert(numtheory:-factorset(n), `+`))
fi
end proc:
f(1):= 2:
map(f, [$1..100]); # Robert Israel, Mar 30 2020
MATHEMATICA
s[n_] := DivisorSum[n, # &, PrimeQ[#] &]; a[1] = 2; a[n_] := Length[ FixedPointList[s, n]] - 1; Array[a, 60, 0] (* Amiram Eldar, Dec 12 2018 *)
PROG
(C++)
int Sum(int x){int acum=0, i=-1; for(; primes[++i]<=x; )if(!(x%primes[i])) acum+=primes[i]; return acum; }
int a(int n){int cn=0, last=n; while(1){cn++; n=Sum(n); if(n==last)break; last=n; } return cn; }
(PARI) a(n)={my(k=1); while(n&&!isprime(n), k++; n=vecsum(factor(n)[, 1])); k} \\ Andrew Howroyd, Dec 12 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Wilmer Emiro Castrillon Calderon, Dec 12 2018
STATUS
approved