Notons rq l'exposant de chaque facteur premier q dans la décomposition de f.
Soient p un facteur premier de n et dq l'ordre multiplicatif de aq modulo p. Alors, dq divise n – 1 mais pas (n – 1)/q, donc (n – 1)/dq est un entier non divisible par q. Or n – 1 est divisible par qrq.
Par conséquent, qrq divise dq et (a fortiori) p – 1. Le produit f des qrq divise donc aussi p – 1
↑(en) H. C. Pocklington(de), « The determination of the prime or composite nature of large numbers by Fermat's theorem », Proc. Cambridge Philos. Soc., vol. 18, 1914-1916, p. 29-30.