login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A341624
a(n) = 0 if n is a deficient number, otherwise a(n) is the number of nondeficient divisors of the last number in the iteration x -> A003961(x) (starting from x=n) for which that count (A341620) is nonzero.
4
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 2, 0, 1, 0, 0, 0, 3, 0, 0, 0, 1, 0, 2, 0, 0, 0, 0, 0, 4, 0, 0, 0, 2, 0, 2, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 3, 0, 2, 0, 0, 0, 5, 0, 0, 0, 0, 0, 2, 0, 0, 0, 1, 0, 6, 0, 0, 0, 0, 0, 2, 0, 3, 0, 0, 0, 5, 0, 0, 0, 1, 0, 4, 0, 0, 0, 0, 0, 5, 0, 0, 0, 2, 0, 2, 0, 1, 0, 0, 0, 6, 0, 0, 0, 3, 0, 2, 0, 0, 0, 0, 0, 1
OFFSET
1,12
PROG
(PARI)
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A341620(n) = sumdiv(n, d, (sigma(d)>=(2*d)));
A341624(n) = { my(t, u=0); while((t=A341620(n))>0, u=t; n = A003961(n)); (u); };
CROSSREFS
Cf. A005100 (positions of zeros).
Differs from A341620 for the first time at n=120, where a(120)=1, while A341620(120)=8.
Cf. also A341508, A341618.
Sequence in context: A143840 A028649 A305558 * A341620 A097798 A065205
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 22 2021
STATUS
approved