login
Complement of A138929: positive integers not of the form 2*p^k, k >= 0, p a prime (also 2).
5

%I #18 Aug 30 2024 16:32:37

%S 1,3,5,7,9,11,12,13,15,17,19,20,21,23,24,25,27,28,29,30,31,33,35,36,

%T 37,39,40,41,42,43,44,45,47,48,49,51,52,53,55,56,57,59,60,61,63,65,66,

%U 67,68,69,70,71,72,73,75,76,77,78,79,80,81,83,84,85,87,88,89,90,91,92,93,95,96,97,99,100

%N Complement of A138929: positive integers not of the form 2*p^k, k >= 0, p a prime (also 2).

%C The complement relative to the positive integers of the present sequence is A138929.

%C The sequence includes all positive integers of the forms (i) odd, (ii) 2^k*p, p an odd prime and k>=2, and (iii) 2^e0*p1^e1*p2^e2 *** pk^ek, k >= 2, with odd primes p1, ..., pk, and each exponent from {e0, ..., ek} is >= 1.

%C For a(n) > 1 a regular a(n)-gon, with length ratio (smallest diagonal)/side rho(a(n)) = 2*cos(Pi/a(n)), the inverse of rho(a(n)), which is an element of the algebraic number field Q(rho(a(n))), is in fact a Q(rho(a(n)))-integer. For a(1)=1 rho(1) = -2, and the inverse is not a Q-integer.

%H Paolo Xausa, <a href="/A230078/b230078.txt">Table of n, a(n) for n = 1..10000</a>

%e Even members a(n) of the form (ii) 2^k*p, p an odd prime and k>=2 are: 12, 20, 24, 28, 36, 40, 44, 48, 52, 56, 68, 72, 76, 80, 88, 92, 96, 100,...

%e Even members a(n) of the form (iii), given above, include 30, 42, 60, 66, 70, 78, 84, 90, ...

%e For the regular 5-gon (pentagon) rho(5) = tau = (1 + sqrt(5))/ 2 (the golden section). The number field is Q(rho(5)), and for the inverse one has 1/rho(5) = -1*1 + 1* rho(5) (in the power basis <1, rho(5)>, in which Q(rho(5))-integers have integer coefficients).

%e For the regular 7-gon rho(7) = 2*cos(Pi/7), (approximately 1.801937736) is of degree 3, and 1/rho(7) = 2*1 + 1*rho(7) - 1*rho(7)^2, (approximately 0.5549581320), hence a Q(rho(7))- integer.

%e For Gauss' regular 17-gon rho(17) = 2*cos(Pi/17) (approximately 1.965946199) is of degree 8 and 1/rho(17) = -4*1+ 10*rho(17)^1 + 10*rho(17)^2 - 15*rho(17)^3 -6*rho(17)^4 + 7*rho(17)^5 + 1*rho(17)^6 -1*rho(17)^7, (approximately 0.5086610), hence this is a Q(rho(17))- integer.

%t With[{upto = 100}, Complement[Range[upto], 2*Join[{1}, Select[Range[upto/2], PrimePowerQ]]]] (* _Paolo Xausa_, Aug 30 2024 *)

%o (Python)

%o from sympy import primepi, integer_nthroot

%o def A230078(n):

%o if n == 1: return 1

%o def f(x): return int(n+1+sum(primepi(integer_nthroot(x>>1,k)[0]) for k in range(1,(x>>1).bit_length())))

%o kmin, kmax = 0,1

%o while f(kmax) > kmax:

%o kmax <<= 1

%o while kmax-kmin > 1:

%o kmid = kmax+kmin>>1

%o if f(kmid) <= kmid:

%o kmax = kmid

%o else:

%o kmin = kmid

%o return kmax # _Chai Wah Wu_, Aug 29 2024

%Y Cf. A138929 (complement), 2*A020513, A230079 (1/rho(a(n))).

%K nonn,easy

%O 1,2

%A _Wolfdieter Lang_, Nov 02 2013