OFFSET
1,1
COMMENTS
EXAMPLE
23 is in the sequence because sigma(24) - sigma(22) = 60 - 36 = 24 = sigma(23).
PROG
(Magma) [n: n in [1..1000000] | IsPrime(n) and SumOfDivisors(n) eq ((SumOfDivisors(n+1)) - (SumOfDivisors(n-1)))];
(Magma) [n: n in [A076530(n)] | IsPrime(n)];
(PARI) is_ok(index)=my(p=prime(index)); p+1==sigma(p+1)-sigma(p-1);
main(size)=my(v=vector(size), index=1); for(i=1, size, while(!is_ok(index), index++); v[i]=prime(index); index++); v \\ Anders Hellström, Jul 14 2015
(PARI) has(p)=p+1==sigma(p+1)-sigma(p-1)
select(has, primes(1000)) \\ Charles R Greathouse IV, Jul 22 2015
CROSSREFS
KEYWORD
nonn,more,hard
AUTHOR
Jaroslav Krizek, Jul 14 2015
STATUS
approved
