OFFSET
1,1
EXAMPLE
a(3) = 13 is a term because 13, the only prime dividing 13, also divides A114707(12) = 13.
MAPLE
R:= 1: v:= 1:
for n from 2 to 10^6 do
v:= v + nops(select(p -> v mod p <> 0, numtheory:-factorset(n)));
R:= R, v;
od:
select(t -> R[t]=R[t-1], [$2..10^6]);
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Jun 27 2024
STATUS
approved