|
%I #24 Sep 08 2022 08:46:16
%S 3,5,9,11,17,33,65,129,231,257,513,1025,2049,4097,8193,16385,32769,
%T 65537,119831,131073,262145,524289,1048577,2097153,4194305,8388609,
%U 16777217,33554433,67108865,134217729,268435457,536870913,1073741825,2147483649,4294967297
%N Numbers n such that sigma(n-1) - phi(n-1) = (3n-5)/2.
%C Numbers n such that A051612(n-1) = (3n-5)/2.
%C Numbers of the form 2^n + 1 for n >= 1 from A000051 are terms.
%C Prime terms are in A270778.
%H Giovanni Resta, <a href="/A270836/b270836.txt">Table of n, a(n) for n = 1..46</a> (terms < 10^13)
%e 17 is a term because sigma(16) - phi(16) = 31-8 = 23 = (3*17-5)/2.
%t Select[Range[10^6], 2 (DivisorSigma[1, # - 1] - EulerPhi[# - 1]) == 3 # - 5 &] (* _Michael De Vlieger_, Mar 24 2016 *)
%o (Magma) [n: n in[1..10^7] | 2*(SumOfDivisors(n-1) - EulerPhi(n-1)) eq 3*n-5]
%o (PARI) lista(nn) = {for(n=2, nn, if(sigma(n-1) - eulerphi(n-1) == (3*n-5)/2, print1(n, ", "))); } \\ _Altug Alkan_, Mar 23 2016
%Y Cf. A000010, A000051, A000203, A051612, A270778, A270837.
%K nonn
%O 1,1
%A _Jaroslav Krizek_, Mar 23 2016
%E a(30)-a(32) from _Michel Marcus_, Apr 05 2016
%E a(33)-a(35) from _Giovanni Resta_, Apr 11 2016
|
