Earlier this evening (morning), I posted a question about showing that a finite number of dyadic squares can fill up an arbitrary proportion of the unit disk. I'm sure there are better ways to prove this, but the way I came up with is to look at each square in the first quadrant or, more specifically, each square's upper-right-hand corner, which has coordinates of the form $\left( \frac{p}{2^k}, \frac{q}{2^k}\right)$. Since this is the unit disk, I know that
$$\sqrt{\left(\frac{p}{2^k} \right)^2 + \left( \frac{q}{2^k} \right)^2 }\leq 1$$
or, equivalently, that $p^2 + q^2 \leq 4^k$.
What I want to do is count the number of solutions there are in $\mathbb{Z}^+\times \mathbb{Z}^+$, but I don't have any good idea on how to go about doing that. I have a lower bound, $\frac{\left(2^k - 1 \right)\left(2^k\right)}{2}$, since $(1,1), (1,2), (2,1), \ldots (1, 2^k - 1), (2, 2^k - 2), \ldots (2^k-1, 1)$ are all solutions. But this isn't nearly strong enough for what I need. Is there a closed form for the size of the solution set for any given $k$? I'm thinking that it would be useful if I had something I could compare to $4^k$ via a limit, though if there are other options, I'd certainly be interested in them.