OFFSET
1,2
COMMENTS
The old name was A196445(n)/2.
From Wolfdieter Lang, Nov 12 2019: (Start)
These are the odd numbers m for which the degree of the algebraic number sin(Pi/(2*m)) (the degree of its minimal polynomial), given by A055035(2*m), is odd. Because A055035(1) = 1, there is just this other instance with odd A055035.
This sequence {a(n)} consists for n >= 2 of all powers >= 1 of each prime p == 3 (mod 4) from A002145, sorted into increasing order. This follows from the factorization of odd m >= 3, and that phi(m)/2 has to be odd.
For a(n), with n >= 2, there is exactly one pair of solutions x = +1 and -1 (the trivial solution) of the congruence x^2 == +1 (mod a(n)), and there is no solution of the congruence x^2 == -1 (mod a(n)). The proof starts with showing this for p == 3 (mod 4). It can be shown that the square 1 appears only for x = 1 if x runs through 1, ..., (p-1)/2. The other x range (p+1)/2, ..., p-1, which has the same squares (mod p), can, by reading backwards, be interpreted as the -x partners of x from the first range. The Legendre symbol (-1, p) = -1 shows the second claim. Then one applies the lifting theorem for powers of primes (see Apostol, Theorem 5.30, p. 121), where only its part (a) is needed, and the step by step lifting to each prime power is unique.
For a(1) = 1 there is just one solution x = 0 of the congruence x^2 == +1 (mod 1), and also of x^2 == -1 (mod 1).
The complementary sequence with odd m >= 3 and even phi(m)/2 is given in A327922.
(End)
REFERENCES
T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, pp. 121-122.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(1) = 1, and for n >= 2 the set of all positive powers of each prime p == 3 (mod 4) (A002145), sorted increasingly.
EXAMPLE
Factorization for n >= 2: 3, 7, 3^2, 11, 19, 23, 3^3, 31, 43, 47, 7^2, 59, 67, 71, 79, 3^4, 83, 103, 107, ... - Wolfdieter Lang, Nov 12 2019
MATHEMATICA
a[n_] := If[n == 2, 1, EulerPhi[n]/{1, 1, 2, 1}[[Mod[n, 4] + 1]]]; aa = {}; Do[If[OddQ[a[n]], AppendTo[aa, n/2]], {n, 2, 1000}]; aa
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Artur Jasinski, Oct 15 2011
EXTENSIONS
Name changed by Wolfdieter Lang, Nov 12 2019
Name edited by Jon E. Schoenfield, Jan 21 2023
STATUS
approved