OFFSET
1,1
COMMENTS
Compare to A095050: Numbers such that all ten digits are needed to write all positive divisors in decimal representation.
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
EXAMPLE
The divisors of 918 are 1, 2, 3, 6, 9, 17, 18, 27, 34, 51, 54, 102, 153, 306, 459, and 918. Every digit appears at least twice. Thus, 918 is in this sequence.
MAPLE
q:= n-> (p-> is(min(seq(coeff(p, x, j), j=0..9))>1))(add(x^i, i=
map(d-> convert(d, base, 10)[], [numtheory[divisors](n)[]]))):
select(q, [$10..2600])[]; # Alois P. Heinz, Apr 21 2022
MATHEMATICA
Select[Range[3000], Length[Transpose[Tally[Flatten[IntegerDigits[Divisors[#]]]]][[2]]] == 10 && Min[Transpose[Tally[Flatten[IntegerDigits[Divisors[#]]]]][[2]]] > 1 &]
PROG
(Python)
from sympy import divisors
def ok(n):
digits_used = {d:0 for d in "0123456789"}
for div in divisors(n, generator=True):
for d in str(div): digits_used[d] += 1
if all(digits_used[d] > 1 for d in "0123456789"): return True
return False
print([k for k in range(2521) if ok(k)]) # Michael S. Branicky, Jun 25 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Tanya Khovanova, Jun 17 2021
STATUS
approved