Found with PFGW Version 3.5.2

2042^3331-3331^2042 is 2-PRP! (25.0703s+0.0050s)
2042^3331-3331^2042 is 3-PRP! (29.3625s+0.0127s)
2042^3331-3331^2042 is 5-PRP! (27.4101s+0.0062s)
2042^3331-3331^2042 is 7-PRP! (26.6940s+0.0066s)
2042^3331-3331^2042 is 11-PRP! (26.4577s+0.0049s)

Primality testing 2042^3331-3331^2042 [N-1/N+1, Brillhart-Lehmer-Selfridge]
Running N-1 test using base 3
Running N+1 test using discriminant 7, base 6+sqrt(7)
Calling N+1 BLS with factored part 0.02% and helper 0.01% (0.08% proof)
2042^3331-3331^2042 is Fermat and Lucas PRP! (130.4173s+0.0157s)

To be completed...