%I #43 Jan 01 2024 08:00:40
%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,2,0,0,0,0,0,1,0,0,
%U 0,1,0,1,0,0,0,0,0,1,0,1,0,0,0,2,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 a(n) is the number of primitive nondeficient numbers (A006039) dividing n.
%C As a simple consequence of the definition of a primitive nondeficient number, a(n) is nonzero if and only if n is nondeficient.
%H Antti Karttunen, <a href="/A337690/b337690.txt">Table of n, a(n) for n = 1..65537</a>
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>.
%F a(n) = Sum_{d|n} A341619(d) = Sum_{d|n} [1==A341620(d)]. - Corrected by _Antti Karttunen_, Feb 21 2021
%F a(A005100(n)) = 0.
%F a(A006039(n)) = 1.
%F a(A023196(n)) >= 1.
%F a(A337479(n)) = A337539(n).
%F a(n) <= A341620(n). - _Antti Karttunen_, Feb 22 2021
%F Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{n>=1} 1/A006039(n) = 0.3... (see A006039 for a better estimate of this constant). - _Amiram Eldar_, Jan 01 2024
%e The least nondeficient number, therefore the least primitive nondeficient number is 6. So a(1) = a(2) = a(3) = a(4) = a(5) = 0 as all primitive nondeficient numbers are larger, and therefore not divisors; and a(6) = 1, as only 1 primitive nondeficient number divides 6, namely 6 itself.
%e 60 has the following 12 divisors: 1, 2, 3, 4, 5, 6, 10, 12, 15, 20, 30, 60. Of these, only 6 and 20 are in A006039, thus a(60) = 2.
%o (PARI)
%o A341619(n) = if(sigma(n) < (2*n), 0, fordiv(n, d, if((d<n)&&(sigma(d) >= 2*d), return(0))); (1)); \\ After code in A071395
%o A337690(n) = sumdiv(n, d, A341619(d));
%Y A006039 (or equivalently, its characteristic function, A341619) is used to define this sequence.
%Y See A000203 and A023196 for definitions of deficient and nondeficient.
%Y Sequences with similar definitions: A080224, A294927, A337539, A341620.
%Y Cf. A302110, A341604.
%Y Positions of 0's: A005100.
%Y Positions of numbers >= k: A023196 (k=1), A337688 (k=2), A337689 (k=3).
%Y Positions of first appearances are given in A337691.
%Y Differs from its derived sequence A341618 for the first time at n=120, where a(120)=2, while A341618(120)=1.
%K nonn
%O 1,60
%A _Antti Karttunen_ and _Peter Munn_, Sep 15 2020
%E Data section extended to 120 terms by _Antti Karttunen_, Feb 21 2021