|
%I #44 Oct 06 2020 02:35:43
%S 0,0,0,0,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,1,0,0,0,1,0,0,0,1,0,1,0,0,0,0,
%T 0,1,0,0,0,1,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,1,0,0,0,1,0,0,0,0,0,1,0,0,
%U 0,1,0,1,0,0,0,0,0,1,0,1,0,0,0,1,0,0,0,1,0,1,0,0,0,0,0,1,0,0,0,1,0,1,0,1,0,0,0,1,0,0,0,1,0,1,0,0,0,0,0,2
%N Number of iterations of x -> A003961(x) needed before the result is deficient (sigma(x) < 2x), when starting from x=n.
%C It holds that a(n) <= A336836(n) for all n, because sigma(n) <= A003961(n) for all n (see A286385 for a proof).
%C The first 3 occurs at n = 19399380, the first 4 at n = 195534950863140268380. See A336389.
%C If x and y are relatively prime (i.e., gcd(x,y) = 1), then a(x*y) >= max(a(x),a(y)). Compare to a similar comment in A336915.
%H Antti Karttunen, <a href="/A336835/b336835.txt">Table of n, a(n) for n = 1..65537</a>
%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>
%F If A294934(n) = 1, a(n) = 0, otherwise a(n) = 1 + a(A003961(n)).
%F From _Antti Karttunen_, Aug 21-Sep 01 2020: (Start)
%F For all n >= 1,
%F a(A046523(n)) >= a(n).
%F a(A071364(n)) >= a(n).
%F a(A108951(n)) = A337474(n).
%F a(A025487(n)) = A337475(n).
%F (End)
%e For n = 120, sigma(120) = 360 >= 2*120, thus 120 is not deficient, and we get the next number by applying the prime shift, A003961(120) = 945, and sigma(945) = 1920 >= 945*2, so neither 945 is deficient, so we prime shift once again, and A003961(945) = 9625, which is deficient, as sigma(9625) = 14976 < 2*9625. Thus after two iteration steps we encounter a deficient number, and therefore a(120) = 2.
%t Array[-1 + Length@ NestWhileList[If[# == 1, 1, Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}]] &, #, DivisorSigma[1, #] >= 2 # &] &, 120] (* _Michael De Vlieger_, Aug 27 2020 *)
%o (PARI)
%o A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
%o A336835(n) = { my(i=0); while(sigma(n) >= (n+n), i++; n = A003961(n)); (i); };
%Y Cf. A003961, A005100, A005940, A023196, A046523, A047802, A071364, A108951, A286385, A294934, A336834, A337474, A337475.
%Y Cf. A336389 (position of the first occurrence of a term >= n).
%Y Differs from A294936 for the first time at n=120.
%Y Cf. also A246271, A252459, A336836 and A336915 for similar iterations.
%K nonn
%O 1,120
%A _Antti Karttunen_, Aug 07 2020
|
