OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := (-1)^e*e; q[1] = True; q[n_] := Plus @@ f @@@ FactorInteger[n] == 0; Select[Range[10^4], q[#] && q[# + 1] &]
PROG
(PARI) is(n) = {my(f = factor(n)); sum(i = 1, #f~, (-1)^f[i, 2]*f[i, 2]) == 0};
is1 = is(1); for(k = 2, 10^4, is2 = is(k); if(is1 && is2, print1(k-1, ", ")); is1 = is2);
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Aug 06 2022
STATUS
approved