login
Numbers n such that phi(n) = phi(n+7), with Euler's totient function phi = A000010.
7

%I #35 Sep 08 2022 08:45:54

%S 5,7,21,45,75,105,285,488,585,765,1148,1275,1358,1785,2528,3465,4088,

%T 6825,9405,12375,14348,15345,16208,16988,23648,25905,25935,42698,

%U 50018,52845,54615,61448,62865,68445,78195,80025,82005,88328,93555,98475

%N Numbers n such that phi(n) = phi(n+7), with Euler's totient function phi = A000010.

%C There are 40 terms below 10^5, 81 terms below 10^6 and 162 terms below 10^7. There are 6606 terms below 10^12. [_Jud McCranie_, Feb 13 2012]

%C _Farideh Firoozbakht_ asks whether there is some a(n+1) = a(n)+7, cf. link.

%C For n < 10^13, the only n such that phi(n-7) = phi(n) = phi(n+7) is 30057431145. - _Giovanni Resta_, Feb 27 2014

%H Jud McCranie, <a href="/A179189/b179189.txt">Table of n, a(n) for n = 1..6606</a> (terms < 10^12)

%H F. Firoozbakht, <a href="http://www.primepuzzles.net/puzzles/puzz_466.htm">Puzzle 466. phi(n-1)=phi(n)=phi(n+1)</a>, in C. Rivera's Primepuzzles.

%H Kevin Ford, <a href="https://arxiv.org/abs/2002.12155">Solutions of phi(n) = phi(n+k) and sigma(n) = sigma(n + k)</a>, arXiv:2002.12155 [math.NT], 2020.

%F A000010(a(n)) = A000010(a(n)+7).

%t Select[Range[100000], EulerPhi[#] == EulerPhi[# + 7] &] (* _Vincenzo Librandi_, Sep 08 2016 *)

%o (PARI) {op=vector(N=7); for( n=1, 1e5, if( op[n%N+1]+0==op[n%N+1]=eulerphi(n), print1(n-N, ", ")))}

%o (Magma) [n: n in [1..100000] | EulerPhi(n) eq EulerPhi(n+7)]; // _Vincenzo Librandi_, Sep 08 2016

%Y Cf. A000010, A001274, A001494, A179186, A179187, A179188, A007015.

%K nonn

%O 1,1

%A _M. F. Hasler_, Jan 05 2011