Theorem pockthi 16536

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 16535 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 16307	. . . . . 6 ⊢ (𝑃 ∈ ℙ → 𝑃 ∈ ℕ)
4	2, 3	ax-mp 5	. . . . 5 ⊢ 𝑃 ∈ ℕ
5		pockthi.e	. . . . . 6 ⊢ 𝐸 ∈ ℕ
6	5	nnnn0i 12171	. . . . 5 ⊢ 𝐸 ∈ ℕ0
7		nnexpcl 13723	. . . . 5 ⊢ ((𝑃 ∈ ℕ ∧ 𝐸 ∈ ℕ0) → (𝑃↑𝐸) ∈ ℕ)
8	4, 6, 7	mp2an 688	. . . 4 ⊢ (𝑃↑𝐸) ∈ ℕ
9	8	a1i 11	. . 3 ⊢ (𝐷 ∈ ℕ → (𝑃↑𝐸) ∈ ℕ)
10		id 22	. . 3 ⊢ (𝐷 ∈ ℕ → 𝐷 ∈ ℕ)
11		pockthi.gt	. . . 4 ⊢ 𝐷 < (𝑃↑𝐸)
12	11	a1i 11	. . 3 ⊢ (𝐷 ∈ ℕ → 𝐷 < (𝑃↑𝐸))
13		pockthi.n	. . . . 5 ⊢ 𝑁 = (𝑀 + 1)
14		pockthi.fac	. . . . . . 7 ⊢ 𝑀 = (𝐷 · (𝑃↑𝐸))
15	1	nncni 11913	. . . . . . . 8 ⊢ 𝐷 ∈ ℂ
16	8	nncni 11913	. . . . . . . 8 ⊢ (𝑃↑𝐸) ∈ ℂ
17	15, 16	mulcomi 10914	. . . . . . 7 ⊢ (𝐷 · (𝑃↑𝐸)) = ((𝑃↑𝐸) · 𝐷)
18	14, 17	eqtri 2766	. . . . . 6 ⊢ 𝑀 = ((𝑃↑𝐸) · 𝐷)
19	18	oveq1i 7265	. . . . 5 ⊢ (𝑀 + 1) = (((𝑃↑𝐸) · 𝐷) + 1)
20	13, 19	eqtri 2766	. . . 4 ⊢ 𝑁 = (((𝑃↑𝐸) · 𝐷) + 1)
21	20	a1i 11	. . 3 ⊢ (𝐷 ∈ ℕ → 𝑁 = (((𝑃↑𝐸) · 𝐷) + 1))
22		prmdvdsexpb 16349	. . . . . . 7 ⊢ ((𝑥 ∈ ℙ ∧ 𝑃 ∈ ℙ ∧ 𝐸 ∈ ℕ) → (𝑥 ∥ (𝑃↑𝐸) ↔ 𝑥 = 𝑃))
23	2, 5, 22	mp3an23 1451	. . . . . 6 ⊢ (𝑥 ∈ ℙ → (𝑥 ∥ (𝑃↑𝐸) ↔ 𝑥 = 𝑃))
24		pockthi.m	. . . . . . . . . . . . 13 ⊢ 𝑀 = (𝐺 · 𝑃)
25		pockthi.g	. . . . . . . . . . . . . 14 ⊢ 𝐺 ∈ ℕ
26	25, 4	nnmulcli 11928	. . . . . . . . . . . . 13 ⊢ (𝐺 · 𝑃) ∈ ℕ
27	24, 26	eqeltri 2835	. . . . . . . . . . . 12 ⊢ 𝑀 ∈ ℕ
28	27	nncni 11913	. . . . . . . . . . 11 ⊢ 𝑀 ∈ ℂ
29		ax-1cn 10860	. . . . . . . . . . 11 ⊢ 1 ∈ ℂ
30	28, 29, 13	mvrraddi 11168	. . . . . . . . . 10 ⊢ (𝑁 − 1) = 𝑀
31	30	oveq2i 7266	. . . . . . . . 9 ⊢ (𝐴↑(𝑁 − 1)) = (𝐴↑𝑀)
32	31	oveq1i 7265	. . . . . . . 8 ⊢ ((𝐴↑(𝑁 − 1)) mod 𝑁) = ((𝐴↑𝑀) mod 𝑁)
33		pockthi.mod	. . . . . . . . 9 ⊢ ((𝐴↑𝑀) mod 𝑁) = (1 mod 𝑁)
34		peano2nn 11915	. . . . . . . . . . . . 13 ⊢ (𝑀 ∈ ℕ → (𝑀 + 1) ∈ ℕ)
35	27, 34	ax-mp 5	. . . . . . . . . . . 12 ⊢ (𝑀 + 1) ∈ ℕ
36	13, 35	eqeltri 2835	. . . . . . . . . . 11 ⊢ 𝑁 ∈ ℕ
37	36	nnrei 11912	. . . . . . . . . 10 ⊢ 𝑁 ∈ ℝ
38	27	nngt0i 11942	. . . . . . . . . . . 12 ⊢ 0 < 𝑀
39	27	nnrei 11912	. . . . . . . . . . . . 13 ⊢ 𝑀 ∈ ℝ
40		1re 10906	. . . . . . . . . . . . 13 ⊢ 1 ∈ ℝ
41		ltaddpos2 11396	. . . . . . . . . . . . 13 ⊢ ((𝑀 ∈ ℝ ∧ 1 ∈ ℝ) → (0 < 𝑀 ↔ 1 < (𝑀 + 1)))
42	39, 40, 41	mp2an 688	. . . . . . . . . . . 12 ⊢ (0 < 𝑀 ↔ 1 < (𝑀 + 1))
43	38, 42	mpbi 229	. . . . . . . . . . 11 ⊢ 1 < (𝑀 + 1)
44	43, 13	breqtrri 5097	. . . . . . . . . 10 ⊢ 1 < 𝑁
45		1mod 13551	. . . . . . . . . 10 ⊢ ((𝑁 ∈ ℝ ∧ 1 < 𝑁) → (1 mod 𝑁) = 1)
46	37, 44, 45	mp2an 688	. . . . . . . . 9 ⊢ (1 mod 𝑁) = 1
47	33, 46	eqtri 2766	. . . . . . . 8 ⊢ ((𝐴↑𝑀) mod 𝑁) = 1
48	32, 47	eqtri 2766	. . . . . . 7 ⊢ ((𝐴↑(𝑁 − 1)) mod 𝑁) = 1
49		oveq2 7263	. . . . . . . . . . . 12 ⊢ (𝑥 = 𝑃 → ((𝑁 − 1) / 𝑥) = ((𝑁 − 1) / 𝑃))
50	25	nncni 11913	. . . . . . . . . . . . . . 15 ⊢ 𝐺 ∈ ℂ
51	4	nncni 11913	. . . . . . . . . . . . . . 15 ⊢ 𝑃 ∈ ℂ
52	50, 51	mulcomi 10914	. . . . . . . . . . . . . 14 ⊢ (𝐺 · 𝑃) = (𝑃 · 𝐺)
53	30, 24, 52	3eqtrri 2771	. . . . . . . . . . . . 13 ⊢ (𝑃 · 𝐺) = (𝑁 − 1)
54	36	nncni 11913	. . . . . . . . . . . . . . 15 ⊢ 𝑁 ∈ ℂ
55	54, 29	subcli 11227	. . . . . . . . . . . . . 14 ⊢ (𝑁 − 1) ∈ ℂ
56	4	nnne0i 11943	. . . . . . . . . . . . . 14 ⊢ 𝑃 ≠ 0
57	55, 51, 50, 56	divmuli 11659	. . . . . . . . . . . . 13 ⊢ (((𝑁 − 1) / 𝑃) = 𝐺 ↔ (𝑃 · 𝐺) = (𝑁 − 1))
58	53, 57	mpbir 230	. . . . . . . . . . . 12 ⊢ ((𝑁 − 1) / 𝑃) = 𝐺
59	49, 58	eqtrdi 2795	. . . . . . . . . . 11 ⊢ (𝑥 = 𝑃 → ((𝑁 − 1) / 𝑥) = 𝐺)
60	59	oveq2d 7271	. . . . . . . . . 10 ⊢ (𝑥 = 𝑃 → (𝐴↑((𝑁 − 1) / 𝑥)) = (𝐴↑𝐺))
61	60	oveq1d 7270	. . . . . . . . 9 ⊢ (𝑥 = 𝑃 → ((𝐴↑((𝑁 − 1) / 𝑥)) − 1) = ((𝐴↑𝐺) − 1))
62	61	oveq1d 7270	. . . . . . . 8 ⊢ (𝑥 = 𝑃 → (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = (((𝐴↑𝐺) − 1) gcd 𝑁))
63		pockthi.gcd	. . . . . . . 8 ⊢ (((𝐴↑𝐺) − 1) gcd 𝑁) = 1
64	62, 63	eqtrdi 2795	. . . . . . 7 ⊢ (𝑥 = 𝑃 → (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)
65		pockthi.a	. . . . . . . . 9 ⊢ 𝐴 ∈ ℕ
66	65	nnzi 12274	. . . . . . . 8 ⊢ 𝐴 ∈ ℤ
67		oveq1 7262	. . . . . . . . . . . 12 ⊢ (𝑦 = 𝐴 → (𝑦↑(𝑁 − 1)) = (𝐴↑(𝑁 − 1)))
68	67	oveq1d 7270	. . . . . . . . . . 11 ⊢ (𝑦 = 𝐴 → ((𝑦↑(𝑁 − 1)) mod 𝑁) = ((𝐴↑(𝑁 − 1)) mod 𝑁))
69	68	eqeq1d 2740	. . . . . . . . . 10 ⊢ (𝑦 = 𝐴 → (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ↔ ((𝐴↑(𝑁 − 1)) mod 𝑁) = 1))
70		oveq1 7262	. . . . . . . . . . . . 13 ⊢ (𝑦 = 𝐴 → (𝑦↑((𝑁 − 1) / 𝑥)) = (𝐴↑((𝑁 − 1) / 𝑥)))
71	70	oveq1d 7270	. . . . . . . . . . . 12 ⊢ (𝑦 = 𝐴 → ((𝑦↑((𝑁 − 1) / 𝑥)) − 1) = ((𝐴↑((𝑁 − 1) / 𝑥)) − 1))
72	71	oveq1d 7270	. . . . . . . . . . 11 ⊢ (𝑦 = 𝐴 → (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁))
73	72	eqeq1d 2740	. . . . . . . . . 10 ⊢ (𝑦 = 𝐴 → ((((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1 ↔ (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
74	69, 73	anbi12d 630	. . . . . . . . 9 ⊢ (𝑦 = 𝐴 → ((((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1) ↔ (((𝐴↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)))
75	74	rspcev 3552	. . . . . . . 8 ⊢ ((𝐴 ∈ ℤ ∧ (((𝐴↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
76	66, 75	mpan 686	. . . . . . 7 ⊢ ((((𝐴↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝐴↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
77	48, 64, 76	sylancr 586	. . . . . 6 ⊢ (𝑥 = 𝑃 → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
78	23, 77	syl6bi 252	. . . . 5 ⊢ (𝑥 ∈ ℙ → (𝑥 ∥ (𝑃↑𝐸) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)))
79	78	rgen 3073	. . . 4 ⊢ ∀𝑥 ∈ ℙ (𝑥 ∥ (𝑃↑𝐸) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1))
80	79	a1i 11	. . 3 ⊢ (𝐷 ∈ ℕ → ∀𝑥 ∈ ℙ (𝑥 ∥ (𝑃↑𝐸) → ∃𝑦 ∈ ℤ (((𝑦↑(𝑁 − 1)) mod 𝑁) = 1 ∧ (((𝑦↑((𝑁 − 1) / 𝑥)) − 1) gcd 𝑁) = 1)))
81	9, 10, 12, 21, 80	pockthg 16535	. 2 ⊢ (𝐷 ∈ ℕ → 𝑁 ∈ ℙ)
82	1, 81	ax-mp 5	1 ⊢ 𝑁 ∈ ℙ

Syntax hints:   → wi 4   ↔ wb 205   ∧ wa 395   = wceq 1539   ∈ wcel 2108  ∀wral 3063  ∃wrex 3064   class class class wbr 5070  (class class class)co 7255  ℝcr 10801  0cc0 10802  1c1 10803   + caddc 10805   · cmul 10807   < clt 10940   − cmin 11135   / cdiv 11562  ℕcn 11903  ℕ₀cn0 12163  ℤcz 12249   mod cmo 13517  ↑cexp 13710   ∥ cdvds 15891   gcd cgcd 16129  ℙcprime 16304

This theorem was proved from axioms:  ax-mp 5  ax-1 6  ax-2 7  ax-3 8  ax-gen 1799  ax-4 1813  ax-5 1914  ax-6 1972  ax-7 2012  ax-8 2110  ax-9 2118  ax-10 2139  ax-11 2156  ax-12 2173  ax-ext 2709  ax-rep 5205  ax-sep 5218  ax-nul 5225  ax-pow 5283  ax-pr 5347  ax-un 7566  ax-cnex 10858  ax-resscn 10859  ax-1cn 10860  ax-icn 10861  ax-addcl 10862  ax-addrcl 10863  ax-mulcl 10864  ax-mulrcl 10865  ax-mulcom 10866  ax-addass 10867  ax-mulass 10868  ax-distr 10869  ax-i2m1 10870  ax-1ne0 10871  ax-1rid 10872  ax-rnegex 10873  ax-rrecex 10874  ax-cnre 10875  ax-pre-lttri 10876  ax-pre-lttrn 10877  ax-pre-ltadd 10878  ax-pre-mulgt0 10879  ax-pre-sup 10880

This theorem depends on definitions:  df-bi 206  df-an 396  df-or 844  df-3or 1086  df-3an 1087  df-tru 1542  df-fal 1552  df-ex 1784  df-nf 1788  df-sb 2069  df-mo 2540  df-eu 2569  df-clab 2716  df-cleq 2730  df-clel 2817  df-nfc 2888  df-ne 2943  df-nel 3049  df-ral 3068  df-rex 3069  df-reu 3070  df-rmo 3071  df-rab 3072  df-v 3424  df-sbc 3712  df-csb 3829  df-dif 3886  df-un 3888  df-in 3890  df-ss 3900  df-pss 3902  df-nul 4254  df-if 4457  df-pw 4532  df-sn 4559  df-pr 4561  df-tp 4563  df-op 4565  df-uni 4837  df-int 4877  df-iun 4923  df-br 5071  df-opab 5133  df-mpt 5154  df-tr 5188  df-id 5480  df-eprel 5486  df-po 5494  df-so 5495  df-fr 5535  df-we 5537  df-xp 5586  df-rel 5587  df-cnv 5588  df-co 5589  df-dm 5590  df-rn 5591  df-res 5592  df-ima 5593  df-pred 6191  df-ord 6254  df-on 6255  df-lim 6256  df-suc 6257  df-iota 6376  df-fun 6420  df-fn 6421  df-f 6422  df-f1 6423  df-fo 6424  df-f1o 6425  df-fv 6426  df-riota 7212  df-ov 7258  df-oprab 7259  df-mpo 7260  df-om 7688  df-1st 7804  df-2nd 7805  df-frecs 8068  df-wrecs 8099  df-recs 8173  df-rdg 8212  df-1o 8267  df-2o 8268  df-oadd 8271  df-er 8456  df-map 8575  df-en 8692  df-dom 8693  df-sdom 8694  df-fin 8695  df-sup 9131  df-inf 9132  df-dju 9590  df-card 9628  df-pnf 10942  df-mnf 10943  df-xr 10944  df-ltxr 10945  df-le 10946  df-sub 11137  df-neg 11138  df-div 11563  df-nn 11904  df-2 11966  df-3 11967  df-n0 12164  df-xnn0 12236  df-z 12250  df-uz 12512  df-q 12618  df-rp 12660  df-fz 13169  df-fzo 13312  df-fl 13440  df-mod 13518  df-seq 13650  df-exp 13711  df-hash 13973  df-cj 14738  df-re 14739  df-im 14740  df-sqrt 14874  df-abs 14875  df-dvds 15892  df-gcd 16130  df-prm 16305  df-odz 16394  df-phi 16395  df-pc 16466

This theorem is referenced by:  1259prm  16765  2503prm  16769  4001prm  16774