Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A001597
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A001597 Perfect powers: m^k where m is an integer and k >= 2.
(Formerly M3326 N1336)
+0
189
1, 4, 8, 9, 16, 25, 27, 32, 36, 49, 64, 81, 100, 121, 125, 128, 144, 169, 196, 216, 225, 243, 256, 289, 324, 343, 361, 400, 441, 484, 512, 529, 576, 625, 676, 729, 784, 841, 900, 961, 1000, 1024, 1089, 1156, 1225, 1296, 1331, 1369, 1444, 1521, 1600, 1681, 1728, 1764 (list; graph; listen)
OFFSET

1,2

COMMENT

Catalan's conjecture (now a theorem) is that 1 occurs just once as a difference, between 8 and 9.

Goldbach showed that Sum 1/(a(n)-1) = 1.

REFERENCES

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).

H. W. Gould, Problem H-170, Fib. Quart., 8 (1970), 268.

R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 66.

D. J. Newman, A Problem Seminar, Springer; see Problem #72.

R. Schoof, Catalan's Conjecture, Springer-Verlag, 2008, p. 1.

LINKS

David W. Wilson, Table of n, a(n) for n = 1..10000

A. Dendane, Power (Exponential) Calculator

Serhat Sevki Dincer, Powers up to 2^50

Alf van der Poorten, Remarks on the sequence of 'perfect' powers

M. Waldschmidt, Open Diophantine problems

Eric Weisstein's World of Mathematics, Perfect Power

FORMULA

Formulae from postings to the Number Theory List by various authors, 2002:

Sum_{i=2}^{infty} sum_{j=2}^{infty} 1/i^j =1;

Sum_{k=1}^infty 1/(a_k-1)=1;

Sum_{k=1}^infty 1/(a_k+1)= pi^2 / 3 - 5/2;

Sum_{k=1}^infty 1/a_k = sum_{n=2}^infty mu(n)(1- zeta(n)) approx = .87446436840494...

For asymptotics see Newman.

EXAMPLE

x + 4*x^2 + 8*x^3 + 9*x^4 + 16*x^5 + 25*x^6 + 27*x^7 + 32*x^8 + 36*x^9 + ...

MATHEMATICA

Union[ Join[{1}, Flatten[ Table[ n^i, {n, 2, Sqrt[1800]}, {i, 2, Log[n, 1800]}]]]]

Join[{1}, Select[Range@1848, GCD @@ Last /@ FactorInteger@# > 1 &]] (* or *)

PROGRAM

(MAGMA) [1] cat [n : n in [2..1000] | IsPower(n) ];

(PARI) {a(n) = local(m, c); if( n<2, n==1, c=1; m=1; while( c<n, m++; if( ispower(m), c++)); m)} /* Michael Somos Aug 05 2009 */

CROSSREFS

Cf. A023055, A023057, A070428, A074981, A025478.

Cf. A089579, A089580 (number of exact powers < 10^n).

Complement of A007916.

Sequence in context: A080366 A001694 A157985 this_sequence A072777 A076292 A090516

Adjacent sequences: A001594 A001595 A001596 this_sequence A001598 A001599 A001600

KEYWORD

nonn,easy,nice

AUTHOR

N. J. A. Sloane (njas(AT)research.att.com).

page 1

Search completed in 0.003 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 25 20:09 EST 2009. Contains 167514 sequences.


AT&T Labs Research