Preprint
Article

This version is not peer-reviewed.

From the Pythagorean Dream to the Fermatian Obstruction: The Unified Chain

Submitted:

05 March 2026

Posted:

06 March 2026

You are already at the latest version

Abstract

For every prime p and every integer a, the backward finite difference δp(a) := aᵖ − (a − 1)ᵖ equals the cyclotomic binary form Φp(a, a − 1) and hence the norm N_Q(ζp)/Q(a − ζp(a − 1)). For p = 3 this specialises to δ3(a) = N_Z[ω](a − ω(a − 1)), connecting the individual cubic finite difference — obtained by differencing the classical sum formula of Nicomachus of Gerasa (∼100 CE) — with the Eisenstein norm that appears in Euler's factorisation of a³ + b³.Starting from the historical identity S3(n) = Tₙ² where Tₙ = n(n + 1)/2, and applying the backward finite difference operator ∇f(n) := f(n) − f(n − 1) — formalised by Taylor (1715) and systematised by Boole (1860) — the Cubic Identity is derived: n³ = (n²/4)[(n + 1)² − (n − 1)²] = Tₙ² − Tₙ₋₁².This identity is extended to all p ≥ 1 via the Universal Faulhaber–Bernoulli Identity (UFBI): nᵖ = 1/(p+1) Σⱼ₌₀ᵖ C(p+1,j) Bⱼ⁺ δp+1−j(n), δm(n) := nᵐ − (n−1)ᵐ.The central contribution of this work is the Unified Chain Formula: ∇Tₙ² = δ3(a) = N_Z[ω](a − ω(a−1)) = Φ3(a, a−1) = N_Q(ζ3)/Q(a − ζ3(a−1)), which connects, in a single proved identity, five centuries of mathematics: Nicomachus (1st century), Boole (19th century), Euler/Eisenstein (18th century), and Gauss/cyclotomic theory (19th–20th centuries). This chain is not present as such in the existing literature; its originality lies in the explicit articulation of these connections, not in the individual equalities, each of which follows from classical results.Beyond the Unified Chain, the following new elements are introduced: (i) the Tower of Norms a³ = Σₖ₌₁ᵃ N(αk), making explicit how each perfect cube is a stack of hexagonal norms; (ii) the Cyclotomic Compatibility Index ICC(n, p), which quantifies the arithmetic obstruction to hᵖ = aᵖ + bᵖ having integer solutions; (iii) the Window Incompatibility Theorem, formalising why the hexagonal windows {a−1, a, a+1} and {b−1, b, b+1} can never merge into a single window {h−1, h, h+1} in Z[ω] for a, b ≥ 2; (iv) the Order Theorem for δm(n), providing a complete characterisation of prime divisibility of finite differences via multiplicative orders; (v) the Extreme Reduction Theorem (ERT), showing that the Order Filter eliminates every pair (a, b) with a ≥ 2 from the equation a³ + b³ = c³, reducing the problem to the case a = 1; (vi) the Fermatian Rigidity Index R(p), a quantitative measure of how far (aᵖ + bᵖ)^(1/p) is from an integer. All results are illustrated throughout by the single running example a = 6, b = 10, and the key number 91 = 7 × 13. Verified over 179,700 pairs with 50-digit precision: zero exceptions. This work does not claim to prove Fermat's Last Theorem, definitively established by Wiles [1].

Keywords: 
;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  

1. Introduction

1.1. Historical Context

The Pythagorean school, founded in Croton (Magna Graecia) around the 6th century BC, established one of the first philosophical-mathematical syntheses in Western history: the belief that “everything is number” ( π α ν τ α α ρ ι θ μ o ς ). For the Pythagoreans, numbers were not mere calculation tools but ontological principles structuring cosmic reality. This worldview found its purest expression in figurate numbers: geometric representations of integers by means of regular point configurations.
The Pythagorean theorem crystallised the harmony between geometry and arithmetic:
h 2 = a 2 + b 2 ,
where h denotes the hypotenuse of a right triangle and a, b its legs. Euclid provided the complete parameterisation of primitive Pythagorean triples:
a = m 2 n 2 , b = 2 m n , h = m 2 + n 2 ,
with m > n > 0 , gcd ( m , n ) = 1 , m n ( mod 2 ) .
Nicomachus of Gerasa (∼60–120 AD) systematised in his Introductio Arithmetica the classical Greek arithmetic knowledge. He established the cumulative sum of cubes:
S 3 ( n ) = k = 1 n k 3 = T n 2 , T n = n ( n + 1 ) 2 .
It is essential to note that Nicomachus did not formulate nor suggest any identity involving symmetric differences ( n + 1 ) 2 ( n 1 ) 2 ; his contribution was limited exclusively to the cumulative sum formula. Attributing the Cubic Identity (derived in Section 3) to him would be an inadmissible historiographic anachronism.
Johann Faulhaber (1580–1635) computed S p ( n ) = k = 1 n k p for p up to at least 17 in his Academia Algebrae (1631) [14], discovering the polynomial nature of these sums and observing that S 3 ( n ) = T n 2 is the only case with a perfect quadratic factorisation.
Jakob Bernoulli (1654–1705) published posthumously in Ars Conjectandi (1713) [8] the general formula for S p ( n ) in terms of the coefficients that today bear his name. The Bernoulli numbers B j vanish for all odd j 3 ; this is a fundamental fact for the Universal Identity.
Pierre de Fermat (1607–1665) wrote in 1637, in the margin of his copy of Diophantus’s Arithmetica: “Cubum autem in duos cubos, aut quadratoquadratum in duos quadratoquadratos, et generaliter nullam in infinitum ultra quadratum potestatem in duos eiusdem nominis fas est dividere cuius rei demonstrationem mirabilem sane detexi. Hanc marginis exiguitas non caperet.” Most historians consider that his “marvellous proof” was incomplete.
Leonhard Euler (1707–1783) provided the first proof of FLT for p = 3 (∼1770) by infinite descent in the ring Z [ ω ] of Eisenstein integers, where ω = e 2 π i / 3 . The key factorisation x 3 + y 3 = ( x + y ) ( x + ω y ) ( x + ω 2 y ) = z 3 requires unique factorisation in Z [ ω ] , a property Euler assumed tacitly and that was rigorously established by Gauss and Kummer [2].
Sophie Germain (1776–1831) developed the first systematic framework for attacking FLT via modular congruences, verifying conditions for all primes p < 100 . She adopted the masculine pseudonym “M. Le Blanc” to correspond with Lagrange and Gauss. Gauss wrote in 1807: “When a person of the sex which, according to our customs and prejudices, must encounter infinitely more difficulties than men to familiarise herself with these thorny researches, succeeds nevertheless in surmounting these obstacles and penetrating their most obscure parts, she must without doubt have the noblest courage, quite extraordinary talents, and superior genius.”
Ernst Kummer (1810–1893) introduced the revolutionary concept of ideal numbers and proved FLT for all regular primes.
Kenneth Ribet proved in 1986 that Frey’s elliptic curve (associated with a hypothetical counterexample to FLT) would violate the Taniyama–Shimura–Weil conjecture, reducing FLT to a problem in the theory of modular forms.
Andrew Wiles (1953–) proved FLT in full generality in 1994 (published in 1995) [1] via the Modularity Theorem for elliptic curves. Wiles’ proof unifies algebraic geometry, representation theory, complex analysis, and number theory.
Table 1 summarises the chronological development.

1.2. The Two Traditions and the Bridge

The sum formula of Nicomachus,
S 3 ( n ) : = k = 1 n k 3 = T n 2 , T n = n ( n + 1 ) 2 ,
is one of the oldest identities in number theory [3]. Applying the backward-difference operator ∇ yields the individual cubic difference δ 3 ( a ) = 3 a 2 3 a + 1 , a formula within reach of any secondary-school student.
More than sixteen centuries later, Euler’s factorisation of a 3 + b 3 over the ring of Eisenstein integers Z [ ω ] employs the norm of the element a ω ( a 1 ) [2]. The central observation of this work is that these two objects are the same:
δ 3 ( a ) = N Z [ ω ] a ω ( a 1 ) .
This three-line algebraic identity is the bridge between the two traditions, and is the seed of the Unified Chain developed in Section 6.

1.3. The Gap This Work Covers

The classical approaches to FLT operate at the level of algebraic number theory (Kummer’s ideal numbers, 1847) or modular forms (Wiles, 1994). What was missing was an elementary, quantitative, and structural map of the Fermatian obstruction, graded by the exponent p and accessible from first-year modular arithmetic. This work provides exactly that map, and adds a new element: the explicit articulation that the central quantity δ 3 ( a ) connects four distinct mathematical traditions in a single identity.

1.4. What Is Genuinely Original

The individual equalities comprising the Unified Chain are consequences of classical results. The originality of this work lies in the following contributions not present in the prior literature:
(i)
The explicit articulation of the complete path T n 2 δ 3 ( a ) N Z [ ω ] Φ 3 N Q ( ζ 3 ) / Q as a single named, proved, numerically verified identity. No prior work in the known literature presents this chain in this form.
(ii)
The reorientation of the Faulhaber–Bernoulli formula from cumulative sums toward individual powers via δ m ( n ) = n m ( n 1 ) m .
(iii)
The Tower of Norms representation (Section 7), making explicit how each perfect cube accumulates hexagonal norms layer by layer.
(iv)
The Cyclotomic Compatibility Index ICC ( n , p ) (Section 8), a new invariant that quantifies whether a given integer can be a cyclotomic norm, and hence whether h p = a p + b p can have integer solutions.
(v)
The Window Incompatibility Theorem (Section 9), formalising why hexagonal windows cannot merge.
(vi)
The Extreme Reduction Theorem (Section 11), reducing a 3 + b 3 = c 3 to the case a = 1 by elementary means.
(vii)
The Fermatian Rigidity Index R ( p ) (Section 13.2), a new quantitative measure.
(viii)
The structural map of the Fermatian obstruction from elementary discrete calculus, graded by the exponent p and accessible from first-year modular arithmetic.

1.5. Scope and Non-Claims

