Theorem pockthi 12496

Description: Pocklington's theorem, which gives a sufficient criterion for a number N to be prime. This is the preferred method for verifying large primes, being much more efficient to compute than trial division. This form has been optimized for application to specific large primes; see pockthg 12495 for a more general closed-form version. (Contributed by Mario Carneiro, 2-Mar-2014.)

Hypotheses

Ref	Expression
pockthi.p	|- P e. Prime
pockthi.g	|- G e. NN
pockthi.m	|- M = ( G x. P )
pockthi.n	|- N = ( M + 1 )
pockthi.d	|- D e. NN
pockthi.e	|- E e. NN
pockthi.a	|- A e. NN
pockthi.fac	|- M = ( D x. ( P ^ E ) )
pockthi.gt	|- D < ( P ^ E )
pockthi.mod	|- ( ( A ^ M ) mod N ) = ( 1 mod N )
pockthi.gcd	|- ( ( ( A ^ G ) - 1 ) gcd N ) = 1

Assertion

Ref	Expression
pockthi	|- N e. Prime

Proof of Theorem pockthi

Dummy variables x y are mutually distinct and distinct from all other variables.

Step	Hyp	Ref	Expression
1		pockthi.d	. 2 |- D e. NN
2		pockthi.p	. . . . . 6 |- P e. Prime
3		prmnn 12248	. . . . . 6 |- ( P e. Prime -> P e. NN )
4	2, 3	ax-mp 5	. . . . 5 |- P e. NN
5		pockthi.e	. . . . . 6 |- E e. NN
6	5	nnnn0i 9248	. . . . 5 |- E e. NN0
7		nnexpcl 10623	. . . . 5 |- ( ( P e. NN /\ E e. NN0 ) -> ( P ^ E ) e. NN )
8	4, 6, 7	mp2an 426	. . . 4 |- ( P ^ E ) e. NN
9	8	a1i 9	. . 3 |- ( D e. NN -> ( P ^ E ) e. NN )
10		id 19	. . 3 |- ( D e. NN -> D e. NN )
11		pockthi.gt	. . . 4 |- D < ( P ^ E )
12	11	a1i 9	. . 3 |- ( D e. NN -> D < ( P ^ E ) )
13		pockthi.n	. . . . 5 |- N = ( M + 1 )
14		pockthi.fac	. . . . . . 7 |- M = ( D x. ( P ^ E ) )
15	1	nncni 8992	. . . . . . . 8 |- D e. CC
16	8	nncni 8992	. . . . . . . 8 |- ( P ^ E ) e. CC
17	15, 16	mulcomi 8025	. . . . . . 7 |- ( D x. ( P ^ E ) ) = ( ( P ^ E ) x. D )
18	14, 17	eqtri 2214	. . . . . 6 |- M = ( ( P ^ E ) x. D )
19	18	oveq1i 5928	. . . . 5 |- ( M + 1 ) = ( ( ( P ^ E ) x. D ) + 1 )
20	13, 19	eqtri 2214	. . . 4 |- N = ( ( ( P ^ E ) x. D ) + 1 )
21	20	a1i 9	. . 3 |- ( D e. NN -> N = ( ( ( P ^ E ) x. D ) + 1 ) )
22		prmdvdsexpb 12287	. . . . . . 7 |- ( ( x e. Prime /\ P e. Prime /\ E e. NN ) -> ( x || ( P ^ E ) <-> x = P ) )
23	2, 5, 22	mp3an23 1340	. . . . . 6 |- ( x e. Prime -> ( x || ( P ^ E ) <-> x = P ) )
24		pockthi.m	. . . . . . . . . . . . 13 |- M = ( G x. P )
25		pockthi.g	. . . . . . . . . . . . . 14 |- G e. NN
26	25, 4	nnmulcli 9004	. . . . . . . . . . . . 13 |- ( G x. P ) e. NN
27	24, 26	eqeltri 2266	. . . . . . . . . . . 12 |- M e. NN
28	27	nncni 8992	. . . . . . . . . . 11 |- M e. CC
29		ax-1cn 7965	. . . . . . . . . . 11 |- 1 e. CC
30	28, 29, 13	mvrraddi 8236	. . . . . . . . . 10 |- ( N - 1 ) = M
31	30	oveq2i 5929	. . . . . . . . 9 |- ( A ^ ( N - 1 ) ) = ( A ^ M )
32	31	oveq1i 5928	. . . . . . . 8 |- ( ( A ^ ( N - 1 ) ) mod N ) = ( ( A ^ M ) mod N )
33		pockthi.mod	. . . . . . . . 9 |- ( ( A ^ M ) mod N ) = ( 1 mod N )
34		peano2nn 8994	. . . . . . . . . . . . 13 |- ( M e. NN -> ( M + 1 ) e. NN )
35	27, 34	ax-mp 5	. . . . . . . . . . . 12 |- ( M + 1 ) e. NN
36	13, 35	eqeltri 2266	. . . . . . . . . . 11 |- N e. NN
37		nnq 9698	. . . . . . . . . . 11 |- ( N e. NN -> N e. QQ )
38	36, 37	ax-mp 5	. . . . . . . . . 10 |- N e. QQ
39	27	nngt0i 9012	. . . . . . . . . . . 12 |- 0 < M
40	27	nnrei 8991	. . . . . . . . . . . . 13 |- M e. RR
41		1re 8018	. . . . . . . . . . . . 13 |- 1 e. RR
42		ltaddpos2 8472	. . . . . . . . . . . . 13 |- ( ( M e. RR /\ 1 e. RR ) -> ( 0 < M <-> 1 < ( M + 1 ) ) )
43	40, 41, 42	mp2an 426	. . . . . . . . . . . 12 |- ( 0 < M <-> 1 < ( M + 1 ) )
44	39, 43	mpbi 145	. . . . . . . . . . 11 |- 1 < ( M + 1 )
45	44, 13	breqtrri 4056	. . . . . . . . . 10 |- 1 < N
46		q1mod 10427	. . . . . . . . . 10 |- ( ( N e. QQ /\ 1 < N ) -> ( 1 mod N ) = 1 )
47	38, 45, 46	mp2an 426	. . . . . . . . 9 |- ( 1 mod N ) = 1
48	33, 47	eqtri 2214	. . . . . . . 8 |- ( ( A ^ M ) mod N ) = 1
49	32, 48	eqtri 2214	. . . . . . 7 |- ( ( A ^ ( N - 1 ) ) mod N ) = 1
50		oveq2 5926	. . . . . . . . . . . 12 |- ( x = P -> ( ( N - 1 ) / x ) = ( ( N - 1 ) / P ) )
51	25	nncni 8992	. . . . . . . . . . . . . . 15 |- G e. CC
52	4	nncni 8992	. . . . . . . . . . . . . . 15 |- P e. CC
53	51, 52	mulcomi 8025	. . . . . . . . . . . . . 14 |- ( G x. P ) = ( P x. G )
54	30, 24, 53	3eqtrri 2219	. . . . . . . . . . . . 13 |- ( P x. G ) = ( N - 1 )
55	36	nncni 8992	. . . . . . . . . . . . . . 15 |- N e. CC
56	55, 29	subcli 8295	. . . . . . . . . . . . . 14 |- ( N - 1 ) e. CC
57	4	nnap0i 9013	. . . . . . . . . . . . . 14 |- P # 0
58	56, 52, 51, 57	divmulapi 8785	. . . . . . . . . . . . 13 |- ( ( ( N - 1 ) / P ) = G <-> ( P x. G ) = ( N - 1 ) )
59	54, 58	mpbir 146	. . . . . . . . . . . 12 |- ( ( N - 1 ) / P ) = G
60	50, 59	eqtrdi 2242	. . . . . . . . . . 11 |- ( x = P -> ( ( N - 1 ) / x ) = G )
61	60	oveq2d 5934	. . . . . . . . . 10 |- ( x = P -> ( A ^ ( ( N - 1 ) / x ) ) = ( A ^ G ) )
62	61	oveq1d 5933	. . . . . . . . 9 |- ( x = P -> ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) = ( ( A ^ G ) - 1 ) )
63	62	oveq1d 5933	. . . . . . . 8 |- ( x = P -> ( ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = ( ( ( A ^ G ) - 1 ) gcd N ) )
64		pockthi.gcd	. . . . . . . 8 |- ( ( ( A ^ G ) - 1 ) gcd N ) = 1
65	63, 64	eqtrdi 2242	. . . . . . 7 |- ( x = P -> ( ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 )
66		pockthi.a	. . . . . . . . 9 |- A e. NN
67	66	nnzi 9338	. . . . . . . 8 |- A e. ZZ
68		oveq1 5925	. . . . . . . . . . . 12 |- ( y = A -> ( y ^ ( N - 1 ) ) = ( A ^ ( N - 1 ) ) )
69	68	oveq1d 5933	. . . . . . . . . . 11 |- ( y = A -> ( ( y ^ ( N - 1 ) ) mod N ) = ( ( A ^ ( N - 1 ) ) mod N ) )
70	69	eqeq1d 2202	. . . . . . . . . 10 |- ( y = A -> ( ( ( y ^ ( N - 1 ) ) mod N ) = 1 <-> ( ( A ^ ( N - 1 ) ) mod N ) = 1 ) )
71		oveq1 5925	. . . . . . . . . . . . 13 |- ( y = A -> ( y ^ ( ( N - 1 ) / x ) ) = ( A ^ ( ( N - 1 ) / x ) ) )
72	71	oveq1d 5933	. . . . . . . . . . . 12 |- ( y = A -> ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) = ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) )
73	72	oveq1d 5933	. . . . . . . . . . 11 |- ( y = A -> ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = ( ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) )
74	73	eqeq1d 2202	. . . . . . . . . 10 |- ( y = A -> ( ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 <-> ( ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) )
75	70, 74	anbi12d 473	. . . . . . . . 9 |- ( y = A -> ( ( ( ( y ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) <-> ( ( ( A ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) ) )
76	75	rspcev 2864	. . . . . . . 8 |- ( ( A e. ZZ /\ ( ( ( A ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) ) -> E. y e. ZZ ( ( ( y ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) )
77	67, 76	mpan 424	. . . . . . 7 |- ( ( ( ( A ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( A ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) -> E. y e. ZZ ( ( ( y ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) )
78	49, 65, 77	sylancr 414	. . . . . 6 |- ( x = P -> E. y e. ZZ ( ( ( y ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) )
79	23, 78	biimtrdi 163	. . . . 5 |- ( x e. Prime -> ( x || ( P ^ E ) -> E. y e. ZZ ( ( ( y ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) ) )
80	79	rgen 2547	. . . 4 |- A. x e. Prime ( x || ( P ^ E ) -> E. y e. ZZ ( ( ( y ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) )
81	80	a1i 9	. . 3 |- ( D e. NN -> A. x e. Prime ( x || ( P ^ E ) -> E. y e. ZZ ( ( ( y ^ ( N - 1 ) ) mod N ) = 1 /\ ( ( ( y ^ ( ( N - 1 ) / x ) ) - 1 ) gcd N ) = 1 ) ) )
82	9, 10, 12, 21, 81	pockthg 12495	. 2 |- ( D e. NN -> N e. Prime )
83	1, 82	ax-mp 5	1 |- N e. Prime

Syntax hints:   -> wi 4   /\ wa 104   <-> wb 105   = wceq 1364   e. wcel 2164   A.wral 2472   E.wrex 2473   class class class wbr 4029  (class class class)co 5918   RRcr 7871   0cc0 7872   1c1 7873   + caddc 7875   x. cmul 7877   < clt 8054   - cmin 8190   / cdiv 8691   NNcn 8982   NN0cn0 9240   ZZcz 9317   QQcq 9684   mod cmo 10393   ^cexp 10609   || cdvds 11930   gcd cgcd 12079   Primecprime 12245

This theorem was proved from axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-ia1 106  ax-ia2 107  ax-ia3 108  ax-in1 615  ax-in2 616  ax-io 710  ax-5 1458  ax-7 1459  ax-gen 1460  ax-ie1 1504  ax-ie2 1505  ax-8 1515  ax-10 1516  ax-11 1517  ax-i12 1518  ax-bndl 1520  ax-4 1521  ax-17 1537  ax-i9 1541  ax-ial 1545  ax-i5r 1546  ax-13 2166  ax-14 2167  ax-ext 2175  ax-coll 4144  ax-sep 4147  ax-nul 4155  ax-pow 4203  ax-pr 4238  ax-un 4464  ax-setind 4569  ax-iinf 4620  ax-cnex 7963  ax-resscn 7964  ax-1cn 7965  ax-1re 7966  ax-icn 7967  ax-addcl 7968  ax-addrcl 7969  ax-mulcl 7970  ax-mulrcl 7971  ax-addcom 7972  ax-mulcom 7973  ax-addass 7974  ax-mulass 7975  ax-distr 7976  ax-i2m1 7977  ax-0lt1 7978  ax-1rid 7979  ax-0id 7980  ax-rnegex 7981  ax-precex 7982  ax-cnre 7983  ax-pre-ltirr 7984  ax-pre-ltwlin 7985  ax-pre-lttrn 7986  ax-pre-apti 7987  ax-pre-ltadd 7988  ax-pre-mulgt0 7989  ax-pre-mulext 7990  ax-arch 7991  ax-caucvg 7992

This theorem depends on definitions:  df-bi 117  df-stab 832  df-dc 836  df-3or 981  df-3an 982  df-tru 1367  df-fal 1370  df-nf 1472  df-sb 1774  df-eu 2045  df-mo 2046  df-clab 2180  df-cleq 2186  df-clel 2189  df-nfc 2325  df-ne 2365  df-nel 2460  df-ral 2477  df-rex 2478  df-reu 2479  df-rmo 2480  df-rab 2481  df-v 2762  df-sbc 2986  df-csb 3081  df-dif 3155  df-un 3157  df-in 3159  df-ss 3166  df-nul 3447  df-if 3558  df-pw 3603  df-sn 3624  df-pr 3625  df-op 3627  df-uni 3836  df-int 3871  df-iun 3914  df-br 4030  df-opab 4091  df-mpt 4092  df-tr 4128  df-id 4324  df-po 4327  df-iso 4328  df-iord 4397  df-on 4399  df-ilim 4400  df-suc 4402  df-iom 4623  df-xp 4665  df-rel 4666  df-cnv 4667  df-co 4668  df-dm 4669  df-rn 4670  df-res 4671  df-ima 4672  df-iota 5215  df-fun 5256  df-fn 5257  df-f 5258  df-f1 5259  df-fo 5260  df-f1o 5261  df-fv 5262  df-isom 5263  df-riota 5873  df-ov 5921  df-oprab 5922  df-mpo 5923  df-1st 6193  df-2nd 6194  df-recs 6358  df-irdg 6423  df-frec 6444  df-1o 6469  df-2o 6470  df-oadd 6473  df-er 6587  df-en 6795  df-dom 6796  df-fin 6797  df-sup 7043  df-inf 7044  df-pnf 8056  df-mnf 8057  df-xr 8058  df-ltxr 8059  df-le 8060  df-sub 8192  df-neg 8193  df-reap 8594  df-ap 8601  df-div 8692  df-inn 8983  df-2 9041  df-3 9042  df-4 9043  df-n0 9241  df-xnn0 9304  df-z 9318  df-uz 9593  df-q 9685  df-rp 9720  df-fz 10075  df-fzo 10209  df-fl 10339  df-mod 10394  df-seqfrec 10519  df-exp 10610  df-ihash 10847  df-cj 10986  df-re 10987  df-im 10988  df-rsqrt 11142  df-abs 11143  df-clim 11422  df-proddc 11694  df-dvds 11931  df-gcd 12080  df-prm 12246  df-odz 12348  df-phi 12349  df-pc 12423

This theorem is referenced by: (None)