OFFSET
1,1
COMMENTS
Contains 6*p*q if p and q are distinct primes > 3. The first term not of this form is 770. - Robert Israel, Jan 09 2025
a(43) = 2002 is the only term coprime to 15. - Charles R Greathouse IV, Jan 13 2025
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
a(n) == 2 (mod 4).
a(n) ~ (1/6)*n log n/log log n. - Charles R Greathouse IV, Jan 13 2025
EXAMPLE
210 is a term because 210=2*3*5*7 is the product of four distinct primes and it is smaller than the sum of its proper divisors 366.
1155 is not a term because 1155=3*5*7*11 is the product of four distinct primes and it is larger than the sum of its proper divisors 1149.
MAPLE
filter:= proc(n) local F, t;
F:= ifactors(n)[2];
F[.., 2] = [1, 1, 1, 1] and mul(t[1]+1, t = F) > 2*n
end proc:
select(filter, [seq(i, i=2..3000, 4)]); # Robert Israel, Jan 09 2025
MATHEMATICA
q[n_] := Module[{f = FactorInteger[n]}, f[[;; , 2]] == {1, 1, 1, 1} && Times @@ (1 + 1/f[[;; , 1]]) > 2]; Select[Range[3000], q] (* Amiram Eldar, Jan 09 2025 *)
PROG
(PARI) list(lim)=my(v=List(select(k->k<=lim, [1430, 1870, 2002, 2090, 2210, 2470, 2530, 2990, 3190, 3230, 3410, 3770, 4030, 4070, 4510, 4730, 5170, 5830]))); forprime(p=5, sqrtint(lim\6), my(t=6*p); forprime(q=p+2, lim\t, listput(v, t*q))); forprime(p=11, lim\70, listput(v, 70*p)); Set(v) \\ Charles R Greathouse IV, Jan 13 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Massimo Kofler, Jan 09 2025
STATUS
approved