OFFSET
1,1
COMMENTS
Let x(0)x(1)...x(q-1)x(q) denote the decimal expansion of a number n. The sequence lists the numbers n such that the prefix p = x(0)x(1)...x(q-1) and the suffix x(1)...x(q-1)x(q) have the same prime distinct divisors.
EXAMPLE
3748 is in the sequence because 374 and 748 have the same prime divisors: {2,11,17}.
MAPLE
with(numtheory):
for n from 10 to 10000 do:
x:=convert(n, base, 10):n1:=nops(x):
s1 := n mod 10^ilog10(n):
s2:=(n-irem(n, 10))/10:
x1:=factorset(s1):x2:=factorset(s2):
if x1 = x2 and x1 <>{}
then
printf(`%d, `, n):
else
fi:
od:
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Jun 27 2014
STATUS
approved