/-
Copyright 2026 The Formal Conjectures Authors.
Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at
https://www.apache.org/licenses/LICENSE-2.0
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
-/
import FormalConjecturesUtilErdős Problem 937
[BBC24] Bajpai, P., Bennett, M. A. and Chan, T. H.,
namespace Erdos937
open NatThe four numbers $a, a+d, a+2d, a+3d$ form a four-term arithmetic progression ($d > 0$) of pairwise coprime powerful numbers.
def IsCoprimePowerfulAP4 (a d : ℕ) : Prop :=
0 < d ∧
a.Powerful ∧ (a + d).Powerful ∧ (a + 2 * d).Powerful ∧ (a + 3 * d).Powerful ∧
a.Coprime (a + d) ∧ a.Coprime (a + 2 * d) ∧ a.Coprime (a + 3 * d) ∧
(a + d).Coprime (a + 2 * d) ∧ (a + d).Coprime (a + 3 * d) ∧ (a + 2 * d).Coprime (a + 3 * d)
Are there infinitely many four-term arithmetic progressions of coprime powerful numbers?
(A number $n$ is Nat.Powerful.)
Erdős [Er76d] asked this; the answer is yes: Bajpai, Bennett and Chan [BBC24] proved that
there are infinitely many four-term arithmetic progressions of pairwise coprime powerful numbers.
(Without coprimality this is easy, and by a theorem of Fermat there are no four
@[category research solved, AMS 11]
theorem erdos_937 :
answer(True) ↔ {p : ℕ × ℕ | IsCoprimePowerfulAP4 p.1 p.2}.Infinite := ⊢ True ↔ {p | IsCoprimePowerfulAP4 p.1 p.2}.Infinite
All goals completed! 🐙
Sanity check for IsCoprimePowerfulAP4: the progression $0, 1, 2, 3$ is not a valid
example, since $2$ is not powerful.
@[category test, AMS 11]
theorem not_isCoprimePowerfulAP4_zero_one : ¬ IsCoprimePowerfulAP4 0 1 := ⊢ ¬IsCoprimePowerfulAP4 0 1
h:IsCoprimePowerfulAP4 0 1⊢ False
h2:(0 + 2 * 1).Powerful⊢ False
exact Nat.not_full_of_prime_mod_prime_sq 2 1 Nat.prime_two (h2:(0 + 2 * 1).Powerful⊢ 2 % 2 ^ (1 + 1) = 2 All goals completed! 🐙) h2
end Erdos937