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