OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..375
PROG
(PARI) sumprime(n, s, fac, i) = fac=factor(n); for(i=1, matsize(fac)[1], s=s+fac[i, 1]); return(s); z(n)=sumprime(n-numdiv(n)); d(n)=numdiv(sumprime(n)+n); for(n=1, 10^6, if(d(n)==z(n), print(n)))
(PARI) sumprime(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]); return(s) } z(n)= { sumprime(n - numdiv(n)) } d(n)= { numdiv(sumprime(n) + n) } { n=0; for (m=1, 10^9, if(d(m)==z(m), write("b063961.txt", n++, " ", m); if (n==375, break)) ) } \\ Harry J. Smith, Sep 04 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Sep 04 2001
STATUS
approved