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

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A057655 The circle problem: number of points (x,y) in square lattice with x^2+y^2 <= n. 13
1, 5, 9, 9, 13, 21, 21, 21, 25, 29, 37, 37, 37, 45, 45, 45, 49, 57, 61, 61, 69, 69, 69, 69, 69, 81, 89, 89, 89, 97, 97, 97, 101, 101, 109, 109, 113, 121, 121, 121, 129, 137, 137, 137, 137, 145, 145, 145, 145, 149, 161, 161, 169, 177, 177, 177 (list; graph; refs; listen; history; internal format)
OFFSET

0,2

REFERENCES

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 106.

F. Fricker, Einfuehrung in die Gitterpunktlehre, Birkhaeuser, Boston, 1982.

P. de la Harpe, Topics in Geometric Group Theory, Univ. Chicago Press, 2000, p. 5.

E. Kraetzel, lattice Points, Kluwer, Dordrecht, 1988.

C. D. Olds, A. Lax and G. P. Davidoff, The Geometry of Numbers, Math. Assoc. Amer., 2000, p. 51.

W. Sierpinski, Elementary Theory of Numbers, Elsevier, North-Holland, 1988.

LINKS

T. D. Noe, Table of n, a(n) for n = 0..1000

F. Richman, Count Gaussian integers of norm less than and equal to n

W. Sierpi\'{n}ski, Elementary Theory of Numbers, Warszawa 1964.

F. Richman, Counting Gaussian integers in a disk

FORMULA

a(n) = 1 + 4*{ [n/1] - [n/3] + [n/5] - [n/7] + ... }. - Gauss

a(n) = 1 + 4*Sum_{ k = 0 .. [sqrt(n)] } [ sqrt(n-k^2) ]. - Liouville (?)

a(n) - Pi*n = O(sqrt(n)) (Gauss). a(n) - Pi*n = O(n^c), c = 23/73 + epsilon ~ 0.3151 (Huxley). If a(n) - Pi*n = O(n^c) then c > 1/4 (Landau, Hardy). It is conjectured that a(n) - Pi*n = O(n^(1/4 + epsilon)) for all epsilon >0.

a(n) = A122510(2,n). [From R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Apr 21 2010]

MATHEMATICA

f[n_] := 1 + 4Sum[ Floor@ Sqrt[n - k^2], {k, 0, Sqrt[n]}]; Table[ f[n], {n, 0, 60}] (from Robert G. Wilson v (rgwv(at)rgwv.com), Jun 16 2006)

PROG

(PARI) a(n)=sum(x=-n, n, sum(y=-n, n, if((sign(x^2+y^2-n)+1)*sign(x^2+y^2-n), 0, 1)))

(Haskell)

a057655 n = length [(x, y) | x <- [-n..n], y <- [-n..n], x^2 + y^2 <= n]

-- Reinhard Zumkeller, Jan 23 2012

CROSSREFS

Partial sums of A004018. Cf. A057656, A057961, A057962. For another version see A000328.

A014198(n) + 1.

Sequence in context: A175363 A073168 A127500 * A175374 A141124 A046255

Adjacent sequences:  A057652 A057653 A057654 * A057656 A057657 A057658

KEYWORD

nonn,easy,nice

AUTHOR

N. J. A. Sloane (njas(AT)research.att.com), Oct 15 2000

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

Content is available under The OEIS End-User License Agreement .

Last modified February 16 15:27 EST 2012. Contains 205930 sequences.