OFFSET
1,2
COMMENTS
Inspired by A328095. To avoid all primes being in the sequence the divisors of k includes k itself.
Contains 10^k, 5*10^k and 6*10^k for all k, 3*10^k, 4*10^k, 7*10^k and 9*10^k for all odd k. - Robert Israel, Nov 11 2019
LINKS
Scott R. Shannon, Table of n, a(n) for n = 1..1000. Note when searching for these numbers one needs to use arbitrary precision packages; the product for 24570000 has 1486 digits.
EXAMPLE
16 is in the sequence as the divisors of 16 are 1,2,4,8,16, and 16*(1*2*4*8*16) = 16*1024 = 16384, and '16384' contains '16' as a substring.
30 is in the sequence as the divisors of 30 are 1,2,3,5,6,10,15,30, and 30*(1*2*3*5*6*10*15*30) = 30*810000 = 24300000, and '24300000' contains '30' as a substring.
MATHEMATICA
f[n_] := n^(1+DivisorSigma[0, n]/2); aQ[n_] := SequenceCount[IntegerDigits[f[n]], IntegerDigits[n]] > 0; Select[Range[3000], aQ] (* Amiram Eldar, Nov 10 2019 *)
PROG
(Magma) a:=[]; for k in [1..3000] do t:=IntegerToString(k*(&*Divisors(k))); s:=IntegerToString(k); if s in t then Append(~a, k); end if; end for; a; // Marius A. Burtea, Nov 10 2019
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Scott R. Shannon, Nov 10 2019
STATUS
approved