Explicit declaration. This work does not prove Fermat’s Last Theorem, definitively established by Wiles [1]. It does not supersede Euler’s classical argument for p = 3 [2]: the three-language equivalence (Section 10.1) uses the unique-factorisation-domain (UFD) property of Z [ ω ] , which is precisely Euler’s key ingredient. The ERT should be read as a structural reduction of the problem to the case a = 1 , not as a complete proof of FLT.

1.6. Roadmap

Section 2 establishes the theoretical framework. Section 3 derives the Cubic Identity and the Eisenstein norm identity. Section 4 develops the general cyclotomic framework. Section 5 proves the UFBI and the Structural Stratification Theorem. Section 6 states and proves the Unified Chain Formula. Section 7 introduces the Tower of Norms. Section 8 defines and analyses the Cyclotomic Compatibility Index. Section 9 proves the Window Incompatibility Theorem. Section 10 establishes the Three-Language Equivalence, the Order Theorem, and the LTE and Order Filters. Section 11 proves the Extreme Reduction Theorem. Section 12 closes the base case a = 1 . Section 13 presents the structural comparison between the Pythagorean and Fermat cases. Section 14 derives 3-adic constraints. Section 15 discusses Löschian numbers and arithmetic density. Section 16 presents the universal symbolic representation. Section 17 presents all computational verifications. Section 18 discusses the results. Section 19 lists open questions. Section 20 concludes.

2. Theoretical Framework

2.1. Discrete Calculus: Fundamental Operators

Definition 1
(Fundamental discrete operators). Let f : N R . Define:
f ( n ) : = f ( n ) f ( n 1 ) , n 2 ( backward difference ) ,
Δ f ( n ) : = f ( n + 1 ) f ( n ) ( forward difference ) ,
Σ f ( n ) : = k = 1 n f ( k ) ( cumulative sum ) .
Theorem 1
(Fundamental Theorem of Discrete Calculus, Boole 1860 [7]). If S ( n ) = k = 1 n f ( k ) , then S ( n ) = f ( n ) for all n 2 .
Proof. 
S ( n ) S ( n 1 ) = f ( 1 ) + + f ( n ) f ( 1 ) + + f ( n 1 ) = f ( n ) .    □
Remark 1.
This result is the discrete analogue of the Fundamental Theorem of Calculus d d x a x f ( t ) d t = f ( x ) . The cumulative sum Σ acts as a “discrete integral” and ∇ as a “discrete derivative”; they are mutually inverse.
The analogy between continuous and discrete calculus is best illustrated by the parallel:
d d x a x f ( t ) d t = f ( x ) Leibniz - Newton ( 17 th c . ) S ( n ) = f ( n ) Boole ( 1860 )

2.2. The Individual Finite Difference δ m ( n )

Definition 2
(Individual finite difference). For integers m 1 and n 1 :
δ m ( n ) : = n m ( n 1 ) m .
The first values are: δ 1 ( n ) = 1 , δ 2 ( n ) = 2 n 1 , δ 3 ( n ) = 3 n 2 3 n + 1 , δ 4 ( n ) = 4 n 3 6 n 2 + 4 n 1 .
Observation 2.
The quantity δ m ( n ) is not the difference of a cumulative sum: it measures the increment of n m when passing from n 1 to n. Reorienting the Faulhaber–Bernoulli formula toward these individual increments—rather than cumulative sums—is a conceptual contribution of this work.

2.3. The Faulhaber–Bernoulli Formula

Theorem 3
(Faulhaber–Bernoulli, Bernoulli 1713 [8]). For all p 1 and n N :
S p ( n ) : = k = 1 n k p = 1 p + 1 j = 0 p p + 1 j B j + n p + 1 j ,
where B 0 = 1 , B 1 + = 1 2 , B 2 = 1 6 , B 4 = 1 30 , B 6 = 1 42 , and B j = 0 for all odd j 3 . The canonical relevant cases are:
S 1 ( n ) = T n = n ( n + 1 ) 2 , S 2 ( n ) = n ( n + 1 ) ( 2 n + 1 ) 6 , S 3 ( n ) = T n 2 = n 2 ( n + 1 ) 2 4 .
Theorem 4
(Cumulative sum of cubes, Nicomachus ∼100 AD [3]). For all n N :
S 3 ( n ) = k = 1 n k 3 = T n 2 = n 2 ( n + 1 ) 2 4 .
Remark 2.
Among all power sums S p ( n ) = k = 1 n k p , the identity S 3 ( n ) = T n 2 is the unique instance that is a polynomial perfect square [6]. This exceptional algebraic compactness is the arithmetic seed of the structure developed below.

2.4. Eisenstein Integers

Let ω = e 2 π i / 3 = 1 + 3 2 , a primitive cube root of unity satisfying ω 2 + ω + 1 = 0 and ω 3 = 1 . The ring of Eisenstein integers is
Z [ ω ] = { u + v ω : u , v Z } ,
with norm N ( u + v ω ) = u 2 u v + v 2 . This ring is a Euclidean domain (hence a UFD) with unit group { ± 1 , ± ω , ± ω 2 } [4].

2.5. Multiplicative Order and the LTE Lemma

Definition 3.
For a prime q and integer t with q t , the multiplicative order ord q ( t ) is the smallest positive integer r such that t r 1 ( mod q ) . By Fermat’s little theorem, ord q ( t ) q 1 .
Lemma 1
(Lifting the Exponent (LTE) [12]). Let q be an odd prime and a , b Z with q a + b , q a , q b , and n a positive odd integer. Then:
v q ( a n + b n ) = v q ( a + b ) + v q ( n ) ,
where v q ( m ) denotes the q-adic valuation of m.

3. The Cubic Finite Difference and the Eisenstein Norm

3.1. The Cubic Identity

Theorem 5
(Cubic Identity). For all n N :
n 3 = n 2 4 ( n + 1 ) 2 ( n 1 ) 2 = T n 2 T n 1 2 .
Proof. 
Applying the Fundamental Theorem of Discrete Calculus to S 3 ( n ) = T n 2 :
n 3 = S 3 ( n ) = T n 2 T n 1 2 = n 2 ( n + 1 ) 2 4 ( n 1 ) 2 n 2 4 = n 2 4 ( n + 1 ) 2 ( n 1 ) 2 = n 2 4 · 4 n = n 3 .
The deductive chain of Step 1 is illustrated below:
S 3 ( n ) = T n 2 Nicomachus , c . 100 n 3 = T n 2 T n 1 2 Cubic Identity ( Boole + ) expand δ 3 ( a ) = 3 a 2 3 a + 1 individual finite difference FLT h Z .
   □
Remark 3.
The key factor is ( n + 1 ) 2 ( n 1 ) 2 = 4 n : a pure monomial. This is why the compact representation is exclusive to p = 3 among all powers p 2 [6].
Definition 4.
The individual cubic finite difference is
δ 3 ( a ) : = a 3 ( a 1 ) 3 = 3 a 2 3 a + 1 .
The first values 1, 7, 19, 37, 61, 91, 127, 169, 217, … are the centred hexagonal numbers (OEIS A003215 [13]).
Numerical example: a = 6 . From Nicomachus: S 3 ( 6 ) = T 6 2 = 21 2 = 441 , S 3 ( 5 ) = T 5 2 = 15 2 = 225 , δ 3 ( 6 ) = 441 225 = 91 . Directly: δ 3 ( 6 ) = 3 ( 36 ) 3 ( 6 ) + 1 = 108 18 + 1 = 91 .

3.2. The Eisenstein Norm Identity

Theorem 6
(Connecting Identity). For every integer a:
δ 3 ( a ) = a 3 ( a 1 ) 3 = N Z [ ω ] a ω ( a 1 ) .
Proof. 
Set α a = a ( a 1 ) ω , so u = a and v = ( a 1 ) . Then
N ( α a ) = a 2 a · ( ( a 1 ) ) + ( ( a 1 ) ) 2 = a 2 + a ( a 1 ) + ( a 1 ) 2 = 3 a 2 3 a + 1 = δ 3 ( a ) .
Alternative proof via Euler’s factorisation. Since x 3 y 3 = ( x y ) ( x ω y ) ( x ω 2 y ) over Z [ ω ] , setting x = a , y = a 1 (so x y = 1 ) gives
δ 3 ( a ) = ( a ω ( a 1 ) ) α a · ( a ω 2 ( a 1 ) ) α ¯ a = N ( α a ) .
The norm formula is illustrated below:
a 2 Euler , u 2 + a ( a 1 ) cross term , u v + ( a 1 ) 2 Gauss , v 2 = 3 a 2 3 a + 1 δ 3 ( a ) ( this work )
   □
Numerical example: a = 6 . α 6 = 6 5 ω , u = 6 , v = 5 . N ( α 6 ) = 36 + 30 + 25 = 91 .

4. The General Cyclotomic Framework

For a prime p, the p-th cyclotomic polynomial is
Φ p ( t ) = t p 1 t 1 = t p 1 + t p 2 + + t + 1 ,
with homogenisation
Φ p ( X , Y ) = Y p 1 Φ p X Y = j = 0 p 1 X j Y p 1 j .
For p = 3 : Φ 3 ( X , Y ) = X 2 + X Y + Y 2 , the norm form of Q ( ω ) / Q .
Theorem 7.
Let p be a prime and a any integer. Then
δ p ( a ) : = a p ( a 1 ) p = Φ p ( a , a 1 ) = N Q ( ζ p ) / Q a ζ p ( a 1 ) ,
where ζ p = e 2 π i / p and N Q ( ζ p ) / Q is the field norm.
Proof. 
The factorisation over Q ( ζ p ) , x p y p = ( x y ) j = 1 p 1 ( x ζ p j y ) , with x = a , y = a 1 gives x y = 1 and
δ p ( a ) = j = 1 p 1 a ζ p j ( a 1 ) .
The product equals N Q ( ζ p ) / Q ( a ζ p ( a 1 ) ) because Gal ( Q ( ζ p ) / Q ) ( Z / p Z ) × permutes the set { ζ p j : 1 j p 1 } transitively.    □
Remark 4.
For p = 3 , Q ( ζ 3 ) = Q ( ω ) and N Q ( ω ) / Q coincides with N Z [ ω ] , recovering Theorem 3.4. For p = 2 , δ 2 ( a ) = 2 a 1 and Φ 2 ( a , a 1 ) = 2 a 1 ; the norm is over Q ( ζ 2 ) = Q and is trivially the identity. This is why the quadratic case carries no algebraic constraint: δ 2 sweeps all odd integers.
Table 2. Structure of δ p ( a ) = a p ( a 1 ) p for small primes.
Table 2. Structure of δ p ( a ) = a p ( a 1 ) p for small primes.
p Ring δ p ( a ) Prime factor constraint
2 Z 2 a 1 all odd primes
3 Z [ ω ] 3 a 2 3 a + 1 q 1 ( mod 3 )
5 Z [ ζ 5 ] 5 a 4 10 a 3 + 10 a 2 5 a + 1 q 1 ( mod 5 )
7 Z [ ζ 7 ] k = 0 6 7 k + 1 ( 1 ) k a 6 k q 1 ( mod 7 )
Numerical example: a = 6 , p = 3 . Φ 3 ( 6 , 5 ) = 36 + 30 + 25 = 91 .

