login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A064439 Numbers n such that sigma(n) - n - 1 = pi(n). 0
4, 55, 65, 95, 125, 145, 155, 185, 205, 2779, 2863, 55297, 174691, 174779, 487903, 1301989, 1302457, 5254751, 6383483, 23140961, 48267437, 59651051, 70111213, 70111247, 92514491, 199445641, 212210443, 514269523, 514269599, 21881358361, 1602278990111 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

a(32) > 3*10^12. - Giovanni Resta, Mar 31 2017

LINKS

Table of n, a(n) for n=1..31.

MATHEMATICA

Select[Range[10^5], DivisorSigma[1, #] - # - 1 == PrimePi[#] &] (* Giovanni Resta, Mar 31 2017 *)

PROG

(PARI) sig(n) = sigma(n)-n-1; pi(x, c=0) = forprime(p=2, x, c++); c for(n=1, 10^8, if(sig(n)==pi(n), print(n)))

(PARI) { n=0; for (m=1, 10^9, if ((sigma(m) - m - 1)==primepi(m), write("b064439.txt", n++, " ", m); if (n==100, break)) ) } \\ Harry J. Smith, Sep 14 2009

CROSSREFS

Cf. A000720, A048050.

Sequence in context: A048371 A077658 A217124 * A133218 A190441 A151576

Adjacent sequences:  A064436 A064437 A064438 * A064440 A064441 A064442

KEYWORD

nonn

AUTHOR

Jason Earls, Oct 01 2001

EXTENSIONS

More terms from Klaus Brockhaus, Oct 05 2001. No further term < 800000.

a(16)-a(17) from Harry J. Smith, Sep 14 2009

a(18)-a(29) from Donovan Johnson, Feb 09 2013

a(30)-a(31) from Giovanni Resta, Mar 31 2017

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 20 01:17 EDT 2019. Contains 325168 sequences. (Running on oeis4.)