OFFSET
1,2
COMMENTS
All primes containing 1 are members.
Sequence is a subsequence of A011531. The first nonprime terms of the sequence are 1, 121, 125, 1207, 1255, 1379, 10201, 10379, 11009, 11209, 12419, 12709, 12755, ... - R. J. Mathar, Jul 26 2007
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
131 is a term. 143 is not a term as the divisor 11 contains two 1's.
MAPLE
isA092911 := proc(n) local digs, digsleft, divs, d, i, j ; digs := convert(n, base, 10) ; divs := numtheory[divisors](n) ; for i from 1 to nops(divs) do digsleft := digs ; d := convert(op(i, divs), base, 10) ; for j in d do if member(j, digsleft, 'jposit') then digsleft := subsop(jposit=NULL, digsleft) ; else RETURN(false) ; fi ; od ; od ; RETURN(true) ; end: for n from 1 to 600 do if isA092911(n) then printf("%d, ", n) ; fi ; od ; # R. J. Mathar, Jul 26 2007
MATHEMATICA
subQ[s1_, s2_] := AllTrue[Count[s1, #] & /@ (First /@ (t = Tally[s2])) - Last /@ t, # >= 0 &]; digQ[n1_, n2_] := subQ[IntegerDigits[n1], IntegerDigits[n2]]; seqQ[n_] := AllTrue[Most@Divisors[n], digQ[n, #] &]; Select[Range[600], seqQ] (* Amiram Eldar, Nov 12 2020 *)
PROG
(Python)
from sympy import divisors
from collections import Counter
def ok(n):
ncounts = Counter(str(n))
for d in divisors(n)[:-1]:
divcounts = Counter(str(d))
if any(ncounts[c] < divcounts[c] for c in divcounts): return False
return True
print(list(filter(ok, range(1, 630)))) # Michael S. Branicky, May 08 2021
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Amarnath Murthy, Mar 14 2004
EXTENSIONS
Corrected and extended by R. J. Mathar, Jul 26 2007
STATUS
approved