OFFSET
1,1
EXAMPLE
a(20) = 1227818 = 2 * 19 * 79 * 409 and
(1227818 - 20) /(2 + 20) = 55809;
(1227818 - 20) /(19 + 20) = 31482;
(1227818 - 20) /(79 + 20) = 12402;
(1227818 - 20) /(409 + 20) = 2862.
MAPLE
f:= proc(n) local k, F;
for k from n+1 do
F:= ifactors(k)[2];
if nops(F)=1 or F[.., 2] <> [1$nops(F)] then next fi;
if andmap(p -> (k - n) mod (p+n) = 0, F[.., 1]) then return k fi
od
end proc:
map(f, [$1..40]); # Robert Israel, Feb 17 2026
PROG
(PARI) isok(k, n) = if (!issquarefree(k) || isprime(k), return(0)); my(f=factor(k)[, 1]); for (i=1, #f, if ((k-n) % (f[i]+n), return(0)); ); return(1);
a(n) = my(k=n+1); while (!isok(k, n), k++); k; \\ Michel Marcus, Mar 30 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Mar 29 2025
EXTENSIONS
More terms from Michel Marcus, Mar 30 2025
STATUS
approved
