OFFSET
1,2
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..1000 (first 100 terms from Harry J. Smith)
PROG
(PARI) sopf(n, s, fac, i)=fac=factor(n); for(i=1, matsize(fac)[1], s=s+fac[i, 1]); return(s);
for(n=1, 10^6, if(sopf(n)==sopf(n^2-1), print(n)))
(PARI) sopf(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]); return(s) }
{ n=0; for (m=1, 10^9, if (sopf(m)==sopf(m^2 - 1), write("b064019.txt", n++, " ", m); if (n==100, break)) ) } \\ Harry J. Smith, Sep 06 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Sep 07 2001
EXTENSIONS
a(21)-a(35) from Donovan Johnson, Jun 15 2009
STATUS
approved