login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A029912 Start with n; repeatedly sum prime factors (counted with multiplicity) and add 1, until reach 1, 6 or a prime. 1

%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_

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 21 12:33 EST 2024. Contains 370235 sequences. (Running on oeis4.)