

A246655


Prime powers: Numbers of the form p^k where p is a prime and k >= 1.


45



2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 37, 41, 43, 47, 49, 53, 59, 61, 64, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 121, 125, 127, 128, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,1


COMMENTS

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 cototients A051953(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


LINKS

Robert Israel, Table of n, a(n) for n = 1..10000
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.
Eric Weisstein's World of Mathematics, Prime Power
Eric Weisstein's World of Mathematics, Projective Plane
Wikipedia, Prime power
Index entries for "core" sequences


FORMULA

a(n) is characterized by A001221(a(n)) = 1.
a(n) is characterized by A014963(a(n)) != 1.
Euler's A000010(a(n)) = a(n)*(1  1/A014963(a(n)).
All three relations above are not true for A000961(n) instead of a(n).


MAPLE

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


MATHEMATICA

Select[Range[222], PrimePowerQ]


PROG

(Sage)
filter(lambda n: sloane.A001221(n) == 1, (1..222))
(PARI)
[p p < [1..222], isprimepower(p)]
(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


CROSSREFS

Cf. A000961, A001221, A014963, A069513, A246547.
Sequence in context: A059046 A144711 A036116 * A000961 A128603 A195943
Adjacent sequences: A246652 A246653 A246654 * A246656 A246657 A246658


KEYWORD

nonn,nice,core,easy


AUTHOR

Peter Luschny and Franklin T. AdamsWatters, Sep 01 2014


STATUS

approved



