Found with PFGW Version 3.7.5

8042^10601-10601^8042 is 2-PRP! (36.6798s+0.0013s)
8042^10601-10601^8042 is 3-PRP! (36.7734s+0.0015s)
8042^10601-10601^8042 is 5-PRP! (36.4995s+0.0013s)
8042^10601-10601^8042 is 7-PRP! (41.0234s+0.0029s)
8042^10601-10601^8042 is 11-PRP! (42.8115s+0.0018s)

Primality testing 8042^10601-10601^8042 [N-1/N+1, Brillhart-Lehmer-Selfridge]
Running N-1 test using base 3
Running N-1 test using base 7
Running N+1 test using discriminant 13, base 2+sqrt(13)
Calling N+1 BLS with factored part 0.01% and helper 0.01% (0.02% proof)
8042^10601-10601^8042 is Fermat and Lucas PRP! (213.1701s+0.0133s)

To be completed...