OFFSET
1,1
COMMENTS
For the purpose here, if a number has repeated prime factors, those are written repeatedly. For example, the factorization of 27 is expressed as (3, 3, 3) rather than (3^3). - Alonso del Arte, Jul 05 2011
LINKS
Klaus Brockhaus, Table of n, a(n) for n = 1..10000
EXAMPLE
Since the prime factorization of 95 is (5, 19), and both 9 and 5 occur in (5, 19), the number 95 is on the list.
Since the prime factorization of 1255 is (5, 251), and 1, 2, and both 5s occur in (5, 251), the number 1255 is on the list.
22 is not on the list because its prime factorization is (2, 11) and that does not have enough 2s. Nor is 25 on the list because for this sequence we express its factorization as (5, 5) rather than (5^2).
MATHEMATICA
Select[Range[2, 5000], Not[PrimeQ[#]] && Sort[DigitCount[FromDigits[Flatten[IntegerDigits/@Flatten[Table[#1, {#2}]&@@@FactorInteger[#]]]]] - DigitCount[#]][[1]] >= 0 &] (* Alonso del Arte, Jun 28 2011, based on HomePrimeStep function by Eric W. Weisstein *)
PROG
(Magma) S:=[]; for n in [1..10000] do if not IsPrime(n) then u:=Intseq(n); f:=Factorization(n); v:=&cat[ [ f[j, 1]: i in [1..f[j, 2]] ]: j in [1..#f] ]; w:=&cat[ Intseq(p): p in v ]; if forall{ a: a in [0..9] | Multiplicity(SequenceToMultiset(u), a) le Multiplicity(SequenceToMultiset(w), a) } then Append(~S, n); end if; end if; end for; S; // Klaus Brockhaus, Jul 09 2011
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Gil Broussard, Jun 24 2011
STATUS
approved