OFFSET
1,3
COMMENTS
EXAMPLE
486486 is in the sequence because the prime divisors of 486486, 864864 and 648648 are 2,3,7,11 and 13.
MAPLE
with(numtheory): T:=array(1..10):U:=array(1..10):
for n from 11 to 10^6 do:
c:=1:x:=convert(n, base, 10):n1:=nops(x):si:=factorset(n):
for i from 1 to n1 do:T[i]:=x[n1-i+1]:od:
for j from 1 to n1-1 do:
for k from 1 to n1-1 do:
U[k]:=T[k+1]:
od:
U[n1]:=T[1]:s:=sum('U[n1-p+1]*10^(p-1)', 'p'=1..n1):
if factorset(s)=si
then
c:=c+1:
else
fi:
for l from 1 to n1 do:
T[l]:=U[l]:
od:
if c=n1
then
printf(`%d, `, n):
else
fi:
od:
od:
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Jun 29 2014
EXTENSIONS
Single-digit numbers added by N. J. A. Sloane, Jul 23 2014
STATUS
approved