1,1

The elements are called prime powers in contrast to the powers of primes which are the numbers of the same form but with k >= 0, cf.A000961.

Every nonzero integer is the product of elements of this sequence which are relatively prime and an element of {-1, 1}. This product is up to a rearrangement of the factors unique. (This statement is the fundamental theorem of arithmetic.)

These numbers are the numbers such that the von Mangoldt function is nonzero.

These numbers are the numbers of elements in finite fields. -Franz Vrabec,Aug 11 2004

A positive integer n is a prime power if and only if nZ is a primary ideal of Z. -John Cremona,Sep 02 2014

Also, numbers n divisible by their cototientsA051953(n). -Ivan Neretin,May 29 2016

Numbers n such that (theta_3(q) - theta_3(q^n)) / 2 is the g.f. of a multiplicative sequence. -Michael Somos,Oct 17 2016

Numbers that are evenly divisible by exactly one prime number. -Lee A. Newberg,May 07 2018

Ram proved that these are precisely the numbers n such that the binomial coefficients n!/(m!(n-m)!) for m = 1..n-1 have a common factor greater than 1 (which is the unique prime dividing n). See Joris, Oestreicher & Steinig for a generalization. -Charles R Greathouse IV,Apr 24 2019

Blagojević & Ziegler prove that for these n and for any convex polygon in the plane, the polygon can be partitioned into n polygons with equal area and equal perimeter. The result is conjectured (by Nandakumar & Rao, who proved the case n = 2) to hold for all n. -Charles R Greathouse IV,Apr 24 2019

Numbers n such thatA367064(n) < 0. -Chai Wah Wu,Nov 06 2023

Robert Israel,Table of n, a(n) for n = 1..10000

Brady Haran and Günter Ziegler,Cannons and Sparrows,Numberphile video (2018)

H. Joris, C. Oestreicher and J. Steinig,The greatest common divisor of certain sets of binomial coefficients,Journal of Number Theory 21 (1985), pp. 101-119.

Pavle V. M. Blagojević and Günter M. Ziegler,Convex equipartitions via equivariant obstruction theory,arXiv:1202.5504 [math.AT], 2012-2013; Israel Journal of Mathematics 200:1 (June 2014), pp 49-77.

Hans Montanus and Ron Westdijk,Cellular Automation and Binomials,Green Blue Mathematics (2022), p. 90.

Laurentiu Panaitopol,Some of the properties of the sequence of powers of prime numbers,Rocky Mountain Journal of Mathematics, Volume 31, Number 4, Winter 2001.

Balak Ram,Common factors of n!/(m!(n-m)!), (m = 1, 2,... n-1),Journal of the Indian Mathematical Club (Madras) 1 (1909), pp. 39-43.

Eric Weisstein's World of Mathematics,Prime Power

Eric Weisstein's World of Mathematics,Projective Plane

Chai Wah Wu,Algorithms for complementary sequences,arXiv:2409.05844 [math.NT], 2024.

Wikipedia,Prime power

a(n) is characterized byA001221(a(n)) = 1.

a(n) is characterized byA014963(a(n))!= 1.

Euler'sA000010(a(n)) = a(n)*(1 - 1/A014963(a(n))).

All three relations above are not true forA000961(n) instead of a(n).

Sum_{k=1..n} 1/a(k) ~ log(log(a(n))) +A077761+A136141.-François Huppé,Jul 31 2024

select(t -> nops(numtheory:-factorset(t))=1, [$1..1000]); #Robert Israel,Sep 01 2014

A246655:= proc(n)A000961(n+1) end proc: #R. J. Mathar,Jan 09 2017

isprimepower:= n -> nops(NumberTheory:-PrimeFactors(n)) = 1: #Peter Luschny,Oct 09 2022

Select[Range[222], PrimePowerQ]

(Sage)

[n for n in (1..222) if sloane.A001221(n) == 1]

(PARI)

[p| p <- [1..222], isprimepower(p)]

(PARI) list(lim)=my(v=List(primes([2, lim\=1]))); for(e=2, logint(lim, 2), forprime(p=2, sqrtnint(lim, e), listput(v, p^e))); Set(v) \\Charles R Greathouse IV,Feb 03 2023

(Python)

from sympy import primerange

m = 10**5

A246655= []

for p in primerange(1, m):

pe = p

while pe < m:

A246655.append(pe)

pe *= p

A246655= sorted(A246655) #Chai Wah Wu,Sep 04 2014

(Python)

from sympy import primepi, integer_nthroot

defA246655(n):

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

kmin, kmax = 1, 2

while f(kmax) >= kmax:

kmax <<= 1

while True:

kmid = kmax+kmin>>1

if f(kmid) < kmid:

kmax = kmid

else:

kmin = kmid

if kmax-kmin <= 1:

break

return kmax #Chai Wah Wu,Aug 20 2024

There are four different sequences which may legitimately be called "prime powers":A000961(p^k, k >= 0),A246655(p^k, k >= 1),A246547(p^k, k >= 2),A025475(p^k, k=0 and k >= 2). When you refer to "prime powers", be sure to specify which of these you mean. AlsoA001597is the sequence of nontrivial powers n^k, n >= 1, k >= 2. -N. J. A. Sloane,Mar 24 2018

Cf.A001221,A014963,A069513,A367064.

Partial sums ofA275120.

Sequence in context:A144711A337935A036116*A000961A128603A195943

Adjacent sequences:A246652A246653A246654*A246656A246657A246658

nonn,nice,core,easy

Peter LuschnyandFranklin T. Adams-Watters,Sep 01 2014

approved