login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A336835 Number of iterations of x -> A003961(x) needed before the result is deficient (sigma(x) < 2x), when starting from x=n. 14
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, 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, 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 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,120

COMMENTS

It holds that a(n) <= A336836(n) for all n, because sigma(n) <= A003961(n) for all n (see A286385 for a proof).

The first 3 occurs at n = 19399380, the first 4 at n = 195534950863140268380. See A336389.

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.

LINKS

Antti Karttunen, Table of n, a(n) for n = 1..65537

Index entries for sequences computed from indices in prime factorization

Index entries for sequences related to sigma(n)

FORMULA

If A294934(n) = 1, a(n) = 0, otherwise a(n) = 1 + a(A003961(n)).

From Antti Karttunen, Aug 21-Sep 01 2020: (Start)

For all n >= 1,

a(A046523(n)) >= a(n).

a(A071364(n)) >= a(n).

a(A108951(n)) = A337474(n).

a(A025487(n)) = A337475(n).

(End)

EXAMPLE

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.

MATHEMATICA

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 *)

PROG

(PARI)

A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961

A336835(n) = { my(i=0); while(sigma(n) >= (n+n), i++; n = A003961(n)); (i); };

CROSSREFS

Cf. A003961, A005100, A005940, A023196, A046523, A047802, A071364, A108951, A286385, A294934, A336834, A337474, A337475.

Cf. A336389 (position of the first occurrence of a term >= n).

Differs from A294936 for the first time at n=120.

Cf. also A246271, A252459, A336836 and A336915 for similar iterations.

Sequence in context: A322358 A322437 A330174 * A037281 A143241 A308064

Adjacent sequences:  A336832 A336833 A336834 * A336836 A336837 A336838

KEYWORD

nonn

AUTHOR

Antti Karttunen, Aug 07 2020

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 17 16:50 EDT 2021. Contains 347487 sequences. (Running on oeis4.)