OFFSET
1,1
COMMENTS
Subsequence of A001358.
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000
EXAMPLE
1469 = 13*113 is in the sequence because reverse(1469)=reverse(13)*reverse(113) => 9641 = 31*311 where 31 and 311 are prime numbers.
MAPLE
with(numtheory):lst:={}:T1:=array(1..300):T2:=array(1..300):k:=0:
for n from 1 to 1000 do:
p:=ithprime(n):xp:=convert(p, base, 10):
np:=nops(xp):sp:=sum('xp[np-i+1]*10^(i-1)', 'i'=1..np):
if type(sp, prime)=true
then
k:=k+1:T1[k]:=p:T2[k]:=sp:
else
fi:
od:
for i from 1 to k do:
for j from i to k do:
x:=T1[i]*T1[j]:y:=convert(x, base, 10):n2:=nops(y):
s:=sum('y[n2-i+1]*10^(i-1)', 'i'=1..n2):
if T2[i]*T2[j]=s
then
lst:=lst union {x}:
else
fi:
od:
od:
print(lst):
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Mar 15 2014
STATUS
approved