OFFSET
1,1
COMMENTS
Most terms end in 0, since 2*5*prime will work if prime contains 2 and 5 as substrings. The other terms are listed in A059402.
There must be at least two maximal prime powers dividing each term. - Harvey P. Dale, Dec 04 2016
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
1197 = 9 * 7 * 19 and all of these are substrings.
MATHEMATICA
psmppQ[n_]:=Module[{pp=#[[1]]^#[[2]]&/@FactorInteger[n], idn= IntegerDigits[ n]}, Length[pp]>1&&And@@Table[ SequenceCount[ idn, IntegerDigits[pp[[i]]]]>0, {i, Length[pp]}]]; Select[Range[ 33000], psmppQ] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Dec 04 2016 *)
PROG
(Haskell)
import Data.List (isInfixOf)
a059401 n = a059401_list !! (n-1)
a059401_list = filter (\x -> a010055 x == 0 &&
all (`isInfixOf` show x) (map show $ a141809_row x)) [1..]
-- Reinhard Zumkeller, Dec 16 2013
CROSSREFS
KEYWORD
base,nice,nonn
AUTHOR
Erich Friedman, Jan 29 2001
EXTENSIONS
Offset corrected by Reinhard Zumkeller, Dec 16 2013
Edited by Peter Munn, Sep 01 2022
STATUS
approved