A passive scalar is advected by a velocity field, with a nonuniform spatial source that maintains concentration inhomogeneities. For example, the scalar could be temperature with a source consisting of hot and cold spots, such that the mean temperature is constant. Which source distributions are best mixed by this velocity field? This question has a straightforward yet rich answer that is relevant to real mixing problems. We use a multiscale measure of steady-state enhancement to mixing and optimize it by a variational approach. We then solve the resulting Euler-Lagrange equation for a perturbed uniform flow and for simple cellular flows. The optimal source distributions have many broad features that are as expected: they avoid stagnation points, favor regions of fast flow, and their contours are aligned such that the flow blows hot spots onto cold and vice versa. However, the detailed structure varies widely with diffusivity and other problem parameters. Though these are model problems, the optimization procedure is simple enough to be adapted to more complex situations.