|
%I #32 Mar 11 2019 08:22:13
%S 4,6,8,9,10,12,14,20,22,24,25,26,28,30,32,34,38,42,44,46,48,49,52,54,
%T 58,60,62,66,68,74,76,78,80,81,82,86,90,92,94,102,106,112,114,116,118,
%U 121,122,124,134,138,140,142,146,148,150,158,160,164,166,168,169,172,174
%N Numbers that are divisible by the number of their nontrivial divisors.
%C We may define the number of divisors of a number n in four ways:
%C (1) A070824(n) = number of nontrivial or real divisors: 1 < d < n;
%C (2) variant of A032741(n) = number of small divisors: 1 and real divisors;
%C (3) A032741(n) = number of big or proper divisors: real divisors and n;
%C (4) A000005(n) = number of all divisors of n: 1, n and real divisors.
%C The case (1), divisibility through the number of nontrivial divisors, defines this sequence.
%D T. Szimeonov, A számok [The numbers], Budapest, 2019, VVMA, 124 p.
%e 1 and the prime numbers do not have any nontrivial divisors; A070824(n) is 0 for n=1 or a prime, and so they are not terms.
%e The only nontrivial divisor of 4 is 2, so A070824(4) = 1; 4 is divisible by 1, so 4 is a term.
%e A070824(15) = 2, and 15 is not divisible by 2, so 15 is not a term.
%t seqQ[n_] := (nd = DivisorSigma[0, n] - 2) > 0 && Divisible[n, nd]; Select[Range[200], seqQ] (* _Amiram Eldar_, Mar 11 2019 *)
%o (PARI) f(n) = if (n==1, 0, numdiv(n)-2); \\ A070824
%o isok(n) = f(n) && !frac(n/f(n)); \\ _Michel Marcus_, Feb 17 2019
%Y Cf. A070824, A054010, A033950.
%K nonn
%O 1,1
%A _Todor Szimeonov_, Feb 05 2019
%E More terms from _Michel Marcus_, Feb 17 2019
|
