A176509 Composite numbers m for which A064380(m) = A000010(m).

8, 27, 125, 128, 343, 1331, 2187, 2197, 4913, 6859, 12167, 24389, 29791, 32768, 50653, 68921, 78125, 79507, 103823, 148877, 205379, 226981, 300763, 357911, 389017, 493039, 571787, 704969, 823543, 912673, 1030301, 1092727, 1225043, 1295029, 1442897, 2048383, 2248091

Offset

1,1

Comments

Theorem. A064380(m) = A000010(m) iff m has the form m=p^(2^k-1), k>=1, p a prime. Eliminating the primes (k=1), the terms of the sequence have this form for k>1. All terms of A030078 (k=2) and A092759 (k=3) and prime powers of A010803 (k=4) are in the sequence, for example.

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Formula

a(n) ~ n^3 log^3 n. - Charles R Greathouse IV, Feb 19 2013

Sum_{n>=1} 1/a(n) = Sum_{k>=2} 1/P(2^k-1) = 0.183077059924063305405..., where P(s) is the prime zeta function. - Amiram Eldar, Jul 11 2024

Mathematica

seq[max_] := Module[{ps = Select[Range[Floor[Surd[max, 3]]], PrimeQ], e, k, s = {}}, Do[e = Floor[Log[ps[[i]], max]]; k = Floor[Log2[e + 1]]; s = Join[s, ps[[i]]^(2^Range[2, k] - 1)], {i, 1, Length[ps]}]; Sort[s]]; seq[3*10^6] (* Amiram Eldar, Mar 26 2023 *)

Prog

(PARI) is(n)=my(e=isprimepower(n)); e>2 && 2^valuation(e+1, 2)==e+1 \\ Charles R Greathouse IV, Feb 19 2013

Crossrefs

Cf. A000010, A010803, A030078, A050376, A056824, A064380, A092759, A176472.

Sequence in context: A331285 A183316 A045971 * A179125 A030160 A030162

Adjacent sequences: A176506 A176507 A176508 * A176510 A176511 A176512

Keyword

nonn

Author

Vladimir Shevelev, Apr 19 2010

Extensions

128 inserted, 1024 deleted, 2187 inserted, 32768 inserted, etc. - R. J. Mathar, Nov 21 2010

More terms from Amiram Eldar, Mar 26 2023

Status

approved