Theorem pockthi 13049

Description: Pocklington's theorem, which gives a sufficient criterion for a number 𝑁 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 13048 for a more general closed-form version. (Contributed by Mario Carneiro, 2-Mar-2014.)

Hypotheses

Ref	Expression
pockthi.p	⊢ 𝑃 ∈ ℙ
pockthi.g	⊢ 𝐺 ∈ ℕ
pockthi.m	⊢ 𝑀 = (𝐺 · 𝑃)
pockthi.n	⊢ 𝑁 = (𝑀 + 1)
pockthi.d	⊢ 𝐷 ∈ ℕ
pockthi.e	⊢ 𝐸 ∈ ℕ
pockthi.a	⊢ 𝐴 ∈ ℕ
pockthi.fac	⊢ 𝑀 = (𝐷 · (𝑃↑𝐸))
pockthi.gt	⊢ 𝐷 < (𝑃↑𝐸)
pockthi.mod	⊢ ((𝐴↑𝑀) mod 𝑁) = (1 mod 𝑁)
pockthi.gcd	⊢ (((𝐴↑𝐺) − 1) gcd 𝑁) = 1

Assertion

Ref	Expression
pockthi	⊢ 𝑁 ∈ ℙ

Proof of Theorem pockthi

Dummy variables 𝑥 𝑦 are mutually distinct and distinct from all other variables.