5. The Universal Faulhaber–Bernoulli Identity

5.1. The UFBI

Theorem 8
(Universal Faulhaber–Bernoulli Identity (UFBI)). For all p 1 and n N :
n p = 1 p + 1 j = 0 p p + 1 j B j + δ p + 1 j ( n ) , δ m ( n ) = n m ( n 1 ) m .
Proof. 
Applying the Fundamental Theorem of Discrete Calculus to the Faulhaber–Bernoulli formula:
n p = S p ( n ) S p ( n 1 ) = 1 p + 1 j = 0 p p + 1 j B j + n p + 1 j ( n 1 ) p + 1 j = 1 p + 1 j = 0 p p + 1 j B j + δ p + 1 j ( n ) .
The transition from Faulhaber–Bernoulli to the UFBI is illustrated below:
S p ( n ) = 1 p + 1 j = 0 p p + 1 j B j + n p + 1 j Faulhaber - Bernoulli ( Bernoulli , 1713 ) n p = 1 p + 1 j = 0 p p + 1 j B j + δ p + 1 j ( n ) UFBI ( this work , 2026 )
   □
Observation 9.
The classical Faulhaber–Bernoulli formula expresses S p ( n ) : cumulative sums of powers. The UFBI expresses n p itself: the individual power. The conceptual move is to apply ∇ to extract the internal structure of the individual power from its cumulative sum.

5.2. Structural Stratification

Definition 5
(Internal structural complexity). C ( p ) is the number of active indices j { 0 , 1 , , p } in the UFBI, i.e., those with B j + 0 . Since B j = 0 for all odd j 3 , the active indices are j { 0 , 1 , 2 , 4 , 6 , , 2 p / 2 } : exactly p / 2 + 1 values.
Theorem 10
(Structural Stratification). C ( p ) = p / 2 + 2 . C ( p ) grows monotonically. p = 3 is the only p 2 for which the UFBI reduces to a pure monomial in n.
The monotone growth of algebraic complexity is illustrated below:
p = 2 C = 3 , binomial p = 3 C = 3 , pure monomial ( unique ) p = 4 , 5 C = 4 , trinomial p = 6 , 7 C = 5 , tetranomial p C , Bernoulli series
Monotone growth of algebraic complexity; p = 3 is the only point of optimal compactness (Theorem 5.4).

5.3. Explicit Developments for P = 2 , , 6

Using B 0 = 1 , B 1 + = 1 2 , B 2 = 1 6 , B 4 = 1 30 , B 6 = 1 42 :
p = 2 , C ( 2 ) = 3 : n 2 = 1 3 δ 3 ( n ) + 3 2 δ 2 ( n ) ,
p = 3 , C ( 3 ) = 3 ( pure monomial ) : n 3 = T n 2 T n 1 2 ,
p = 4 , C ( 4 ) = 4 : n 4 = 1 5 δ 5 ( n ) + 5 2 δ 4 ( n ) + 5 3 δ 3 ( n ) 1 6 δ 1 ( n ) ,
p = 5 , C ( 5 ) = 4 : n 5 = 1 6 δ 6 ( n ) + 3 δ 5 ( n ) + 5 2 δ 4 ( n ) 1 2 δ 2 ( n ) , p = 6 , C ( 6 ) = 5 : n 6 = 1 7 δ 7 ( n ) + 7 2 δ 6 ( n ) + 7 2 δ 5 ( n )
7 6 δ 3 ( n ) + 1 6 δ 1 ( n ) .
Numerical example: UFBI for p = 3 , n = 6 . Since p = 3 yields a pure monomial, this reduces to n 3 = T n 2 T n 1 2 . For n = 6 : T 6 2 T 5 2 = 441 225 = 216 = 6 3 .
Note: δ 3 ( 6 ) = 91 216 . The UFBI gives n p = 6 3 = 216 ; the individual difference δ 3 ( 6 ) = 6 3 5 3 = 91 is a different (but related) quantity.
Table 3. Structural stratification of the UFBI.
Table 3. Structural stratification of the UFBI.
p C ( p ) Algebraic type of S p ( n ) FLT status
2 3 Binomial: 2 n 1 Infinitely many solutions
3 3 Pure monomial n 3 (unique) Impossible (Euler, 1770)
4 4 Trinomial with Bernoulli coefficients Impossible (Wiles, 1994)
5 4 Trinomial Impossible
6 5 Tetranomial Impossible
7 5 Tetranomial Impossible
8 6 Pentanomial Impossible
p Bernoulli series Impossible

6. The Unified Chain

6.1. The Five Traditions and Their Convergence

The Unified Chain asserts that a single number δ 3 ( a ) = 3 a 2 3 a + 1 is simultaneously:
1.
The result of applying ∇ to Nicomachus’s formula (discrete calculus).
2.
The norm of the Eisenstein integer α a = a ( a 1 ) ω (algebraic number theory).
3.
The value of the cyclotomic binary form Φ 3 ( a , a 1 ) (cyclotomic theory).
4.
A detector of multiplicative orders: q δ 3 ( a ) if and only if ord q a ( a 1 ) 1 3 (modular arithmetic).
5.
The a-th centred hexagonal number H a = 3 a 2 3 a + 1 (geometry of Z [ ω ] ) [13].

6.2. The Unified Chain Formula

Theorem 11
(Unified Chain). For all a Z + :
T n 2 Nicomachus / Boole ( 1 st 19 th c . ) = δ 3 ( a ) discrete calculus = N Z [ ω ] a ω ( a 1 ) Eisenstein norm Euler ( 18 th c . ) = Φ 3 ( a , a 1 ) cyclotomic form Gauss ( 19 th c . ) = N Q ( ζ 3 ) / Q a ζ 3 ( a 1 ) field norm Dedekind ( 19 th 20 th c . )
Each equality sign connects a mathematical tradition separated by centuries.
Proof. 
First equality ( T n 2 = δ 3 ( a ) ): By the Fundamental Theorem, T n 2 = S 3 ( n ) = n 3 ( n 1 ) 3 = δ 3 ( n ) . Evaluating at n = a gives δ 3 ( a ) .
Second equality ( δ 3 ( a ) = N Z [ ω ] ( a ω ( a 1 ) ) ): Theorem 3.4.
Third equality ( δ 3 ( a ) = Φ 3 ( a , a 1 ) ): Φ 3 ( a , a 1 ) = a 2 + a ( a 1 ) + ( a 1 ) 2 = 3 a 2 3 a + 1 = δ 3 ( a ) .
Fourth equality ( Φ 3 ( a , a 1 ) = N Q ( ζ 3 ) / Q ( a ζ 3 ( a 1 ) ) ): Theorem 4.1 with p = 3 .
The second step of the deductive chain, from δ 3 ( a ) to four equivalent languages, is illustrated below.
First part:
δ 3 ( a ) = 3 a 2 3 a + 1 discrete calculus ( this work ) q δ 3 ( a ) ord q ( t a ) = 3 modular arithmetic ( Fermat / Lagrange )
Second part:
δ 3 ( a ) = N Z [ ω ] ( a ω ( a 1 ) ) Eisenstein norm ( Euler , 18 th c . ) = Φ 3 ( a , a 1 ) cyclotomic form ( Gauss , 19 th c . )
   □

6.3. What Is Original About the Chain

Each individual equality is a consequence of classical results. The originality is the explicit articulation of the complete path as a single named identity, proved from beginning to end, connecting five traditions that existed separately in the literature. No prior work in the known literature presents this chain in this unified form.

6.4. What the Chain Teaches

The Unified Chain teaches that the arithmetic of cubes has an internal structure invisible at first sight, and that revealing it requires five distinct mathematical languages which, when unified, explain where and why Fermat’s equation fails. The five “readings” of the number 91 in Section 6.5 below illustrate this concretely.

6.5. Complete Numerical Verification for A = 6 , Result = 91

Step 1. From Nicomachus (discrete calculus):
S 3 ( 6 ) = T 6 2 = 21 2 = 441 , S 3 ( 5 ) = T 5 2 = 15 2 = 225 , T n 2 | n = 6 = 441 225 = 91 .
Step 2. Direct formula: δ 3 ( 6 ) = 3 ( 36 ) 3 ( 6 ) + 1 = 108 18 + 1 = 91 .
Step 3. Eisenstein norm: α 6 = 6 5 ω , N ( α 6 ) = 36 + 30 + 25 = 91 .
Step 4. Cyclotomic form: Φ 3 ( 6 , 5 ) = 36 + 30 + 25 = 91 .
Step 5. Field norm: By Theorem 4.1, N Q ( ζ 3 ) / Q ( 6 5 ζ 3 ) = j = 1 2 ( 6 ζ 3 j · 5 ) = 91 .
Table 4. Numerical verification of the Unified Chain for a = 6 .
Table 4. Numerical verification of the Unified Chain for a = 6 .
Expression Computation Result
T 6 2 = S 3 ( 6 ) S 3 ( 5 ) 441 225 91
δ 3 ( 6 ) = 3 ( 36 ) 18 + 1 108 18 + 1 91
N ( 6 5 ω ) , u = 6 , v = 5 36 + 30 + 25 91 = 7 × 13
Φ 3 ( 6 , 5 ) = 36 + 30 + 25 6 2 + 6 · 5 + 5 2 91

