OFFSET
1,1
COMMENTS
If k is in the sequence then so is k*m for positive m. - David A. Corneth, Aug 19 2022
Numbers that are divisible by at least one of 7, 8, 9, 10, 12, 15, 22, 33, 39, 52, 55, 68, 102, 114, 138. For proof, see link. - Robert Israel, Sep 02 2022
The asymptotic density of this sequence is 17819629/37182145 = 0.479252... . - Amiram Eldar, Aug 08 2023
LINKS
EXAMPLE
10 is in the sequence since 10 = 2+2+2+1+1+1+1, where each summand divides 10.
MAPLE
filter:= n -> ormap(t -> n mod t = 0, [7, 8, 9, 10, 12, 15, 22, 33, 39, 52, 55, 68, 102, 114, 138]):
select(filter, [$1..200]); # Robert Israel, Sep 02 2022
MATHEMATICA
q[n_, k_] := AnyTrue[Tuples[Divisors[n], k], Total[#] == n &]; Select[Range[130], q[#, 7] &] (* Amiram Eldar, Aug 19 2022 *)
PROG
(PARI) isok(k) = my(d=divisors(k)); forpart(p=k, if (setintersect(d, Set(p)) == Set(p), return(1)), , [7, 7]); \\ Michel Marcus, Aug 19 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Aug 18 2022
STATUS
approved