Step	Hyp	Ref	Expression
1		pockthi.d	. 2 ⊢ 𝐷 ∈ ℕ
2		pockthi.p	. . . . . 6 ⊢ 𝑃 ∈ ℙ
3		prmnn 12800	. . . . . 6 ⊢ (𝑃 ∈ ℙ → 𝑃 ∈ ℕ)
4	2, 3	ax-mp 5	. . . . 5 ⊢ 𝑃 ∈ ℕ
5		pockthi.e	. . . . . 6 ⊢ 𝐸 ∈ ℕ
6	5	nnnn0i 9500	. . . . 5 ⊢ 𝐸 ∈ ℕ0
7		nnexpcl 10910	. . . . 5 ⊢ ((𝑃 ∈ ℕ ∧ 𝐸 ∈ ℕ0) → (𝑃↑𝐸) ∈ ℕ)
8	4, 6, 7	mp2an 426	. . . 4 ⊢ (𝑃↑𝐸) ∈ ℕ
9	8	a1i 9	. . 3 ⊢ (𝐷 ∈ ℕ → (𝑃↑𝐸) ∈ ℕ)
10		id 19	. . 3 ⊢ (𝐷 ∈ ℕ → 𝐷 ∈ ℕ)
11		pockthi.gt	. . . 4 ⊢ 𝐷 < (𝑃↑𝐸)
12	11	a1i 9	. . 3 ⊢ (𝐷 ∈ ℕ → 𝐷 < (𝑃↑𝐸))
13		pockthi.n	. . . . 5 ⊢ 𝑁 = (𝑀 + 1)
14		pockthi.fac	. . . . . . 7 ⊢ 𝑀 = (𝐷 · (𝑃↑𝐸))
15	1	nncni 9243	. . . . . . . 8 ⊢ 𝐷 ∈ ℂ
16	8	nncni 9243	. . . . . . . 8 ⊢ (𝑃↑𝐸) ∈ ℂ
17	15, 16	mulcomi 8276	. . . . . . 7 ⊢ (𝐷 · (𝑃↑𝐸)) = ((𝑃↑𝐸) · 𝐷)
18	14, 17	eqtri 2253	. . . . . 6 ⊢ 𝑀 = ((𝑃↑𝐸) · 𝐷)
19	18	oveq1i 6059	. . . . 5 ⊢ (𝑀 + 1) = (((𝑃↑𝐸) · 𝐷) + 1)
20	13, 19	eqtri 2253	. . . 4 ⊢ 𝑁 = (((𝑃↑𝐸) · 𝐷) + 1)
21	20	a1i 9	. . 3 ⊢ (𝐷 ∈ ℕ → 𝑁 = (((𝑃↑𝐸) · 𝐷) + 1))
22		prmdvdsexpb 12839	. . . . . . 7 ⊢ ((𝑥 ∈ ℙ ∧ 𝑃 ∈ ℙ ∧ 𝐸 ∈ ℕ) → (𝑥 ∥ (𝑃↑𝐸) ↔ 𝑥 = 𝑃))
23	2, 5, 22	mp3an23 1366	. . . . . 6 ⊢ (𝑥 ∈ ℙ → (𝑥 ∥ (𝑃↑𝐸) ↔ 𝑥 = 𝑃))
24		pockthi.m	. . . . . . . . . . . . 13 ⊢ 𝑀 = (𝐺 · 𝑃)
25		pockthi.g	. . . . . . . . . . . . . 14 ⊢ 𝐺 ∈ ℕ
26	25, 4	nnmulcli 9255	. . . . . . . . . . . . 13 ⊢ (𝐺 · 𝑃) ∈ ℕ
27	24, 26	eqeltri 2305	. . . . . . . . . . . 12 ⊢ 𝑀 ∈ ℕ
28	27	nncni 9243	. . . . . . . . . . 11 ⊢ 𝑀 ∈ ℂ
29		ax-1cn 8216	. . . . . . . . . . 11 ⊢ 1 ∈ ℂ
30	28, 29, 13	mvrraddi 8486	. . . . . . . . . 10 ⊢ (𝑁 − 1) = 𝑀
31	30	oveq2i 6060	. . . . . . . . 9 ⊢ (𝐴↑(𝑁 − 1)) = (𝐴↑𝑀)
32	31	oveq1i 6059	. . . . . . . 8 ⊢ ((𝐴↑(𝑁 − 1)) mod 𝑁) = ((𝐴↑𝑀) mod 𝑁)
33		pockthi.mod	. . . . . . . . 9 ⊢ ((𝐴↑𝑀) mod 𝑁) = (1 mod 𝑁)
34		peano2nn 9245	. . . . . . . . . . . . 13 ⊢ (𝑀 ∈ ℕ → (𝑀 + 1) ∈ ℕ)
35	27, 34	ax-mp 5	. . . . . . . . . . . 12 ⊢ (𝑀 + 1) ∈ ℕ
36	13, 35	eqeltri 2305	. . . . . . . . . . 11 ⊢ 𝑁 ∈ ℕ
37		nnq 9961	. . . . . . . . . . 11 ⊢ (𝑁 ∈ ℕ → 𝑁 ∈ ℚ)
38	36, 37	ax-mp 5	. . . . . . . . . 10 ⊢ 𝑁 ∈ ℚ
39	27	nngt0i 9263	. . . . . . . . . . . 12 ⊢ 0 < 𝑀
40	27	nnrei 9242	. . . . . . . . . . . . 13 ⊢ 𝑀 ∈ ℝ
41		1re 8269	. . . . . . . . . . . . 13 ⊢ 1 ∈ ℝ
42		ltaddpos2 8723	. . . . . . . . . . . . 13 ⊢ ((𝑀 ∈ ℝ ∧ 1 ∈ ℝ) → (0 < 𝑀 ↔ 1 < (𝑀 + 1)))
43	40, 41, 42	mp2an 426	. . . . . . . . . . . 12 ⊢ (0 < 𝑀 ↔ 1 < (𝑀 + 1))
44	39, 43	mpbi 145	. . . . . . . . . . 11 ⊢ 1 < (𝑀 + 1)
45	44, 13	breqtrri 4135	. . . . . . . . . 10 ⊢ 1 < 𝑁
46		q1mod 10714	. . . . . . . . . 10 ⊢ ((𝑁 ∈ ℚ ∧ 1 < 𝑁) → (1 mod 𝑁) = 1)
47	38, 45, 46	mp2an 426	. . . . . . . . 9 ⊢ (1 mod 𝑁) = 1
48	33, 47	eqtri 2253	. . . . . . . 8 ⊢ ((𝐴↑𝑀) mod 𝑁) = 1
49	32, 48	eqtri 2253	. . . . . . 7 ⊢ ((𝐴↑(𝑁 − 1)) mod 𝑁) = 1
50		oveq2 6057	. . . . . . . . . . . 12 ⊢ (𝑥 = 𝑃 → ((𝑁 − 1) / 𝑥) = ((𝑁 − 1) / 𝑃))
51	25	nncni 9243	. . . . . . . . . . . . . . 15 ⊢ 𝐺 ∈ ℂ
52	4	nncni 9243	. . . . . . . . . . . . . . 15 ⊢ 𝑃 ∈ ℂ
53	51, 52	mulcomi 8276	. . . . . . . . . . . . . 14 ⊢ (𝐺 · 𝑃) = (𝑃 · 𝐺)
54	30, 24, 53	3eqtrri 2258	. . . . . . . . . . . . 13 ⊢ (𝑃 · 𝐺) = (𝑁 − 1)
55	36	nncni 9243	. . . . . . . . . . . . . . 15 ⊢ 𝑁 ∈ ℂ
56	55, 29	subcli 8545	. . . . . . . . . . . . . 14 ⊢ (𝑁 − 1) ∈ ℂ
57	4	nnap0i 9264	. . . . . . . . . . . . . 14 ⊢ 𝑃 # 0
58	56, 52, 51, 57	divmulapi 9036	. . . . . . . . . . . . 13 ⊢ (((𝑁 − 1) / 𝑃) = 𝐺 ↔ (𝑃 · 𝐺) = (𝑁 − 1))
59	54, 58	mpbir 146	. . . . . . . . . . . 12 ⊢ ((𝑁 − 1) / 𝑃) = 𝐺
60	50, 59	eqtrdi 2281	. . . . . . . . . . 11 ⊢ (𝑥 = 𝑃 → ((𝑁 − 1) / 𝑥) = 𝐺)
61	60	oveq2d 6065	. . . . . . . . . 10 ⊢ (𝑥 = 𝑃 → (𝐴↑((𝑁 − 1) / 𝑥)) = (𝐴↑𝐺))
62	61	oveq1d 6064	. . . . . . . . 9 ⊢ (𝑥 = 𝑃 → ((𝐴↑((𝑁 − 1) / 𝑥)) − 1) = ((𝐴↑𝐺) − 1))
63	62	oveq1d 6064	. . . . . . . 8 ⊢ (𝑥 = 𝑃 → (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = (((𝐴↑𝐺) − 1) gcd 𝑁))
64		pockthi.gcd	. . . . . . . 8 ⊢ (((𝐴↑𝐺) − 1) gcd 𝑁) = 1
65	63, 64	eqtrdi 2281	. . . . . . 7 ⊢ (𝑥 = 𝑃 → (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)
66		pockthi.a	. . . . . . . . 9 ⊢ 𝐴 ∈ ℕ
67	66	nnzi 9594	. . . . . . . 8 ⊢ 𝐴 ∈ ℤ
68		oveq1 6056	. . . . . . . . . . . 12 ⊢ (𝑦 = 𝐴 → (𝑦↑(𝑁 − 1)) = (𝐴↑(𝑁 − 1)))
69	68	oveq1d 6064	. . . . . . . . . . 11 ⊢ (𝑦 = 𝐴 → ((𝑦↑(𝑁 − 1)) mod 𝑁) = ((𝐴↑(𝑁 − 1)) mod 𝑁))
70	69	eqeq1d 2241	. . . . . . . . . 10 ⊢ (𝑦 = 𝐴 → (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ↔ ((𝐴↑(𝑁 − 1)) mod 𝑁) = 1))
71		oveq1 6056	. . . . . . . . . . . . 13 ⊢ (𝑦 = 𝐴 → (𝑦↑((𝑁 − 1) / 𝑥)) = (𝐴↑((𝑁 − 1) / 𝑥)))
72	71	oveq1d 6064	. . . . . . . . . . . 12 ⊢ (𝑦 = 𝐴 → ((𝑦↑((𝑁 − 1) / 𝑥)) − 1) = ((𝐴↑((𝑁 − 1) / 𝑥)) − 1))
73	72	oveq1d 6064	. . . . . . . . . . 11 ⊢ (𝑦 = 𝐴 → (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁))
74	73	eqeq1d 2241	. . . . . . . . . 10 ⊢ (𝑦 = 𝐴 → ((((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1 ↔ (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
75	70, 74	anbi12d 473	. . . . . . . . 9 ⊢ (𝑦 = 𝐴 → ((((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1) ↔ (((𝐴↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)))
76	75	rspcev 2920	. . . . . . . 8 ⊢ ((𝐴 ∈ ℤ ∧ (((𝐴↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
77	67, 76	mpan 424	. . . . . . 7 ⊢ ((((𝐴↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
78	49, 65, 77	sylancr 414	. . . . . 6 ⊢ (𝑥 = 𝑃 → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
79	23, 78	biimtrdi 163	. . . . 5 ⊢ (𝑥 ∈ ℙ → (𝑥 ∥ (𝑃↑𝐸) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)))
80	79	rgen 2595	. . . 4 ⊢ ∀𝑥 ∈ ℙ (𝑥 ∥ (𝑃↑𝐸) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
81	80	a1i 9	. . 3 ⊢ (𝐷 ∈ ℕ → ∀𝑥 ∈ ℙ (𝑥 ∥ (𝑃↑𝐸) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)))
82	9, 10, 12, 21, 81	pockthg 13048	. 2 ⊢ (𝐷 ∈ ℕ → 𝑁 ∈ ℙ)
83	1, 82	ax-mp 5	1 ⊢ 𝑁 ∈ ℙ

Syntax hints:   → wi 4   ∧ wa 104   ↔ wb 105   = wceq 1398   ∈ wcel 2203  ∀wral 2520  ∃wrex 2521   class class class wbr 4108  (class class class)co 6049  ℝcr 8122  0cc0 8123  1c1 8124   + caddc 8126   · cmul 8128   < clt 8304   − cmin 8440   / cdiv 8942  ℕcn 9233  ℕ₀cn0 9492  ℤcz 9573  ℚcq 9947   mod cmo 10680  ↑cexp 10896   ∥ cdvds 12466   gcd cgcd 12642  ℙcprime 12797

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 619  ax-in2 620  ax-io 717  ax-5 1496  ax-7 1497  ax-gen 1498  ax-ie1 1542  ax-ie2 1543  ax-8 1553  ax-10 1554  ax-11 1555  ax-i12 1556  ax-bndl 1558  ax-4 1559  ax-17 1575  ax-i9 1579  ax-ial 1583  ax-i5r 1584  ax-13 2205  ax-14 2206  ax-ext 2214  ax-coll 4224  ax-sep 4227  ax-nul 4235  ax-pow 4286  ax-pr 4321  ax-un 4553  ax-setind 4658  ax-iinf 4709  ax-cnex 8214  ax-resscn 8215  ax-1cn 8216  ax-1re 8217  ax-icn 8218  ax-addcl 8219  ax-addrcl 8220  ax-mulcl 8221  ax-mulrcl 8222  ax-addcom 8223  ax-mulcom 8224  ax-addass 8225  ax-mulass 8226  ax-distr 8227  ax-i2m1 8228  ax-0lt1 8229  ax-1rid 8230  ax-0id 8231  ax-rnegex 8232  ax-precex 8233  ax-cnre 8234  ax-pre-ltirr 8235  ax-pre-ltwlin 8236  ax-pre-lttrn 8237  ax-pre-apti 8238  ax-pre-ltadd 8239  ax-pre-mulgt0 8240  ax-pre-mulext 8241  ax-arch 8242  ax-caucvg 8243

This theorem depends on definitions:  df-bi 117  df-stab 839  df-dc 843  df-3or 1006  df-3an 1007  df-tru 1401  df-fal 1404  df-nf 1510  df-sb 1812  df-eu 2083  df-mo 2084  df-clab 2219  df-cleq 2225  df-clel 2228  df-nfc 2373  df-ne 2413  df-nel 2508  df-ral 2525  df-rex 2526  df-reu 2527  df-rmo 2528  df-rab 2529  df-v 2814  df-sbc 3042  df-csb 3138  df-dif 3212  df-un 3214  df-in 3216  df-ss 3223  df-nul 3508  df-if 3620  df-pw 3670  df-sn 3694  df-pr 3695  df-op 3697  df-uni 3914  df-int 3949  df-iun 3992  df-br 4109  df-opab 4171  df-mpt 4172  df-tr 4208  df-id 4413  df-po 4416  df-iso 4417  df-iord 4486  df-on 4488  df-ilim 4489  df-suc 4491  df-iom 4712  df-xp 4754  df-rel 4755  df-cnv 4756  df-co 4757  df-dm 4758  df-rn 4759  df-res 4760  df-ima 4761  df-iota 5311  df-fun 5353  df-fn 5354  df-f 5355  df-f1 5356  df-fo 5357  df-f1o 5358  df-fv 5359  df-isom 5360  df-riota 6002  df-ov 6052  df-oprab 6053  df-mpo 6054  df-1st 6333  df-2nd 6334  df-recs 6535  df-irdg 6600  df-frec 6621  df-1o 6646  df-2o 6647  df-oadd 6650  df-er 6766  df-en 6975  df-dom 6976  df-fin 6977  df-sup 7274  df-inf 7275  df-pnf 8306  df-mnf 8307  df-xr 8308  df-ltxr 8309  df-le 8310  df-sub 8442  df-neg 8443  df-reap 8845  df-ap 8852  df-div 8943  df-inn 9234  df-2 9292  df-3 9293  df-4 9294  df-n0 9493  df-xnn0 9560  df-z 9574  df-uz 9850  df-q 9948  df-rp 9983  df-fz 10339  df-fzo 10473  df-fl 10626  df-mod 10681  df-seqfrec 10806  df-exp 10897  df-ihash 11134  df-cj 11520  df-re 11521  df-im 11522  df-rsqrt 11676  df-abs 11677  df-clim 11957  df-proddc 12230  df-dvds 12467  df-gcd 12643  df-prm 12798  df-odz 12900  df-phi 12901  df-pc 12976

This theorem is referenced by: (None)