|
%I #11 Aug 21 2022 11:47:31
%S 1,3,5,5,6,6,7,7,7,7,7,7,7,7,7,7,7,7,7,7,11,7,7,7,11,7,7,7,11,11,11,
%T 11,7,7,13,11,7,7,17,7,13,13,7,7,7,7,7,7,7,13,11,7,7,7,17,7,23,11,13,
%U 13,7,7,7,13,19,17,7,7,7,7,13,13,7,7,7,7,19,19,7,7,13,7,7,7,23,7
%N Start with n; repeatedly sum prime factors (counted with multiplicity) and add 1, until reach 1, 6 or a prime.
%C If p is in A023200 then a(3*p) = p+4. It appears that all n > 35 such that a(n) > n/3 are 3*p for p in A023200. - _Robert Israel_, Dec 18 2019
%H Robert Israel, <a href="/A029912/b029912.txt">Table of n, a(n) for n = 1..10000</a>
%e 20 -> 2+2+5+1 = 10 -> 2+5+1 = 8 -> 2+2+2+1 = 7 so a(20)=7.
%p f:= proc(n) option remember;
%p local v;
%p v:= add(t[1]*t[2],t=ifactors(n)[2])+1;
%p if v = 1 or v = 6 or isprime(v) then return v fi;
%p procname(v)
%p end proc:
%p map(f, [$1..100]); # _Robert Israel_, Dec 18 2019
%t a[n_] := a[n] = If[n==1, 1, Module[{v}, v = Sum[t[[1]]*t[[2]], {t, FactorInteger[n]}]+1; If[v==1 || v==6 || PrimeQ[v], Return[v]]; a[v]]];
%t a /@ Range[100] (* _Jean-François Alcover_, Aug 21 2022, after _Robert Israel_ *)
%K nonn
%O 1,2
%A _Dann Toliver_
|
