Found with PFGW Version 3.6.0

8731^10002-10002^8731 is 2-PRP! (86.1949s+0.0047s)
8731^10002-10002^8731 is 3-PRP! (85.8009s+0.0059s)
8731^10002-10002^8731 is 5-PRP! (88.1250s+0.0059s)
8731^10002-10002^8731 is 7-PRP! (86.0426s+0.0046s)
8731^10002-10002^8731 is 11-PRP! (85.8206s+0.0040s)

Primality testing 8731^10002-10002^8731 [N-1/N+1, Brillhart-Lehmer-Selfridge]
Running N-1 test using base 5
Running N+1 test using discriminant 11, base 1+sqrt(11)
Calling N-1 BLS with factored part 0.03% and helper 0.01% (0.10% proof)
8731^10002-10002^8731 is Fermat and Lucas PRP! (387.7135s+0.0045s)

To be completed...