OFFSET
1,2
COMMENTS
Numbers where all the digits of the prime factors are either 1 or 0. Multiplicative closure of A020449.
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
FORMULA
{k | all prime factors of k are in A020449}. - Michael S. Branicky, Jun 26 2025
EXAMPLE
121 = 11 * 11 is in the sequence, since its only prime factor 11 does not have any digits larger than 1.
1001 = 7 * 13 * 11 is not in the sequence since 7 and 13 has digits larger than 1.
MATHEMATICA
Select[Range[2*10^6], AllTrue[IntegerDigits/@First/@FactorInteger[#]//Flatten, #<2&]&] (* James C. McMahon, Jun 28 2025 *)
PROG
(Python)
from sympy import primefactors
def ok(n): return all(set(str(p)) <= set("01") for p in primefactors(n))
print([k for k in range(1, 2*10**6) if ok(k)]) # Michael S. Branicky, Jun 26 2025
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Jens Ahlström, Jun 26 2025
STATUS
approved