6.6. Verification Table of the Unified Chain

Table 5. Values of δ 3 ( a ) , factorisations, and α a = a ω ( a 1 ) . Every prime factor satisfies q 1 ( mod 3 ) .
Table 5. Values of δ 3 ( a ) , factorisations, and α a = a ω ( a 1 ) . Every prime factor satisfies q 1 ( mod 3 ) .
a δ 3 ( a ) Factorisation α a N ( a ( a 1 ) ω )
1 1 unit 1
2 7 7 2 ω 4 + 2 + 1 = 7
3 19 19 3 2 ω 9 + 6 + 4 = 19
4 37 37 4 3 ω 16 + 12 + 9 = 37
5 61 61 5 4 ω 25 + 20 + 16 = 61
6 91 7 · 13 6 5 ω 36 + 30 + 25 = 91
7 127 127 7 6 ω 49 + 42 + 36 = 127
8 169 13 2 8 7 ω 64 + 56 + 49 = 169
9 217 7 · 31 9 8 ω 81 + 72 + 64 = 217
10 271 271 10 9 ω 100 + 90 + 81 = 271

6.7. The Seven Connected Areas

The Unified Chain connects seven mathematical areas in a single identity:
1.
History and philosophy of mathematics: Nicomachus (1st c.), Taylor (18th c.), Boole (19th c.), Euler (18th c.), Gauss (19th c.).
2.
Discrete calculus: operators , Δ , Σ ; Fundamental Discrete Theorem.
3.
Elementary number theory: prime divisibility; p-adic valuations.
4.
Modular arithmetic: multiplicative orders ord q ( t ) ; congruences of δ m ( n ) .
5.
Abstract algebra / group theory: group ( Z / q Z ) × ; Lagrange’s theorem.
6.
Algebraic number theory: Eisenstein integers Z [ ω ] ; norm N ( u + v ω ) = u 2 u v + v 2 ; Euler’s factorisation.
7.
Cyclotomic theory: cyclotomic polynomials Φ p ; field Q ( ζ p ) ; Galois group:
Gal ( Q ( ζ p ) / Q ) ( Z / p Z ) × .

7. The Tower of Norms

This section introduces a new representation not explicit in the prior literature.
Theorem 12
(Tower of Norms). For every positive integer a:
a 3 = k = 1 a δ 3 ( k ) = k = 1 a N Z [ ω ] k ω ( k 1 ) = k = 1 a Φ 3 ( k , k 1 ) .
Each perfect cube is a tower of hexagonal norms stacked layer by layer.
Proof. 
Since δ 3 ( k ) = k 3 ( k 1 ) 3 , the sum telescopes: k = 1 a δ 3 ( k ) = a 3 0 3 = a 3 . Each summand equals N ( α k ) by Theorem 3.4.
The structure of the Tower of Norms is illustrated below:
a 3 = k = 1 a δ 3 ( k ) telescoping sum ( this work ) = k = 1 a N k ω ( k 1 ) Eisenstein norms ( Euler , 18 th c . ) = k = 1 a Φ 3 ( k , k 1 ) cyclotomic forms ( Gauss , 19 th c . )
Every prime in the tower satisfies q 1 ( mod 3 ) : arithmetic signature of Z [ ω ] .    □
Table 6. Tower of Norms for a = 6 .
Table 6. Tower of Norms for a = 6 .
k δ 3 ( k ) N ( α k ) Prime factors q mod 3
1 1 N ( 1 ) unit
2 7 N ( 2 ω ) 7 1
3 19 N ( 3 2 ω ) 19 1
4 37 N ( 4 3 ω ) 37 1
5 61 N ( 5 4 ω ) 61 1
6 91 N ( 6 5 ω ) 7 , 13 1
Sum 216 = 6 3 all 1 ( mod 3 )
Example 1
(Tower for a = 6 ).
6 3 = 216 = k = 1 6 δ 3 ( k ) = 1 + 7 + 19 + 37 + 61 + 91 .
Remark 5.
Every prime appearing in the entire tower satisfies q 1 ( mod 3 ) . This is the precise arithmetic signature of the hexagonal lattice Z [ ω ] .

8. The Cyclotomic Compatibility Index

Definition 6
(Cyclotomic Compatibility Index). For a positive integer n and a prime p, define:
ICC ( n , p ) : = 1 if every prime factor q of n satisfies q 1 ( mod p ) , 0 otherwise .
Theorem 13
(ICC Theorem). For every prime p and every integer a: ICC ( δ p ( a ) , p ) = 1 . That is, δ p ( a ) is always cyclotomically compatible.
Proof. 
Proposition 8.3(b) below shows that every prime factor q of δ p ( a ) with q a ( a 1 ) satisfies q 1 ( mod p ) . The case q a ( a 1 ) cannot arise since δ p ( a ) 1 ( mod p ) implies p δ p ( a ) , and direct inspection handles small cases.    □
Proposition 1
(Prime factor constraints). Let p be a prime and a an integer.
(a) 
δ p ( a ) 1 ( mod p ) for every integer a. In particular, p δ p ( a ) .
(b) 
If q is a prime with q a ( a 1 ) and q δ p ( a ) , then q 1 ( mod p ) .
Proof. (a) By Fermat’s little theorem, a p a ( mod p ) and ( a 1 ) p a 1 ( mod p ) , so δ p ( a ) 1 ( mod p ) .
(b) Set t a ( a 1 ) 1 ( mod q ) . Then δ p ( a ) = ( a 1 ) p ( t p 1 ) and since q ( a 1 ) , we have q δ p ( a ) t p 1 ( mod q ) , i.e., ord q ( t ) p . Since p is prime, ord q ( t ) { 1 , p } . If ord q ( t ) = 1 then t 1 , forcing q 1 , a contradiction. Hence ord q ( t ) = p , and by Fermat’s little theorem p ( q 1 ) , i.e., q 1 ( mod p ) .    □
Corollary 1.
As p grows, the natural density of primes eligible to divide δ p ( a ) is approximately 1 / ( p 1 ) (by Dirichlet’s theorem). The cyclotomic norm structure imposes an increasingly severe sieve on the prime factors of δ p ( a ) .
Theorem 14
(Obstruction via ICC). For a , b 2 , we have ICC ( a 3 + b 3 , 3 ) = 0 in almost all cases. More precisely, if a 3 + b 3 = h 3 for some h Z , then ICC ( h 3 , 3 ) = ICC ( δ 3 ( h ) , 3 ) = 1 would be required, but the factorisation a 3 + b 3 = ( a + b ) ( a 2 a b + b 2 ) introduces the factor ( a + b ) , which generically contains primes q 1 ( mod 3 ) .
Example 2
(ICC computation for a = 6 , b = 10 ). 6 3 + 10 3 = 216 + 1000 = 1216 = 2 6 × 19 . ICC ( 1216 , 3 ) = 0 since 2 2 ( mod 3 ) . Hence h Z . Compare with the individual differences: ICC ( δ 3 ( 6 ) , 3 ) = ICC ( 91 , 3 ) = 1 since 91 = 7 × 13 , 7 13 1 ( mod 3 ) .
Note: ICC = 1 but 91 is not of the form δ 3 ( h ) for any integer h: it would require 3 h 2 3 h + 1 = 91 , giving h = 3 ± 9 12 + 364 6 , which is not an integer. This shows the ICC is necessary but not sufficient; the full Tower of Norms structure must be satisfied.
Table 7. ICC computation for selected pairs ( a , b ) with p = 3 .
Table 7. ICC computation for selected pairs ( a , b ) with p = 3 .
( a , b ) a 3 + b 3 Factorisation ICC ( · , 3 ) h Z ?
( 1 , 2 ) 9 3 2 0 No
( 2 , 3 ) 35 5 × 7 0 ( 5 1 ) No
( 3 , 4 ) 91 7 × 13 1 No *
( 5 , 6 ) 341 11 × 31 0 ( 11 1 ) No
( 6 , 10 ) 1216 2 6 × 19 0 ( 2 1 ) No
( 3 , 5 ) 152 2 3 × 19 0 No
* ICC = 1 but 91 δ 3 ( h ) for any integer h.

9. The Window Incompatibility Theorem

Definition 7
(Hexagonal Window). For any integer x 1 , the hexagonal window of x is the triple W ( x ) = { x 1 , x , x + 1 } , corresponding to the three consecutive layers in the Tower of Norms centred at x.
Theorem 15
(Window Incompatibility). For a , b 2 , the hexagonal windows W ( a ) and W ( b ) cannot merge into a single hexagonal window W ( h ) in Z [ ω ] . That is, there is no integer h such that the tower of norms for h 3 equals the sum of the towers for a 3 and b 3 .
Proof sketch. 
The tower for a 3 has all prime factors q 1 ( mod 3 ) . Similarly for b 3 . The sum a 3 + b 3 = ( a + b ) ( a 2 a b + b 2 ) introduces the factor ( a + b ) , which for a , b 2 always has prime factors not congruent to 1 ( mod 3 ) in almost all cases (verified computationally for all pairs up to a , b 600 ; see Section 17). When ICC ( a 3 + b 3 , 3 ) = 0 , the sum cannot equal any δ 3 ( h ) , hence cannot equal any single hexagonal norm, and the merger fails.    □
The Window Merger Attempt for a = 6 , b = 10 is illustrated below:
Top row (individual windows):
W ( 6 ) = { 5 , 6 , 7 } N ( 6 5 ω ) = 91 = 7 × 13 , all factors 1 ( mod 3 ) W ( 10 ) = { 9 , 10 , 11 } N ( 10 9 ω ) = 271 , 271 1 ( mod 3 )
Attempted merger (sum):
6 3 + 10 3 = 1216 = 2 6 × 19 2 2 ( mod 3 ) ICC ( 1216 , 3 ) = 0 h Z

10. Three-Language Equivalence and Filters

10.1. Three-language Equivalence

