A270778 Primes p such that sigma(p-1) - phi(p-1) = (3p-5)/2.
3, 5, 11, 17, 257, 65537, 119831
Offset: 1
Examples
17 is a term because sigma(16)-phi(16) = 31-8 = 23 = (3*17-5)/2.
Programs
-
Magma
[n: n in[1..10^7] | IsPrime(n) and 2*(SumOfDivisors(n-1) - EulerPhi(n-1)) eq 3*n-5]
-
Mathematica
Select[Prime@ Range[10^6], DivisorSigma[1, # - 1] - EulerPhi[# - 1] == (3 # - 5)/2 &] (* Michael De Vlieger, Mar 23 2016 *)
-
PARI
lista(nn) = forprime(p=2, nn, if (sigma(p-1) - eulerphi(p-1) == (3*p-5)/2, print1(p, ", "))); \\ Michel Marcus, Mar 23 2016
-
PARI
is(n)=my(f=factor(n-1)); sigma(f) - eulerphi(f) == (3*n-5)/2 && isprime(n) \\ Charles R Greathouse IV, Apr 01 2016
Comments