Abstract. Let f be a nonconstant squarefree polynomial. Which of the values f (c + 1), f (c + 2), . . . , f (c + H) are locally square at all small primes? This paper presents an algorithm that answers this question in time H/M 2+o(1) for an average small c as H → ∞, where M = H 1/log 2 log H . In contrast, the usual method takes time H/M 1+o(1) . This paper also presents the results of two record-setting computations: an enumeration of locally square integers up to 24 · 2 64 , and an enumeration of locally square values of x 3 + y 7 for small x and y.