OFFSET
1,1
LINKS
Harry J. Smith and Donovan Johnson, Table of n, a(n) for n = 1..1000 (first 80 terms from Harry J. Smith)
MATHEMATICA
snpd[n_]:=Module[{fi=FactorInteger[n]}, {Total[Flatten[Table[#[[1]], {#[[2]]}]&/@fi]], Total[Transpose[fi][[2]]]}]; Flatten[Position[Partition[ Array[ snpd, 5600000], 2, 1], _?(#[[1]]==#[[2]]&), {1}, Heads->False]] (* Harvey P. Dale, Nov 29 2013 *)
PROG
(PARI) sopfr(n)= { local(c, f, s=0); G=0; f=factor(n); for(i=1, matsize(f)[1], c=f[i, 2]; G+=c; s+=f[i, 1]*c); return(s) } { n=c=s=0; for (m=1, 10^9, us=sopfr(m + 1); if(c==G && s==us, write("b064351.txt", n++, " ", m); if (n==100, break)); c=G; s=us ) } \\ Harry J. Smith, Sep 12 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Oct 15 2001
EXTENSIONS
a(25) - a(27) from Harry J. Smith, Sep 12 2009
STATUS
approved