OFFSET
1,3
COMMENTS
For n > 0, 2*a(n) + 1 is a term of A020884. This is because 2*a(n) + 1 is odd and every odd number is the difference of the squares of two consecutive numbers and hence are coprime.
For n > 0, (2*a(n) + 1) * (2*a(n) + 2) is a term of A024364. This is because (2*a(n) + 1) * (2*a(n) + 2) = 2*((a(n) + 1)^2 + (a(n) + 1) * a(n)) and gcd((a(n) + 1), a(n)) = 1.
For n > 0, a(n) is congruent to 1 or 4 mod 6.
2*a(n) + 1 is congruent to 1 or 3 mod 6 and is a term of A047241.
2*a(n) + 2 is congruent to 2 or 4 mod 6 and is a term of A047235.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..5416 (calculated using the b-file at A001274)
FORMULA
a(n) = (A299535(n) - 2) / 2.
EXAMPLE
0 is a term because phi(1) = phi(2) = 1.
1 is a term because phi(3) = phi(4) = 2.
7 is a term because phi(15) = phi(16) = 8.
MATHEMATICA
Select[Range[0, 9999], EulerPhi[2# + 1] == EulerPhi[2# + 2] &] (* Alonso del Arte, Jul 05 2019 *)
Select[(#-1)/2&/@SequencePosition[EulerPhi[Range[900000]], {x_, x_}][[All, 1]], IntegerQ] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Aug 24 2019 *)
PROG
(PARI) lista(nn) = for(n=0, nn, if(eulerphi(2*n+1) == eulerphi(2*n+2), print1(n, ", ")));
lista(430000)
CROSSREFS
KEYWORD
nonn
AUTHOR
Torlach Rush, May 24 2019
STATUS
approved