login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A246655 Prime powers: Numbers of the form p^k where p is a prime and k >= 1. 47
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. Adams-Watters, Sep 01 2014

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy .

Last modified July 24 16:44 EDT 2017. Contains 289775 sequences.