Yeah, I did it this way at first, looping through the primes and testing for square-ness. But then it seemed better to do it the other way, though it probably doesn't make much difference.

By the way, I tried to use math.plt, but it wouldn't compile for v3.99:

setup-plt: number-theory.ss:1578:21: compile: unbound variable in module in: reverse!