|
%I #12 Jul 29 2024 06:18:02
%S 1,7,34,115,298,677,1431,2853,5401,9813,17289,29687,49901,82315,
%T 133617,213801,337742,527231,814218,1244962,1886069,2833392,4223121,
%U 6247850,9179762,13399951,19440903,28043882,40233399,57423953,81557903
%N a(0) = 1; a(n) = half of the a(n-1)-th even nontotient number.
%e a(5) = 677 = 1354/2, where 1354 is the a(4) = 298th term in A005277. [Corrected by _Sean A. Irvine_, Jul 28 2024]
%t a = Table[0, {10^6}]; Do[b = EulerPhi[n]/2; If[b < 10^6 + 1, a[[b]] = 1], {n, 1, 2 10^7}]; b = Select[ Range[10^6], a[[ # ]] == 0 &]; c = 1; Do[ Print[c]; c = b[[c]], {n, 0, 18}]
%Y Cf. A005277, A072416.
%K nonn
%O 0,2
%A _Robert G. Wilson v_, Jun 20 2002
%E a(19)-a(30) from _Donovan Johnson_, May 30 2010
|
