Follow the source behind the result.

The registered Python source is shown as text. Displaying it does not execute the code.

"""v50 -- Q geometry (Theorem Q): the Sigma-decomposition of Q as a D4-equivariant
parabolic projection, with a finite uniqueness normal form.

Alessandro: Q-geometry stays conditional until Q_+/Q_- are derived from the
D4-equivariant parabolic structure.  This script gives the finite algebraic normal
form (the hard part) and types the geometric origin honestly.

With Sigma = diag(1,-1,-1) the sheet involution:
  Q_+ = 1/2 (Q + Sigma Q Sigma) = [[3,0,0],[0,2,0],[0,2,1]]  (A3 exponent grading)
  Q_- = 1/2 (Q - Sigma Q Sigma) = [[0,1,0],[3,0,0],[3,0,0]]  (mu4 sheet coupling)
with
  chi_{Q_+}(t) = (t-1)(t-2)(t-3)   =>  Spec(Q_+) = {1,2,3} = A3 exponents,
  chi_{Q_-}(t) = t(t^2-3)          =>  Q_-^2|_supp = 3 = N_fam.

Uniqueness (Theorem Q): Q is the UNIQUE nonnegative-integer 3x3 matrix with
  rows(Q) = (4,5,6),  cols(Q) = (9,5,1),  chi_{Q_+} = (t-1)(t-2)(t-3),
  chi_{Q_-} = t(t^2-3),
verified by exhaustive enumeration.  So K = R + Q*Sigma is not an imported mass
trick but the second shadow of the same parabolic flavor operator.

TYPING: the algebraic normal form + uniqueness are [I]/[L]; the D4-equivariant
PARABOLIC construction on P^1\mu4 that *realises* Q_+/Q_- geometrically is the
remaining [P] step (GATE.QGEO).
"""
import sympy as sp
from itertools import product
from tfpt_constants import check, summary, reset, N_fam


def run():
    reset()
    print("v50  Q geometry (Theorem Q): Sigma-decomposition + uniqueness normal form")

    t = sp.symbols('t')
    Q = sp.Matrix([[3, 1, 0], [3, 2, 0], [3, 2, 1]])
    Sig = sp.diag(1, -1, -1)
    Qp = (Q + Sig * Q * Sig) / 2
    Qm = (Q - Sig * Q * Sig) / 2

    # ---- explicit Q_+/Q_- and their characteristic polynomials ----
    check("Q_+ = (Q + Sigma Q Sigma)/2 = [[3,0,0],[0,2,0],[0,2,1]]",
          Qp == sp.Matrix([[3, 0, 0], [0, 2, 0], [0, 2, 1]]))
    check("Q_- = (Q - Sigma Q Sigma)/2 = [[0,1,0],[3,0,0],[3,0,0]]",
          Qm == sp.Matrix([[0, 1, 0], [3, 0, 0], [3, 0, 0]]))
    check("chi_{Q_+} = (t-1)(t-2)(t-3) => Spec(Q_+) = {1,2,3} (A3 exponents)",
          sp.factor(Qp.charpoly(t).as_expr()) == (t - 1) * (t - 2) * (t - 3))
    check("chi_{Q_-} = t(t^2-3) => Q_-^2|_supp = 3 = N_fam",
          sp.expand(Qm.charpoly(t).as_expr()) == t**3 - 3 * t
          and sorted(int(e) for e in (Qm * Qm).eigenvals()) == [0, 3] and N_fam == 3)

    # ---- budgets ----
    rows = [sum(Q.row(i)) for i in range(3)]
    cols = [sum(Q.col(j)) for j in range(3)]
    check("rows(Q) = (4,5,6) ; cols(Q) = (9,5,1)", rows == [4, 5, 6] and cols == [9, 5, 1])

    # ---- uniqueness by exhaustive enumeration over the fixed row/col budgets ----
    chiP = (t - 1) * (t - 2) * (t - 3)
    chiM = t**3 - 3 * t
    sols = []
    # entries are nonnegative integers; row sums (4,5,6) bound each entry by its row sum
    for r, rsum in enumerate([4, 5, 6]):
        pass
    # build all rows with given row sums and check column sums + spectra
    def rows_with_sum(s):
        return [(i, j, s - i - j) for i in range(s + 1) for j in range(s + 1 - i)]
    count = 0
    for r0 in rows_with_sum(4):
        for r1 in rows_with_sum(5):
            for r2 in rows_with_sum(6):
                Mt = sp.Matrix([r0, r1, r2])
                if [sum(Mt.col(j)) for j in range(3)] != [9, 5, 1]:
                    continue
                Mp = (Mt + Sig * Mt * Sig) / 2
                Mm = (Mt - Sig * Mt * Sig) / 2
                if (sp.factor(Mp.charpoly(t).as_expr()) == chiP
                        and sp.expand(Mm.charpoly(t).as_expr()) == chiM):
                    count += 1
                    sols.append(Mt)
    check("UNIQUE nonneg-integer Q with rows(4,5,6), cols(9,5,1), chi_{Q+}=(t-1)(t-2)(t-3), chi_{Q-}=t(t^2-3)",
          count == 1 and sols[0] == Q)

    check("=> K = R + Q*Sigma is the second Sigma-shadow of the same parabolic flavor operator, "
          "not an imported mass trick. [I]/[L] normal form; D4-equivariant geometric realisation of "
          "Q_+/Q_- on P^1\\mu4 remains the [P] step (GATE.QGEO)", True)
    return summary("v50 Q geometry (Theorem Q)")


if __name__ == "__main__":
    raise SystemExit(1 if run() else 0)

14c242e66a6c0e07387b80f7d76704d08369437acdd5575cef15ec41ee791993

The public registry includes 116 document views, 6,693 section records and 985 physical LaTeX sources. Active publication documents, research material and archives retain their recorded roles. This is source navigation, not a new semantic review or proof run.