|
Search: id:A014198
|
|
|
| A014198 |
|
Number of integer solutions to x^2+y^2 <= n excluding (0,0). |
|
+0
5
|
|
| 0, 4, 8, 8, 12, 20, 20, 20, 24, 28, 36, 36, 36, 44, 44, 44, 48, 56, 60, 60, 68, 68, 68, 68, 68, 80, 88, 88, 88, 96, 96, 96, 100, 100, 108, 108, 112, 120, 120, 120, 128, 136, 136, 136, 136, 144, 144, 144, 144, 148, 160, 160, 168, 176, 176, 176, 176, 176, 184, 184
(list; graph; listen)
|
|
|
OFFSET
|
0,2
|
|
|
COMMENT
|
a(32)/32 = 100/32 = 3.125; lim a(n)/n = Pi.
|
|
REFERENCES
|
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, th. 339
|
|
LINKS
|
Eric Weisstein's World of Mathematics, Link to a section of The World of Mathematics.
|
|
EXAMPLE
|
For n=2 the 8 solutions are (x,y) = (+-1,0), (0,+-1), (+-1,+-1).
|
|
MAPLE
|
with(numtheory): for n from 0 to 60 do h[ n ] := nops([ thue( abs( x^2+y^2) <= n, [ x, y ] ) ]); od; [ seq(h[ i ], i=0..60) ];
|
|
PROGRAM
|
(PARI) a(n)=local(j); j=sqrtint(n); sum(x=-j, j, sum(y=-j, j, x^2+y^2<=n))-1
|
|
CROSSREFS
|
Cf. A014200. 4 * partial sums of A002654 = -1 + partial sums of A004018.
Sequence in context: A029679 A072541 A098352 this_sequence A083744 A114027 A005877
Adjacent sequences: A014195 A014196 A014197 this_sequence A014199 A014200 A014201
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
njas
|
|
|
Search completed in 0.002 seconds
|
