login
Smallest k such that sopf(n+k) = sopf(k), where sopf = A008472.
5

%I #32 Jul 31 2021 05:53:52

%S 5,2,7,4,114,2,5,8,13,10,25,4,5,2,19,16,85,6,5,5,209,22,25,3,493,26,

%T 31,4,20,2,5,32,7,34,516,12,33,38,10,10,99,6,5,44,57,46,25,6,5,50,49,

%U 52,52,18,855,8,61,58,295,4,261,2,91,64,602,6,5,68,21,10,25,9,7,74,13,76

%N Smallest k such that sopf(n+k) = sopf(k), where sopf = A008472.

%D J. Earls, Mathematical Bliss, Pleroma Publications, 2009, pages 99-100. ASIN: B002ACVZ6O [From _Jason Earls_, Nov 26 2009]

%H Harry J. Smith, <a href="/A065925/b065925.txt">Table of n, a(n) for n=1..1000</a>

%H Carlos Rivera, <a href="http://www.primepuzzles.net/conjectures/conj_025.htm">Conjecture 25. sopf(n) = sopf(n+k)</a>, The Prime Puzzles and Problems Connection.

%e a(6) = 2 because A008472(2) = A008472(6+2) = 2, but A008472(1) = 0 doesn't equal A008472(6+1) = 7.

%t Table[k = 1; While[Total[FactorInteger[n + k][[All, 1]]] != Total[FactorInteger[k][[All, 1]]], k++]; k, {n, 76}] (* _Michael De Vlieger_, Jan 11 2017 *)

%o (PARI)

%o sopf(n) = local(fac, i); fac=factor(n); sum(i=1,matsize(fac)[1],fac[i,1])

%o A065925(m)={local(k,n); for(k=1,m,n=1; while(sopf(n)!=sopf(n+k), n++); print1(n,","))} \\ _Klaus Brockhaus_

%o (Python)

%o from sympy import primefactors

%o from itertools import count, dropwhile

%o def sopf(n): return sum(p for p in primefactors(n))

%o def a(n):

%o k = 1

%o while sopf(n+k) != sopf(k): k += 1

%o return k

%o print([a(n) for n in range(1, 77)]) # _Michael S. Branicky_, May 02 2021

%Y Cf. A008472, A065926, A065927.

%K nonn

%O 1,1

%A _Jason Earls_, Nov 28 2001