OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..500 from Harry J. Smith)
PROG
(PARI) sopf(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]); return(s) } sopfr(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]*f[i, 2]); return(s) } { n=0; for (m=1, 10^9, if (sopf(m)==sopfr(m + 1), write("b064678.txt", n++, " ", m); if (n==500, break)) ) } \\ Harry J. Smith, Sep 22 2009
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jason Earls, Oct 10 2001
STATUS
approved