A368107 Prime powers p^m such that p | m.

4, 16, 27, 64, 256, 729, 1024, 3125, 4096, 16384, 19683, 65536, 262144, 531441, 823543, 1048576, 4194304, 9765625, 14348907, 16777216, 67108864, 268435456, 387420489, 1073741824, 4294967296, 10460353203, 17179869184, 30517578125, 68719476736, 274877906944, 282429536481

Offset

1,1

Comments

Proper subset of A072873, which in turn is a proper subset of A342090.

This sequence represents the prime power block in A072873 and A342090.

A342090 \ {a(n)} is in A126706.

A072873 \ {{1} U {a(n)}} is in A286708, in turn a proper subset of A001694.

Contains A051674.

Links

Michael De Vlieger, Table of n, a(n) for n = 1..3351

Formula

Sum_{n>=1} 1/a(n) = Sum_{n>=1} 1/A088730(n) = 0.372116188498... . - Amiram Eldar, Jan 20 2024

Example

This sequence contains prime powers of the following form:
  2^2, 2^4, i.e., 2^k such that k is even.
  3^3, 3^6, 3^9, i.e., 3^k such that 3 | k.
  5^5, 5^10, 5^15, i.e., 5^k such that 5 | k, etc.

Maple

N:= 10^13: # for terms <= N
R:= NULL:
for i from 1 do
  p:= ithprime(i);
  if p^p > N then break fi;
  R:= R, seq(p^k, k=p..floor(log[p](N)), p);
od:
sort([R]); # Robert Israel, Jan 16 2024

Mathematica

nn = 10^12; i = 1; p = 2; While[p^p <= nn, p = NextPrime[p] ];
MapIndexed[Set[S[First[#2]], #1] &, Prime@ Range@ PrimePi[p] ];
Union@ Reap[
    While[j = S[i];
    While[S[i]^j < nn,
      Sow[S[i]^j]; j += S[i] ]; j > 2,
    i++] ][[-1, 1]]

Prog

(Python)
import heapq
from itertools import islice
from sympy import nextprime
def agen(): # generator of terms
    v, h, m, nextp = 4, [(4, 2)], 4, 3
    while True:
        v, p = heapq.heappop(h)
        yield v
        if v >= m:
            m = nextp**nextp
            heapq.heappush(h, (m, nextp))
            nextp = nextprime(nextp)
        heapq.heappush(h, (v*p**p, p))
print(list(islice(agen(), 31))) # Michael S. Branicky, Jan 16 2024

Crossrefs

Cf. A001694, A051674, A072873, A088730, A126706, A246547, A286708, A342090.

Sequence in context: A072873 A361078 A072653 * A008478 A201009 A111260

Adjacent sequences: A368104 A368105 A368106 * A368108 A368109 A368110

Keyword

nonn,easy

Author

Michael De Vlieger, Jan 15 2024

Status

approved