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