Teoria kompozycji form kwadratowych stanowi jedno z najważniejszych osiągnięć w historii teorii liczb, będąc jednocześnie doskonałym przykładem na to, jak małe obserwacje mogą prowadzić do wielkich odkryć. Początki tej teorii sięgają czasów starożytnych, a jej rozwój miał miejsce na przestrzeni wieków, poczynając od Fermata, przez Eulera, aż po najważniejsze wkłady Lagrange'a i Legendre'a. Właśnie Legendre (1798) położył fundamenty pod formalną analizę kompozycji form kwadratowych, traktując je jako problem związany z dzielnikami kwadratowymi tej samej dyskryminanty.

Legendre w swojej pracy z 1798 roku zaprezentował problem polegający na znalezieniu nowego dzielnika kwadratowego dla dwóch danych form kwadratowych o tej samej formule t2+au2t^2 + au^2. Z perspektywy współczesnej, jest to nic innego jak klasyczny problem komponowania form kwadratowych o tej samej dyskryminancie D=4aD = -4a. Choć Legendre był pionierem w tym zakresie, to dopiero Gauss rozwinął pełną teorię kompozycji form kwadratowych, rozważając przypadek ogólny i wprowadzając pojęcie współczynników kompozycji form.

Gauss zrealizował rozwiązanie w postaci obliczenia współczynników dla form kwadratowych Q1Q2Q_1 \circ Q_2, gdzie Q1Q_1 i Q2Q_2 to formy kwadratowe, a \circ oznacza operację kompozycji. Jednakże pełne rozwiązanie dla kongruencji (90.8)(90.8) i współczynników kompozycji, jak {sj,sj}\{s_j, s'_j\}, zostało podane przez Arndta, co spowodowało, że niektóre elementy teorii kompozycji są tradycyjnie przypisywane właśnie jemu.

Przełomowym momentem w tej teorii było wprowadzenie przez Dirichleta (1851) pojęcia kongruencji, które miało kluczowe znaczenie dla dalszego rozwoju matematyki algebraicznej. W pracy Dirichleta pojawiło się pojęcie grupy Γ\Gamma oraz jej zastosowanie do opisu równoważności form kwadratowych, co pozwoliło na bardziej formalne podejście do zagadnienia.

Kiedy spojrzymy na sposób, w jaki formy kwadratowe są składane, można dostrzec pewną analogię do faktoryzacji liczb pierwszych w dziedzinach algebraicznych. Na przykład, w przypadku pola Q(D)Q(\sqrt{D}), dekompozycja liczb pierwszych w idealnych pierwiastkach może przebiegać na różne sposoby w zależności od wartości symbolu Kroneckera κD(p)\kappa_D(p). Jeśli κD(p)=1\kappa_D(p) = 1, wtedy liczba pierwsza pp rozkłada się na dwa różne ideały pierwsze, jeśli κD(p)=0\kappa_D(p) = 0, to jest to kwadrat ideału pierwszego, a jeśli κD(p)=1\kappa_D(p) = -1, to pp jest ideałem maksymalnym.