OFFSET
1,2
COMMENTS
3^m*10^k for k, m > 0 are terms of this sequence. - Chai Wah Wu, Jun 23 2020
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000
Wikipedia, Palindrome
EXAMPLE
The sequence of rows of A296150 indexed by the terms of this sequence begins: (1), (11), (21), (111), (1111), (321), (11111), (2211), (111111), (3221), (1111111), (4321), (222111), (11111111), (32221), (33211), (111111111), (322221), (332211).
MATHEMATICA
nrmpalQ[n_]:=With[{f=If[n==1, {}, FactorInteger[n]]}, And[PrimePi/@ Sort[First/@f] == Range[ Length[f]], Reverse[Last/@f] == Last/@f]]; Select[Range[100], nrmpalQ]
upto = 10^20; pL[n_] := Block[{p = Prime@Range@n, h = Ceiling[n/2]}, Take[p, h] Reverse@ If[n == 2 h, Take[p, -h], Prepend[ Take[p, 1-h], 1]]]; ric[v_, p_] := If[p == {}, AppendTo[L, v], Block[{w = v}, While[w <= upto, ric[w, Rest@ p]; w *= First@ p]]]; np = 1; L = {1}; While[(b = Times @@ Prime[Range@ np]) <= upto, ric[b, pL[np++]]]; Sort[L] (* Giovanni Resta, Jun 23 2020 *)
PROG
(Python)
from sympy import factorint, primepi
A317087_list = [1]
for n in range(1, 10**5):
d = factorint(n)
k, l = sorted(d.keys()), len(d)
if l > 0 and l == primepi(max(d)):
for i in range(l//2):
if d[k[i]] != d[k[l-i-1]]:
break
else:
A317087_list.append(n) # Chai Wah Wu, Jun 23 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jul 21 2018
STATUS
approved