OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
MATHEMATICA
f[p_, e_] := If[e==1, p+1, 1]; s[n_] := Times @@ f @@@ FactorInteger[n]; s1=0; seq = {}; Do[s2 = s[n]; If[s1 == s2, AppendTo[seq, n-1]]; s1 = s2, {n, 1, 10000}]; seq
PROG
(PARI) s(n)={sumdiv(n, d, d*issquarefree(d)*(gcd(d, n/d) == 1))}
{ for(k=1, 10^6, if(s(k)==s(k+1), print1(k, ", "))) } \\ Andrew Howroyd, Sep 28 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Sep 28 2019
STATUS
approved