OFFSET
1,1
COMMENTS
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
EXAMPLE
a(10) = 14641 = 11^4 (the smallest odd value in this sequence).
a(11) = 21312 = 2^6 * 3^2 * 37.
MAPLE
isA046101 := proc(n) local ifs, f ; ifs := ifactors(n)[2] ; for f in ifs do if op(2, f) >= 4 then RETURN(true) ; fi ; od: RETURN(false) ; end: isA002113 := proc(n) local digs, i ; digs := convert(n, base, 10) ; for i from 1 to nops(digs) do if op(i, digs) <> op(-i, digs) then RETURN(false) ; fi ; od: RETURN(true) ; end: isA133514 := proc(n) isA046101(n) and isA002113(n) ; end: for n from 1 to 100000 do if isA133514(n) then printf("%d, ", n) ; fi ; od: # R. J. Mathar, Jan 12 2008
# second Maple program:
q:= n->StringTools[IsPalindrome](""||n) and max(map(i->i[2], ifactors(n)[2]))>3:
select(q, [$1..70000])[]; # Alois P. Heinz, Sep 27 2023
MATHEMATICA
a = {}; For[n = 2, n < 100000, n++, If[FromDigits[Reverse[IntegerDigits[n]]] == n, b = 0; For[l = 1, l < Length[FactorInteger[n]] + 1, l++, If[FactorInteger[n][[l, 2]] > 3, b = 1]]; If[b == 1, AppendTo[a, n]]]]; a (* Stefan Steinerberger, Dec 26 2007 *)
Select[Range@100000, PalindromeQ@#&&3<Max@Last@Transpose@FactorInteger@#&] (* Hans Rudolf Widmer, Sep 27 2023 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Jonathan Vos Post, Nov 30 2007
EXTENSIONS
More terms from Stefan Steinerberger, Dec 26 2007
More terms from R. J. Mathar, Jan 12 2008
STATUS
approved