OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..1000
EXAMPLE
n=100: sigma(100) - 2*100 = 17, a prime.
PROG
(PARI) j=[]; for(n=1, 50000, x=sigma(n)-2*n; if(x>0 && isprime(x), j=concat(j, n))); j
(PARI) { n=0; for (m=1, 10^9, if (isprime(sigma(m) - 2*m), write("b064271.txt", n++, " ", m); if (n==1000, break)) ) } \\ Harry J. Smith, Sep 11 2009
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jason Earls, Sep 23 2001
STATUS
approved