A260071 Primes p such that sigma(p) = sigma(p+1) - sigma(p-1).
2, 3, 23, 970388922263, 991817878343, 1677028870823
Offset: 1
Examples
23 is in the sequence because sigma(24) - sigma(22) = 60 - 36 = 24 = sigma(23).
Programs
-
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
Comments