OFFSET
1,2
COMMENTS
Multiplicative closure of A036953.
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
FORMULA
{k | all prime factors of k are in A036953}. - Michael S. Branicky, Jun 26 2025
EXAMPLE
202 is in the sequence since the prime factors 2 and 101 both have all digits smaller than 3.
34 is not in the sequence since it has the prime factor 17 that have a digit larger than 2.
MATHEMATICA
A385344Q[k_] := AllTrue[FactorInteger[k][[All, 1]], Max[IntegerDigits[#]] < 3 &];
Select[Range[10000], A385344Q] (* Paolo Xausa, Jun 28 2025 *)
PROG
(Python)
from sympy import primefactors
def ok(n): return all(set(str(f)) <= set("012") for f in primefactors(n))
print([k for k in range(1, 6000) if ok(k)]) # Michael S. Branicky, Jun 26 2025
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Jens Ahlström, Jun 26 2025
STATUS
approved