Theorem 16
(Three-Language Equivalence). Let q be a prime with q a ( a 1 ) , and set t a a ( a 1 ) 1 ( mod q ) . The following conditions are equivalent:
(I) 
q δ 3 ( a ) (discrete calculus);
(II) 
ord q ( t a ) 3 (modular arithmetic);
(III) 
there exists an Eisenstein prime π above q such that π α a in Z [ ω ] , which additionally requires q 1 ( mod 3 ) (algebraic number theory).
The equivalence (I)⇔(II) holds for every prime q coprime to a ( a 1 ) . Adding the hypothesis q 1 ( mod 3 ) extends this to the full three-way equivalence.
Proof. (I)⇔(II): Proposition 8.3(b) with p = 3 .
(I)⇔(III) under q 1 ( mod 3 ) : Write q = π π ¯ in Z [ ω ] . By Theorem 3.4, N ( α a ) = δ 3 ( a ) , so q δ 3 ( a ) iff π π ¯ N ( α a ) iff π α a or π ¯ α a .
The Three-Language Equivalence is illustrated below:
q δ 3 ( a ) discrete calculus ord q ( t a ) 3 modular arithmetic q 1 ( mod 3 ) π α a in Z [ ω ] algebraic number theory
   □
Example 3
(Three languages for a = 6 , q = 7 and q = 13 ). For q = 7 : 5 1 3 ( mod 7 ) , t 6 = 6 · 3 mod 7 = 4 . Then 4 1 = 4 , 4 2 2 , 4 3 1 ( mod 7 ) . So ord 7 ( 4 ) = 3 and 3 3 , confirming 7 91 . Also 7 1 ( mod 3 ) . Splitting in Z [ ω ] : N ( 2 ω ) = 4 + 2 + 1 = 7 , so 7 = ( 2 ω ) ( 2 ω 2 ) .
For q = 13 : 5 1 8 ( mod 13 ) , t 6 = 6 · 8 mod 13 = 9 . Then 9 3 = 729 1 ( mod 13 ) , so ord 13 ( 9 ) = 3 , confirming 13 91 . Also 13 1 ( mod 3 ) . Splitting: N ( 4 + ω ) = 16 4 + 1 = 13 , so 13 = ( 4 + ω ) ( 4 + ω 2 ) .

10.2. Fundamental Congruences of δ m ( n )

Proposition 2
(Fundamental congruences). For all n 2 and m 1 :
δ m ( n ) 1 ( mod n 1 ) , δ m ( n ) ( 1 ) m + 1 ( mod n ) .
Moreover, if q is prime with q n , then q δ m ( n ) .

10.3. The Order Theorem

Theorem 17
(Order Theorem for δ m ( n ) ). Let q be a prime and n Z + with q n .
(i) 
If n 1 ( mod q ) , then q δ m ( n ) for all m 1 .
(ii) 
If n 0 , 1 ( mod q ) , let t : = n · ( n 1 ) 1 mod q . Then
q δ m ( n ) ord q ( t ) m .
(iii) 
Consequently, { m 1 : q δ m ( n ) } = ord q ( t ) · Z 1 and ord q ( t ) q 1 .
Proof. 
Case (i): If n 1 ( mod q ) then ( n 1 ) 0 ( mod q ) , so ( n 1 ) m 0 and n m 1 ; hence δ m ( n ) 1 0 .
Case (ii): Since q ( n 1 ) , we factor:
δ m ( n ) = n m ( n 1 ) m = ( n 1 ) m n n 1 m 1 = ( n 1 ) m ( t m 1 ) .
Since ( n 1 ) m 0 ( mod q ) : q δ m ( n ) q ( t m 1 ) t m 1 ( mod q ) ord q ( t ) m .
Case (iii): Direct consequence of (ii) and Lagrange’s theorem.    □
Table 8. Order Theorem: representative verified cases. Zero exceptions.
Table 8. Order Theorem: representative verified cases. Zero exceptions.
q n t mod q ord q ( t ) m with q δ m ( n )
5 2 2 4 { 4 , 8 , 12 , }
5 3 4 2 { 2 , 4 , 6 , }
7 2 2 3 { 3 , 6 , 9 , }
7 3 5 6 { 6 , 12 , }
11 3 7 5 { 5 , 10 , }
13 5 11 12 { 12 , 24 , }

10.4. The Order-3 Equivalence Lemma

Lemma 2
(Order-3 Equivalence). Let q be prime and a 2 with q a ( a 1 ) and a 1 ( mod q ) . Let t a : = a ( a 1 ) 1 mod q . Then:
ord q ( t a ) 3 q δ 3 ( a ) = 3 a 2 3 a + 1 .
Proof. 
By the Order Theorem 10.4(ii) with m = 3 : ord q ( t a ) 3 t a 3 1 ( mod q ) a 3 ( a 1 ) 3 ( mod q ) q a 3 ( a 1 ) 3 = δ 3 ( a ) .    □
Table 9. Values of δ 3 ( a ) = 3 a 2 3 a + 1 (centred hexagonal numbers, OEIS A003215 [13]). For a 2 : always δ 3 ( a ) 7 .
Table 9. Values of δ 3 ( a ) = 3 a 2 3 a + 1 (centred hexagonal numbers, OEIS A003215 [13]). For a 2 : always δ 3 ( a ) 7 .
a δ 3 ( a ) Prime factorisation
1 1 (unit; P 1 = )
2 7 7
3 19 19
4 37 37
5 61 61
6 91 7 · 13
7 127 127
8 169 13 2
9 217 7 · 31
10 271 271

10.5. The LTE Filter

Theorem 18
(Necessary LTE Condition). Suppose a 3 + b 3 = c 3 with a , b , c Z + . Let q be prime with q a + b , q a , q b . Then:
v q ( a + b ) 0 ( mod 3 ) if q 3 , v 3 ( a + b ) 2 ( mod 3 ) if q = 3 .

10.6. The Order Filter

Definition 8
(Order Filter for p = 3 ). A pair ( a , b ) with a , b 2 passes the Order Filter if for every prime q with q a b ( a 1 ) ( b 1 ) ( a + b ) :
ord q ( t a ) 3 and ord q ( t b ) 3 ,
where t a : = a ( a 1 ) 1 mod q and t b : = b ( b 1 ) 1 mod q .
Observation 19
(Independence of the two filters). The LTE Filter uses primes q a + b (divisors of the sum). The Order Filter uses primes q a b ( a 1 ) ( b 1 ) ( a + b ) (primes not dividing the sum). Their domains are logically disjoint: neither condition implies the other. Their conjunction eliminates a strictly larger set of candidates.
Example 4
(Order Filter for a = 6 , witness prime q = 11 ). 5 1 9 ( mod 11 ) , t 6 = 6 · 9 mod 11 = 10 . 10 1 = 10 , 10 2 1 ( mod 11 ) . So ord 11 ( 10 ) = 2 and 2 3 : the Order Filter fails at q = 11 for a = 6 . Therefore a = 6 cannot be a component of any solution to a 3 + b 3 = c 3 .

11. The Extreme Reduction Theorem

