In this article, we consider the following capacitated covering problem. We are given a set P of n points and a set B of balls from some metric space, and a positive integer U that represents the capacity of each of the balls in B. We would like to compute a subset B ⊆ B of balls and assign each point in P to some ball in B that contains it, such that the number of points assigned to any ball is at most U . The objective function that we would like to minimize is the cardinality of B .We consider this problem in arbitrary metric spaces as well as Euclidean spaces of constant dimension. In the metric setting, even the uncapacitated version of the problem is hard to approximate to within a logarithmic factor. In the Euclidean setting, the best known approximation guarantee in dimensions 3 and higher is logarithmic in the number of points. Thus we focus on obtaining "bi-criteria" approximations. In particular, we are allowed to expand the balls in our solution by some factor, but optimal solutions do not have that flexibility. Our main result is that allowing constant factor expansion of the input balls suffices to obtain constant approximations for these problems. In fact, in the Euclidean setting, only (1 + ) factor expansion is sufficient for any > 0, with the approximation factor being a polynomial in 1/ . We obtain these results using a unified scheme for rounding the natural LP relaxation; this scheme may be useful for other capacitated covering problems. We also complement these bi-criteria approximations by obtaining hardness of approximation results that shed light on our understanding of these problems. complexity theoretic assumptions [14]. The same is true for the MUC, as demonstrated by the following reduction from Set Cover. We take a ball of radius 1 corresponding to each set, and a point corresponding to each element. If an element is in a set, then the distance between the center of the corresponding ball and the point is 1. We consider the metric space induced by the centers and the points. It is easy to see that any solution for this instance of MUC directly gives a solution for the input instance of the general Set Cover, implying that for MUC, it is not possible to get any approximation guarantee better than the O(log n) bound for Set Cover.The MUC in fixed dimensional Euclidean spaces has been extensively studied. One interesting variant is when the allowed set B of balls consists of all unit balls. Hochbaum and Maass [19] gave a polynomial time approximation scheme (PTAS) for this using a grid shifting strategy. When B is an arbitrary finite set of balls, the problem seems to be much harder. An O(1) approximation algorithm in the 2-dimensional Euclidean plane was given by Brönnimann and Goodrich [9]. More recently, a PTAS was obatined by Mustafa and Ray [26]. In dimensions 3 and higher, the best known approximation guarantee is still O(log n). Motivated by this, Har-Peled and Lee [18] gave a PTAS for a bi-criteria version where the algorithm is allowed to expand the input balls by a (1 + ) fac...