login
Values never taken by phi(j)/2 for any j: a(n) = A005277(n)/2.
12

%I #21 Nov 22 2015 13:34:59

%S 7,13,17,19,25,31,34,37,38,43,45,47,49,57,59,61,62,67,71,73,76,77,79,

%T 85,87,91,93,94,97,101,103,107,109,115,117,118,121,122,123,124,127,

%U 129,133,137,139,142,143,145,149,151,152,154,157,159

%N Values never taken by phi(j)/2 for any j: a(n) = A005277(n)/2.

%C Because the degree of the minimal polynomial of cos(2*Pi/k) is phi(k)/2, the degree can never be a number in this sequence. - _Artur Jasinski_, Feb 23 2011

%H Charles R Greathouse IV, <a href="/A079695/b079695.txt">Table of n, a(n) for n = 1..10000</a>

%e A005277(1)=14, therefore a(1)=7.

%t phiQ[m_] := Select[Range[m + 1, 2 m*Product[(1 - 1/(k*Log[k]))^(-1), {k, 2, DivisorSigma[0, m]}]], EulerPhi[#] == m &, 1] != {}; t = Select[Range[2, 320], phiQ]/2; Select[Range@ Max@ t, !MemberQ[t, #] &] (* _Michael De Vlieger_, Mar 22 2015, after _Jean-François Alcover_ at A002180 *)

%o (PARI) is(n)=!istotient(2*n) \\ _Charles R Greathouse IV_, Mar 23 2015

%o (Haskell)

%o import Data.List.Ordered (minus)

%o a079695 n = a079695_list !! (n-1)

%o a079695_list = [1..] `minus` a002180_list

%o -- _Reinhard Zumkeller_, Nov 22 2015

%Y Cf. A005277 (nontotients), A002180 (complementary sequence).

%K nonn

%O 1,1

%A _Jon Perry_, Jan 31 2003