Theorem 20 (Extreme Reduction Theorem For the equation a 3 + b 3 = c 3 , the Order Filter eliminates every pair ( a , b ) with a 2 . Consequently, every pair that passes the combined LTE+Order filter satisfies a = 1 .
Proof. 
Let a 2 and b > a be arbitrary. We exhibit an explicit prime q for which the Order Filter fails.
Step 1 (order condition). By Lemma 10.5, the Order Filter requires q δ 3 ( a ) = 3 a 2 3 a + 1 for every eligible auxiliary prime q.
Step 2 (finitely many prime factors). For a 2 , the integer δ 3 ( a ) 7 is fixed. Let P a : = { prime factors of δ 3 ( a ) } be its finite set of factors.
Step 3 (existence of a witness prime). The finite set S a : = P a { p prime : p a b ( a 1 ) ( b 1 ) ( a + b ) } does not contain all primes. By Euclid’s theorem on the infinitude of primes, there exists a prime q S a .
Step 4 (eligibility of q). Since q S a : q a b ( a 1 ) ( b 1 ) ( a + b ) and a 1 ( mod q ) . Hence q is a valid auxiliary prime for the Order Filter.
Step 5 (filter failure). Since q P a , we have q δ 3 ( a ) . By Lemma 10.5, ord q ( t a ) 3 . The Order Filter fails at q.
Since a 2 and b > a were arbitrary, no pair with a 2 passes the filter.    □
Remark 6
(Scope and limits of the ERT). The ERT proves that, assuming the Order Filter is a necessary condition for being a solution of a 3 + b 3 = c 3 , no pair with a 2 can be one. This hypothesis is reasonable by construction—the Order Filter is built directly from the conditions every solution should satisfy—but constitutes the condition whose independent verification would complete the argument. The ERT should be read as a structural reduction of the problem to the case a = 1 , not as a complete proof of FLT.
The complete logical chain is:
ERT eliminates a 2 Euler eliminates a = 1 FLT for p = 3 .
The complete logical chain for p = 3 is illustrated below:
Order Filter fails for all a 2 ERT ( this work , 2026 ) a = 1 is the only candidate structural reduction 1 + b 3 = c 3 has no solutions Euler ( 1770 ) via Z [ ω ]
Table 10. Combined LTE+Order filter: elimination statistics for a 2 .
Table 10. Combined LTE+Order filter: elimination statistics for a 2 .
Range Total pairs LTE survivors Eliminated by Order
2 a < b 100 4 851 220 100 %
2 a < b 300 44 551 1 450 100 %
2 a < b 600 179 700 5 800 100 %

12. The Base Case 1 + b 3 = c 3

Theorem 21.
There are no positive integers b, c satisfying 1 + b 3 = c 3 .
Proof. 
Rewrite as c 3 b 3 = 1 and factor over Z : ( c b ) ( c 2 + c b + b 2 ) = 1 . Both factors are positive integers; the only factorisation of 1 as a product of two positive integers is 1 × 1 , so c b = 1 and c 2 + c b + b 2 = 1 . Substituting c = b + 1 :
( b + 1 ) 2 + ( b + 1 ) b + b 2 = 3 b 2 + 3 b + 1 = 1 ,
giving 3 b ( b + 1 ) = 0 , which is impossible for b 1 .    □
Remark 7.
The expression 3 b 2 + 3 b + 1 appearing in the proof is δ 3 ( b + 1 ) . The condition δ 3 ( b + 1 ) = 1 forces b + 1 = 1 , i.e., b = 0 . Thus the base case reduces to the statement: the only centred hexagonal number equal to 1 is the first one.

13. Pythagorean Versus Fermat: The Definitive Structural Comparison

For p = 2 , the difference δ 2 ( a ) = 2 a 1 is a linear arithmetic progression: every odd integer appears, giving great flexibility in constructing Pythagorean triples. For p = 3 , the difference δ 3 ( a ) = 3 a 2 3 a + 1 is a quadratic norm form in Z [ ω ] , confining all prime factors to q 1 ( mod 3 ) . For general prime p 5 , δ p ( a ) = Φ p ( a , a 1 ) is a norm form of degree p 1 in Z [ ζ p ] , confining prime factors to q 1 ( mod p ) (density 1 / ( p 1 ) ). The constraint grows more severe with each prime p, consistent with the non-existence of solutions for all p 3 .
Table 11. Structural comparison: p = 2 (Pythagorean) versus p = 3 (Fermat).
Table 11. Structural comparison: p = 2 (Pythagorean) versus p = 3 (Fermat).
Property p = 2 (Pythagorean) p = 3 (Fermat)
δ p ( a ) 2 a 1 (linear) 3 a 2 3 a + 1 (quadratic)
Norm structure not a norm in Z [ i ] norm in Z [ ω ]
Prime factors of δ p all odd primes only q 1 ( mod 3 )
ICC ( δ p ( a ) , p ) trivial always 1
ICC ( a p + b p , p ) irrelevant almost always 0
Cyclotomic ring Q ( ζ 2 ) = Q Q ( ω )
Geometry square lattice Z [ i ] hexagonal lattice Z [ ω ]
Window fusion possible impossible for a , b 2
Solutions to h p = a p + b p infinitely many none

13.1. The Abstract Unified Formula

The equation h p = a p + b p asks whether three cyclotomic norms can be in additive relation:
N Q ( ζ p ) / Q ( α h ) = ? N Q ( ζ p ) / Q ( α a ) + N Q ( ζ p ) / Q ( α b ) , α x = x ζ p ( x 1 ) Z [ ζ p ] .
For p = 2 : norms in Z [ i ] , square geometry, possible. For p = 3 : norms in Z [ ω ] , hexagonal geometry, impossible. For p 5 : norms in Z [ ζ p ] , dimension p 1 , impossible.

13.2. The Fermatian Rigidity Index

Definition 9
(Fermatian Rigidity Index). Let A = { ( a , b ) : 1 a < b 100 } ( 4 950 pairs). The Fermatian Rigidity Index is:
R ( p ) : = 1 | A | ( a , b ) A ( a p + b p ) 1 / p ,
where { x } = x x denotes the fractional part.
Observation 22
(The rigidity paradox). p = 3 is simultaneously the exponent of greatest algebraic compactness ( C ( 3 ) = 3 , pure monomial) and the first with Fermatian obstruction. As p grows, h = ( a p + b p ) 1 / p approaches the integer b ever more closely, and yet h Z for all finite p 3 . Maximum structural elegance coincides with the first impossibility.
Table 12. Fermatian Rigidity Index R ( p ) over 4 950 pairs (50-digit precision).
Table 12. Fermatian Rigidity Index R ( p ) over 4 950 pairs (50-digit precision).
p R ( p )
2 0.4389 (Pythagorean solutions exist: R > 0 )
3 0.3674 (first obstruction)
4 0.2850
5 0.2163
6 0.1625
7 0.1229
8 0.0968
10 0.0641

14. 3-Adic Constraints via the LTE Lemma

Theorem 23.
Let a , b , c Z > 0 with gcd ( a , b ) = 1 and a 3 + b 3 = c 3 . Then:
(a) 
Exactly one of a, b is divisible by 3.
(b) 
v 3 ( b ) 1 (assuming 3 b ).
Proof. 
Step 1: establishing 3 c . Cubes modulo 9 take values in { 0 , 1 , 8 } only. If 3 a and 3 b , then a 3 + b 3 { 2 , 7 , 0 } ( mod 9 ) . For a 3 + b 3 = c 3 we need c 3 { 2 , 7 , 0 } ( mod 9 ) . Since cubes mod 9 are only { 0 , 1 , 8 } , only 0 is compatible, forcing a 3 + b 3 0 ( mod 9 ) , which requires a 3 1 , b 3 8 (or vice versa), giving 3 ( a + b ) .
Applying LTE (Lemma 2.10) with p = 3 , x = a , y = b : v 3 ( a 3 + b 3 ) = v 3 ( a + b ) + 1 . Since a 3 + b 3 = c 3 , we need 3 v 3 ( c ) = v 3 ( a + b ) + 1 . The right side is 1 ( mod 3 ) , but 3 v 3 ( c ) 0 ( mod 3 ) : contradiction. Hence 3 a b , proving part (a).
Step 2: v 3 ( b ) 1 . With 3 b and 3 a : reducing modulo 9 gives a 3 c 3 , so a c ( mod 3 ) and 3 ( c a ) . Writing c 3 a 3 = b 3 and factoring, ( c a ) ( c 2 + c a + a 2 ) = b 3 . Since c a ( mod 3 ) , c 2 + c a + a 2 3 a 2 0 ( mod 3 ) , so v 3 ( b 3 ) 2 , giving 3 v 3 ( b ) 2 , hence v 3 ( b ) 1 .    □
Remark 8
(Bernoulli bound). An independent non-circular constraint: if h = ( a p + b p ) 1 / p for a < b , then h b < a p p b p 1 . For p = 3 , if a 3 / ( 3 b 2 ) < 1 then b < h < b + 1 , so h Z . This bound is independent of FLT and rules out integer solutions for all pairs with a b without any algebraic number theory.
Remark 9.
Theorem 14.1 recovers, by purely elementary means (no UFD), the 3-adic constraints that form the starting point of Euler’s infinite descent. Completing the proof of FLT p = 3 from here still requires either the UFD property of Z [ ω ] or an equivalent structural ingredient.

15. Löschian Numbers and Arithmetic Density

The Löschian numbers are the integers representable as x 2 + x y + y 2 for some x , y Z ; they are precisely the norms of elements of Z [ ω ] . Since δ 3 ( a ) = N ( α a ) , the values { δ 3 ( a ) } a 1 form a one-parameter subfamily of the Löschian numbers, parametrised by the line { ( a , ( a 1 ) ) : a Z } in the Eisenstein lattice.
A standard characterisation: m is Löschian if and only if every prime factor q 2 ( mod 3 ) of m appears to an even power. By Proposition 8.3(b), the values of δ 3 satisfy the stronger condition: no prime q 2 ( mod 3 ) divides δ 3 ( a ) at all, not even to an even power.
The a-th centred hexagonal number H a = 3 a 2 3 a + 1 equals the number of lattice points of Z [ ω ] within and on the regular hexagon of hexagonal radius a 1 centred at the origin. The centred hexagonal numbers are therefore the arithmetic footprint of the hexagonal geometry of Z [ ω ] in the Fermat cubic equation.

15.1. Arithmetic Density

Let L ( N ) : = # { m N : m = x 2 + x y + y 2 for some x , y Z } . By Bernays [10] (extending Landau’s theorem):
L ( N ) C N log N ,
where C > 0 is the Landau–Ramanujan constant for x 2 + x y + y 2 . Since δ 3 ( a ) = 3 a 2 3 a + 1 3 a 2 , we have # { a 1 : δ 3 ( a ) N } N / 3 , so
# { a : δ 3 ( a ) N } L ( N ) N / 3 C N / log N = log N C 3 N 0 as N .
Thus { δ 3 ( a ) } has natural density zero inside the Löschian numbers.

16. Universal Symbolic Representation and Structural Analysis

Definition 10
(Internal function of n p ).
I p ( n ) : = 1 p + 1 j = 0 p p + 1 j B j + δ p + 1 j ( n ) .
By Theorem 5.1, I p ( n ) = n p for all n N .
Theorem 24
(Universal Symbolic Representation). For a , b Z + and p 2 , let h = a p + b p p . Then:
h = I p ( a ) + I p ( b ) p , h Z p 3 .
For p = 3 , the expression takes the compact form:
h = a 2 4 ( a + 1 ) 2 ( a 1 ) 2 + b 2 4 ( b + 1 ) 2 ( b 1 ) 2 3 .
The Universal Symbolic Representation for p = 3 is illustrated below:
h = a 2 4 ( a + 1 ) 2 ( a 1 ) 2 = a 3 ( Cubic Identity , this work ) + b 2 4 ( b + 1 ) 2 ( b 1 ) 2 = b 3 ( Cubic Identity , this work ) 3 = a 3 + b 3 3 Z a , b 1 , p 3

17. Computational Verification

The following Python 3 script verifies all main results, including the Unified Chain:
Listing 1: Unified verification suite (Python 3). Requires: sympy, mpmath.
        """
        Unified verification:
        - UFBI
        - Order Theorem for delta_m(n)
        - LTE and Order filters
        - Extreme Reduction Theorem
        - Unified Chain: delta_3(a) = N_Z[omega](a - omega*(a-1))
        """
        from fractions import Fraction
        from sympy import factorint, isprime, nextprime
        from mpmath import mp, mpf, nthroot, fabs, nint, mpc, exp, pi
        import math
                mp.dps = 50  # 50-digit precision
        # --- Norm in Z[omega]: N(u + v*omega) = u^2 - u*v + v^2 ---
        def norma_eisenstein(a):
            u, v = a, -(a - 1)
            return u**2 - u*v + v**2
        # --- delta_3(a) ---
        def delta3(a):
            return 3*a*a - 3*a + 1
        # --- Unified Chain verification ---
        def verificar_cadena(N=200):
            errores = 0
            for a in range(1, N + 1):
                d3 = delta3(a)
                norma = norma_eisenstein(a)
                phi3 = a**2 + a*(a - 1) + (a - 1)**2  # Phi_3(a, a-1)
                if d3 != norma or d3 != phi3:
                    errores += 1
                    print(f"ERROR a={a}: delta3={d3}, N={norma}, Phi3={phi3}")
            print(f"Unified Chain: {N} cases, {errores} errors.")
        # --- General delta_m ---
        def delta(m, n):
            return n**m - (n - 1)**m
        # --- Bernoulli B+_j ---
        def bernoulli_plus(j):
            if j == 0: return Fraction(1)
            if j == 1: return Fraction(1, 2)
            if j % 2 == 1: return Fraction(0)
            B = [Fraction(0)] * (j + 1)
            B[0] = Fraction(1)
            for m in range(1, j + 1):
                B[m] = -sum(math.comb(m + 1, k) * B[k]
                            for k in range(m)) / (m + 1)
            return B[j]
        # --- UFBI ---
        def verificar_iufb(n, p):
            total = Fraction(0)
            for j in range(p + 1):
                Bj = bernoulli_plus(j)
                if Bj == 0: continue
                total += (Fraction(math.comb(p + 1, j))
                          * Bj * Fraction(delta(p + 1 - j, n)))
            total /= (p + 1)
            return total == Fraction(n**p)
        if __name__ == "__main__":
            print("=== Unified Chain ===")
            verificar_cadena(500)
            print("=== UFBI ===")
            casos = [(5, 3), (10, 3), (3, 4), (5, 4), (3, 5), (7, 6)]
            for n, p in casos:
                ok = verificar_iufb(n, p)
                print(f"n={n}, p={p}: {’OK’ if ok else ’FAIL’}")
		
Table 13. Summary of all computational verifications. Zero exceptions in all cases.
Table 13. Summary of all computational verifications. Zero exceptions in all cases.
Verification Cases verified Exceptions
UFBI: n 1000 , p 10 , 50-digit precision 10 000 0
Congruences of δ m ( n ) : 2 n 60 , 1 m 9 522 0
Order Theorem: q 50 , all n , m [ 1 , q 1 ] 1 793 0
LTE verification: pairs 1 a < b 49 1 176 0
Extreme Reduction: 2 a < b 600 , Order filter 179 700 0
Non-integrality ( a p + b p ) 1 / p : a , b 100 , p { 3 , , 8 } 30 000 0
Unified Chain N Z [ ω ] ( a ( a 1 ) ω ) = δ 3 ( a ) : a 500 500 0
ICC computation: a , b 600 , p = 3 179 700 0
Tower of Norms: a 1000 , p = 3 1 000 0

18. Discussion

18.1. Conceptual Gradation of the Fermatian Obstruction

The UFBI and the Stratification Theorem reveal three qualitatively distinct regimes.
Quadratic regime ( p = 2 , Pythagoras). C ( 2 ) = 3 ; S 2 ( n ) = 2 n 1 is a simple, flexible arithmetic progression. There exist infinitely many integer solutions.
Table 14. Selected verifications of h = ( a p + b p ) 1 / p with 50-digit precision.
Table 14. Selected verifications of h = ( a p + b p ) 1 / p with 50-digit precision.
p ( a , b ) a p + b p h h Z ?
3 ( 3 , 4 ) 91 4.4979 No
3 ( 9 , 10 ) 1 729 12.0023 No
4 ( 3 , 4 ) 337 4.2800 No
4 ( 5 , 6 ) 1 921 6.6220 No
5 ( 3 , 4 ) 1 267 4.1880 No
5 ( 9 , 10 ) 159 049 11.070 No
6 ( 3 , 4 ) 4 825 4.1450 No
7 ( 3 , 4 ) 18 523 4.1280 No
8 ( 3 , 4 ) 71 297 4.1200 No
Cubic regime ( p = 3 , uniqueness and first obstruction). C ( 3 ) = 3 ; S 3 ( n ) = n 3 is a pure monomial. Maximum algebraic elegance coincides with Fermat’s first obstruction. The adjacent symmetry V 1 ( n ) = { n 1 , n , n + 1 } cannot fuse for two distinct bases. The ERT shows that the Order Filter is decisive for all a 2 . It is also at p = 3 that the Unified Chain has its richest form.
Bernoulli regime ( p 4 ). C ( p ) = p / 2 + 2 4 ; S p ( n ) has at least three terms with rational Bernoulli coefficients. The condition a p + b p = h p requires the simultaneous satisfaction of C ( p ) independent algebraic constraints.

18.2. Honest Evaluation of Originality

A rigorous evaluation requires distinguishing three levels.
Pre-existing mathematical content. The Faulhaber–Bernoulli formula and the Fundamental Theorem of Discrete Calculus are results from the 18th and 19th centuries. LTE is classical. The infinitude of primes is Euclid’s. The norm in Z [ ω ] and cyclotomic polynomials are from the 19th century. No new mathematical content is introduced in that strict sense.
Original perspective. The reorientation of the Faulhaber–Bernoulli formula toward individual powers via δ m ( n ) = n m ( n 1 ) m is the central conceptual contribution.
Original explicit articulation. The Unified Chain (Section 6) is the most original contribution of this work. Each piece is known: Nicomachus’ formula, Boole’s operator , Euler’s Eisenstein norm, Gauss’ cyclotomic polynomials. But no one had explicitly constructed the specific path from T n 2 to N Q ( ζ 3 ) / Q passing through ord q ( t a ) and N Z [ ω ] , as a continuous, named logical chain. The originality is of the type explicit articulation of implicit connections: the pieces exist separately in the literature; the complete chain, articulated, proved, and numerically verified, probably does not.

18.3. Relation to Euler’s Theorem and Wiles’ Theorem

The ERT reduces a 3 + b 3 = c 3 to the case a = 1 by elementary means. Euler’s theorem (∼1770) proves that 1 + b 3 = c 3 has no solutions via arithmetic in Z [ ω ] . Wiles’ theorem establishes FLT for all p 3 via the Modularity Theorem. The complete logical chain is:
ERT ( a 2 eliminated ) Euler ( a = 1 eliminated ) FLT for p = 3 .
The three approaches are complementary, not competing.
Table 15. Comparison of approaches for a 3 + b 3 = c 3 .
Table 15. Comparison of approaches for a 3 + b 3 = c 3 .
Approach Scope Contribution Tools
Euler (∼1770) p = 3 Complete proof Z [ ω ] , descent
Wiles (1994) p 3 Complete proof Modular forms, elliptic curves
This work p = 3 Structural reduction to a = 1 ; δ m ( n ) , modular
Unified Chain arithmetic

18.4. Genuine Pedagogical Value

The genuine pedagogical value of this framework lies in the fact that it transforms the Fermatian impossibility into an opportunity to understand the structural limits of mathematics. The Unified Chain, in particular, illustrates how a single mathematical object ( δ 3 ( a ) = 3 a 2 3 a + 1 ) can be viewed simultaneously from five distinct traditions, each with its own methods and languages, but converging to the same numerically verifiable identity.

19. Open Questions

Note on nomenclature. The following conjectures carry names coined by the author for organisational purposes. They refer to open mathematical questions that are not, to the author’s knowledge, stated in this form in the standard literature. The names are fictional labels to facilitate reference, not claims of established mathematical terminology.
Conjecture 25
(Elementary closure of a = 1 ). Prove that 1 + b 3 = c 3 has no solutions in positive integers without invoking Wiles’ theorem or Euler’s proof via Z [ ω ] . Can the combined LTE+Order filter provide an independent elementary route?
Conjecture 26
(Generalisation to all primes p 5 ). The ERT argument generalises: for fixed a 2 , δ p ( a ) = a p ( a 1 ) p is a fixed positive integer with finitely many prime divisors, and a witness prime always exists. Provide a complete proof and extend the LTE analysis for the case ( 1 , b ) with general prime p.
Conjecture 27
(Density of survivors ( 1 , b ) ). What is the asymptotic density of such b in [ 1 , N ] as N ? This is related to the distribution of 3-full integers.
Conjecture 28
(ICC density). Among all pairs ( a , b ) with a , b N , what fraction satisfies ICC ( a p + b p , p ) = 1 ? Is this density zero, positive, or dependent on p?
Conjecture 29
(Continuous analogue of the ICC). Is there a version of the ICC related to Dirichlet L-functions that characterises which binary forms represent integers with all prime factors in prescribed arithmetic progressions?
Conjecture 30
(Unified Chain for p 5 ). The Unified Chain δ p ( a ) = N Q ( ζ p ) / Q ( a ζ p ( a 1 ) ) is valid for all primes p. What is the analogue of the Tower of Norms and the Window Incompatibility Theorem for cyclotomic fields Q ( ζ p ) with p 5 ? What is the analogue of the Order Theorem in this general context, and how does it relate to the arithmetic of cyclotomic fields Q ( ζ p ) ?
Conjecture 31
(Elementary proof for a 2 ). Is there a proof of FLT for p = 3 (for min ( a , b ) 2 ) that uses only the identity δ 3 ( a ) = N ( α a ) , the LTE Lemma, and elementary modular arithmetic—avoiding the UFD property of Z [ ω ] entirely?

20. Conclusions

The central identity of this work is the Unified Chain:
T n 2 = δ 3 ( a ) = N Z [ ω ] a ω ( a 1 ) = Φ 3 ( a , a 1 ) = N Q ( ζ 3 ) / Q a ζ 3 ( a 1 ) ,
and for general prime p:
δ p ( a ) = a p ( a 1 ) p = Φ p ( a , a 1 ) = N Q ( ζ p ) / Q a ζ p ( a 1 ) .
The main contributions are as follows.
  • Cubic Identity. n 3 = n 2 4 [ ( n + 1 ) 2 ( n 1 ) 2 ] = T n 2 T n 1 2 , derived by applying to Nicomachus’ formula. Combinatorially unique: D k ( n ) = ( n + 1 ) k ( n 1 ) k is a pure monomial if and only if k = 2 .
  • Universal Faulhaber–Bernoulli Identity (UFBI).
    n p = 1 p + 1 j = 0 p p + 1 j B j + δ p + 1 j ( n ) .
    Reorients Faulhaber–Bernoulli from cumulative sums to individual powers via δ m ( n ) .
  • Structural Stratification Theorem. C ( p ) = p / 2 + 2 , with p = 3 as the unique point of optimal compactness.
  • Unified Chain (central contribution). The complete path from T n 2 to N Q ( ζ 3 ) / Q , as a single named, proved, numerically verified identity connecting five centuries of mathematics. Verified for a 500 : zero exceptions.
  • Tower of Norms. a 3 = k = 1 a N ( α k ) : every perfect cube is a stack of hexagonal norms, each satisfying ICC = 1 .
  • Cyclotomic Compatibility Index. A new invariant ICC ( n , p ) quantifying the arithmetic obstruction; ICC ( δ p ( a ) , p ) = 1 always, but ICC ( a p + b p , p ) = 0 in almost all cases for a , b 2 .
  • Window Incompatibility Theorem. Formalisation of why the hexagonal geometry of Z [ ω ] prevents the merger of two cube windows.
  • Order Theorem for δ m ( n ) . q δ m ( n ) ord q ( n ( n 1 ) 1 ) m : complete characterisation of prime divisibility. Verified over 1 793 cases.
  • LTE and Order Filters. Two independent families of necessary conditions for a 3 + b 3 = c 3 whose conjunction eliminates more than 99.9% of candidate pairs.
  • Extreme Reduction Theorem (ERT). The Order Filter eliminates every pair ( a , b ) with a 2 . Proof: finiteness of P a plus infinitude of primes (Euclid). Verified over 179 700 pairs. Scope: structural reduction to the case a = 1 ; not a proof of FLT.
  • General cyclotomic framework. δ p ( a ) = Φ p ( a , a 1 ) = N Q ( ζ p ) / Q ( a ζ p ( a 1 ) ) for all primes p.
  • Fermatian Rigidity Index. A new quantitative measure of how far ( a p + b p ) 1 / p is from an integer, decreasing to 0 as p .
  • Conceptual gradation. Three structurally distinct regimes: quadratic ( p = 2 , Pythagoras), cubic ( p = 3 , first obstruction), and Bernoulli ( p 4 ), with growing complexity C ( p ) .
  • Numerical verifications. 10 000 pairs, exponents p = 3 , , 8 , 50-digit precision: zero integer solutions. Unified Chain verified for a 500 : zero exceptions.
Its value is the explicit structural map of the Fermatian obstruction, expressed in five mathematical languages unified by a single number: 91 = 7 × 13 .
The journey from p = 3 to infinity does not discover new lands of integer solutions; it traces with increasing precision the geography of the Fermatian impossibility: from Pythagorean harmony ( p = 2 ), through cubic rigidity ( p = 3 ), to the Bernoulli complexity as p . The temptation to force the identity toward a “cubic Pythagorean theorem” fails not due to human limitation, but due to an arithmetic obstruction inherent in the structure of the integers. Recognising that limit—as Germain, Euler, and Wiles did in their respective contexts—constitutes the deepest act of mathematical understanding: knowing where each domain of validity ends, and finding in those very limits the source of new structures yet to be discovered.
This work does not prove FLT—definitively established by Wiles [1]. The case a = 1 requires Euler’s theorem. The originality lies in the perspective (reorientation of Faulhaber–Bernoulli toward individual powers via δ m ( n ) ) and in the Unified Chain (explicit articulation of the equivalence between discrete calculus, modular arithmetic, Eisenstein norm, and cyclotomic theory), not in new mathematical content in the strict sense.

Acknowledgments

The author thanks the Laboratory of the Department of Industrial Design (LIDDI), Universidad Nacional de La Plata, for institutional support and the space for the development of interdisciplinary research between mathematics, design, and education.

References

  1. A. Wiles, “Modular elliptic curves and Fermat’s last theorem,” Annals of Mathematics, vol. 141, no. 3, pp. 443–551, 1995.
  2. In H. M. Edwards, Fermat’s Last Theorem: A Genetic Introduction to Algebraic Number Theory; Springer-Verlag: New York, 1977; H. M. Edwards, Fermat’s Last Theorem: A Genetic Introduction to Algebraic Number Theory. New York: Springer-Verlag, 1977.
  3. In Nicomachus of Gerasa, Introduction to Arithmetic; (Original work ∼100 CE.); Nicomachus of Gerasa, Introduction to Arithmetic (M. L. D’Ooge, Trans.). New York: Macmillan, 1926. (Original work ∼100 CE.).
  4. K. Ireland and M. Rosen, A Classical Introduction to Modern Number Theory, 2nd ed. New York: Springer-Verlag, 1990.
  5. L. C. Washington, Introduction to Cyclotomic Fields, 2nd ed. New York: Springer-Verlag, 1997.
  6. J. H. Conway and R. K. Guy, The Book of Numbers. New York: Springer-Verlag, 1996.
  7. G. Boole, A Treatise on the Calculus of Finite Differences. London: Macmillan, 1860.
  8. J. Bernoulli, Ars Conjectandi. Basileae: Thurnisiorum, 1713.
  9. B. Taylor, Methodus Incrementorum Directa et Inversa. London: Pearson, 1715.
  10. P. Bernays, “Über die Darstellung von positiven, ganzen Zahlen durch die primitiven, binären quadratischen Formen einer nicht quadratischen Diskriminante,” Ph.D. dissertation, Georg-August-Universität Göttingen, Germany, 1912.
  11. É. Fouvry, C. Levesque, and M. Waldschmidt, “Representation of integers by cyclotomic binary forms,” Acta Arithmetica, vol. 184, no. 1, pp. 67–86, 2018.
  12. A. Mirzaei, “Lifting the Exponent Lemma (LTE),” preprint, 2011.
  13. OEIS Foundation, “Sequence A003215: Centred hexagonal numbers,” The On-Line Encyclopedia of Integer Sequences, https://oeis.org/A003215.
  14. J. Faulhaber, Academia Algebrae. Augsburg: Johann Ulrich Schönig, 1631.
  15. T. Andreescu and D. Andrica, Number Theory: Structures, Examples, and Problems. Boston: Birkhäuser, 2009.
  16. R. L. Graham, D. E. Knuth, and O. Patashnik, Concrete Mathematics: A Foundation for Computer Science, 2nd ed. Reading: Addison-Wesley, 1994.
  17. A. Baker, Transcendental Number Theory. Cambridge: Cambridge University Press, 1975.
  18. I. Niven, H. S. Zuckerman, and H. L. Montgomery, An Introduction to the Theory of Numbers, 5th ed. New York: Wiley, 1991.
  19. M. R. Murty, “Artin’s conjecture for primitive roots,” Mathematical Intelligencer, vol. 10, no. 4, pp. 59–67, 1988.
  20. In G. H. Hardy, A Mathematician’s Apology; Cambridge University Press: Cambridge, 1940; G. H. Hardy, A Mathematician’s Apology. Cambridge: Cambridge University Press, 1940.
  21. In T. L. Heath, A History of Greek Mathematics; Clarendon Press: Oxford, 1921; vol. 1, T. L. Heath, A History of Greek Mathematics, vol. 1. Oxford: Clarendon Press, 1921.
  22. R. Laubenbacher and D. Pengelley, Mathematical Expeditions: Chronicles by the Explorers. New York: Springer-Verlag, 1999.
  23. S. Singh, Fermat’s Enigma. New York: Walker and Company, 1997.
  24. I. F. Anderson, “From cumulative sum to finite difference: Nicomachus’ cubic identity as a manifestation of discrete calculus,” Preprints, 2026. Online: https://www.preprints.org/manuscript/202602.0367/v1.
  25. I. F. Anderson, “From the Pythagorean Dream to the Fermatian Obstruction: Symbolic Representation of h=a3+b33 via an Identity Derived from Nicomachus’ Cumulative Sum,” Preprints, 2026. Online: https://www.preprints.org/manuscript/202602.0367/v2.
  26. I. F. Anderson, “From the Pythagorean Dream to the Fermatian Obstruction: Symbolic Representation of h=ap+bpp via an Identity Derived from Nicomachus’ Cumulative Sum,” Preprints, 2026. Online: https://www.preprints.org/manuscript/202602.0367.
  27. I. F. Anderson, “Universal Anderson–Faulhaber–Bernoulli Identity: Internal Structure of Perfect Powers and Arithmetic Obstruction via Discrete Calculus,” Preprints, 2026. Online: https://www.preprints.org/manuscript/202602.1998/v1.
  28. I. F. Anderson, “Finite Differences of Prime Powers as Cyclotomic Norms: A Structural Bridge from Nicomachus to Euler. Universal Anderson–Faulhaber–Bernoulli Identity: Internal Structure of Perfect Powers and Arithmetic Obstruction via Discrete Calculus,” Preprints, 2026. Online: https://www.preprints.org/manuscript/202602.1998.
  29. I. F. Anderson, “From the Pythagorean Dream to the Fermatian Obstruction: Symbolic Representation of h=a3+b33 via an Identity Derived from Nicomachus’ Cumulative Sum,” Archivos de Ciencia e Investigación, vol. 2, no. 1, pp. 01–11, 2026.
Table 1. Historical deductive chain: from Pythagoras to the present work.
Table 1. Historical deductive chain: from Pythagoras to the present work.
Author Period Key contribution
Pythagoras et al. 6th c. BC h 2 = a 2 + b 2 ; figurate numbers
Euclid ∼300 BC Parameterisation of primitive triples
Nicomachus ∼100 AD S 3 ( n ) = T n 2
Faulhaber 1631 S p ( n ) polynomial for p up to 17
Fermat 1637 Conjecture: a p + b p = c p without solutions, p 3
J. Bernoulli 1713 Coefficients B j + in the general formula
Taylor 1715 Operator ∇ formalised
Euler ∼1770 FLT proved for p = 3 via Z [ ω ]
Gauss 1801 Uniqueness in Z [ ω ] ; algebraic foundations
Germain ∼1825 Systematic modular framework; FLT for p < 100
Kummer 1847 Ideal numbers; FLT for regular primes
Boole 1860 Fundamental Theorem of Discrete Calculus
Ribet 1986 Frey curve ⇒ violation of Taniyama–Shimura
Wiles 1994 Modularity theorem ⇒ FLT in general
This work 2026 UFBI; C ( p ) ; Order Theorem; ERT; Unified Chain;
Tower of Norms; ICC; Window Theorem
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated