on veut seulement être sûr que c’est elle, ne pas se tromper sur l’identité […] cet extrait algébrique... — Albertine disparue
Taking a cue from Marcel Proust, I propose an algebraic approach to partition identities, via regular grammars. Many partition identities take the form ‘A ( n ) = B ( n ) A(n)=B(n) A ( n ) = B ( n ) ’, where A ( n ) A(n) A ( n ) and B ( n ) B(n) B ( n ) count the partitions of n n n satisfying constraints given by regular expressions. I give in Theorem A a method that (if it doesn’t run forever) produces algebraically verifiable certificates of such identities. As an application, I prove in Theorem B the three-row, level-two case of Primc’s conjecture.
A teaser: Euler's identity
Here is an illustration of how it works. In Euler’s identity, the archetypal example, A ( n ) A(n) A ( n ) counts the partitions of n n n into odd parts and B ( n ) B(n) B ( n ) those into distinct parts. Draw the Ferrers diagram of λ ⊢ n \lambda\vdash n λ ⊢ n with λ i \lambda_i λ i boxes in row i i i , the longest row on top, and read its lower right boundary from bottom left to top right, writing R R R for a unit step to the right and U U U for a unit step up. Nonempty partitions thus correspond to words in R { R , U } ∗ U R\{R,U\}^{*}U R { R , U } ∗ U ; we ignore the empty partition unless stated otherwise. The odd and distinct constraints correspond to the regular languages
L o = R ( U + R R ) ∗ U , L d = R ( R + U R ) ∗ U L_o=R(U+RR)^{*}U,\qquad L_d=R(R+UR)^{*}U L o = R ( U + R R ) ∗ U , L d = R ( R + U R ) ∗ U
whose generating functions are the series
H o = R ( 1 − U − R 2 ) − 1 U , H d = R ( 1 − R − U R ) − 1 U H_o=R(1-U-R^2)^{-1}U,\qquad H_d=R(1-R-UR)^{-1}U H o = R ( 1 − U − R 2 ) − 1 U , H d = R ( 1 − R − U R ) − 1 U
in the skew field of fractions of the quantum plane A = Q [ q ] ⟨ U , R ⟩ / ( R U − q U R ) \mathbb{A}=\mathbb{Q}[q]\langle U,R\rangle/(RU-qUR) A = Q [ q ] ⟨ U , R ⟩ / ( R U − q U R ) . Expanded as sums over words, these rational expressions define elements of the q q q -adic completion A ^ \widehat{\mathbb{A}} A of A \mathbb{A} A : the outer letters R R R and U U U make the sums converge, see §3.5. This ring is convenient because the boundary word of λ = ( λ 1 ≥ ⋯ ≥ λ ℓ > 0 ) \lambda=(\lambda_1\geq\cdots\geq\lambda_\ell>0) λ = ( λ 1 ≥ ⋯ ≥ λ ℓ > 0 ) satisfies
w ( λ ) = R λ ℓ U R λ ℓ − 1 − λ ℓ U ⋯ R λ 1 − λ 2 U = q n U R λ 1 in A . (1) w(\lambda)=R^{\lambda_\ell}UR^{\lambda_{\ell-1}-\lambda_\ell}U\cdots R^{\lambda_1-\lambda_2}U=q^{n}UR^{\lambda_1}
\quad\text{in }\mathbb{A}. \tag*{(1)} w ( λ ) = R λ ℓ U R λ ℓ − 1 − λ ℓ U ⋯ R λ 1 − λ 2 U = q n U R λ 1 in A . (1)
every pair formed by a letter R R R and a later letter U U U contributes a factor q q q , and the U U U recording a part j j j has exactly j j j such predecessors.
Every s ∈ A ^ s\in\widehat{\mathbb{A}} s ∈ A has a unique normal form ∑ i , j s i , j U i R j \sum_{i,j}s_{i,j}U^iR^j ∑ i , j s i , j U i R j with s i , j ∈ Q [ [ q ] ] s_{i,j}\in\mathbb{Q}[[q]] s i , j ∈ Q [[ q ]] . Its normal-ordered evaluation
ε ( s ) = ∑ i , j s i , j ∈ Q [ [ q ] ] \varepsilon(s)=\sum_{i,j}s_{i,j}\in\mathbb{Q}[[q]] ε ( s ) = i , j ∑ s i , j ∈ Q [[ q ]]
is linear but not multiplicative: ε ( R U ) = q ≠ ε ( R ) ε ( U ) \varepsilon(RU)=q\ne\varepsilon(R)\varepsilon(U) ε ( R U ) = q = ε ( R ) ε ( U ) . The generating functions of odd and of distinct partitions are ε ( H o ) \varepsilon(H_o) ε ( H o ) and ε ( H d ) \varepsilon(H_d) ε ( H d ) , and Euler’s identity is ε ( H o − H d ) = 0 \varepsilon(H_o-H_d)=0 ε ( H o − H d ) = 0 . A certificate of it is a decomposition
H o − H d = ( 1 − U ) a + b ( 1 − R ) , a , b ∈ A ^ . H_o-H_d=(1-U)a+b(1-R),\qquad a,b\in\widehat{\mathcal A}. H o − H d = ( 1 − U ) a + b ( 1 − R ) , a , b ∈ A .
Left multiplication by U U U and right multiplication by R R R shift normal exponents without changing coefficients, so ε ( ( 1 − U ) a ) = ε ( b ( 1 − R ) ) = 0 \varepsilon((1-U)a)=\varepsilon(b(1-R))=0 ε (( 1 − U ) a ) = ε ( b ( 1 − R )) = 0 by telescoping, and such a decomposition proves the identity.
The series a , b a,b a , b are typically not rational, but they can be chosen holonomic , that is, defined by finitely many q q q -difference equations. The normal symbol of s s s is the commutative series ∑ i , j s i , j x i R j \sum_{i,j}s_{i,j}x^iR^j ∑ i , j s i , j x i R j , and its marginal is obtained by setting x = 1 x=1 x = 1 . Since U i R j ⋅ U = q j U i + 1 R j U^iR^j\cdot U=q^jU^{i+1}R^j U i R j ⋅ U = q j U i + 1 R j , right multiplication by U U U acts on symbols by g ↦ x Υ g g\mapsto x\,\Upsilon g g ↦ x Υ g , where
( Υ g ) ( x , R ) = g ( x , q R ) , (\Upsilon g)(x,R)=g(x,qR), ( Υ g ) ( x , R ) = g ( x , q R ) ,
and right multiplication by R R R acts by g ↦ R g g\mapsto Rg g ↦ R g ; on marginals, U U U acts by Υ \Upsilon Υ alone. Since Υ R = q R Υ \Upsilon R=qR\Upsilon Υ R = q R Υ , the assignment U ↦ Υ U\mapsto\Upsilon U ↦ Υ , R ↦ R R\mapsto R R ↦ R embeds A \mathcal A A into the opposite of the algebra of q q q -difference operators in R R R . Consequently, a series defined by a linear system over A \mathcal A A , such as H o H_o H o or H d H_d H d , has a symbol defined by a linear system of q q q -difference equations in R R R , with x x x as a parameter. A certificate then consists of such systems for a a a and b b b , finitely many values selecting their solutions, and finite data guaranteeing that their coefficients tend q q q -adically to 0. For instance, write ν \nu ν for the q q q -adic valuation, and call a series q q q -integral if it has no negative powers of q q q . If b = q Υ b ′ b=q\Upsilon b' b = q Υ b ′ with b ′ b' b ′ q q q -integral, that is, b ′ ∈ Q [ [ q , R ] ] b'\in\mathbb{Q}[[q,R]] b ′ ∈ Q [[ q , R ]] , then [ R n ] b = q n + 1 [ R n ] b ′ [R^n]b=q^{n+1}[R^n]b' [ R n ] b = q n + 1 [ R n ] b ′ has valuation at least n + 1 n+1 n + 1 , so b ∈ A ^ b\in\widehat{\mathcal A} b ∈ A . A holonomic kernel for b ′ b' b ′ , together with finite data certifying its integrality, is such a guarantee.
For Euler’s identity, moving the outer letters inwards gives H o = q U ( 1 − q U − q 2 R 2 ) − 1 R H_o=qU(1-qU-q^2R^2)^{-1}R H o = q U ( 1 − q U − q 2 R 2 ) − 1 R and H d = q U ( 1 − q R − q 2 U R ) − 1 R H_d=qU(1-qR-q^2UR)^{-1}R H d = q U ( 1 − q R − q 2 U R ) − 1 R . Hence H o − H d H_o-H_d H o − H d has normal symbol
D ( x , R ) = q R ( m o , x − m d , x ) , { ( 1 − q x Υ − q 2 R 2 ) m o , x = x , ( 1 − q R − q 2 x R Υ ) m d , x = x . D(x,R)=qR(m_{o,x}-m_{d,x}),
\qquad
\begin{cases}
(1-qx\Upsilon-q^2R^2)m_{o,x}=x,\\
(1-qR-q^2xR\Upsilon)m_{d,x}=x.
\end{cases} D ( x , R ) = q R ( m o , x − m d , x ) , { ( 1 − q x Υ − q 2 R 2 ) m o , x = x , ( 1 − q R − q 2 x R Υ ) m d , x = x .
Both operators are congruent to 1 modulo q q q , so these equations have unique solutions in Q [ x , R ] [ [ q ] ] \mathbb{Q}[x,R][[q]] Q [ x , R ] [[ q ]] . The canonical a , b a,b a , b are given by
( 1 − x ) a = D ( x , R ) − D ( 1 , R ) , ( 1 − R ) b = D ( 1 , R ) , (1-x)a=D(x,R)-D(1,R),\qquad(1-R)b=D(1,R), ( 1 − x ) a = D ( x , R ) − D ( 1 , R ) , ( 1 − R ) b = D ( 1 , R ) ,
and are holonomic. It remains to certify their convergence in A ^ \widehat{\mathcal A} A . The integral kernels (§ 3.7 )
u = 1 + R 2 α ( Υ ) u , v = 1 + R 2 β ( Υ ) v , α ( U ) = ( 1 − q U ) − 1 , β ( U ) = ( 1 + q U ) ( 1 + q 2 U ) , u=1+R^2\alpha(\Upsilon)u,\qquad v=1+R^2\beta(\Upsilon)v,\qquad\alpha(U)=(1-qU)^{-1},\qquad\beta(U)=(1+qU)(1+q^2U), u = 1 + R 2 α ( Υ ) u , v = 1 + R 2 β ( Υ ) v , α ( U ) = ( 1 − q U ) − 1 , β ( U ) = ( 1 + q U ) ( 1 + q 2 U ) ,
select u , v ∈ Q [ [ q , R ] ] u,v\in\mathbb{Q}[[q,R]] u , v ∈ Q [[ q , R ]] , and χ ( U ) = 1 − q 2 U \chi(U)=1-q^2U χ ( U ) = 1 − q 2 U is a dilation witness for θ = β / α \theta=\beta/\alpha θ = β / α :
θ ( U ) χ ( U ) = ( 1 − q 2 U 2 ) ( 1 − q 4 U 2 ) = χ ( U 2 ) χ ( q 2 U 2 ) . \theta(U)\chi(U)=(1-q^2U^2)(1-q^4U^2)=\chi(U^2)\chi(q^2U^2). θ ( U ) χ ( U ) = ( 1 − q 2 U 2 ) ( 1 − q 4 U 2 ) = χ ( U 2 ) χ ( q 2 U 2 ) .
Thus u − v = q 2 Υ w u-v=q^2\Upsilon w u − v = q 2 Υ w with w ∈ Q [ [ q , R ] ] w\in\mathbb{Q}[[q,R]] w ∈ Q [[ q , R ]] , see Proposition 3.7 , and the equations of u u u , v v v , m o , 1 m_{o,1} m o , 1 and m d , 1 m_{d,1} m d , 1 give
b = ( 1 + R ( 1 − q Υ ) − 1 ) u − ( 1 + R ( 1 + q Υ ) ) v = q Υ b ′ (2) b=(1+R(1-q\Upsilon)^{-1})u-(1+R(1+q\Upsilon))v=q\Upsilon b' \tag*{(2)} b = ( 1 + R ( 1 − q Υ ) − 1 ) u − ( 1 + R ( 1 + q Υ )) v = q Υ b ′ (2)
with b ′ = q w + R ( 1 − q Υ ) − 1 ( w + Υ v ) ∈ Q [ [ q , R ] ] b'=qw+R(1-q\Upsilon)^{-1}(w+\Upsilon v)\in\mathbb{Q}[[q,R]] b ′ = q w + R ( 1 − q Υ ) − 1 ( w + Υ v ) ∈ Q [[ q , R ]] . Thus b ∈ A ^ b\in\widehat{\mathcal A} b ∈ A ; moreover a ∈ A ^ a\in\widehat{\mathcal A} a ∈ A by Theorem A(2), and the certificate is complete. Corollary 4.2 gives the details, and §4.1 applies the same mechanism to Glaisher’s identity, of which Euler’s is the case λ = 2 \lambda=2 λ = 2 .
Main results
The scheme of §1 applies to any two partition languages in R ⟨ R , U ⟩ ∗ U \mathcal R\langle R,U\rangle^*U R ⟨ R , U ⟩ ∗ U recognized by unambiguous finite automata, and more generally to weighted automata recording colours or other statistics; a boundary word is then counted with the total weight of its accepting paths. Section 3.2 explains how classical constraints are compiled into such automata. Stripping the first R R R and the last U U U yields interior automata with s 1 s_1 s 1 and s 2 s_2 s 2 states. They give a normalized realization
w Q = r , S = w e , Q ≡ I ( m o d q ) , w\mathbf{Q}=r,\qquad S=we,\qquad\mathbf{Q}\equiv I\pmod q, w Q = r , S = w e , Q ≡ I ( mod q ) ,
of dimension s = s 1 + s 2 s=s_1+s_2 s = s 1 + s 2 , of the difference S S S of the two series; see § 3.5 . We write K = Q ( q ) \mathbb{K}=\mathbb{Q}(q) K = Q ( q ) , and use the Pochhammer symbols ( a ; q ) n = ∏ j = 0 n − 1 ( 1 − a q j ) (a;q)_n=\prod_{j=0}^{n-1}(1-aq^j) ( a ; q ) n = ∏ j = 0 n − 1 ( 1 − a q j ) and ( a 1 , … , a t ; q ) n = ∏ i ( a i ; q ) n (a_1,\ldots,a_t;q)_n=\prod_i(a_i;q)_n ( a 1 , … , a t ; q ) n = ∏ i ( a i ; q ) n , also for n = ∞ n=\infty n = ∞ whenever the product converges q q q -adically.
Theorem A. Let H 1 , H 2 H_1,H_2 H 1 , H 2 be the path series (§ 3.1 ) of finite weighted automata with weights in Q ( q ) \mathbb{Q}(q) Q ( q ) regular at q = 0 q=0 q = 0 , all of whose accepted words lie in R { R , U } ∗ U R\{R,U\}^{*}U R { R , U } ∗ U . Let S = H 1 − H 2 S=H_1-H_2 S = H 1 − H 2 have normal symbol D ( x , R ) D(x,R) D ( x , R ) and a normalized realization of dimension s s s , and put
a ( x , R ) = D ( x , R ) − D ( 1 , R ) 1 − x , b ( R ) = D ( 1 , R ) 1 − R = ∑ n ≥ 0 b n R n . a(x,R)=\frac{D(x,R)-D(1,R)}{1-x},\qquad b(R)=\frac{D(1,R)}{1-R}=\sum_{n\geq0}b_nR^n. a ( x , R ) = 1 − x D ( x , R ) − D ( 1 , R ) , b ( R ) = 1 − R D ( 1 , R ) = n ≥ 0 ∑ b n R n . (1) The series a a a and b b b are components of the unique solutions of explicit square inhomogeneous systems of linear q q q -difference equations, of sizes 2 s + 1 2s+1 2 s + 1 and s + 1 s+1 s + 1 , computed from the automata. In particular, a a a is holonomic in ( x , R ) (x,R) ( x , R ) , b b b is holonomic in R R R , and every coefficient b n b_n b n lies in K \mathbb{K} K .
(2) The coefficients of a a a satisfy ν ( a i , j ) ≥ i + j \nu(a_{i,j})\geq i+j ν ( a i , j ) ≥ i + j ; in particular a ∈ A ^ a\in\widehat{\mathcal{A}} a ∈ A .
(3) The following are equivalent:
(a) ε ( H 1 ) = ε ( H 2 ) \varepsilon(H_1)=\varepsilon(H_2) ε ( H 1 ) = ε ( H 2 ) ;
(b) b n → 0 b_n\to0 b n → 0 q q q -adically;
(c) ν ( b n ) ≥ n + 1 \nu(b_n)\geq n+1 ν ( b n ) ≥ n + 1 for all n n n ;
(d) b = q Υ φ b=q\Upsilon\varphi b = q Υ φ for some φ ∈ Q [ [ q , R ] ] \varphi\in\mathbb{Q}[[q,R]] φ ∈ Q [[ q , R ]] ;
(e) b ∈ Q [ R ] [ [ q ] ] b\in\mathbb{Q}[R][[q]] b ∈ Q [ R ] [[ q ]] , so that S = ( 1 − U ) a + b ( 1 − R ) S=(1-U)a+b(1-R) S = ( 1 − U ) a + b ( 1 − R ) in A ^ \widehat{\mathcal{A}} A .
In (d), φ = q − 1 b ( R / q ) \varphi=q^{-1}b(R/q) φ = q − 1 b ( R / q ) , and φ \varphi φ is holonomic. In (e), ( a , b ) (a,b) ( a , b ) is the only pair with S = ( 1 − U ) a + b ( 1 − R ) S=(1-U)a+b(1-R) S = ( 1 − U ) a + b ( 1 − R ) , a ∈ A ^ a\in\widehat{\mathcal{A}} a ∈ A and b ∈ Q [ R ] [ [ q ] ] b\in\mathbb{Q}[R][[q]] b ∈ Q [ R ] [[ q ]] .
Theorem A is proved in § 3.6 . It reduces a partition identity to the integrality of a single holonomic series φ \varphi φ , defined by finite data. The formal boundaries a , b a,b a , b exist whether or not the identity holds, and a false identity is detected after computing finitely many b n b_n b n : since ν ( ε ( S ) − b n ) ≥ n + 1 \nu(\varepsilon(S)-b_n)\geq n+1 ν ( ε ( S ) − b n ) ≥ n + 1 unconditionally (Lemma 3.5 ), the identity is false if and only if ν ( b n ) < n + 1 \nu(b_n)<n+1 ν ( b n ) < n + 1 for some n n n . For a true identity, it remains to prove integrality of φ \varphi φ . A direct certificate (§ 3.7 ) does so by finitely many exact checks on holonomic series; alternatively, any independent proof of the identity yields (3d) (Corollary 3.6 ).
This separates discovery from verification, in the spirit of Wilf–Zeilberger certificates [23 , 24 , 25 ] . Certificates are found by finite linear-algebra searches (§ 3.8 and Appendix § 5 ) for integral kernels, dilation witnesses, finite recurrences, telescoping certificates and quadratic invariants, and then checked by exact algebraic manipulations. Since the quantum plane records area, many combinatorial descriptions — allowed residues, controlled differences between coloured parts, forbidden patterns — compile into rational series to which Theorem A applies (§ 3.2 ).
Many partition identities equate restrictions on differences or multiplicities with restrictions on the residue classes of the parts. Classical examples include Euler’s identity between distinct and odd parts (modulo 2) [8 §326] , the Rogers–Ramanujan identities (modulo 5) [19 ] , Schur’s theorem (modulo 6) [21 ] , and the Göllnitz–Gordon identities (modulo 8) [11 , 13 ] . Later examples include Capparelli’s identity (modulo 12), conjectured in his 1988 thesis [6 ] and first proved by Andrews [2 ] , and Nandi’s three identities (modulo 14), conjectured in 2014 [16 §8.1] and proved by Takigiku and Tsuchioka [22 ] . These examples suggest regular languages as a common framework; Nandi’s difference conditions already involve forbidden patterns of unbounded length. Primc’s coloured partition conjectures [18 ] belong to the same tradition; our main application, Theorem B below, is their three-row case at level two, with modulus 8.
Primc's identities at level two
We prove the three-row, level-two case of Primc’s conjecture [18 Conjecture 2.1] . It equates two classes of coloured partitions: one defined by restrictions on neighbouring multiplicities, the other by congruence classes modulo eight.
For the first class, every positive even part has two colours, called top and bottom , and every odd part has one colour, called middle . Write τ n \tau_n τ n , β n \beta_n β n for the multiplicities of the top and bottom parts 2 n 2n 2 n , and μ n \mu_n μ n for the multiplicity of the part 2 n − 1 2n-1 2 n − 1 . Fix nonnegative integers k 0 k_0 k 0 , k 1 k_1 k 1 , k 2 k_2 k 2 of sum 2, the level , and set τ 0 = k 2 \tau_0=k_2 τ 0 = k 2 , β 0 = k 0 \beta_0=k_0 β 0 = k 0 . The partition is admissible if
max ( τ n − 1 , τ n ) + μ n + max ( β n − 1 , β n ) ≤ 2 ( n ≥ 1 ) . (3) \max(\tau_{n-1},\tau_n)+\mu_n+\max(\beta_{n-1},\beta_n)\le2 \qquad(n\ge1). \tag*{(3)} max ( τ n − 1 , τ n ) + μ n + max ( β n − 1 , β n ) ≤ 2 ( n ≥ 1 ) . (3)
Equivalently, each of the four sums obtained by choosing one of τ n − 1 , τ n \tau_{n-1},\tau_n τ n − 1 , τ n , one of β n − 1 , β n \beta_{n-1},\beta_n β n − 1 , β n , and μ n \mu_n μ n is at most two; in Primc’s description by three-row arrays, these are the sums along the four downward paths through μ n \mu_n μ n (see § 3.2 ). The value k 1 = 2 − k 0 − k 2 k_1=2-k_0-k_2 k 1 = 2 − k 0 − k 2 imposes no further condition. There are six boundary conditions; write Π k 0 , k 1 , k 2 ( q ) \Pi_{k_0,k_1,k_2}(q) Π k 0 , k 1 , k 2 ( q ) for the generating function of the corresponding admissible partitions, including the empty one.
Theorem B. For each N ≥ 0 N\ge0 N ≥ 0 , the number of admissible coloured partitions of N N N with the boundary condition in the first column equals the number of coloured partitions of N N N with the following allowed residues modulo eight:
Each allowed colour has unrestricted multiplicity, and all other residues are forbidden. The cases ( 0 , 0 , 2 ) (0,0,2) ( 0 , 0 , 2 ) and ( 0 , 1 , 1 ) (0,1,1) ( 0 , 1 , 1 ) follow by interchanging the top and bottom colours.
For example, for the boundary ( 2 , 0 , 0 ) (2,0,0) ( 2 , 0 , 0 ) the local rule (3) forbids the part 1 and the top-coloured part 2, and the theorem states
Π 2 , 0 , 0 ( q ) = 1 ( q 2 , q 3 , q 4 , q 4 , q 5 , q 6 ; q 8 ) ∞ . \Pi_{2,0,0}(q)=\frac{1}{(q^2,q^3,q^4,q^4,q^5,q^6;q^8)_\infty}. Π 2 , 0 , 0 ( q ) = ( q 2 , q 3 , q 4 , q 4 , q 5 , q 6 ; q 8 ) ∞ 1 .
Figure 1 shows the two automata for this boundary; the tags on the letters U U U make them unambiguous and are forgotten in A A A (§ 3.2 ).
Figure 1. The two automata for boundary (2, 0, 0). Top: the residue automaton; U 1 , U 2 U_1,U_2 U 1 , U 2 distinguish the colours at 4. Bottom: circles record the preceding even pair τ β \tau\beta τ β , with the number of middle parts as exponent; rectangles record the current even pair. Double borders mark acceptance. Identify the wavy cuts: each outgoing R R R -edge resumes at the same vertical position on the left. The ten plain rectangles also have an R R R -edge to circle 00 (omitted). Repeated labels have different allowed continuations, determined by (3) . All tagged U U U ’s become U U U in A A A .
Primc’s Conjecture 2.1 only asserts that these generating functions are periodic products; the four products of Theorem B are those he found numerically [18 p. 53] . Kanade, Russell, Tsuchioka and Warnaar proposed an explicit product for every width and level [15 Conjecture 4.2] , whose case n = k = 2 n=k=2 n = k = 2 , with the same labelling of the boundary, gives the four products above. They note that the three-row case appears to be open, and they prove functional equations [15 Propositions 4.4 and 4.5] and quadruple-sum expressions [15 Theorem 4.10] for the four generating functions, refined by the number of parts, but not their product evaluations. The level-one case, of modulus 6, follows from the coloured partition theorem [14 Theorem 1.2] , and also from a multisum identity of Russell [20 ] together with Bressoud’s identities; see [15 §4] .
Both partition classes are recognized by small automata: the first needs only to remember the preceding even multiplicities, and the second only the part size modulo eight and its colour. The proof in §4.2 reduces the four-state transfer system of the local conditions to two Göllnitz–Gordon systems, by a quadratic invariant and two rational telescoping certificates, all found by the searches of §3.8. The classical Göllnitz–Gordon identities and Jacobi’s triple product then give the products in the theorem. This proves ϵ ( H 1 ) = ϵ ( H 2 ) \epsilon(H_1)=\epsilon(H_2) ϵ ( H 1 ) = ϵ ( H 2 ) for the path series of the two automata, and Corollary 3.6 turns it into convergent holonomic certificates for them. The new step is the reduction of Primc’s four-state system to these classical identities.
In §4 we also treat Glaisher’s identity, the two Rogers–Ramanujan identities, Bessenrodt’s refinement of Euler’s identity by a global partition statistic, Berkovich–Uncu identities involving part sizes and indices, and an identity between two five-state automata, found by computer search, which fits none of the previous patterns.
Holonomic certificates and their discovery
Partition languages
General background on partitions and q q q -series may be found in Andrews [1 ] . In the boundary word w ( λ ) w(\lambda) w ( λ ) of (1) , R R R raises the current height , the number of letters R R R read so far, and U U U records a part equal to the height; we write ℓ ( λ ) \ell(\lambda) ℓ ( λ ) for the number of parts of λ \lambda λ and ∣ λ ∣ |\lambda| ∣ λ ∣ for their sum.
Regular languages are those described by regular expressions, built from letters by concatenation, union and Kleene star; equivalently, those recognized by finite automata, namely finite directed graphs with an initial state, a set of final states and edges labelled R R R or U U U , whose language consists of the labels of paths from the initial state to a final state. If T T T is the transition matrix of the automaton, whose entries are sums of letters, and r r r , e e e are the indicator vectors of the initial and final states, then the path series r ( I − T ) − 1 e = ∑ n ≥ 0 r T n e r(I-T)^{-1}e = \sum_{n\geq0} rT^n e r ( I − T ) − 1 e = ∑ n ≥ 0 r T n e counts the accepting paths by their labels; for an unambiguous automaton, it is the series of its language. Weighted automata allow other coefficients, for instance multiplicities of colours or markers of statistics.
Compiling constraints into automata
Let f h f_h f h be the multiplicity of the part h h h in λ \lambda λ . Grouping the letters of w ( λ ) w(\lambda) w ( λ ) by height gives
w ( λ ) = R U f 1 R U f 2 ⋯ R U f m , m = λ 1 , f m ≥ 1. (4) w(\lambda)=RU^{f_1}RU^{f_2}\cdots RU^{f_m},\qquad m=\lambda_1,\quad f_m\geq1. \tag*{(4)} w ( λ ) = R U f 1 R U f 2 ⋯ R U f m , m = λ 1 , f m ≥ 1. (4)
Reading a boundary word from left to right therefore scans the multiplicities in increasing order of the parts, and a constraint is regular as soon as it can be checked during this scan with finite memory. We make this precise for the constraints occurring in Rogers–Ramanujan type identities: congruence conditions on the parts, and difference conditions on the multiplicities, possibly of coloured parts.
Colours and tagged letters. Following Primc [18 ] , let N = N 1 ∪ ⋯ ∪ N r \mathcal{N}=\mathcal{N}_1\cup\cdots\cup\mathcal{N}_r N = N 1 ∪ ⋯ ∪ N r be a multiset of positive integers, the elements of N c \mathcal{N}_c N c having colour c c c . A coloured partition on N \mathcal{N} N is given by multiplicities f h , c f_{h,c} f h , c , which vanish unless h ∈ N c h\in\mathcal{N}_c h ∈ N c ; we write f h = ( f h , 1 , … , f h , r ) f_h=(f_{h,1},\ldots,f_{h,r}) f h = ( f h , 1 , … , f h , r ) . Encode it by the word
w ( λ ) = R U 1 f 1 , 1 ⋯ U r f 1 , r R U 1 f 2 , 1 ⋯ U r f 2 , r ⋯ R U 1 f m , 1 ⋯ U r f m , r w(\lambda)=RU_1^{f_{1,1}}\cdots U_r^{f_{1,r}}\,RU_1^{f_{2,1}}\cdots U_r^{f_{2,r}}\cdots RU_1^{f_{m,1}}\cdots U_r^{f_{m,r}} w ( λ ) = R U 1 f 1 , 1 ⋯ U r f 1 , r R U 1 f 2 , 1 ⋯ U r f 2 , r ⋯ R U 1 f m , 1 ⋯ U r f m , r
over the alphabet { R , U 1 , … , U r } \{R,U_1,\ldots,U_r\} { R , U 1 , … , U r } of tagged letters, in which the letters at each height appear in increasing order of colour. This encoding is injective. The tags make the automata below unambiguous; they are forgotten in A A A , where each U c U_c U c becomes U U U , so that path series count coloured partitions, with multiplicity when several of them have the same uncoloured boundary word.
Congruence conditions. Suppose that every N c \mathcal{N}_c N c is a union of residue classes modulo M M M . The encoding above then has regular image: consider the automaton whose states are the pairs ( ρ , j ) (\rho,j) ( ρ , j ) , where ρ ∈ Z / M \rho\in\mathbb{Z}/M ρ ∈ Z / M records the current height modulo M M M and j ∈ { 0 , … , r } j\in\{0,\ldots,r\} j ∈ { 0 , … , r } the last colour read at this height (j = 0 j=0 j = 0 if none), together with an initial state ι \iota ι ; its final states are the ( ρ , j ) (\rho,j) ( ρ , j ) with j ≥ 1 j\geq1 j ≥ 1 , and its transitions are
ι → R ( 1 , 0 ) , ( ρ , j ) → R ( ρ + 1 , 0 ) , ( ρ , j ) → U c ( ρ , c ) if j ≤ c and ρ ∈ N c ( m o d M ) . \iota\xrightarrow{R} (1,0),\qquad(\rho,j)\xrightarrow{R}(\rho+1,0),\qquad(\rho,j)\xrightarrow{U_c}(\rho,c)\ \text{if }j\leq c\text{ and }\rho\in\mathcal{N}_c\pmod M. ι R ( 1 , 0 ) , ( ρ , j ) R ( ρ + 1 , 0 ) , ( ρ , j ) U c ( ρ , c ) if j ≤ c and ρ ∈ N c ( mod M ) .
It is deterministic and accepts the encodings of all nonempty coloured partitions on N \mathcal{N} N , so its generating function is
∏ c = 1 r ∏ h ∈ N c 1 1 − q h − 1 = ∏ t = 1 M 1 ( q t ; q M ) ∞ c t − 1 , \prod_{c=1}^{r}\prod_{h\in\mathcal{N}_c}\frac{1}{1-q^h}-1=\prod_{t=1}^{M}\frac{1}{(q^t;q^M)_\infty^{c_t}}-1, c = 1 ∏ r h ∈ N c ∏ 1 − q h 1 − 1 = t = 1 ∏ M ( q t ; q M ) ∞ c t 1 − 1 ,
where c t c_t c t is the number of colours available at the residue of t t t . Thus every periodic product with nonnegative exponents is, up to the empty partition, the generating function of a regular language; Figure 2 (right) shows a variant, with one colour and a separate final state, for M = 5 M=5 M = 5 . A product with negative exponents, such as Primc’s conjectured ( q 3 ; q 6 ) ∞ / ( q , q 2 , q 4 , q 5 ; q 6 ) ∞ (q^3;q^6)_\infty/(q,q^2,q^4,q^5;q^6)_\infty ( q 3 ; q 6 ) ∞ / ( q , q 2 , q 4 , q 5 ; q 6 ) ∞ for the boundary ( 0 , 1 , 0 ) (0,1,0) ( 0 , 1 , 0 ) , is first cleared of its numerator: multiplying both sides by ( q 3 ; q 6 ) ∞ − 1 (q^3;q^6)_\infty^{-1} ( q 3 ; q 6 ) ∞ − 1 adds, on the side of the difference conditions, one colour available at the parts ≡ 3 ( m o d 6 ) \equiv3 \pmod{6} ≡ 3 ( mod 6 ) and subject to no other condition.
Figure 2. Automata for the Rogers–Ramanujan identities: gap conditions (left) and residue conditions (right). Each automaton reads the complete boundary word, from its unique initial state s s s to its unique final state f f f . Numbers on the right record height modulo 5.
Difference conditions. Difference conditions constrain the multiplicities in a window of consecutive heights. When the multiplicities they involve are bounded, only finitely many multiplicity vectors occur in a window, and the conditions define a regular language:
Lemma 3.1 (Window conditions). Let r r r , M M M , W ≥ 1 W \ge1 W ≥ 1 and K ≥ 0 K \ge0 K ≥ 0 be integers. Let W ⊆ Z / M × { 0 , … , K } r W \mathcal{W} \subseteq\mathbb{Z}/M \times\{0,\ldots,K\}^{rW} W ⊆ Z / M × { 0 , … , K } r W contain every tuple ( ρ , 0 , … , 0 ) (\rho,0,\ldots,0) ( ρ , 0 , … , 0 ) , and fix vectors v 2 − W , … , v 0 ∈ { 0 , … , K } r v_{2-W},\ldots,v_0 \in\{0,\ldots,K\}^r v 2 − W , … , v 0 ∈ { 0 , … , K } r .
The coloured partitions such that
( h m o d M , f h − W + 1 , … , f h ) ∈ W for all h ≥ 1 , (5) (h \bmod M, f_{h-W+1},\ldots,f_h) \in\mathcal{W} \qquad\text{for all } h \ge1, \tag*{(5)} ( h mod M , f h − W + 1 , … , f h ) ∈ W for all h ≥ 1 , (5) where f h = v h f_h=v_h f h = v h for h ≤ 0 h\le0 h ≤ 0 and f h = 0 f_h=0 f h = 0 for h > λ 1 h>\lambda_1 h > λ 1 , are encoded by a regular language, recognized by a deterministic automaton with at most 1 + M ( K + 1 ) r W 1+M(K+1)^{rW} 1 + M ( K + 1 ) r W states.
Proof. After reading a prefix ending at height h h h , the automaton records h h h mod M M M , the vectors f h − W + 1 , … , f h − 1 f_{h-W+1},\ldots,f_{h-1} f h − W + 1 , … , f h − 1 , and the part of f h f_h f h read so far. The first letter R R R leads from the initial state to height 1 1 1 , with recorded vectors v 2 − W , … , v 0 v_{2-W},\ldots,v_0 v 2 − W , … , v 0 . A letter U c U_c U c increases the c c c -th entry of the current vector; it is rejected if a larger colour has already been read at this height, or if the entry would exceed K K K . A letter R R R is rejected unless the completed window lies in W \mathcal{W} W ; otherwise the oldest vector is discarded and the height increases by one. A state reached by a letter U c U_c U c is final if the current window, and the windows at heights h + 1 , … , h + W − 1 h+1,\ldots,h+W-1 h + 1 , … , h + W − 1 obtained from it by appending zero vectors, lie in W \mathcal{W} W ; all later windows are zero. □
Only the information needed for later checks has to be stored, and the automaton can then be minimized. Colours on which W \mathcal{W} W imposes no condition need not be bounded: the automaton reads their letters without recording them. Conditions of different kinds combine, since intersections of regular languages are regular. For instance, Gordon’s conditions [12 ] f h + f h + 1 ≤ k f_h+f_{h+1}\le k f h + f h + 1 ≤ k and f 1 ≤ a f_1\le a f 1 ≤ a , with 0 ≤ a ≤ k 0\le a\le k 0 ≤ a ≤ k , are (5) with r = M = 1 r=M=1 r = M = 1 , W = 2 W=2 W = 2 , K = k K=k K = k and the single initial vector v 0 = k − a v_0=k-a v 0 = k − a : the initial condition becomes a difference condition involving an imaginary multiplicity at height 0. For k = 1 k=1 k = 1 and a = 1 , 0 a=1,0 a = 1 , 0 these are the Rogers–Ramanujan languages L 1 , L 2 \mathcal{L}_1,\mathcal{L}_2 L 1 , L 2 of Figure 2 (left).
Primc’s arrays. Primc’s colours and conditions [18 ] are of this kind. For w = 2 ℓ − 1 w=2\ell-1 w = 2 ℓ − 1 , the multiset N w o d d \mathcal{N}^{\mathrm{odd}}_w N w odd has one colour for each row i = 1 , … , w i=1,\ldots,w i = 1 , … , w of an array, row i i i containing the parts of the parity of i + 1 i+1 i + 1 : the even parts have ℓ \ell ℓ colours, and the odd parts ℓ − 1 \ell-1 ℓ − 1 . Two entries are adjacent if they lie in adjacent rows and their parts differ by one, and a downward path is a sequence of entries, one in each row from top to bottom, any two consecutive ones being adjacent. The level-k k k difference conditions require that the multiplicities along every downward path sum to at most k k k . The initial conditions add imaginary multiplicities: the boundary ( k 0 , … , k ℓ ) (k_0,\ldots,k_\ell) ( k 0 , … , k ℓ ) , with k 0 + ⋯ + k ℓ = k k_0+\cdots+k_\ell=k k 0 + ⋯ + k ℓ = k , is the list ( k 0 , k 1 , 0 , k 2 , 0 , … , 0 , k ℓ − 1 , k ℓ ) (k_0,k_1,0,k_2,0,\ldots,0,k_{\ell-1},k_\ell) ( k 0 , k 1 , 0 , k 2 , 0 , … , 0 , k ℓ − 1 , k ℓ ) , read from the bottom row to the top row, placed at the part 0 0 0 in the odd rows and at the part − 1 -1 − 1 in the even rows.
The parts along a downward path lie in w w w consecutive integers, and every multiplicity lies on some downward path. Hence the admissible coloured partitions are defined by (5), with r = W = w r=W=w r = W = w , M = 2 M=2 M = 2 , K = k K=k K = k , the set W \mathcal{W} W of windows that respect the parities of the rows and in which every downward path has sum at most k k k , and initial vectors v 2 − w , … , v 0 v_{2-w},\ldots,v_0 v 2 − w , … , v 0 that vanish except for v − 1 , v 0 v_{-1},v_0 v − 1 , v 0 , which carry the imaginary multiplicities. By Lemma 3.1, for every width, level and boundary, the admissible coloured partitions form a regular language. Every product formula of [15 Conjecture 4.2] , and every product found numerically in [18 ] , once cleared of negative exponents as above, therefore takes the form ε ( H 1 ) = ε ( H 2 ) \varepsilon(H_1)=\varepsilon(H_2) ε ( H 1 ) = ε ( H 2 ) of Theorem A. For w = 3 w=3 w = 3 , k ≤ 2 k\le2 k ≤ 2 and for w = 5 w=5 w = 5 , k = 1 k=1 k = 1 , the resulting automata reproduce the first terms of all products listed in [18 pp. 53–54] .
For w = 3 w=3 w = 3 , with the multiplicities τ n , μ n , β n \tau_n,\mu_n,\beta_n τ n , μ n , β n used in Theorem B , the array is
τ 0 τ 1 τ 2 τ 3 ⋯ μ 0 μ 1 μ 2 μ 3 μ 4 ⋯ β 0 β 1 β 2 β 3 ⋯ \begin{aligned}
&\tau_0 && \tau_1 && \tau_2 && \tau_3 && \cdots\\
&&\mu_0 && \mu_1 && \mu_2 && \mu_3 && \mu_4 && \cdots\\
&&&\beta_0 && \beta_1 && \beta_2 && \beta_3 && \cdots
\end{aligned} τ 0 μ 0 τ 1 β 0 μ 1 τ 2 β 1 μ 2 τ 3 β 2 μ 3 ⋯ β 3 μ 4 ⋯ ⋯
with columns indexed by the parts − 1 , 0 , 1 , 2 , … -1,0,1,2,\ldots − 1 , 0 , 1 , 2 , … , and imaginary multiplicities τ 0 = k 2 \tau_0=k_2 τ 0 = k 2 , μ 0 = k 1 \mu_0=k_1 μ 0 = k 1 , β 0 = k 0 \beta_0=k_0 β 0 = k 0 . The four downward paths through μ n \mu_n μ n give (3) ; those through μ 0 \mu_0 μ 0 have sum at most k 0 + k 1 + k 2 = 2 k_0+k_1+k_2=2 k 0 + k 1 + k 2 = 2 , so k 1 k_1 k 1 imposes no condition. Once an even height is complete, the automaton only needs to remember the pair ( τ , β ) (\tau,\beta) ( τ , β ) read there; § 4.2 reduces it to four states in this way. Figure 1 shows both automata for the boundary ( 2 , 0 , 0 ) (2,0,0) ( 2 , 0 , 0 ) .
The completion and the verification theorem
Put
A = Q [ q ] ⟨ U , R ⟩ / ( R U − q U R ) , A ^ = lim N A / q N A = Q [ U , R ] [ [ q ] ] , B ^ = Q [ R ] [ [ q ] ] . A=\mathbb{Q}[q]\langle U,R\rangle/(RU-qUR),\qquad\widehat{A}=\lim_{N}A/q^N A=\mathbb{Q}[U,R][[q]],\qquad\widehat{B}=\mathbb{Q}[R][[q]]. A = Q [ q ] ⟨ U , R ⟩ / ( R U − q U R ) , A = N lim A / q N A = Q [ U , R ] [[ q ]] , B = Q [ R ] [[ q ]] .
The notation on the right refers to normal symbols; multiplication in A ^ \widehat{A} A remains the quantum product. Completion is taken before inverting q q q , so the coefficient of each power of q q q is a polynomial in U , R U,R U , R , not an arbitrary series. For s = ∑ k ≥ 0 q k s k ( U , R ) ∈ A ^ s=\sum_{k\geq0}q^k s_k(U,R)\in\widehat{A} s = ∑ k ≥ 0 q k s k ( U , R ) ∈ A , with s k s_k s k in normal form, the marginal and the evaluation of § 1 are
ρ ( s ) = ∑ k ≥ 0 q k s k ( 1 , R ) , ε ( s ) = ∑ k ≥ 0 q k s k ( 1 , 1 ) . \rho(s)=\sum_{k\geq0}q^k s_k(1,R),\qquad\varepsilon(s)=\sum_{k\geq0}q^k s_k(1,1). ρ ( s ) = k ≥ 0 ∑ q k s k ( 1 , R ) , ε ( s ) = k ≥ 0 ∑ q k s k ( 1 , 1 ) .
Both are continuous and linear. A series b = ∑ n ≥ 0 b n R n b=\sum_{n\geq0}b_nR^n b = ∑ n ≥ 0 b n R n with coefficients in Q ( ( q ) ) \mathbb{Q}((q)) Q (( q )) belongs to B ^ [ q − 1 ] \widehat{B}[q^{-1}] B [ q − 1 ] precisely when ν ( b n ) → + ∞ \nu(b_n)\to+\infty ν ( b n ) → + ∞ , and to B ^ \widehat{B} B if moreover all b n b_n b n have nonnegative valuation. Similarly, a normally ordered double series ∑ i , j a i , j U i R j \sum_{i,j}a_{i,j}U^iR^j ∑ i , j a i , j U i R j belongs to A ^ [ q − 1 ] \widehat{A}[q^{-1}] A [ q − 1 ] precisely when ν ( a i , j ) → + ∞ \nu(a_{i,j})\to+\infty ν ( a i , j ) → + ∞ as i + j → ∞ i+j\to\infty i + j → ∞ . In practice we prove bounds of the form
ν ( a i , j ) ≥ δ ( i + j ) − c , ν ( b n ) ≥ δ n − c ( δ > 0 ) . \nu(a_{i,j})\geq\delta(i+j)-c,\qquad\nu(b_n)\geq\delta n-c\qquad(\delta>0). ν ( a i , j ) ≥ δ ( i + j ) − c , ν ( b n ) ≥ δ n − c ( δ > 0 ) .
Theorem 3.2 (Verification). Suppose S ∈ A ^ S\in\widehat{A} S ∈ A and
S = ( 1 − U ) a + b ( 1 − R ) (6) S=(1-U)a+b(1-R) \tag*{(6)} S = ( 1 − U ) a + b ( 1 − R ) (6) with a , b ∈ A ^ [ q − 1 ] a,b\in\widehat{A}[q^{-1}] a , b ∈ A [ q − 1 ] . Then ε ( S ) = 0 \varepsilon(S)=0 ε ( S ) = 0 . Conversely, if ε ( S ) = 0 \varepsilon(S)=0 ε ( S ) = 0 , there is a unique pair satisfying (6) with a ∈ A ^ a\in\widehat{A} a ∈ A and b ∈ B ^ b\in\widehat{B} b ∈ B . If D ( x , R ) D(x,R) D ( x , R ) is the normal symbol of S S S , the normal symbols of this pair are
a ( x , R ) = D ( x , R ) − D ( 1 , R ) 1 − x , b ( R ) = D ( 1 , R ) 1 − R . (7) a(x,R)=\frac{D(x,R)-D(1,R)}{1-x},\qquad b(R)=\frac{D(1,R)}{1-R}. \tag*{(7)} a ( x , R ) = 1 − x D ( x , R ) − D ( 1 , R ) , b ( R ) = 1 − R D ( 1 , R ) . (7)
Proof. Left multiplication by U U U and right multiplication by R R R shift the corresponding normal exponent without changing the coefficient. Consequently ε ( ( 1 − U ) a ) = ε ( b ( 1 − R ) ) = 0 \varepsilon((1-U)a)=\varepsilon(b(1-R))=0 ε (( 1 − U ) a ) = ε ( b ( 1 − R )) = 0 , first for polynomials and then by continuity. This proves the first assertion. Conversely, put C = ρ ( S ) = D ( 1 , R ) C=\rho(S)=D(1,R) C = ρ ( S ) = D ( 1 , R ) . At each q q q -degree, the polynomial D k ( x , R ) − D k ( 1 , R ) D_k(x,R)-D_k(1,R) D k ( x , R ) − D k ( 1 , R ) is divisible by 1 − x 1-x 1 − x , and C k ( 1 ) = 0 C_k(1)=0 C k ( 1 ) = 0 implies that C k ( R ) C_k(R) C k ( R ) is divisible by 1 − R 1-R 1 − R . These divisions construct (7) in the indicated completions, and the symbol of ( 1 − U ) a + b ( 1 − R ) (1-U)a+b(1-R) ( 1 − U ) a + b ( 1 − R ) is ( 1 − x ) a + ( 1 − R ) b = D (1-x)a+(1-R)b=D ( 1 − x ) a + ( 1 − R ) b = D . If b b b is independent of U U U , applying ρ \rho ρ to (6) gives C = b ( 1 − R ) C=b(1-R) C = b ( 1 − R ) , which determines b b b ; then a a a is determined as well, since left multiplication by 1 − U 1-U 1 − U is injective, coefficientwise in q q q .
Certificates
Holonomicity is not needed for the telescoping argument itself, it provides finite descriptions of a a a and b b b . Set
O = K ⟨ R , Y ⟩ / ( Y R − q R Y ) , \mathcal{O}=\mathbb{K}\langle R,\mathcal{Y}\rangle/(\mathcal{Y}R-qR\mathcal{Y}), O = K ⟨ R , Y ⟩ / ( Y R − q R Y ) ,
acting on K [ [ R ] ] \mathbb{K}[[R]] K [[ R ]] by Y g ( R ) = g ( q R ) \mathcal{Y}g(R)=g(qR) Y g ( R ) = g ( q R ) . An operator followed by a polynomial in parentheses denotes the action on that polynomial; thus P ( R ) P(R) P ( R ) is the image of the series R R R under P ∈ O P\in\mathcal{O} P ∈ O . A series in K [ [ R ] ] \mathbb{K}[[R]] K [[ R ]] is holonomic if it is annihilated by a nonzero element of O \mathcal{O} O , equivalently if its q q q -shift span over K ( R ) \mathbb{K}(R) K ( R ) is finite-dimensional. For a normal symbol in x , R x,R x , R we require the analogous finiteness of the joint shift span over K ( x , R ) \mathbb{K}(x,R) K ( x , R ) , under Υ \Upsilon Υ and Ξ \Xi Ξ , where ( Ξ g ) ( x , R ) = g ( q x , R ) (\Xi g)(x,R)=g(qx,R) ( Ξ g ) ( x , R ) = g ( q x , R ) .
A certificate for a partition identity ε ( S ) = 0 \varepsilon(S)=0 ε ( S ) = 0 consists of
(c1) operator kernels defining b b b over O \mathcal{O} O and a a a over its analogue in x , R x,R x , R ;
(c2) initial conditions selecting the intended solutions of those kernels;
(c3) exact algebraic identities verifying (6);
(c4) a proof that ν ( b n ) ≥ n + 1 \nu(b_n) \ge n+1 ν ( b n ) ≥ n + 1 : a direct certificate (Definition 3.9 ) or an independent proof of the identity (Corollary 3.6 ).
No certificate is needed for the convergence of a a a , by Theorem A(2). The selecting data and (c4) are essential: a nonzero annihilating operator does not distinguish between its solutions, and the coefficients of a formal solution need not tend to 0.
Constructing the holonomic kernels
Consider a partition automaton. Let T ( U , R ) T(U,R) T ( U , R ) be the transition matrix of its interior automaton, r 0 r_0 r 0 the row of weights of the interior states reached by the first letter R R R , and e 0 e_0 e 0 the column of weights of the interior states from which the last letter U U U leads to acceptance; thus H = R r 0 ( I − T ( U , R ) ) − 1 e 0 U H=Rr_0(\mathbb{I}-T(U,R))^{-1}e_0U H = R r 0 ( I − T ( U , R ) ) − 1 e 0 U . The monomial identity
R W ( U , R ) U = q U W ( q U , q R ) R , RW(U,R)U=qUW(qU,qR)R, R W ( U , R ) U = q U W ( q U , q R ) R ,
valid for every word W W W , gives a normalized realization
w Q = r , H = w e , Q = I − T ( q U , q R ) , r = q U r 0 , e = e 0 R . (8) wQ=r,\qquad H=we,\qquad Q=\mathbb{I}-T(qU,qR),\qquad r=qUr_0,\qquad e=e_0R. \tag*{(8)} w Q = r , H = w e , Q = I − T ( q U , q R ) , r = q U r 0 , e = e 0 R . (8)
Every interior letter now carries a factor q q q , so Q ≡ I ( m o d q ) Q\equiv\mathbb{I}\pmod q Q ≡ I ( mod q ) , and the selected solution is the geometric series w = r ∑ k ≥ 0 ( I − Q ) k w=r\sum_{k\geq0}(\mathbb{I}-Q)^k w = r ∑ k ≥ 0 ( I − Q ) k , which converges in A ^ 1 × s \widehat{A}^{1\times s} A 1 × s ; this is why the word expansions of §1 converge. For a difference S = H 1 − H 2 S=H_1-H_2 S = H 1 − H 2 we use the block-diagonal realization, with output column ( e 1 , − e 2 ) (e_1,-e_2) ( e 1 , − e 2 ) . More generally, we allow entries of Q , r , e Q,r,e Q , r , e in K ⟨ U , R ⟩ / ( R U − q U R ) K\langle U,R\rangle/(RU-qUR) K ⟨ U , R ⟩ / ( R U − q U R ) that are regular at q = 0 q=0 q = 0 , provided Q ≡ I ( m o d q ) Q\equiv\mathbb{I}\pmod q Q ≡ I ( mod q ) .
Define the anti-homomorphism Φ : A → O \Phi:A\to O Φ : A → O by
Φ ( U i R j ) = R j Υ i . \Phi(U^iR^j)=R^j\Upsilon^i. Φ ( U i R j ) = R j Υ i .
It expresses the right action on marginals: ρ ( u v ) = Φ ( v ) ρ ( u ) \rho(uv)=\Phi(v)\rho(u) ρ ( uv ) = Φ ( v ) ρ ( u ) for normally ordered monomials. Hence
L j i = Φ ( Q i j ) , f = ρ ( r ) T , d i = Φ ( e i ) , L m = f , C = ρ ( S ) = d m , (9) L_{ji}=\Phi(Q_{ij}),\qquad f=\rho(r)^{\mathsf T},\qquad d_i=\Phi(e_i),\qquad Lm=f,\qquad C=\rho(S)=dm, \tag*{(9)} L j i = Φ ( Q ij ) , f = ρ ( r ) T , d i = Φ ( e i ) , L m = f , C = ρ ( S ) = d m , (9)
where m = ρ ( w ) T m=\rho(w)^{\mathsf T} m = ρ ( w ) T . Since L ≡ I ( m o d q ) L\equiv\mathbb{I}\pmod q L ≡ I ( mod q ) , this system has a unique solution in B ^ s \widehat{B}^s B s . Every coefficient of m m m and C C C lies in K K K : writing L = ∑ j R j L j ( Υ ) L=\sum_jR^jL_j(\Upsilon) L = ∑ j R j L j ( Υ ) , the coefficient of R n R^n R n is obtained by solving with the matrix L 0 ( q n ) ≡ I ( m o d q ) L_0(q^n)\equiv\mathbb{I}\pmod q L 0 ( q n ) ≡ I ( mod q ) . Since C C C is linear in f f f , the factor q q q of f = q r 0 T f=qr_0^{\mathsf T} f = q r 0 T may be moved into d d d ; the examples of §4 do so.
Adjoin an auxiliary coordinate u = 1 u=1 u = 1 , which satisfies ( Υ − 1 ) u = 0 (\Upsilon-1)u=0 ( Υ − 1 ) u = 0 . The canonical formal boundary b b b is selected by the finite matrix kernel
G b ( m u b ) = 0 , G b = ( L − f 0 0 Υ − 1 0 − d 0 1 − R ) . (10) G_b
\begin{pmatrix}
m\\
u\\
b
\end{pmatrix}
=0,\qquad
G_b=
\begin{pmatrix}
L & -f & 0\\
0 & \Upsilon-1 & 0\\
-d & 0 & 1-R
\end{pmatrix}. \tag*{(10)} G b m u b = 0 , G b = L 0 − d − f Υ − 1 0 0 0 1 − R . (10)
The last row means ( 1 − R ) b = C (1-R)b=C ( 1 − R ) b = C in K [ [ R ] ] K[[R]] K [[ R ]] , with no convergence assumed. Thus, writing C = ∑ c n R n C=\sum c_nR^n C = ∑ c n R n ,
b = ∑ n ≥ 0 b n R n , b n = ∑ j = 0 n c j . (11) b=\sum_{n\geq0}b_nR^n,\qquad b_n=\sum_{j=0}^{n}c_j. \tag*{(11)} b = n ≥ 0 ∑ b n R n , b n = j = 0 ∑ n c j . (11)
For path series as in Theorem A, the boundary word of a partition contains as many letters R R R as its largest part. Hence c j c_j c j is the difference of the weighted counts of the partitions with largest part j j j , and b n b_n b n that of the partitions with largest part at most n n n .
For a a a , use the full right action on symbols, Φ x ( U i R j ) = x i R j Υ i \Phi_x(U^iR^j)=x^iR^j\Upsilon^i Φ x ( U i R j ) = x i R j Υ i , and compile L x m x = f x L_xm_x=f_x L x m x = f x , D = d x m x D=d_xm_x D = d x m x ; a subscript 1 denotes x = 1 x=1 x = 1 . The normal symbol a a a is selected by
G a ( m x m 1 u a ) = 0 , G a = ( L x 0 − f x 0 0 L 1 − f 1 0 0 0 Υ − 1 0 − d x d 1 0 1 − x ) . (12) G_a
\begin{pmatrix}
m_x\\
m_1\\
u\\
a
\end{pmatrix}
=0,\qquad
G_a=
\begin{pmatrix}
L_x & 0 & -f_x & 0\\
0 & L_1 & -f_1 & 0\\
0 & 0 & \Upsilon-1 & 0\\
-d_x & d_1 & 0 & 1-x
\end{pmatrix}. \tag*{(12)} G a m x m 1 u a = 0 , G a = L x 0 0 − d x 0 L 1 0 d 1 − f x − f 1 Υ − 1 0 0 0 0 1 − x . (12)
The last row is the first formula in (7). These homogeneous kernels have sizes s + 2 s+2 s + 2 and 2 s + 2 2s+2 2 s + 2 ; they are used below to prove holonomicity, and require no scalar elimination. They do not determine their solutions by themselves, since ( Υ − 1 ) u = 0 (\Upsilon-1)u=0 ( Υ − 1 ) u = 0 only forces u u u to be independent of R R R . Substituting u = 1 u=1 u = 1 , that is, deleting the auxiliary row and column, gives square inhomogeneous systems of sizes s + 1 s+1 s + 1 and 2 s + 1 2s+1 2 s + 1 , with right-hand sides formed by f f f , respectively f x f_x f x and f 1 f_1 f 1 . Their solutions in K [ [ R ] ] s + 1 \mathbb{K}[[R]]^{s+1} K [[ R ] ] s + 1 and K [ [ x , R ] ] 2 s + 1 \mathbb{K}[[x,R]]^{2s+1} K [[ x , R ] ] 2 s + 1 are unique: the equations for m m m and m 1 m_1 m 1 are solved coefficientwise as above; those for m x m_x m x are solved degree by degree in x , R x,R x , R , since every letter carries a factor x x x or R R R , so that the part of I − L x \mathrm{I}-L_x I − L x that does not raise the total degree comes from letter-free entries of I − Q \mathrm{I}-Q I − Q and is divisible by q q q ; and 1 − R 1-R 1 − R , 1 − x 1-x 1 − x are injective.
Lemma 3.3 (Finite elimination). Let G G G be a nonsingular n × n n\times n n × n matrix over K ( R ) [ Υ ] \mathbb{K}(R)[\Upsilon] K ( R ) [ Υ ] , of Υ \Upsilon Υ -degree at most δ \delta δ . For each column k k k , a finite linear kernel computation gives a row V V V and a nonzero scalar operator Ω \Omega Ω such that V G = Ω e k T VG=\Omega e_k^{\mathsf T} V G = Ω e k T , with deg V ≤ ( n − 1 ) δ \deg V\le(n-1)\delta deg V ≤ ( n − 1 ) δ .
Proof. Put t = ( n − 1 ) δ t=(n-1)\delta t = ( n − 1 ) δ and seek V j = ∑ i = 0 t v j i ( R ) Υ i V_j=\sum_{i=0}^{t}v_{ji}(R)\Upsilon^i V j = ∑ i = 0 t v j i ( R ) Υ i . Annihilating the other columns gives at most ( n − 1 ) ( t + δ + 1 ) (n-1)(t+\delta+1) ( n − 1 ) ( t + δ + 1 ) linear equations over K ( R ) \mathbb{K}(R) K ( R ) in n ( t + 1 ) n(t+1) n ( t + 1 ) unknowns, one more than the number of equations. A nonzero solution exists, and its remaining entry Ω \Omega Ω cannot vanish because G G G is nonsingular. A common left scalar denominator gives polynomial coefficients in R R R . □
The matrices (10) and (12) are nonsingular over the corresponding Ore skew fields: normalization makes each state matrix invertible in its q q q -adic operator completion, hence nonsingular over the skew field (a left kernel vector, cleared of denominators, would vanish in the completion), and the remaining diagonal blocks Υ − 1 \Upsilon-1 Υ − 1 , 1 − R 1-R 1 − R , 1 − x 1-x 1 − x are nonzero. The lemma therefore proves that all selected components are holonomic. The skew field is used only for elimination; its inverses are never applied to series. For example, Υ − 1 \Upsilon-1 Υ − 1 is invertible there, but annihilates the nonzero series 1 1 1 .
The symbol a a a is jointly holonomic, not merely holonomic with x x x as a parameter. Left multiplication acts on symbols through the homomorphism Υ i R j ↦ x i R j Ξ j \Upsilon^iR^j\mapsto x^iR^j\Xi^j Υ i R j ↦ x i R j Ξ j ; the opposite realization Q v = e Qv=e Q v = e , S = r v S=rv S = r v therefore gives a pure Ξ \Xi Ξ -equation for D D D , and the preceding construction a pure Υ \Upsilon Υ -equation. If their orders are p 1 , p 2 p_1,p_2 p 1 , p 2 , every mixed shift Ξ i Υ j D \Xi^i\Upsilon^jD Ξ i Υ j D reduces to the span of those with i < p 1 i<p_1 i < p 1 , j < p 2 j<p_2 j < p 2 : apply Ξ i − p 1 Υ j \Xi^{i-p_1}\Upsilon^j Ξ i − p 1 Υ j to the first equation to lower i i i without increasing j j j , and Ξ i Υ j − p 2 \Xi^i\Upsilon^{j-p_2} Ξ i Υ j − p 2 to the second to lower j j j without increasing i i i ; shifted leading coefficients remain nonzero in K ( x , R ) \mathbb{K}(x,R) K ( x , R ) . Finite shift rank is preserved by sums and rational multiplication, and C C C is holonomic in R R R and constant in x x x , so (7) proves the assertion for a a a .
The kernels can also be replaced by scalar operators with finite initial data.
Lemma 3.4 (Exceptional indices). Let Ω = ∑ j ≥ j 0 R ℓ j ω j ( Υ ) ∈ O \Omega=\sum_{j\ge j_0}R^{\ell_j}\omega_j(\Upsilon)\in\mathcal{O} Ω = ∑ j ≥ j 0 R ℓ j ω j ( Υ ) ∈ O be nonzero, with ω j 0 ≠ 0 \omega_{j_0}\ne0 ω j 0 = 0 , and put E Ω = { n ≥ 0 : ω j 0 ( q n ) = 0 } E_\Omega=\{n\ge0:\omega_{j_0}(q^n)=0\} E Ω = { n ≥ 0 : ω j 0 ( q n ) = 0 } .
(1) If Ω g = 0 \Omega g=0 Ω g = 0 for g = ∑ g n R n ∈ K [ [ R ] ] g=\sum g_nR^n\in\mathbb{K}[[R]] g = ∑ g n R n ∈ K [[ R ]] , then
ω j 0 ( q n ) g n = − ∑ j > j 0 ω j ( q n + j 0 − j ) g n + j 0 − j , (13) \omega_{j_0}(q^n)g_n=-\sum_{j>j_0}\omega_j(q^{n+j_0-j})g_{n+j_0-j}, \tag*{(13)} ω j 0 ( q n ) g n = − j > j 0 ∑ ω j ( q n + j 0 − j ) g n + j 0 − j , (13) where negative subscripts mean zero. Hence g g g is determined by its coefficients at E Ω E_\Omega E Ω , and g = 0 g=0 g = 0 if and only if these vanish.
(2) Let ω ( Z ) = ∑ i γ i Z i ∈ K [ Z ] \omega(Z)=\sum_i\gamma_iZ^i\in\mathbb{K}[Z] ω ( Z ) = ∑ i γ i Z i ∈ K [ Z ] be nonzero, let i 0 i_0 i 0 be the least index with γ i 0 ≠ 0 \gamma_{i_0}\ne0 γ i 0 = 0 , and let n 0 n_0 n 0 be the maximum of ( ν ( γ i 0 ) − ν ( γ i ) ) / ( i − i 0 ) (\nu(\gamma_{i_0})-\nu(\gamma_i))/(i-i_0) ( ν ( γ i 0 ) − ν ( γ i )) / ( i − i 0 ) over the indices i > i 0 i>i_0 i > i 0 with γ i ≠ 0 \gamma_i\ne0 γ i = 0 (n 0 = − ∞ n_0=-\infty n 0 = − ∞ if there are none). Then ν ( ω ( q n ) ) = ν ( γ i 0 ) + i 0 n \nu(\omega(q^n))=\nu(\gamma_{i_0})+i_0n ν ( ω ( q n )) = ν ( γ i 0 ) + i 0 n for every integer n > n 0 n>n_0 n > n 0 ; in particular ω ( q n ) ≠ 0 \omega(q^n)\ne0 ω ( q n ) = 0 . Applied to ω = ω j 0 \omega=\omega_{j_0} ω = ω j 0 , this shows that E Ω E_\Omega E Ω is contained in { n ≤ n 0 } \{n\le n_0\} { n ≤ n 0 } , hence finite, and E Ω E_\Omega E Ω is computed by evaluating ω j 0 ( q n ) ∈ K \omega_{j_0}(q^n)\in\mathbb{K} ω j 0 ( q n ) ∈ K for these finitely many n n n .
Proof. (1) is the coefficient of R n + j 0 R^{n+j_0} R n + j 0 in Ω g = 0 \Omega g=0 Ω g = 0 . (2) At Z = q n Z=q^n Z = q n , the monomial γ i Z i \gamma_iZ^i γ i Z i becomes γ i q i n \gamma_iq^{in} γ i q in , of valuation ν ( γ i ) + i n \nu(\gamma_i)+in ν ( γ i ) + in . For i > i 0 i>i_0 i > i 0 with γ i ≠ 0 \gamma_i\ne0 γ i = 0 ,
ν ( γ i 0 ) + i 0 n < ν ( γ i ) + i n ⟺ n > ν ( γ i 0 ) − ν ( γ i ) i − i 0 . \nu(\gamma_{i_0})+i_0n<\nu(\gamma_i)+in\Longleftrightarrow n>\frac{\nu(\gamma_{i_0})-\nu(\gamma_i)}{i-i_0}. ν ( γ i 0 ) + i 0 n < ν ( γ i ) + in ⟺ n > i − i 0 ν ( γ i 0 ) − ν ( γ i ) .
Hence, for n > n 0 n>n_0 n > n 0 , the term γ i 0 q i 0 n \gamma_{i_0}q^{i_0n} γ i 0 q i 0 n has strictly smaller valuation than every other term of ω ( q n ) \omega(q^n) ω ( q n ) , so the sum has the valuation of this term.
Values at E Ω E_\Omega E Ω thus select at most one solution, and the normalized construction supplies a consistent one.
Convergence and proof of Theorem A
The kernels produce formal holonomic boundaries whether or not the identity holds: the one-word partition language with series S = R U = q U R S=RU=qUR S = R U = q U R gives b = q R / ( 1 − R ) b=qR/(1-R) b = q R / ( 1 − R ) , so b n = q b_n=q b n = q for every n ≥ 1 n\geq1 n ≥ 1 . For partition languages, a single area estimate gives the precise criterion.
Lemma 3.5 (Area bound). In the setting of Theorem A , write D ( x , R ) = ∑ i , j d i , j x i R j D(x,R)=\sum_{i,j}d_{i,j}x^iR^j D ( x , R ) = ∑ i , j d i , j x i R j and C = D ( 1 , R ) = ∑ j c j R j C=D(1,R)=\sum_j c_jR^j C = D ( 1 , R ) = ∑ j c j R j . Then
ν ( d i , j ) ≥ i + j − 1 , ν ( c j ) ≥ j , ν ( ε ( S ) − b n ) ≥ n + 1 , ν ( a i , j ) ≥ i + j . (14) \nu(d_{i,j})\geq i+j-1,\qquad\nu(c_j)\geq j,\qquad\nu(\varepsilon(S)-b_n)\geq n+1,\qquad\nu(a_{i,j})\geq i+j. \tag*{(14)} ν ( d i , j ) ≥ i + j − 1 , ν ( c j ) ≥ j , ν ( ε ( S ) − b n ) ≥ n + 1 , ν ( a i , j ) ≥ i + j . (14) Proof. A word in R { R , U } ∗ U R\{R,U\}^{*}U R { R , U } ∗ U with i i i letters U U U and j j j letters R R R has area at least i + j − 1 i+j-1 i + j − 1 : count the pairs involving its first R R R or its last U U U . The weights have nonnegative valuation, so ν ( d i , j ) ≥ i + j − 1 \nu(d_{i,j})\geq i+j-1 ν ( d i , j ) ≥ i + j − 1 , and ν ( c j ) ≥ j \nu(c_j)\geq j ν ( c j ) ≥ j because i ≥ 1 i\geq1 i ≥ 1 . Next, ε ( S ) − b n = ∑ j > n c j \varepsilon(S)-b_n=\sum_{j>n}c_j ε ( S ) − b n = ∑ j > n c j . Finally, ( x t − 1 ) / ( 1 − x ) = − ( 1 + x + ⋯ + x t − 1 ) (x^t-1)/(1-x)=-(1+x+\cdots+x^{t-1}) ( x t − 1 ) / ( 1 − x ) = − ( 1 + x + ⋯ + x t − 1 ) gives a i , j = − ∑ t > i d t , j a_{i,j}=-\sum_{t>i}d_{t,j} a i , j = − ∑ t > i d t , j , of valuation at least min t > i ( t + j − 1 ) = i + j \min_{t>i}(t+j-1)=i+j min t > i ( t + j − 1 ) = i + j .
In particular, ( b n ) (b_n) ( b n ) always converges to ε ( S ) \varepsilon(S) ε ( S ) , and
ε ( S ) = 0 ⟺ b n ⟶ 0 ⟺ ν ( b n ) ≥ n + 1 for all n ⟺ b = q Υ ϕ , ϕ ∈ Q [ [ q , R ] ] . (15) \begin{aligned}
\varepsilon(S)=0 &\Longleftrightarrow b_n\longrightarrow0\\
&\Longleftrightarrow\nu(b_n)\geq n+1\text{ for all }n\\
&\Longleftrightarrow b=q\Upsilon\phi,\quad\phi\in\mathbb{Q}[[q,R]]. \tag*{(15)}
\end{aligned} ε ( S ) = 0 ⟺ b n ⟶ 0 ⟺ ν ( b n ) ≥ n + 1 for all n ⟺ b = q Υ ϕ , ϕ ∈ Q [[ q , R ]] . (15)
Indeed, the last condition says exactly that ν ( q − n − 1 b n ) ≥ 0 \nu(q^{-n-1}b_n)\geq0 ν ( q − n − 1 b n ) ≥ 0 . Even arbitrarily slow decay along an unbounded subsequence suffices: it forces ε ( S ) = 0 \varepsilon(S)=0 ε ( S ) = 0 and hence the uniform bound n + 1 n+1 n + 1 . Conversely, if ε ( S ) ≠ 0 \varepsilon(S)\ne0 ε ( S ) = 0 , then ν ( b n ) = ν ( ε ( S ) ) < n + 1 \nu(b_n)=\nu(\varepsilon(S))<n+1 ν ( b n ) = ν ( ε ( S )) < n + 1 for all n ≥ ν ( ε ( S ) ) n\geq\nu(\varepsilon(S)) n ≥ ν ( ε ( S )) , a finite disproof.
Corollary 3.6 (Indirect certificates). If ε ( H 1 ) = ε ( H 2 ) \varepsilon(H_1)=\varepsilon(H_2) ε ( H 1 ) = ε ( H 2 ) is proved by any means, then the canonical pair ( a , b ) (a,b) ( a , b ) of Theorem A, selected by the kernels of §3.5, satisfies S = ( 1 − U ) a + b ( 1 − R ) S=(1-U)a+b(1-R) S = ( 1 − U ) a + b ( 1 − R ) with ν ( a i , j ) ≥ i + j \nu(a_{i,j})\geq i+j ν ( a i , j ) ≥ i + j and ν ( b n ) ≥ n + 1 \nu(b_n)\geq n+1 ν ( b n ) ≥ n + 1 .
Proof. Combine Lemma 3.5, (15) and Theorem 3.2.
Proof of Theorem A. For (1), set u = 1 u=1 u = 1 in (12) and (10); §3.5 proves uniqueness, that the coefficients b n b_n b n lie in K \mathbb{K} K , and holonomicity in the stated variables. Lemma 3.5 proves (2), and (15) the equivalence of (3a)–(3d). A kernel for ϕ = q − 1 b ( R / q ) \phi=q^{-1}b(R/q) ϕ = q − 1 b ( R / q ) follows by replacing R R R by R / q R/q R / q in a kernel for b b b . Condition (3c) implies (3e), and Theorem 3.2 gives (3e)⇒ \Rightarrow ⇒ (3a) together with the uniqueness statement.
Certifying convergence
Items (c1)–(c3) of a certificate are computed from the automata (§3.5), and by (15) item (c4) is equivalent to the identity itself. We supply (c4) in one of two ways. A direct certificate (Definition 3.9) proves that the preimage ϕ = q − 1 b ( R / q ) \phi=q^{-1}b(R/q) ϕ = q − 1 b ( R / q ) is q q q -integral by finitely many exact checks on holonomic series. An indirect certificate is an independent proof of the identity, which Corollary 3.6 converts into convergence.
Let F ⊂ K ( U ) \mathcal{F}\subset\mathbb{K}(U) F ⊂ K ( U ) be the ring of quotients of a polynomial in Q [ q , U ] \mathbb{Q}[q,U] Q [ q , U ] by a polynomial in 1 + q Q [ q , U ] 1+q\mathbb{Q}[q,U] 1 + q Q [ q , U ] . For g ∈ F g\in\mathcal{F} g ∈ F , the operator g ( Υ ) g(\Upsilon) g ( Υ ) multiplies the coefficient of R n R^n R n by g ( q n ) ∈ Q [ [ q ] ] g(q^n)\in\mathbb{Q}[[q]] g ( q n ) ∈ Q [[ q ]] . Hence the integral operators ∑ j R j g j ( Υ ) \sum_jR^jg_j(\Upsilon) ∑ j R j g j ( Υ ) , with g j ∈ F g_j\in\mathcal{F} g j ∈ F , preserve Q [ [ q , R ] ] \mathbb{Q}[[q,R]] Q [[ q , R ]] ; they form a ring, since g ( Υ ) R j = R j g ( q j Υ ) g(\Upsilon)R^j=R^jg(q^j\Upsilon) g ( Υ ) R j = R j g ( q j Υ ) . Call a matrix L L L of integral operators an integral kernel if the R R R -free part of Υ − L \Upsilon-L Υ − L has entries in q F q\mathcal{F} q F . Then Υ − L \Upsilon-L Υ − L raises the total degree in q q q and R R R , so for f ∈ Q [ [ q , R ] ] s f\in\mathbb{Q}[[q,R]]^s f ∈ Q [[ q , R ] ] s the series m = ∑ k ≥ 0 ( Υ − L ) k f m=\sum_{k\geq0}(\Upsilon-L)^kf m = ∑ k ≥ 0 ( Υ − L ) k f converges in Q [ [ q , R ] ] s \mathbb{Q}[[q,R]]^s Q [[ q , R ] ] s and solves L m = f Lm=f L m = f . It is the only solution in Q ( ( q ) ) [ [ R ] ] s \mathbb{Q}((q))[[R]]^s Q (( q )) [[ R ] ] s : the coefficient of R n R^n R n in L m = f Lm=f L m = f determines that of m m m from lower ones, through a matrix congruent to Υ \Upsilon Υ modulo q q q . If, in addition, f ∈ K [ [ R ] ] s f\in\mathbb{K[[R]]}^s f ∈ K [[ R ]] s , the same induction gives m ∈ K [ [ R ] ] s m\in\mathbb{K[[R]]}^s m ∈ K [[ R ]] s .
Integral kernels do not always suffice: in (2), the boundary of Euler’s identity combines two series whose coefficients are q q q -adic units, and its convergence rests on a cancellation between them. The following series account for such cancellations. For ℓ ≥ 1 \ell\ge1 ℓ ≥ 1 and α ∈ F \alpha\in\mathcal{F} α ∈ F , the integral kernel 1 − R ℓ α ( Υ ) 1-R^{\ell}\alpha(\Upsilon) 1 − R ℓ α ( Υ ) selects
u α = ∑ k ≥ 0 ( ∏ j = 0 k − 1 α ( q ℓ j ) ) R ℓ k , ( 1 − R ℓ α ( Υ ) ) u α = 1. u_{\alpha}=\sum_{k\ge0}\left(\prod_{j=0}^{k-1}\alpha\left(q^{\ell^{j}}\right)\right)R^{\ell^{k}},\qquad\left(1-R^{\ell}\alpha(\Upsilon)\right)u_{\alpha}=1. u α = k ≥ 0 ∑ ( j = 0 ∏ k − 1 α ( q ℓ j ) ) R ℓ k , ( 1 − R ℓ α ( Υ ) ) u α = 1.
Proposition 3.7 (First-order series). Fix ℓ ≥ 1 \ell\ge1 ℓ ≥ 1 . Let α , β \alpha,\beta α , β be quotients of polynomials in 1 + q Q [ q , U ] 1+q\mathbb{Q}[q,U] 1 + q Q [ q , U ] , and suppose that θ = β / α \theta=\beta/\alpha θ = β / α satisfies θ − 1 ∈ q c U r F \theta-1\in q^{c}U^{r}\mathcal{F} θ − 1 ∈ q c U r F for some integers c ≥ 0 c\ge0 c ≥ 0 and r ≥ 1 r\ge1 r ≥ 1 .
(1) The following are equivalent:
(a) ∏ j ≥ 0 θ ( q ℓ j ) = 1 \prod_{j\ge0}\theta\left(q^{\ell^{j}}\right)=1 ∏ j ≥ 0 θ ( q ℓ j ) = 1 ;
(b) u α − u β = q c Υ r ψ u_{\alpha}-u_{\beta}=q^{c}\Upsilon^{r}\psi u α − u β = q c Υ r ψ for some ψ ∈ Q [ [ q , R ] ] \psi\in\mathbb{Q}[[q,R]] ψ ∈ Q [[ q , R ]] ;
(c) the coefficients of u α − u β u_{\alpha}-u_{\beta} u α − u β tend to 0 0 0 .
(2) Condition (a) holds if there are an integer λ ≥ 2 \lambda\ge2 λ ≥ 2 and a quotient χ \chi χ of polynomials in 1 + q Q [ q , U ] 1+q\mathbb{Q}[q,U] 1 + q Q [ q , U ] , with χ ( 0 ) = 1 \chi(0)=1 χ ( 0 ) = 1 , such that
θ ( U ) χ ( U ) = ∏ i = 0 λ − 1 χ ( q ℓ i U λ ) . (16) \theta(U)\chi(U)=\prod_{i=0}^{\lambda-1}\chi\left(q^{\ell^{i}}U^{\lambda}\right). \tag*{(16)} θ ( U ) χ ( U ) = i = 0 ∏ λ − 1 χ ( q ℓ i U λ ) . (16)
Proof. The coefficient of R ℓ k R^{\ell k} R ℓ k in u α u_{\alpha} u α is u k = ∏ j < k α ( q ℓ j ) ∈ 1 + q Q [ [ q ] ] u_k=\prod_{j<k}\alpha\left(q^{\ell^{j}}\right)\in1+q\mathbb{Q}[[q]] u k = ∏ j < k α ( q ℓ j ) ∈ 1 + q Q [[ q ]] , that in u β u_{\beta} u β is u k π k u_k\pi_k u k π k with π k = ∏ j < k θ ( q ℓ j ) \pi_k=\prod_{j<k}\theta\left(q^{\ell^{j}}\right) π k = ∏ j < k θ ( q ℓ j ) , and all other coefficients vanish. Since ν ( θ ( q ℓ j ) − 1 ) ≥ c + r ℓ j \nu\left(\theta\left(q^{\ell^{j}}\right)-1\right)\ge c+r\ell j ν ( θ ( q ℓ j ) − 1 ) ≥ c + r ℓ j , the products π k \pi_k π k converge to π ∞ = ∏ j ≥ 0 θ ( q ℓ j ) \pi_{\infty}=\prod_{j\ge0}\theta\left(q^{\ell^{j}}\right) π ∞ = ∏ j ≥ 0 θ ( q ℓ j ) .
(1) If π ∞ = 1 \pi_{\infty}=1 π ∞ = 1 , then π k = ∏ j ≥ k θ ( q ℓ j ) − 1 ≡ 1 ( m o d q c + r ℓ k ) \pi_k=\prod_{j\ge k}\theta\left(q^{\ell^{j}}\right)^{-1}\equiv1\pmod{q^{c+r\ell k}} π k = ∏ j ≥ k θ ( q ℓ j ) − 1 ≡ 1 ( mod q c + r ℓ k ) , so the coefficient u k ( 1 − π k ) u_k(1-\pi_k) u k ( 1 − π k ) of R ℓ k R^{\ell k} R ℓ k in u α − u β u_{\alpha}-u_{\beta} u α − u β has valuation at least c + r ℓ k c+r\ell k c + r ℓ k . This is (b), because [ R n ] q c Υ r ψ = q c + r n [ R n ] ψ [R^n]q^c\Upsilon^r\psi=q^{c+rn}[R^n]\psi [ R n ] q c Υ r ψ = q c + r n [ R n ] ψ . Next, (b) implies (c) because r ≥ 1 r\ge1 r ≥ 1 . Finally, (c) means that u k ( 1 − π k ) → 0 u_k(1-\pi_k)\to0 u k ( 1 − π k ) → 0 ; since the u k u_k u k are units, π k → 1 \pi_k\to1 π k → 1 , so π ∞ = 1 \pi_{\infty}=1 π ∞ = 1 .
(2) Since χ ( 0 ) = 1 \chi(0)=1 χ ( 0 ) = 1 , we have χ − 1 ∈ U F \chi-1\in U\mathcal{F} χ − 1 ∈ U F , so each χ ( q ℓ t ) \chi\left(q^{\ell^t}\right) χ ( q ℓ t ) is a unit congruent to 1 1 1 modulo q ℓ t q^{\ell^t} q ℓ t . Evaluate (16) at U = q ℓ j U=q^{\ell^j} U = q ℓ j and multiply over 0 ≤ j < k 0\le j<k 0 ≤ j < k . As ( i , j ) (i,j) ( i , j ) ranges over { 0 , … , λ − 1 } × { 0 , … , k − 1 } \{0,\ldots,\lambda-1\}\times\{0,\ldots,k-1\} { 0 , … , λ − 1 } × { 0 , … , k − 1 } , the integer i + λ j i+\lambda j i + λj ranges over { 0 , … , λ k − 1 } \{0,\ldots,\lambda k-1\} { 0 , … , λk − 1 } , so
π k = ∏ t = 0 λ k − 1 χ ( q ℓ t ) ∏ t = 0 k − 1 χ ( q ℓ t ) = ∏ t = k λ k − 1 χ ( q ℓ t ) ≡ 1 ( m o d q ℓ k ) . \pi_k=\frac{\prod_{t=0}^{\lambda k-1}\chi\left(q^{\ell^t}\right)}{\prod_{t=0}^{k-1}\chi\left(q^{\ell^t}\right)}=\prod_{t=k}^{\lambda k-1}\chi\left(q^{\ell^t}\right)\equiv1\pmod{q^{\ell^k}}. π k = ∏ t = 0 k − 1 χ ( q ℓ t ) ∏ t = 0 λk − 1 χ ( q ℓ t ) = t = k ∏ λk − 1 χ ( q ℓ t ) ≡ 1 ( mod q ℓ k ) .
Hence π ∞ = 1 \pi_{\infty}=1 π ∞ = 1 . □
The coefficients of u α u_{\alpha} u α and u β u_{\beta} u β are q q q -adic units, so neither series converges at R = 1 R=1 R = 1 . Part (1) reduces the convergence of their difference to the product identity (a), in which θ ( q ℓ j ) \theta\left(q^{\ell^j}\right) θ ( q ℓ j ) is the quotient of the ratios β ( q ℓ j ) \beta\left(q^{\ell^j}\right) β ( q ℓ j ) and α ( q ℓ j ) \alpha\left(q^{\ell^j}\right) α ( q ℓ j ) of consecutive coefficients of u β u_{\beta} u β and u α u_{\alpha} u α ; in Corollary 3.8 below, (a) is the partition identity itself. Part (2) certifies (a) by a single identity in K ( U ) \mathop{\rm K}(U) K ( U ) , which we call a dilation witness . First-order series arise as bounded counts of languages defined by independent conditions on the multiplicities f h f_h f h of (4).
Corollary 3.8 (Periodic product languages). Fix ℓ ≥ 1 \ell\ge1 ℓ ≥ 1 and, for each i ∈ Z / ℓ i\in\mathbb{Z}/\ell i ∈ Z / ℓ , eventually periodic sets M i , 1 , M i , 2 ⊆ Z ≥ 0 M_{i,1},M_{i,2}\subseteq\mathbb{Z}_{\ge0} M i , 1 , M i , 2 ⊆ Z ≥ 0 containing 0 0 0 , with generating functions F i ( x ) = ∑ m ∈ M i , 1 x m F_i(x)=\sum_{m\in M_{i,1}}x^m F i ( x ) = ∑ m ∈ M i , 1 x m and G i ( x ) = ∑ m ∈ M i , 2 x m G_i(x)=\sum_{m\in M_{i,2}}x^m G i ( x ) = ∑ m ∈ M i , 2 x m . Let H 1 H_1 H 1 be the path series of the language of partitions in which, for every h ≥ 1 h\ge1 h ≥ 1 , the multiplicity f h f_h f h of h h h lies in M h m o d ℓ , 1 M_{h\bmod\ell,1} M h mod ℓ , 1 , and define H 2 H_2 H 2 likewise with the sets M i , 2 M_{i,2} M i , 2 . Put
α ( U ) = ∏ i = 1 ℓ F i ( q i U ) , β ( U ) = ∏ i = 1 ℓ G i ( q i U ) . \alpha(U)=\prod_{i=1}^{\ell}F_i\left(q^iU\right),\qquad\beta(U)=\prod_{i=1}^{\ell}G_i\left(q^iU\right). α ( U ) = i = 1 ∏ ℓ F i ( q i U ) , β ( U ) = i = 1 ∏ ℓ G i ( q i U ) . If θ = β / α \theta=\beta/\alpha θ = β / α has a dilation witness, then ϵ ( H 1 ) = ϵ ( H 2 ) \epsilon(H_1)=\epsilon(H_2) ϵ ( H 1 ) = ϵ ( H 2 ) , and the canonical boundaries of H 1 − H 2 H_1-H_2 H 1 − H 2 satisfy ν ( a i , j ) ≥ i + j \nu(a_{i,j})\ge i+j ν ( a i , j ) ≥ i + j and ν ( b n ) ≥ n + 1 \nu(b_n)\ge n+1 ν ( b n ) ≥ n + 1 .
Proof. The languages are regular by (4) : an automaton needs to remember only the height modulo ℓ \ell ℓ and the number of letters U U U read at the current height, up to the sum of a preperiod and a period of the sets M i , 1 , M i , 2 M_{i,1}, M_{i,2} M i , 1 , M i , 2 . Each F i ( x ) F_i(x) F i ( x ) is a quotient of polynomials in 1 + x Q [ x ] 1+x\mathbb{Q}[x] 1 + x Q [ x ] , so α \alpha α , β \beta β and θ \theta θ are quotients of polynomials in 1 + q Q [ q , U ] 1+q\mathbb{Q}[q,U] 1 + q Q [ q , U ] ; moreover θ ( 0 ) = 1 \theta(0)=1 θ ( 0 ) = 1 , so θ − 1 ∈ U F \theta-1\in UF θ − 1 ∈ U F .
By Proposition 3.7 , with c = 0 c=0 c = 0 and r = 1 r=1 r = 1 , the coefficients of u α − u β u_\alpha-u_\beta u α − u β tend to 0 0 0 . Grouping the part sizes ℓ j + 1 , … , ℓ j + ℓ \ell j+1,\ldots,\ell j+\ell ℓ j + 1 , … , ℓ j + ℓ for 0 ≤ j < k 0\leq j<k 0 ≤ j < k , the partitions of the first language with largest part at most ℓ k \ell k ℓ k have generating function ∏ h = 1 ℓ k F h ( q ) = ∏ j < k α ( q ℓ j ) \prod_{h=1}^{\ell k}F_h(q)=\prod_{j<k}\alpha(q^{\ell j}) ∏ h = 1 ℓ k F h ( q ) = ∏ j < k α ( q ℓ j ) , and similarly for the second. By the remark after (11) , b ℓ k b_{\ell k} b ℓ k is therefore the coefficient of R ℓ k R^{\ell k} R ℓ k in u α − u β u_\alpha-u_\beta u α − u β , so b ℓ k → 0 b_{\ell k}\to0 b ℓ k → 0 .
Since decay along a subsequence suffices in (15) , ε ( H 1 ) = ε ( H 2 ) \varepsilon(H_1)=\varepsilon(H_2) ε ( H 1 ) = ε ( H 2 ) and ν ( b n ) ≥ n + 1 \nu(b_n)\geq n+1 ν ( b n ) ≥ n + 1 ; finally, ν ( a i , j ) ≥ i + j \nu(a_{i,j})\geq i+j ν ( a i , j ) ≥ i + j by Theorem A(2). □
Glaisher’s identity (Proposition 4.1) is the basic example.
Definition 3.9. A direct certificate for the convergence of a canonical boundary b b b consists of
finitely many integral inputs : components of solutions of L m = f Lm=f L m = f , with L L L an integral kernel and f ∈ Q [ q ] [ R ] s f\in\mathbb{Q}[q][R]^s f ∈ Q [ q ] [ R ] s , and series ψ \psi ψ as in Proposition 3.7 (1b), each with a dilation witness;
a series ϕ ′ \phi' ϕ ′ obtained from the integral inputs and from Q [ q ] [ R ] \mathbb{Q}[q][R] Q [ q ] [ R ] by sums and integral operators;
operator identities and selecting values proving b = q Υ ϕ ′ b=q^\Upsilon\phi' b = q Υ ϕ ′ .
By construction ϕ ′ ∈ Q [ [ q , R ] ] \phi'\in\mathbb{Q}[[q,R]] ϕ ′ ∈ Q [[ q , R ]] , so a direct certificate proves Theorem A(3d), and hence all of Theorem A(3). An identity between selected solutions is checked either by explicit operator multipliers or by a common annihilator and the finitely many exceptional coefficients of Lemma 3.4.
Euler’s identity (Corollary 4.2) and the five-state identity of § 4.7 have direct certificates. The other identities of §4 have indirect ones: Glaisher’s identity by Corollary 3.8, and those of §§§ 4.2 –§ 4.6 by finite telescoping identities, a quadratic invariant (for Primc) and a recurrence with an explicit solution (for Berkovich–Uncu), together with explicit valuation estimates; the proof for Primc also uses the Göllnitz–Gordon identities and Jacobi’s triple product. Appendix A describes further finite tests, not needed for these examples, which bound the valuations of b n b_n b n instead of exhibiting an integral preimage.
Finite searches
The searches were developed to certify the vanishing differences between the 222 354 minimal deterministic interior automata with at most four states, enumerated with one initial state and partial transitions: grouped by their series to order q 50 q^{50} q 50 , they form 19 556 non-singleton families. For bounded operator degrees, the binding identities of a direct certificate, and the telescoping certificates and quadratic invariants below, are solutions of finite linear systems; first-order data and dilation witnesses of bounded degree can be enumerated and checked by (16) ; and rational shift quotients for a proposed bounded sum are checked by clearing denominators, as in §4.5. Our JULIA implementation (§4.7) uses such bounded searches, and Appendix A describes further ones.
Telescoping. Let ( f n ) n ≥ 0 (f_n)_{n\geq0} ( f n ) n ≥ 0 satisfy f n + 1 = r ( q n ) f n f_{n+1}=r(q^n)f_n f n + 1 = r ( q n ) f n for a rational function r r r , in which further parameters, such as v v v in §4.2, may occur. A telescoping certificate for a residual ω \omega ω is a rational function Ψ \Psi Ψ such that
r ( ζ ) Ψ ( q ζ ) − Ψ ( ζ ) = ω ( ζ ) . (17) r(\zeta)\Psi(q\zeta)-\Psi(\zeta)=\omega(\zeta). \tag*{(17)} r ( ζ ) Ψ ( q ζ ) − Ψ ( ζ ) = ω ( ζ ) . (17)
This is an identity between rational functions of an indeterminate ζ \zeta ζ , checked by clearing denominators; for a fixed denominator and a degree bound, Ψ \Psi Ψ is found by solving a linear system. Evaluating at ζ = q n \zeta=q^n ζ = q n , multiplying by f n f_n f n and summing over 0 ≤ n ≤ N 0\leq n\leq N 0 ≤ n ≤ N gives
∑ n = 0 N f n ω ( q n ) = f N + 1 Ψ ( q N + 1 ) − f 0 Ψ ( 1 ) , \sum_{n=0}^{N} f_n\omega(q^n)=f_{N+1}\Psi(q^{N+1})-f_0\Psi(1), n = 0 ∑ N f n ω ( q n ) = f N + 1 Ψ ( q N + 1 ) − f 0 Ψ ( 1 ) ,
provided that no denominator vanishes at the points q n q^n q n ; by Lemma 3.4 (2), this takes finitely many evaluations. For an infinite sum, the upper flux f N + 1 Ψ ( q N + 1 ) f_{N+1}\Psi(q^{N+1}) f N + 1 Ψ ( q N + 1 ) must also tend to 0. This is checked symbolically, treating N N N and q N q^N q N as independent symbols. By Lemma 3.4 (2), applied to numerators and denominators, ν ( r ( q n ) ) \nu(r(q^n)) ν ( r ( q n )) and ν ( Ψ ( q n ) ) \nu(\Psi(q^n)) ν ( Ψ ( q n )) are affine functions of n n n for large n n n , read off from the terms of lowest degree in ζ \zeta ζ . Hence
ν ( f N + 1 Ψ ( q N + 1 ) ) = ν ( f 0 ) + ∑ n = 0 N ν ( r ( q n ) ) + ν ( Ψ ( q N + 1 ) ) \nu\left(f_{N+1}\Psi(q^{N+1})\right)=\nu(f_0)+\sum_{n=0}^{N}\nu(r(q^n))+\nu(\Psi(q^{N+1})) ν ( f N + 1 Ψ ( q N + 1 ) ) = ν ( f 0 ) + n = 0 ∑ N ν ( r ( q n )) + ν ( Ψ ( q N + 1 ))
is a quadratic polynomial in N N N for large N N N , and the flux tends to 0 when this polynomial tends to + ∞ +\infty + ∞ , typically because it is at least c ′ N 2 − c c'N^2-c c ′ N 2 − c with c ′ > 0 c'>0 c ′ > 0 . The same count applies to other valuations, such as the v v v -adic one in § 4.2 . All checks are thus manipulations of rational functions in ζ \zeta ζ and of polynomials in N N N , which any computer algebra system performs. Several summation variables are handled by sums of such differences.
Quadratic invariants. For a system F ( z ) = M ( z ) F ( q h z ) F(z)=M(z)F(q^hz) F ( z ) = M ( z ) F ( q h z ) with F F F formal at z = 0 z=0 z = 0 , a quadratic invariant is a polynomial symmetric matrix J ( z ) J(z) J ( z ) such that
M ( z ) T J ( z ) M ( z ) = λ ( z ) J ( q h z ) M(z)^{\mathsf{T}}J(z)M(z)=\lambda(z)J(q^hz) M ( z ) T J ( z ) M ( z ) = λ ( z ) J ( q h z )
for some polynomial λ ( z ) \lambda(z) λ ( z ) . For fixed λ \lambda λ , this equation is linear in the coefficients of J J J . If λ ( 0 ) = 0 \lambda(0)=0 λ ( 0 ) = 0 , it forces F T J F = 0 F^{\mathsf{T}}JF=0 F T J F = 0 , because this series g ( z ) g(z) g ( z ) satisfies g ( z ) = λ ( z ) g ( q h z ) g(z)=\lambda(z)g(q^hz) g ( z ) = λ ( z ) g ( q h z ) and so has no term of lowest order. Such a relation can reveal a factorization. In § 4.2 , it is det H = 0 \det H=0 det H = 0 for a 2 × 2 2\times2 2 × 2 matrix H H H formed from the four states, on which the transfer acts by H ( v ) = Γ + ( v ) H ( q v ) Γ − T H(v)=\Gamma_+(v)H(qv)\Gamma_-^{\mathsf{T}} H ( v ) = Γ + ( v ) H ( q v ) Γ − T , that is, through the tensor product of two two-state transfers. The solution is the rank-one matrix u ( v ) u ( − v ) T u(v)u(-v)^{\mathsf{T}} u ( v ) u ( − v ) T , and the components of u u u are found by telescoping.
Partition identities
Glaisher's identity
For an integer λ ≥ 2 \lambda\ge2 λ ≥ 2 , let H 1 H_1 H 1 and H 2 H_2 H 2 be the path series of the languages of partitions with no part divisible by λ \lambda λ , respectively in which no part occurs λ \lambda λ or more times. They are as in Corollary 3.8 , with ℓ = λ \ell=\lambda ℓ = λ , M 0 , 1 = { 0 } M_{0,1}=\{0\} M 0 , 1 = { 0 } , M i , 1 = Z ≥ 0 M_{i,1}=\mathbb{Z}_{\ge0} M i , 1 = Z ≥ 0 for i ≢ 0 i\not\equiv0 i ≡ 0 , and M i , 2 = { 0 , … , λ − 1 } M_{i,2}=\{0,\ldots,\lambda-1\} M i , 2 = { 0 , … , λ − 1 } ; thus
F i ( x ) = 1 1 − x ( i ≢ 0 ) , F 0 ( x ) = 1 , G i ( x ) = 1 + x + ⋯ + x λ − 1 . F_i(x)=\frac{1}{1-x}\quad(i\not\equiv0),\qquad F_0(x)=1,\qquad G_i(x)=1+x+\cdots+x^{\lambda-1}. F i ( x ) = 1 − x 1 ( i ≡ 0 ) , F 0 ( x ) = 1 , G i ( x ) = 1 + x + ⋯ + x λ − 1 .
For λ = 2 \lambda=2 λ = 2 , they are the series H o H_o H o and H d H_d H d of §1.
Proposition 4.1 (following Glaisher [10 ] ). Let λ ≥ 2 \lambda\ge2 λ ≥ 2 be an integer. For every n n n , the partitions of n n n with no part divisible by λ \lambda λ are as many as those in which no part occurs λ \lambda λ or more times; for λ = 2 \lambda=2 λ = 2 , this is Euler’s identity. The canonical boundaries of H 1 − H 2 H_1-H_2 H 1 − H 2 satisfy ν ( a i , j ) ≥ i + j \nu(a_{i,j})\ge i+j ν ( a i , j ) ≥ i + j and ν ( b n ) ≥ n + 1 \nu(b_n)\ge n+1 ν ( b n ) ≥ n + 1 .
Proof. In the notation of Corollary 3.8 , since G i ( x ) = ( 1 − x λ ) / ( 1 − x ) G_i(x)=(1-x^\lambda)/(1-x) G i ( x ) = ( 1 − x λ ) / ( 1 − x ) ,
α ( U ) = ∏ i = 1 λ − 1 1 1 − q i U , β ( U ) = ∏ i = 1 λ 1 − q λ i U λ 1 − q i U . \alpha(U)=\prod_{i=1}^{\lambda-1}\frac{1}{1-q^iU},\qquad\beta(U)=\prod_{i=1}^{\lambda}\frac{1-q^{\lambda i}U^\lambda}{1-q^{i}U}. α ( U ) = i = 1 ∏ λ − 1 1 − q i U 1 , β ( U ) = i = 1 ∏ λ 1 − q i U 1 − q λi U λ .
The factors 1 − q i U 1-q^iU 1 − q i U with i < λ i<\lambda i < λ cancel in θ = β / α \theta=\beta/\alpha θ = β / α , so
θ ( U ) = ∏ i = 1 λ − 1 ( 1 − q λ i U λ ) 1 − q λ U . \theta(U)=\frac{\prod_{i=1}^{\lambda-1}(1-q^{\lambda i}U^\lambda)}{1-q^{\lambda}U}. θ ( U ) = 1 − q λ U ∏ i = 1 λ − 1 ( 1 − q λi U λ ) .
The polynomial χ ( U ) = 1 − q λ U \chi(U)=1-q^\lambda U χ ( U ) = 1 − q λ U is a dilation witness: both sides of (16) equal the numerator of θ \theta θ . Corollary 3.8 concludes.
The dilation witness encodes the classical argument 1 + x + ⋯ + x λ − 1 = ( 1 − x λ ) / ( 1 − x ) 1+x+\cdots+x^{\lambda-1}=(1-x^\lambda)/(1-x) 1 + x + ⋯ + x λ − 1 = ( 1 − x λ ) / ( 1 − x ) . The proof identifies b k b_k b k through its combinatorial meaning. For λ = 2 \lambda=2 λ = 2 , the same first-order data give the direct certificate of the teaser (Definition 3.9 ), whose checks use only the kernels.
Corollary 4.2 (Euler’s certificate). In the notation of §1, the canonical boundary b b b of H o − H d H_o-H_d H o − H d satisfies (2), and b = q Υ b ′ b=q\Upsilon b' b = q Υ b ′ with b ′ = q w + R ( 1 − q Υ ) − 1 ( w + Υ v ) ∈ Q [ [ q , R ] ] b'=qw+R(1-q\Upsilon)^{-1}(w+\Upsilon v)\in\mathbb{Q}[[q,R]] b ′ = q w + R ( 1 − q Υ ) − 1 ( w + Υ v ) ∈ Q [[ q , R ]] .
Proof. For λ = 2 \lambda=2 λ = 2 , the data in the proof of Proposition 4.1 are
α ( U ) = 1 1 − q U , β ( U ) = B ( U ) = ( 1 + q U ) ( 1 + q 2 U ) , χ ( U ) = 1 − q 2 U . \alpha(U)=\frac{1}{1-qU},\qquad\beta(U)=B(U)=(1+qU)(1+q^2U),\qquad\chi(U)=1-q^2U. α ( U ) = 1 − q U 1 , β ( U ) = B ( U ) = ( 1 + q U ) ( 1 + q 2 U ) , χ ( U ) = 1 − q 2 U .
and θ ( U ) = ( 1 − q U ) ( 1 + q U ) ( 1 + q 2 U ) \theta(U)=(1-qU)(1+qU)(1+q^2U) θ ( U ) = ( 1 − q U ) ( 1 + q U ) ( 1 + q 2 U ) . Thus u = u α u=u_\alpha u = u α and v = u β v=u_\beta v = u β are the series of §1, and Proposition 3.7}, with ℓ = 2 \ell=2 ℓ = 2 , c = 2 c=2 c = 2 and r = 1 r=1 r = 1 , gives again u − v = q 2 Υ w u-v=q^2\Upsilon w u − v = q 2 Υ w with w ∈ Q [ [ q , R ] ] w\in\mathbb{Q}[[q,R]] w ∈ Q [[ q , R ]] .
In the computations below, g ( Υ ) R ℓ = R ℓ g ( q ℓ Υ ) g(\Upsilon)R^\ell=R^\ell g(q^\ell\Upsilon) g ( Υ ) R ℓ = R ℓ g ( q ℓ Υ ) .
Put B o = ( 1 + R ( 1 − q Υ ) − 1 ) u B_o=(1+R(1-q\Upsilon)^{-1})u B o = ( 1 + R ( 1 − q Υ ) − 1 ) u and B d = ( 1 + R ( 1 + q Υ ) ) v B_d=(1+R(1+q\Upsilon))v B d = ( 1 + R ( 1 + q Υ )) v . By the equations of u u u and v v v ,
( 1 − q Υ − q 2 R 2 ) Υ ( 1 − q Υ ) − 1 u = Υ ( u − R 2 ( 1 − q Υ ) − 1 u ) = 1 , (1-q\Upsilon-q^2R^2)\Upsilon(1-q\Upsilon)^{-1}u=\Upsilon(u-R^2(1-q\Upsilon)^{-1}u)=1, ( 1 − q Υ − q 2 R 2 ) Υ ( 1 − q Υ ) − 1 u = Υ ( u − R 2 ( 1 − q Υ ) − 1 u ) = 1 ,
( 1 − q R − q 2 R Υ ) Υ B d = Υ ( 1 − R ( 1 + q Υ ) ) B d = Υ ( 1 − R 2 B ( Υ ) ) v = 1. (1-qR-q^2R\Upsilon)\Upsilon B_d=\Upsilon(1-R(1+q\Upsilon))B_d=\Upsilon(1-R^2B(\Upsilon))v=1. ( 1 − q R − q 2 R Υ ) Υ B d = Υ ( 1 − R ( 1 + q Υ )) B d = Υ ( 1 − R 2 B ( Υ )) v = 1.
By uniqueness of normalized solutions, m o , 1 = Υ ( 1 − q Υ ) − 1 u m_{o,1}=\Upsilon(1-q\Upsilon)^{-1}u m o , 1 = Υ ( 1 − q Υ ) − 1 u and m d , 1 = Υ B d m_{d,1}=\Upsilon B_d m d , 1 = Υ B d . The same equations give
( 1 − R ) B o = 1 + R ( ( 1 − q Υ ) − 1 − 1 ) u = 1 + q R m o , 1 , (1-R)B_o=1+R((1-q\Upsilon)^{-1}-1)u=1+qRm_{o,1}, ( 1 − R ) B o = 1 + R (( 1 − q Υ ) − 1 − 1 ) u = 1 + q R m o , 1 ,
( 1 − R ) B d = 1 + q R Υ B d = 1 + q R m d , 1 . (1-R)B_d=1+qR\Upsilon B_d=1+qRm_{d,1}. ( 1 − R ) B d = 1 + q R Υ B d = 1 + q R m d , 1 .
Hence ( 1 − R ) ( B o − B d ) = D ( 1 , R ) (1-R)(B_o-B_d)=D(1,R) ( 1 − R ) ( B o − B d ) = D ( 1 , R ) , so b = B o − B d b=B_o-B_d b = B o − B d , which is (2). Finally, u = v + q 2 Υ w u=v+q^2\Upsilon w u = v + q 2 Υ w ,
( 1 − q Υ ) − 1 − ( 1 + q Υ ) = q 2 ( 1 − q Υ ) − 1 (1-q\Upsilon)^{-1}-(1+q\Upsilon)=q^2(1-q\Upsilon)^{-1} ( 1 − q Υ ) − 1 − ( 1 + q Υ ) = q 2 ( 1 − q Υ ) − 1
and q 2 R Υ = q Υ R q^2R\Upsilon=q\Upsilon R q 2 R Υ = q Υ R give
B o − B d = q 2 Υ w + q 2 R Υ ( 1 − q Υ ) − 1 ( w + Υ v ) = q Υ b ′ . B_o-B_d=q^2\Upsilon w+q^2R\Upsilon(1-q\Upsilon)^{-1}(w+\Upsilon v)=q\Upsilon b'. B o − B d = q 2 Υ w + q 2 R Υ ( 1 − q Υ ) − 1 ( w + Υ v ) = q Υ b ′ .
Primc’s three-row identities at level two
We now prove Theorem B , following the strategy of §§1–3. Both partition classes are recognized by automata (Figure 1 ). In A, left multiplication turns the automaton of the admissible partitions into a linear system of q q q -difference equations in the variable x x x of the normal symbols, which marks the number of parts; at R = 1 R=1 R = 1 , it reduces to a four-state system. A quadratic invariant factors this system into two two-state systems, and two telescoping certificates solve these in closed form; both were found by the searches of §3.8. At x = 1 x=1 x = 1 , the Göllnitz–Gordon identities and Jacobi’s triple product then give the products of Theorem B , and Corollary 3.6 converts the resulting equalities into convergent holonomic certificates for the automata.
From automata to four states
The admissible partitions are recognized by the window automaton of Lemma 3.1 , with r = W = 3 r=W=3 r = W = 3 , M = 2 M=2 M = 2 and K = 2 K=2 K = 2 , and the residue side by the congruence automaton with M = 8 M=8 M = 8 (§3.2); Figure 1 shows both for the boundary ( 2 , 0 , 0 ) (2,0,0) ( 2 , 0 , 0 ) . We use the multiplicities τ n \tau_n τ n , β n \beta_n β n , μ n \mu_n μ n and the boundary values defined before Theorem B . Once an even height is complete, the window automaton only needs to remember the pair ( τ , β ) (\tau,\beta) ( τ , β ) read there, and the boundary only fixes the initial pair ( τ 0 , β 0 ) = ( k 2 , k 0 ) (\tau_0,\beta_0)=(k_2,k_0) ( τ 0 , β 0 ) = ( k 2 , k 0 ) . For a pair ( τ , β ) (\tau,\beta) ( τ , β ) , let F τ , β ( z ) F_{\tau,\beta}(z) F τ , β ( z ) be the generating function of the admissible tails: the finitely supported sequences ( μ n , τ n , β n ) n ≥ 1 (\mu_n,\tau_n,\beta_n)_{n\geq1} ( μ n , τ n , β n ) n ≥ 1 satisfying (3) with ( τ 0 , β 0 ) = ( τ , β ) (\tau_0,\beta_0)=(\tau,\beta) ( τ 0 , β 0 ) = ( τ , β ) , with z z z marking the number of parts and q q q their total size, measured from the even height at which ( τ , β ) (\tau,\beta) ( τ , β ) was read; the zero tail contributes 1. The next odd and even heights carry μ \mu μ parts of size 1 and τ ′ + β ′ \tau'+\beta' τ ′ + β ′ parts of size 2, and shift every later part by 2; this gives the factor ( z q ) μ ( z q 2 ) τ ′ + β ′ (zq)^\mu(zq^2)^{\tau'+\beta'} ( z q ) μ ( z q 2 ) τ ′ + β ′ and the substitution z ↦ q 2 z z\mapsto q^2z z ↦ q 2 z , while (3) allows 0 ≤ μ ≤ 2 − max ( τ , τ ′ ) − max ( β , β ′ ) 0\leq\mu\leq2-\max(\tau,\tau')-\max(\beta,\beta') 0 ≤ μ ≤ 2 − max ( τ , τ ′ ) − max ( β , β ′ ) . Hence
F τ , β ( z ) = ∑ τ ′ , β ′ ≥ 0 ( z q 2 ) τ ′ + β ′ ∑ μ = 0 2 − max ( τ , τ ′ ) − max ( β , β ′ ) ( z q ) μ F τ ′ , β ′ ( q 2 z ) . F_{\tau,\beta}(z)=\sum_{\tau',\beta'\geq0}(zq^2)^{\tau'+\beta'}\sum_{\mu=0}^{2-\max(\tau,\tau')-\max(\beta,\beta')}(zq)^\mu F_{\tau',\beta'}(q^2z). F τ , β ( z ) = τ ′ , β ′ ≥ 0 ∑ ( z q 2 ) τ ′ + β ′ μ = 0 ∑ 2 − m a x ( τ , τ ′ ) − m a x ( β , β ′ ) ( z q ) μ F τ ′ , β ′ ( q 2 z ) .
where a negative upper bound contributes zero. This is the opposite realization of §3.5, before normalization, with the letters of two consecutive heights grouped together. Since F τ , β = F β , τ F_{\tau,\beta}=F_{\beta,\tau} F τ , β = F β , τ by top–bottom reflection, the four states 00, 01, 02, 11 remain. For instance, after the pair ( 0 , 0 ) (0,0) ( 0 , 0 ) , the next pair ( 0 , 0 ) (0,0) ( 0 , 0 ) allows μ ≤ 2 \mu\leq2 μ ≤ 2 , the pairs ( 0 , 1 ) (0,1) ( 0 , 1 ) and ( 1 , 0 ) (1,0) ( 1 , 0 ) allow $\mu\leq1, and the pairs ( 0 , 2 ) (0,2) ( 0 , 2 ) , ( 2 , 0 ) (2,0) ( 2 , 0 ) and ( 1 , 1 ) (1,1) ( 1 , 1 ) allow only μ = 0 \mu=0 μ = 0 ; this is the first row of M below. Set F = ( F 00 , F 01 , F 02 , F 11 ) T F=(F_{00},F_{01},F_{02},F_{11})^{\mathsf T} F = ( F 00 , F 01 , F 02 , F 11 ) T , s = q z s=qz s = q z and t = q 2 z t=q^2z t = q 2 z . Then
F ( z ) = M ( z ) F ( q 2 z ) , M ( z ) = ( 1 + s + s 2 2 t ( 1 + s ) 2 t 2 t 2 1 + s t ( 2 + s ) t 2 t 2 1 t t 2 0 1 2 t 0 t 2 ) . (18) F(z)=\mathrm{M}(z)F(q^2z),\qquad\mathrm{M}(z)=\begin{pmatrix}1+s+s^2&2t(1+s)&2t^2&t^2\\1+s&t(2+s)&t^2&t^2\\1&t&t^2&0\\1&2t&0&t^2\end{pmatrix}. \tag*{(18)} F ( z ) = M ( z ) F ( q 2 z ) , M ( z ) = 1 + s + s 2 1 + s 1 1 2 t ( 1 + s ) t ( 2 + s ) t 2 t 2 t 2 t 2 t 2 0 t 2 t 2 0 t 2 . (18)
The empty tail gives F ( 0 ) = ( 1 , 1 , 1 , 1 ) T F(0)=(1,1,1,1)^{\mathsf T} F ( 0 ) = ( 1 , 1 , 1 , 1 ) T . Since M ( 0 ) \mathrm{M}(0) M ( 0 ) is an idempotent of rank one, I − q 2 n M ( 0 ) \mathbb{I}-q^{2n}\mathrm{M}(0) I − q 2 n M ( 0 ) is invertible for n ≥ 1 n\ge1 n ≥ 1 , and coefficient induction proves uniqueness in Q ( q ) [ [ z ] ] 4 \mathbb{Q}(q)[[z]]^4 Q ( q ) [[ z ] ] 4 .
The four generating functions of Theorem B are, in order, Π 2 , 0 , 0 = F 02 ( 1 ) \Pi_{2,0,0}=F_{02}(1) Π 2 , 0 , 0 = F 02 ( 1 ) , Π 1 , 1 , 0 = F 01 ( 1 ) \Pi_{1,1,0}=F_{01}(1) Π 1 , 1 , 0 = F 01 ( 1 ) , Π 1 , 0 , 1 = F 11 ( 1 ) \Pi_{1,0,1}=F_{11}(1) Π 1 , 0 , 1 = F 11 ( 1 ) and Π 0 , 2 , 0 = F 00 ( 1 ) \Pi_{0,2,0}=F_{00}(1) Π 0 , 2 , 0 = F 00 ( 1 ) . With z z z marking the number of parts, F 02 F_{02} F 02 , F 00 F_{00} F 00 , F 01 F_{01} F 01 , F 11 F_{11} F 11 are the series D 2 Λ 0 ( 2 ) D^{(2)}_{2\Lambda_0} D 2 Λ 0 ( 2 ) , D 2 Λ 0 ( 2 ) D^{(2)}_{2\Lambda_0} D 2 Λ 0 ( 2 ) , D Λ 0 + Λ 1 ( 2 ) D^{(2)}_{\Lambda_0+\Lambda_1} D Λ 0 + Λ 1 ( 2 ) , D Λ 0 + Λ 2 ( 2 ) D^{(2)}_{\Lambda_0+\Lambda_2} D Λ 0 + Λ 2 ( 2 ) of [15 §4] .
Factoring the four-state system
The search for quadratic invariants of §3.8, applied to (18), finds
J ( z ) = ( 0 0 0 1 0 − 2 0 0 0 0 2 q z 0 1 0 0 0 ) , λ ( z ) = − q 5 z 3 ( 1 − q z ) , J(z)=\begin{pmatrix}0&0&0&1\\0&-2&0&0\\0&0&2qz&0\\1&0&0&0\end{pmatrix},\qquad\lambda(z)=-q^5z^3(1-qz), J ( z ) = 0 0 0 1 0 − 2 0 0 0 0 2 q z 0 1 0 0 0 , λ ( z ) = − q 5 z 3 ( 1 − q z ) ,
and polynomial arithmetic checks that M ( z ) T J ( z ) M ( z ) = λ ( z ) J ( q 2 z ) \mathrm{M}(z)^{\mathsf T}J(z)\mathrm{M}(z)=\lambda(z)J(q^2z) M ( z ) T J ( z ) M ( z ) = λ ( z ) J ( q 2 z ) . The search takes J J J symmetric of degree at most one in z z z , with 20 unknown coefficients; comparing determinants, with det M ( z ) = q 11 z 6 ( 1 − q z 2 ) \det\mathrm{M}(z)=q^{11}z^6(1-qz^2) det M ( z ) = q 11 z 6 ( 1 − q z 2 ) , suggests λ = ± q 5 z 3 ( 1 − q z ) \lambda=\pm q^5z^3(1-qz) λ = ± q 5 z 3 ( 1 − q z ) . The ten entries of M T J M − λ J ( q 2 z ) \mathrm{M}^{\mathsf T}J\mathrm{M}-\lambda J(q^2z) M T J M − λ J ( q 2 z ) have degree at most five in z z z , so this is a 60 × 20 60\times20 60 × 20 linear system over Q ( q ) \mathbb{Q}(q) Q ( q ) ; its kernel is spanned by J J J for the minus sign and is zero for the plus sign. Since λ ( 0 ) = 0 \lambda(0)=0 λ ( 0 ) = 0 , this forces F T J F = 0 F^{\mathsf T}JF=0 F T J F = 0 , that is,
F 00 ( z ) F 11 ( z ) = F 01 ( z ) 2 − q z F 02 ( z ) 2 (19) F_{00}(z)F_{11}(z)=F_{01}(z)^2-qzF_{02}(z)^2 \tag*{(19)} F 00 ( z ) F 11 ( z ) = F 01 ( z ) 2 − q z F 02 ( z ) 2 (19)
This relation suggests assembling the four series into a matrix of determinant zero. Introduce v v v by v 2 = q z v^2=qz v 2 = q z and define
H ( v ) = ( F 00 ( v 2 / q ) F 01 ( v 2 / q ) + v F 02 ( v 2 / q ) F 01 ( v 2 / q ) − v F 02 ( v 2 / q ) F 11 ( v 2 / q ) ) H(v)=\begin{pmatrix}F_{00}(v^2/q)&F_{01}(v^2/q)+vF_{02}(v^2/q)\\F_{01}(v^2/q)-vF_{02}(v^2/q)&F_{11}(v^2/q)\end{pmatrix} H ( v ) = ( F 00 ( v 2 / q ) F 01 ( v 2 / q ) − v F 02 ( v 2 / q ) F 01 ( v 2 / q ) + v F 02 ( v 2 / q ) F 11 ( v 2 / q ) )
and
Γ ± ( v ) = ( 1 ± v + v 2 q v 2 1 q v 2 ) . \Gamma_{\pm}(v)=\begin{pmatrix}1\pm v+v^2&qv^2\\1&qv^2\end{pmatrix}. Γ ± ( v ) = ( 1 ± v + v 2 1 q v 2 q v 2 ) .
By (18), direct multiplication gives the exact identity
H ( v ) = Γ + ( v ) H ( q v ) Γ − ( v ) T ; (20) H(v)=\Gamma_{+}(v)H(qv)\Gamma_{-}(v)^{\mathsf T}; \tag*{(20)} H ( v ) = Γ + ( v ) H ( q v ) Γ − ( v ) T ; (20)
thus the transfer acts on H H H through the tensor product of two two-state transfers, and (19) reads det H = 0 \det H=0 det H = 0 . Equation (20) has a unique solution with the all-ones constant matrix. Indeed, put Γ 0 = Γ ± ( 0 ) \Gamma_0=\Gamma_{\pm}(0) Γ 0 = Γ ± ( 0 ) and P ( Z ) = Γ 0 Z Γ 0 T \mathcal{P}(Z)=\Gamma_0Z\Gamma_0^{\mathsf T} P ( Z ) = Γ 0 Z Γ 0 T for 2 × 2 2\times2 2 × 2 matrices Z Z Z , and write H ( v ) = ∑ n H n v n H(v)=\sum_n H_nv^n H ( v ) = ∑ n H n v n . The coefficient of v n v^n v n in (20) determines ( I − q n P ) H n (\mathbb{I}-q^n\mathcal{P})H_n ( I − q n P ) H n from H 0 , … , H n − 1 H_0,\ldots,H_{n-1} H 0 , … , H n − 1 ; since P 2 = P \mathcal{P}^2=\mathcal{P} P 2 = P , the operator I − q n P \mathbb{I}-q^n\mathcal{P} I − q n P has, for n > 0 n>0 n > 0 , the explicit inverse
( I − q n P ) − 1 = I + q n 1 − q n P , (\mathbb{I}-q^n\mathcal{P})^{-1}=\mathbb{I}+\frac{q^n}{1-q^n}\mathcal{P}, ( I − q n P ) − 1 = I + 1 − q n q n P ,
and for n = 0 n=0 n = 0 the image of P \mathcal{P} P consists of the multiples of the all-ones matrix.
Let u ( v ) = ( g ( v ) , h ( v ) ) T u(v)=(g(v),h(v))^{\mathsf T} u ( v ) = ( g ( v ) , h ( v ) ) T solve
u ( v ) = Γ + ( v ) u ( q v ) , u ( 0 ) = ( 1 , 1 ) T . (21) u(v)=\Gamma_{+}(v)u(qv),\qquad u(0)=(1,1)^{\mathsf T}. \tag*{(21)} u ( v ) = Γ + ( v ) u ( q v ) , u ( 0 ) = ( 1 , 1 ) T . (21)
The same argument proves uniqueness. Since Γ − ( v ) = Γ + ( − v ) \Gamma_{-}(v)=\Gamma_{+}(-v) Γ − ( v ) = Γ + ( − v ) , the matrix u ( v ) u ( − v ) T u(v)u(-v)^{\mathsf T} u ( v ) u ( − v ) T solves (20) with the correct constant matrix. Consequently
F 00 ( v 2 / q ) = g ( v ) g ( − v ) , F 11 ( v 2 / q ) = h ( v ) h ( − v ) , (22) F_{00}(v^2/q)=g(v)g(-v),\qquad F_{11}(v^2/q)=h(v)h(-v), \tag*{(22)} F 00 ( v 2 / q ) = g ( v ) g ( − v ) , F 11 ( v 2 / q ) = h ( v ) h ( − v ) , (22)
F 01 ( v 2 / q ) = g ( v ) h ( − v ) + h ( v ) g ( − v ) 2 , F 02 ( v 2 / q ) = g ( v ) h ( − v ) − h ( v ) g ( − v ) 2 v . (23) F_{01}(v^2/q)=\frac{g(v)h(-v)+h(v)g(-v)}{2},\qquad F_{02}(v^2/q)=\frac{g(v)h(-v)-h(v)g(-v)}{2v}. \tag*{(23)} F 01 ( v 2 / q ) = 2 g ( v ) h ( − v ) + h ( v ) g ( − v ) , F 02 ( v 2 / q ) = 2 v g ( v ) h ( − v ) − h ( v ) g ( − v ) . (23)
Two finite telescoping certificates
The solution of (21) is
g ( v ) = ∑ n ≥ 0 q n ( n − 1 ) / 2 v n ( − v ; q ) n ( q ; q ) n , h ( v ) = ∑ n ≥ 0 q n ( n + 1 ) / 2 v n ( − v ; q ) n ( q ; q ) n . (24) g(v)=\sum_{n\ge0}\frac{q^{n(n-1)/2}v^n(-v;q)_n}{(q;q)_n},\qquad h(v)=\sum_{n\ge0}\frac{q^{n(n+1)/2}v^n(-v;q)_n}{(q;q)_n}. \tag*{(24)} g ( v ) = n ≥ 0 ∑ ( q ; q ) n q n ( n − 1 ) /2 v n ( − v ; q ) n , h ( v ) = n ≥ 0 ∑ ( q ; q ) n q n ( n + 1 ) /2 v n ( − v ; q ) n . (24)
To verify this, set
f n ( v ) = q n ( n − 1 ) / 2 v n ( − v ; q ) n ( q ; q ) n , ζ = q n , f_n(v)=\frac{q^{n(n-1)/2}v^n(-v;q)_n}{(q;q)_n},\qquad\zeta=q^n, f n ( v ) = ( q ; q ) n q n ( n − 1 ) /2 v n ( − v ; q ) n , ζ = q n ,
so that h = ∑ n ζ f n h=\sum_n\zeta f_n h = ∑ n ζ f n and
f n + 1 ( v ) f n ( v ) = r ( ζ ) : = v ζ ( 1 + v ζ ) 1 − q ζ , f n ( q v ) f n ( v ) = ζ ( 1 + v ζ ) 1 + v . \frac{f_{n+1}(v)}{f_n(v)}=r(\zeta):=\frac{v\zeta(1+v\zeta)}{1-q\zeta},\qquad\frac{f_n(qv)}{f_n(v)}=\frac{\zeta(1+v\zeta)}{1+v}. f n ( v ) f n + 1 ( v ) = r ( ζ ) := 1 − q ζ v ζ ( 1 + v ζ ) , f n ( v ) f n ( q v ) = 1 + v ζ ( 1 + v ζ ) .
The two telescoping certificates
Ψ 0 ( ζ ) = ζ − 1 , Ψ 1 ( ζ ) = − v ζ ( 1 − ζ ) 1 + v \Psi_0(\zeta)=\zeta-1,\qquad\Psi_1(\zeta)=-\frac{v\zeta(1-\zeta)}{1+v} Ψ 0 ( ζ ) = ζ − 1 , Ψ 1 ( ζ ) = − 1 + v v ζ ( 1 − ζ )
satisfy (17) with the following residuals:
r ( ζ ) Ψ 0 ( q ζ ) − Ψ 0 ( ζ ) = 1 − ζ − v ζ ( 1 + v ζ ) , (25) r(\zeta)\Psi_0(q\zeta)-\Psi_0(\zeta)=1-\zeta-v\zeta(1+v\zeta), \tag*{(25)} r ( ζ ) Ψ 0 ( q ζ ) − Ψ 0 ( ζ ) = 1 − ζ − v ζ ( 1 + v ζ ) , (25)
r ( ζ ) Ψ 1 ( q ζ ) − Ψ 1 ( ζ ) = v ζ ( 1 − ζ ) − q v 2 ζ 2 ( 1 + v ζ ) 1 + v . (26) r(\zeta)\Psi_1(q\zeta)-\Psi_1(\zeta)=\frac{v\zeta(1-\zeta)-qv^2\zeta^2(1+v\zeta)}{1+v}. \tag*{(26)} r ( ζ ) Ψ 1 ( q ζ ) − Ψ 1 ( ζ ) = 1 + v v ζ ( 1 − ζ ) − q v 2 ζ 2 ( 1 + v ζ ) . (26)
Multiply by f n ( v ) f_n(v) f n ( v ) and sum from 0 to N N N ; the left sides telescope. Both lower fluxes vanish because Ψ i ( 1 ) = 0 \Psi_i(1)=0 Ψ i ( 1 ) = 0 . The upper fluxes are f N + 1 ( v ) Ψ i ( q N + 1 ) f_{N+1}(v)\Psi_i(q^{N+1}) f N + 1 ( v ) Ψ i ( q N + 1 ) , of v v v -valuation at least N + 1 N+1 N + 1 , so they tend to zero in Q [ [ v ] ] \mathbb{Q}[[v]] Q [[ v ]] . Equations (25)–(26) therefore give
g ( v ) − h ( v ) = v ( 1 + v ) g ( q v ) , h ( v ) = g ( q v ) + q v 2 h ( q v ) . g(v)-h(v)=v(1+v)g(qv),\qquad h(v)=g(qv)+qv^2h(qv). g ( v ) − h ( v ) = v ( 1 + v ) g ( q v ) , h ( v ) = g ( q v ) + q v 2 h ( q v ) .
These are precisely (21), and uniqueness proves (24). Both certificates were found by the search of §3.8: with Ψ \Psi Ψ a polynomial in ζ \zeta ζ of degree at most one, respectively two, and an unknown scalar multiple of the residual, clearing the denominator 1 − q ζ 1-q\zeta 1 − q ζ of r r r turns (17) into linear systems over Q ( q , v ) \mathbb{Q}(q,v) Q ( q , v ) of sizes 4 × 3 4\times3 4 × 3 and 5 × 4 5\times4 5 × 4 , each with one-dimensional kernel.
Göllnitz–Gordon specialization
Put q = r 2 q=r^2 q = r 2 and v = r v=r v = r . The two sums in (24) become the classical Göllnitz–Gordon series:
g ( r ; r 2 ) = ∑ n ≥ 0 r n 2 ( − r ; r 2 ) n ( r 2 ; r 2 ) n = 1 ( r , r 4 , r 7 ; r 8 ) ∞ , (27) g(r;r^2)=\sum_{n\ge0}\frac{r^{n^2}(-r;r^2)_n}{(r^2;r^2)_n}=\frac{1}{(r,r^4,r^7;r^8)_\infty}, \tag*{(27)} g ( r ; r 2 ) = n ≥ 0 ∑ ( r 2 ; r 2 ) n r n 2 ( − r ; r 2 ) n = ( r , r 4 , r 7 ; r 8 ) ∞ 1 , (27)
h ( r ; r 2 ) = ∑ n ≥ 0 r n 2 + 2 n ( − r ; r 2 ) n ( r 2 ; r 2 ) n = 1 ( r 3 , r 4 , r 5 ; r 8 ) ∞ . (28) h(r;r^2)=\sum_{n\ge0}\frac{r^{n^2+2n}(-r;r^2)_n}{(r^2;r^2)_n}=\frac{1}{(r^3,r^4,r^5;r^8)_\infty}. \tag*{(28)} h ( r ; r 2 ) = n ≥ 0 ∑ ( r 2 ; r 2 ) n r n 2 + 2 n ( − r ; r 2 ) n = ( r 3 , r 4 , r 5 ; r 8 ) ∞ 1 . (28)
The final equalities are the Göllnitz–Gordon identities [11 , 13 ] . Since n 2 n^2 n 2 and n 2 + 2 n n^2+2n n 2 + 2 n have the parity of n n n , replacing r r r by − r -r − r in (27)–(28) gives the products for g ( − r ; r 2 ) g(-r;r^2) g ( − r ; r 2 ) and h ( − r ; r 2 ) h(-r;r^2) h ( − r ; r 2 ) .
The series g g g , h h h and F τ , β ( v 2 / q ) F_{\tau,\beta}(v^2/q) F τ , β ( v 2 / q ) lie in Q [ [ q , v ] ] \mathbb{Q}[[q,v]] Q [[ q , v ]] (a tail with k k k parts has size at least k k k , and z k = q − k v 2 k z^k=q^{-k}v^{2k} z k = q − k v 2 k ), and q ↦ r 2 q\mapsto r^2 q ↦ r 2 , v ↦ r v\mapsto r v ↦ r defines a continuous homomorphism Q [ [ q , v ] ] → Q [ [ r ] ] \mathbb{Q}[[q,v]]\to\mathbb{Q}[[r]] Q [[ q , v ]] → Q [[ r ]] ; so the identities above survive this specialization, which corresponds to z = 1 z=1 z = 1 .
Pairing the odd factors in the products at r r r and − r -r − r gives
F 00 ( 1 ; q ) = g ( r ; r 2 ) g ( − r ; r 2 ) = 1 ( q , q 7 ; q 8 ) ∞ ( q 2 ; q 4 ) ∞ 2 , F_{00}(1;q)=g(r;r^2)g(-r;r^2)=\frac{1}{(q,q^7;q^8)_\infty(q^2;q^4)_\infty^2}, F 00 ( 1 ; q ) = g ( r ; r 2 ) g ( − r ; r 2 ) = ( q , q 7 ; q 8 ) ∞ ( q 2 ; q 4 ) ∞ 2 1 ,
F 11 ( 1 ; q ) = h ( r ; r 2 ) h ( − r ; r 2 ) = 1 ( q 3 , q 5 ; q 8 ) ∞ ( q 2 ; q 4 ) ∞ 2 . F_{11}(1;q)=h(r;r^2)h(-r;r^2)=\frac{1}{(q^3,q^5;q^8)_\infty(q^2;q^4)_\infty^2}. F 11 ( 1 ; q ) = h ( r ; r 2 ) h ( − r ; r 2 ) = ( q 3 , q 5 ; q 8 ) ∞ ( q 2 ; q 4 ) ∞ 2 1 .
These are the products for the last and third rows of Theorem B: residues 1, 7, respectively 3, 5, have one colour, and residues 2, 6 have two colours.
The two cross-products
Write θ ( y ; p ) = ( y , p / y ; p ) ∞ \theta(y;p)=(y,p/y;p)_\infty θ ( y ; p ) = ( y , p / y ; p ) ∞ . Jacobi’s triple product [1 Theorem 2.8] gives, in Z [ y 1 ± 1 , y 2 ± 1 ] [ [ p ] ] \mathbb{Z}[y_1^{\pm1},y_2^{\pm1}][[p]] Z [ y 1 ± 1 , y 2 ± 1 ] [[ p ]] ,
θ ( − y 1 ; p ) θ ( y 2 ; p ) + θ ( y 1 ; p ) θ ( − y 2 ; p ) = 2 ( p 2 ; p 2 ) ∞ 2 ( p ; p ) ∞ 2 θ ( y 1 y 2 ; p 2 ) θ ( p y 1 / y 2 ; p 2 ) . (29) \theta(-y_1;p)\theta(y_2;p)+\theta(y_1;p)\theta(-y_2;p)=2\frac{(p^2;p^2)_\infty^2}{(p;p)_\infty^2}\theta(y_1y_2;p^2)\theta(py_1/y_2;p^2). \tag*{(29)} θ ( − y 1 ; p ) θ ( y 2 ; p ) + θ ( y 1 ; p ) θ ( − y 2 ; p ) = 2 ( p ; p ) ∞ 2 ( p 2 ; p 2 ) ∞ 2 θ ( y 1 y 2 ; p 2 ) θ ( p y 1 / y 2 ; p 2 ) . (29)
θ ( − y 1 ; p ) θ ( y 2 ; p ) − θ ( y 1 ; p ) θ ( − y 2 ; p ) = 2 y 1 ( p 2 ; p 2 ) ∞ 2 ( p ; p ) ∞ 2 θ ( y 2 / y 1 ; p 2 ) θ ( p y 1 y 2 ; p 2 ) . (30) \theta(-y_1;p)\theta(y_2;p)-\theta(y_1;p)\theta(-y_2;p)=2y_1\frac{(p^2;p^2)_\infty^2}{(p;p)_\infty^2}\theta(y_2/y_1;p^2)\theta(py_1y_2;p^2). \tag*{(30)} θ ( − y 1 ; p ) θ ( y 2 ; p ) − θ ( y 1 ; p ) θ ( − y 2 ; p ) = 2 y 1 ( p ; p ) ∞ 2 ( p 2 ; p 2 ) ∞ 2 θ ( y 2 / y 1 ; p 2 ) θ ( p y 1 y 2 ; p 2 ) . (30)
These two dissections follow by multiplying the left sides by ( p ; p ) ∞ 2 (p;p)_\infty^2 ( p ; p ) ∞ 2 and expanding the two bilateral series; the sign factor is ( − 1 ) n ± ( − 1 ) m (-1)^n\pm(-1)^m ( − 1 ) n ± ( − 1 ) m . For the plus sign put m = k + l m=k+l m = k + l , n = k − l n=k-l n = k − l : the exponent becomes k 2 + l 2 − k k^2+l^2-k k 2 + l 2 − k and the monomial ( y 1 y 2 ) k ( y 1 / y 2 ) l (y_1y_2)^k(y_1/y_2)^l ( y 1 y 2 ) k ( y 1 / y 2 ) l . For the minus sign put m = k + l + 1 m=k+l+1 m = k + l + 1 , n = k − l n=k-l n = k − l : the exponent becomes k 2 + l 2 + l k^2+l^2+l k 2 + l 2 + l and the monomial y 1 ( y 1 y 2 ) k ( y 1 / y 2 ) l y_1(y_1y_2)^k(y_1/y_2)^l y 1 ( y 1 y 2 ) k ( y 1 / y 2 ) l . In each case the two independent sums are triple products.
The products g ( r ) h ( − r ) g(r)h(-r) g ( r ) h ( − r ) and h ( r ) g ( − r ) h(r)g(-r) h ( r ) g ( − r ) have common denominator
D = ( r 4 , r 8 ) ∞ 2 ( r 2 , r 6 , r 10 , r 14 ; r 16 ) ∞ , D=(r^4,r^8)_\infty^2(r^2,r^6,r^{10},r^{14};r^{16})_\infty, D = ( r 4 , r 8 ) ∞ 2 ( r 2 , r 6 , r 10 , r 14 ; r 16 ) ∞ ,
and their sum and difference have numerators θ ( − r ; r 8 ) θ ( r 3 ; r 8 ) ± θ ( r ; r 8 ) θ ( − r 3 ; r 8 ) \theta(-r;r^8)\theta(r^3;r^8)\pm\theta(r;r^8)\theta(-r^3;r^8) θ ( − r ; r 8 ) θ ( r 3 ; r 8 ) ± θ ( r ; r 8 ) θ ( − r 3 ; r 8 ) . Take y 1 = r y_1=r y 1 = r , y 2 = r 3 y_2=r^3 y 2 = r 3 , p = r 8 p=r^8 p = r 8 in the two dissections. Since ( r 16 ; r 16 ) ∞ 2 / ( r 8 ; r 8 ) ∞ 2 = ( r 8 ; r 16 ) ∞ − 2 (r^{16};r^{16})_\infty^2/(r^8;r^8)_\infty^2=(r^8;r^{16})_\infty^{-2} ( r 16 ; r 16 ) ∞ 2 / ( r 8 ; r 8 ) ∞ 2 = ( r 8 ; r 16 ) ∞ − 2 , these numerators are
2 ( r 4 , r 6 , r 10 , r 12 ; r 16 ) ∞ ( r 8 ; r 16 ) ∞ 2 , 2 r ( r 2 , r 4 , r 12 , r 14 ; r 16 ) ∞ ( r 8 ; r 16 ) ∞ 2 . \frac{2(r^4,r^6,r^{10},r^{12};r^{16})_\infty}{(r^8;r^{16})_\infty^2},
\qquad
\frac{2r(r^2,r^4,r^{12},r^{14};r^{16})_\infty}{(r^8;r^{16})_\infty^2}. ( r 8 ; r 16 ) ∞ 2 2 ( r 4 , r 6 , r 10 , r 12 ; r 16 ) ∞ , ( r 8 ; r 16 ) ∞ 2 2 r ( r 2 , r 4 , r 12 , r 14 ; r 16 ) ∞ .
Substituting into (23) and cancelling factors gives
F 01 ( 1 ; r 2 ) = 1 ( r 2 , r 4 , r 8 , r 8 , r 12 , r 14 ; r 16 ) ∞ , F_{01}(1;r^2)=\frac{1}{(r^2,r^4,r^8,r^8,r^{12},r^{14};r^{16})_\infty}, F 01 ( 1 ; r 2 ) = ( r 2 , r 4 , r 8 , r 8 , r 12 , r 14 ; r 16 ) ∞ 1 ,
F 02 ( 1 ; r 2 ) = 1 ( r 4 , r 6 , r 8 , r 8 , r 10 , r 12 ; r 16 ) ∞ , F_{02}(1;r^2)=\frac{1}{(r^4,r^6,r^8,r^8,r^{10},r^{12};r^{16})_\infty}, F 02 ( 1 ; r 2 ) = ( r 4 , r 6 , r 8 , r 8 , r 10 , r 12 ; r 16 ) ∞ 1 ,
giving the second and first rows of Theorem B. In both products residue 4 has two colours; the other four listed residues have one colour each. The bilateral series converge r r r -adically, since the triple-product summands have positive quadratic valuation in both summation indices. This completes the proof of Theorem B.
Certificates for the automata
Let S S S be the difference of the two automata of Figure 1 . For the boundary ( 2 , 0 , 0 ) (2,0,0) ( 2 , 0 , 0 ) , they have 29 and 9 interior states, so s = 38 s=38 s = 38 , and the homogeneous kernels G a , G b G_a,G_b G a , G b of (12) and (10) have sizes 78 and 40. The proof above establishes ε ( S ) = 0 \varepsilon(S)=0 ε ( S ) = 0 , so Corollary 3.6 supplies convergent holonomic certificates for the original automata, with ν ( a i , j ) ≥ i + j \nu(a_{i,j})\ge i+j ν ( a i , j ) ≥ i + j and ν ( b n ) ≥ n + 1 \nu(b_n)\ge n+1 ν ( b n ) ≥ n + 1 , without eliminating their matrices to scalar operators. The implementation (§ 4.8 ) lists both automata and computes the sparse operator entries of these kernels.
The Rogers–Ramanujan identities
This classical case shows that the whole proof, product side included, reduces to finite rational identities of the kind found by the searches of § 3.8 . Write ( q ) n = ( q ; q ) n (q)_n=(q;q)_n ( q ) n = ( q ; q ) n . For i = 1 , 2 i=1,2 i = 1 , 2 , the language
L i = R i R ∗ U ( R 2 R ∗ U ) ∗ \mathcal{L}_i=R^iR^*U(R^2R^*U)^* L i = R i R ∗ U ( R 2 R ∗ U ) ∗
encodes nonempty partitions with successive parts differing by at least two and smallest part at least i i i : successive blocks of R R R ’s record the smallest part and the successive gaps. Subtracting the staircase i , i + 2 , … , i + 2 ( n − 1 ) i,i+2,\ldots,i+2(n-1) i , i + 2 , … , i + 2 ( n − 1 ) gives
1 + ε ( H L i ) = ∑ n ≥ 0 q n 2 + ( i − 1 ) n ( q ) n . 1+\varepsilon(H_{\mathcal{L}_i})=\sum_{n\ge0}\frac{q^{n^2+(i-1)n}}{(q)_n}. 1 + ε ( H L i ) = n ≥ 0 ∑ ( q ) n q n 2 + ( i − 1 ) n .
Let P 1 , P 2 \mathcal{P}_1,\mathcal{P}_2 P 1 , P 2 be the regular languages of partitions into parts with residues in Σ 1 = { 1 , 4 } \Sigma_1=\{1,4\} Σ 1 = { 1 , 4 } , respectively Σ 2 = { 2 , 3 } \Sigma_2=\{2,3\} Σ 2 = { 2 , 3 } , modulo five; they have five height states, with a U U U -loop at each allowed residue, see Figure 2 . Both pairs of automata are instances of §3.2. The Rogers–Ramanujan identities [19 ] assert
1 + ε ( H L 1 ) = 1 ( q , q 4 ; q 5 ) ∞ , 1 + ε ( H L 2 ) = 1 ( q 2 , q 3 ; q 5 ) ∞ . (31) 1+\varepsilon(H_{\mathcal{L}_1})=\frac{1}{(q,q^4;q^5)_\infty},\qquad1+\varepsilon(H_{\mathcal{L}_2})=\frac{1}{(q^2,q^3;q^5)_\infty}. \tag*{(31)} 1 + ε ( H L 1 ) = ( q , q 4 ; q 5 ) ∞ 1 , 1 + ε ( H L 2 ) = ( q 2 , q 3 ; q 5 ) ∞ 1 . (31)
We derive (31) from finite identities, by two Bailey steps and Jacobi’s triple product.
For a parameter γ \gamma γ , call ( α , β ) (\alpha,\beta) ( α , β ) a Bailey pair if
β m = ∑ r = 0 m α r ( q ) m − r ( γ q ; q ) m + r . (32) \beta_m=\sum_{r=0}^{m}\frac{\alpha_r}{(q)_{m-r}(\gamma q;q)_{m+r}}. \tag*{(32)} β m = r = 0 ∑ m ( q ) m − r ( γ q ; q ) m + r α r . (32)
The unit pairs needed here have β m = δ m , 0 \beta_m=\delta_{m,0} β m = δ m , 0 and
γ = 1 : α 0 = 1 , α r = ( − 1 ) r q r ( r − 1 ) / 2 ( 1 + q r ) ( r > 0 ) , (33) \gamma=1:\quad\alpha_0=1,\qquad\alpha_r=(-1)^r q^{r(r-1)/2}(1+q^r)\quad(r>0), \tag*{(33)} γ = 1 : α 0 = 1 , α r = ( − 1 ) r q r ( r − 1 ) /2 ( 1 + q r ) ( r > 0 ) , (33)
γ = q : α r = ( − 1 ) r q r ( r − 1 ) / 2 ( 1 − q 2 r + 1 ) 1 − q ( r ≥ 0 ) . (34) \gamma=q:\quad\alpha_r=\frac{(-1)^r q^{r(r-1)/2}(1-q^{2r+1})}{1-q}\quad(r\ge0). \tag*{(34)} γ = q : α r = 1 − q ( − 1 ) r q r ( r − 1 ) /2 ( 1 − q 2 r + 1 ) ( r ≥ 0 ) . (34)
Both are finite identities. For (33), pairing the indices r , − r r,-r r , − r turns the right side of (32) into
∑ r = − m m ( − 1 ) r q r ( r − 1 ) / 2 ( q ) m − r ( q ) m + r = ( − 1 ) m q m ( m + 1 ) / 2 ( q ) 2 m ( q − m ; q ) 2 m . \sum_{r=-m}^{m}\frac{(-1)^r q^{r(r-1)/2}}{(q)_{m-r}(q)_{m+r}}=\frac{(-1)^m q^{m(m+1)/2}}{(q)_{2m}}(q^{-m};q)_{2m}. r = − m ∑ m ( q ) m − r ( q ) m + r ( − 1 ) r q r ( r − 1 ) /2 = ( q ) 2 m ( − 1 ) m q m ( m + 1 ) /2 ( q − m ; q ) 2 m .
For (34), pair r , − r − 1 r,-r-1 r , − r − 1 instead, obtaining
∑ r = − m − 1 m ( − 1 ) r q r ( r − 1 ) / 2 ( q ) m − r ( q ) m + r + 1 = ( − 1 ) m + 1 q ( m + 1 ) ( m + 2 ) / 2 ( q ) 2 m + 1 ( q − m − 1 ; q ) 2 m + 1 . \sum_{r=-m-1}^{m}\frac{(-1)^r q^{r(r-1)/2}}{(q)_{m-r}(q)_{m+r+1}}=\frac{(-1)^{m+1}q^{(m+1)(m+2)/2}}{(q)_{2m+1}}(q^{-m-1};q)_{2m+1}. r = − m − 1 ∑ m ( q ) m − r ( q ) m + r + 1 ( − 1 ) r q r ( r − 1 ) /2 = ( q ) 2 m + 1 ( − 1 ) m + 1 q ( m + 1 ) ( m + 2 ) /2 ( q − m − 1 ; q ) 2 m + 1 .
The right-hand sides vanish for m > 0 m > 0 m > 0 and equal 1 for m = 0 m = 0 m = 0 . These evaluations are instances of the finite q q q -binomial identity
( y ; q ) N = ∑ j = 0 N ( − 1 ) j q j ( j − 1 ) / 2 [ N j ] q y j , [ N j ] q = ( q ) N ( q ) j ( q ) N − j . (y;q)_N=\sum_{j=0}^{N}(-1)^j q^{j(j-1)/2}\begin{bmatrix}N\\j\end{bmatrix}_q y^j,\qquad\begin{bmatrix}N\\j\end{bmatrix}_q=\frac{(q)_N}{(q)_j(q)_{N-j}}. ( y ; q ) N = j = 0 ∑ N ( − 1 ) j q j ( j − 1 ) /2 [ N j ] q y j , [ N j ] q = ( q ) j ( q ) N − j ( q ) N .
which follows coefficientwise from ( y ; q ) N + 1 = ( 1 − y q N ) ( y ; q ) N (y;q)_{N+1}=(1-yq^N)(y;q)_N ( y ; q ) N + 1 = ( 1 − y q N ) ( y ; q ) N and the Gaussian Pascal recurrence.
The Bailey step has a rational telescoping certificate. Put
h ( m , k ) = κ k q k 2 ( q ) m − k ( q ) k ( κ q ; q ) k , ζ = q m , ξ = q k , h(m,k)=\frac{\kappa^kq^{k^2}}{(q)_{m-k}(q)_k(\kappa q;q)_k},\qquad\zeta=q^m,\qquad\xi=q^k, h ( m , k ) = ( q ) m − k ( q ) k ( κ q ; q ) k κ k q k 2 , ζ = q m , ξ = q k ,
and
r ( ζ , ξ ) = κ q ξ 2 ( 1 − ζ / ξ ) ( 1 − q ξ ) ( 1 − κ q ξ ) , Ψ ( ζ , ξ ) = − ζ ( 1 − ξ ) ( 1 − κ ξ ) ξ , r(\zeta,\xi)=\frac{\kappa q\xi^2(1-\zeta/\xi)}{(1-q\xi)(1-\kappa q\xi)},\qquad\Psi(\zeta,\xi)=-\frac{\zeta(1-\xi)(1-\kappa\xi)}{\xi}, r ( ζ , ξ ) = ( 1 − q ξ ) ( 1 − κ q ξ ) κ q ξ 2 ( 1 − ζ / ξ ) , Ψ ( ζ , ξ ) = − ξ ζ ( 1 − ξ ) ( 1 − κ ξ ) ,
so that r ( ζ , ξ ) = h ( m , k + 1 ) / h ( m , k ) r(\zeta,\xi)=h(m,k+1)/h(m,k) r ( ζ , ξ ) = h ( m , k + 1 ) / h ( m , k ) . The rational identity
r ( ζ , ξ ) Ψ ( ζ , q ξ ) − Ψ ( ζ , ξ ) = ( 1 − ζ ) ( 1 − κ ζ ) − ( 1 − ζ / ξ ) (35) r(\zeta,\xi)\Psi(\zeta,q\xi)-\Psi(\zeta,\xi)=(1-\zeta)(1-\kappa\zeta)-(1-\zeta/\xi) \tag*{(35)} r ( ζ , ξ ) Ψ ( ζ , q ξ ) − Ψ ( ζ , ξ ) = ( 1 − ζ ) ( 1 − κ ζ ) − ( 1 − ζ / ξ ) (35)
is checked by clearing denominators. Extend h ( m , k ) h(m,k) h ( m , k ) by zero outside 0 ≤ k ≤ m 0\leq k\leq m 0 ≤ k ≤ m . Multiplying by h ( m , k ) h(m,k) h ( m , k ) and summing over k k k gives ( 1 − q m ) ( 1 − κ q m ) F m = F m − 1 (1-q^m)(1-\kappa q^m)F_m=F_{m-1} ( 1 − q m ) ( 1 − κ q m ) F m = F m − 1 for F m = ∑ k h ( m , k ) F_m=\sum_k h(m,k) F m = ∑ k h ( m , k ) , because ( 1 − q m − k ) h ( m , k ) = h ( m − 1 , k ) (1-q^{m-k})h(m,k)=h(m-1,k) ( 1 − q m − k ) h ( m , k ) = h ( m − 1 , k ) ; the lower flux vanishes at ξ = 1 \xi=1 ξ = 1 , the upper flux vanishes because h ( m , m + 1 ) = 0 h(m,m+1)=0 h ( m , m + 1 ) = 0 , and F 0 = 1 F_0=1 F 0 = 1 . Thus
∑ k = 0 m h ( m , k ) = 1 ( q ) m ( κ q ; q ) m . (36) \sum_{k=0}^{m}h(m,k)=\frac{1}{(q)_m(\kappa q;q)_m}. \tag*{(36)} k = 0 ∑ m h ( m , k ) = ( q ) m ( κ q ; q ) m 1 . (36)
The certificate Ψ \Psi Ψ is found by the same search, with denominator ξ \xi ξ , a numerator of degree at most two in ξ \xi ξ over Q ( q , κ , ζ ) \mathbb{Q}(q,\kappa,\zeta) Q ( q , κ , ζ ) , and an unknown multiple of the right side of (35) : clearing ξ ( 1 − q ξ ) ( 1 − κ q ξ ) \xi(1-q\xi)(1-\kappa q\xi) ξ ( 1 − q ξ ) ( 1 − κ q ξ ) gives a 5 × 4 5\times4 5 × 4 linear system with one-dimensional kernel. Substituting (32) into
α r ′ = γ r q r 2 α r , β m ′ = ∑ j = 0 m γ j q j 2 β j ( q ) m − j , \alpha'_r=\gamma^r q^{r^2}\alpha_r,\qquad\beta'_m=\sum_{j=0}^{m}\frac{\gamma^j q^{j^2}\beta_j}{(q)_{m-j}}, α r ′ = γ r q r 2 α r , β m ′ = j = 0 ∑ m ( q ) m − j γ j q j 2 β j ,
and applying (36) with m ↦ m − r m\mapsto m-r m ↦ m − r , k = j − r k=j-r k = j − r and κ = γ q 2 r \kappa=\gamma q^{2r} κ = γ q 2 r , proves that ( α ′ , β ′ ) (\alpha',\beta') ( α ′ , β ′ ) is again a Bailey pair.
Apply this finite transformation twice to either unit pair, and multiply the relation at index M M M by ( q ) M (q)_M ( q ) M . The result is the finite equality
∑ n = 0 M γ n q n 2 ( q ) n ( q ) M ( q ) M − n = ∑ r = 0 M γ 2 r q 2 r 2 α r ( q ) M ( q ) M − r ( γ q ; q ) M + r . (37) \sum_{n=0}^{M}\frac{\gamma^n q^{n^2}}{(q)_n}\frac{(q)_M}{(q)_{M-n}}=\sum_{r=0}^{M}\gamma^{2r}q^{2r^2}\alpha_r\frac{(q)_M}{(q)_{M-r}(\gamma q;q)_{M+r}}. \tag*{(37)} n = 0 ∑ M ( q ) n γ n q n 2 ( q ) M − n ( q ) M = r = 0 ∑ M γ 2 r q 2 r 2 α r ( q ) M − r ( γ q ; q ) M + r ( q ) M . (37)
Put γ = q δ \gamma=q^\delta γ = q δ with δ ∈ { 0 , 1 } \delta\in\{0,1\} δ ∈ { 0 , 1 } . Both sides converge as M → ∞ M\to\infty M → ∞ , with errors of valuation at least M + 1 + δ M+1+\delta M + 1 + δ . On the left, a nonconstant factorial ratio contributes at least M − n + 1 M-n+1 M − n + 1 , and
n 2 + δ n + M − n + 1 ≥ M + 1 + δ ( n ≥ 1 ) . n^2+\delta n+M-n+1\geq M+1+\delta\qquad(n\geq1). n 2 + δ n + M − n + 1 ≥ M + 1 + δ ( n ≥ 1 ) .
On the right the numerator has valuation w δ ( r ) = ( 5 r 2 + ( 4 δ − 1 ) r ) / 2 w_\delta(r)=(5r^2+(4\delta-1)r)/2 w δ ( r ) = ( 5 r 2 + ( 4 δ − 1 ) r ) /2 for r > 0 r>0 r > 0 , and w δ ( r ) + M − r + 1 ≥ M + 1 + δ w_\delta(r)+M-r+1\geq M+1+\delta w δ ( r ) + M − r + 1 ≥ M + 1 + δ . The denominator replacement at r = 0 r=0 r = 0 begins at M + 1 + δ M+1+\delta M + 1 + δ ; the remaining denominator errors and omitted tails have larger valuation. Hence (37) gives
∑ n ≥ 0 q n 2 ( q ) n = 1 ( q ) ∞ ∑ r ∈ Z ( − 1 ) r q ( 5 r 2 − r ) / 2 , \sum_{n\geq0}\frac{q^{n^2}}{(q)_n}=\frac{1}{(q)_\infty}\sum_{r\in\mathbb{Z}}(-1)^r q^{(5r^2-r)/2}, n ≥ 0 ∑ ( q ) n q n 2 = ( q ) ∞ 1 r ∈ Z ∑ ( − 1 ) r q ( 5 r 2 − r ) /2 ,
∑ n ≥ 0 q n 2 + n ( q ) n = 1 ( q ) ∞ ∑ r ∈ Z ( − 1 ) r q ( 5 r 2 + 3 r ) / 2 . \sum_{n\geq0}\frac{q^{n^2+n}}{(q)_n}=\frac{1}{(q)_\infty}\sum_{r\in\mathbb{Z}}(-1)^r q^{(5r^2+3r)/2}. n ≥ 0 ∑ ( q ) n q n 2 + n = ( q ) ∞ 1 r ∈ Z ∑ ( − 1 ) r q ( 5 r 2 + 3 r ) /2 .
The second equality pairs r ≥ 0 r\geq0 r ≥ 0 with − r − 1 -r-1 − r − 1 , and cancels the factor 1 − q 1-q 1 − q against ( q 2 ; q ) ∞ (q^2;q)_\infty ( q 2 ; q ) ∞ .
Jacobi’s triple product follows from a finite identity:
( y ; p ) N ( p / y ; p ) N = ∑ k = − N N ( − 1 ) k p k ( k − 1 ) / 2 y k ( 2 N N + k ) p . (y;p)_N(p/y;p)_N=\sum_{k=-N}^{N}(-1)^k p^{k(k-1)/2}y^k\binom{2N}{N+k}_p. ( y ; p ) N ( p / y ; p ) N = k = − N ∑ N ( − 1 ) k p k ( k − 1 ) /2 y k ( N + k 2 N ) p .
This is the preceding finite q q q -binomial theorem, after a shift of the summation index. At ( p , y ) = ( q 5 , q 2 ) (p,y)=(q^5,q^2) ( p , y ) = ( q 5 , q 2 ) , respectively ( q 5 , q 4 ) (q^5,q^4) ( q 5 , q 4 ) , its two sides converge with errors of valuation at least 5 N + 2 5N+2 5 N + 2 , respectively 5 N + 1 5N+1 5 N + 1 : the Gaussian coefficient differs from 1 / ( q 5 ; q 5 ) ∞ 1/(q^5;q^5)_\infty 1/ ( q 5 ; q 5 ) ∞ by a series of valuation at least 5 ( N − ∣ k ∣ + 1 ) 5(N-|k|+1) 5 ( N − ∣ k ∣ + 1 ) , and combining this with the quadratic theta exponent gives the stated bounds, including the omitted tail. The resulting products are ( q 2 , q 3 , q 5 ; q 5 ) ∞ (q^2,q^3,q^5;q^5)_\infty ( q 2 , q 3 , q 5 ; q 5 ) ∞ and ( q , q 4 , q 5 ; q 5 ) ∞ (q,q^4,q^5;q^5)_\infty ( q , q 4 , q 5 ; q 5 ) ∞ . Cancellation with ( q ) ∞ (q)_\infty ( q ) ∞ proves (31).
To bind this calculation to the quantum-plane certificates, put S i = H L i − H P i S_i=H_{L_i}-H_{P_i} S i = H L i − H P i . Each gap language has three interior states and each residue language five, so the normalized difference system has eight states. Let E E E be the 5 × 5 5\times5 5 × 5 cyclic matrix with E j , j − 1 = 1 E_{j,j-1}=1 E j , j − 1 = 1 , indices read modulo five; let I Σ I_\Sigma I Σ be the diagonal matrix with entries 1 j ∈ Σ 1_{j\in\Sigma} 1 j ∈ Σ , and χ Σ \chi_\Sigma χ Σ its diagonal column. Set
L g a p = ( 1 − q R 0 − q R − q Υ 1 0 0 − q R 1 ) , L Σ = I 5 − q R E − q Υ I Σ . L_{\mathrm{gap}}=\begin{pmatrix}1-qR&0&-qR\\-q\Upsilon&1&0\\0&-qR&1\end{pmatrix},\qquad L_\Sigma=I_5-qRE-q\Upsilon I_\Sigma. L gap = 1 − q R − q Υ 0 0 1 − q R − q R 0 1 , L Σ = I 5 − q R E − q Υ I Σ .
The data for S i S_i S i are
L i = diag ( L g a p , L Σ i ) , f i = ( e 2 i − 1 e 1 ) , d i = q R ( e 1 T − χ Σ i T ) . L_i=\operatorname{diag}(L_{\mathrm{gap}},L_{\Sigma_i}),\qquad f_i=\begin{pmatrix}e_{2i-1}\\e_1\end{pmatrix},\qquad d_i=qR\begin{pmatrix}e_1^{\mathsf T}&-\chi_{\Sigma_i}^{\mathsf T}\end{pmatrix}. L i = diag ( L gap , L Σ i ) , f i = ( e 2 i − 1 e 1 ) , d i = q R ( e 1 T − χ Σ i T ) .
where the first basis vector in f i f_i f i belongs to Q 3 \mathbb{Q}^3 Q 3 and the second to Q 5 \mathbb{Q}^5 Q 5 . For the marked system use L i , x = L i ∣ Υ ↦ x Υ L_{i,x}=L_i|_{\Upsilon\mapsto x\Upsilon} L i , x = L i ∣ Υ ↦ x Υ , f i , x = x f i f_{i,x}=xf_i f i , x = x f i , d i , x = d i d_{i,x}=d_i d i , x = d i . The general construction gives G b G_b G b of size 10 and G a G_a G a of size 18; with their constant coordinate set to 1 and their state vectors the normalized solutions, their last coordinates are the required a , b a,b a , b . Since the identities above give ε ( S i ) = 0 \varepsilon(S_i)=0 ε ( S i ) = 0 , Corollary 3.6 yields ν ( a i , j ) ≥ i + j \nu(a_{i,j})\ge i+j ν ( a i , j ) ≥ i + j and ν ( b n ) ≥ n + 1 \nu(b_n)\ge n+1 ν ( b n ) ≥ n + 1 , which completes both Rogers–Ramanujan certificates.
A global constraint: alternating sums
For a partition λ = ( λ 1 ≥ ⋯ ≥ λ ℓ > 0 ) \lambda=(\lambda_1\ge\cdots\ge\lambda_\ell>0) λ = ( λ 1 ≥ ⋯ ≥ λ ℓ > 0 ) , put
α ( λ ) = λ 1 − λ 2 + λ 3 − ⋯ + ( − 1 ) ℓ − 1 λ ℓ . \alpha(\lambda)=\lambda_1-\lambda_2+\lambda_3-\cdots+(-1)^{\ell-1}\lambda_\ell. α ( λ ) = λ 1 − λ 2 + λ 3 − ⋯ + ( − 1 ) ℓ − 1 λ ℓ .
Bessenrodt’s refinement of Euler’s identity, obtained from Sylvester’s bijection [5 , 7 ] , states that for every n , k ≥ 0 n,k\ge0 n , k ≥ 0 the number of distinct-part partitions of n n n with α ( λ ) = k \alpha(\lambda)=k α ( λ ) = k equals the number of odd-part partitions of n n n with exactly k k k parts. This is a global constraint on the partition, rather than a restriction on an individual part or on an individual adjacent difference. Its weighted form is
∑ λ d i s t i n c t z α ( λ ) q ∣ λ ∣ = ∑ μ o d d z ℓ ( μ ) q ∣ μ ∣ = 1 ( z q ; q 2 ) ∞ . (38) \sum_{\lambda\ \mathrm{distinct}}z^{\alpha(\lambda)}q^{|\lambda|}
=\sum_{\mu\ \mathrm{odd}}z^{\ell(\mu)}q^{|\mu|}
=\frac{1}{(zq;q^2)_\infty}. \tag*{(38)} λ distinct ∑ z α ( λ ) q ∣ λ ∣ = μ odd ∑ z ℓ ( μ ) q ∣ μ ∣ = ( z q ; q 2 ) ∞ 1 . (38)
The refinement is still recognized by a finite weighted automaton. In the boundary word of a distinct partition, give a letter R R R an additional weight z z z precisely when an odd number of letters U U U follow it; the total weight is z α ( λ ) z^{\alpha(\lambda)} z α ( λ ) . Writing
B 0 = R ( 1 − R ) − 1 U , B 1 = z R ( 1 − z R ) − 1 U , B_0=R(1-R)^{-1}U,\qquad B_1=zR(1-zR)^{-1}U, B 0 = R ( 1 − R ) − 1 U , B 1 = z R ( 1 − z R ) − 1 U ,
the two weighted regular series, including the empty partition, are
H d e = ( 1 + B 1 ) ( 1 − B 0 B 1 ) − 1 , H d o = 1 + R ( 1 − z U − R 2 ) − 1 z U . (39) H_d^{\mathrm{e}}=(1+B_1)(1-B_0B_1)^{-1},\qquad H_d^{\mathrm{o}}=1+R(1-zU-R^2)^{-1}zU. \tag*{(39)} H d e = ( 1 + B 1 ) ( 1 − B 0 B 1 ) − 1 , H d o = 1 + R ( 1 − z U − R 2 ) − 1 z U . (39)
Indeed, even-length distinct partitions have block pattern ( B 0 B 1 ) ∗ (B_0B_1)^* ( B 0 B 1 ) ∗ , and odd-length ones have pattern B 1 ( B 0 B 1 ) ∗ B_1(B_0B_1)^* B 1 ( B 0 B 1 ) ∗ ; the second language allows U U U only at odd heights, each weighted z z z . Extracting [ z k ] [z^k] [ z k ] is also a finite-state operation: record the accumulated exponent up to k k k and reject paths that exceed it. Thus the identities at each fixed value of the statistic are also between regular languages.
We prove (38) by telescoping. For a distinct partition of length m m m , set δ i = λ i − λ i + 1 ≥ 1 \delta_i=\lambda_i-\lambda_{i+1}\ge1 δ i = λ i − λ i + 1 ≥ 1 , where λ m + 1 = 0 \lambda_{m+1}=0 λ m + 1 = 0 . Then
∣ λ ∣ = ∑ i = 1 m i δ i , α ( λ ) = ∑ 1 ≤ i ≤ m i o d d δ i . |\lambda|=\sum_{i=1}^{m}i\delta_i,\qquad\alpha(\lambda)=\sum_{\substack{1\le i\le m\\i\ \mathrm{odd}}}\delta_i. ∣ λ ∣ = i = 1 ∑ m i δ i , α ( λ ) = 1 ≤ i ≤ m i odd ∑ δ i .
Pairing lengths 2 j − 1 2j-1 2 j − 1 and 2 j 2j 2 j therefore gives
z j q 2 j 2 − j ( z q ; q 2 ) j ( q 2 ; q 2 ) j − 1 + z j q 2 j 2 + j ( z q ; q 2 ) j ( q 2 ; q 2 ) j = z j q 2 j 2 − j ( z q ; q 2 ) j ( q 2 ; q 2 ) j . \frac{z^{j}q^{2j^{2}-j}}{(zq;q^{2})_{j}(q^{2};q^{2})_{j-1}}+\frac{z^{j}q^{2j^{2}+j}}{(zq;q^{2})_{j}(q^{2};q^{2})_{j}}=\frac{z^{j}q^{2j^{2}-j}}{(zq;q^{2})_{j}(q^{2};q^{2})_{j}}. ( z q ; q 2 ) j ( q 2 ; q 2 ) j − 1 z j q 2 j 2 − j + ( z q ; q 2 ) j ( q 2 ; q 2 ) j z j q 2 j 2 + j = ( z q ; q 2 ) j ( q 2 ; q 2 ) j z j q 2 j 2 − j .
Including the empty partition and putting η = z q \eta=zq η = z q , Q = q 2 Q=q^{2} Q = q 2 , the left side of (38) is
V ( η ) = ∑ j ≥ 0 f j ( η ) , f j ( η ) = η j Q j ( j − 1 ) ( Q ; Q ) j ( η ; Q ) j . V(\eta)=\sum_{j\ge0}f_j(\eta),\qquad f_j(\eta)=\frac{\eta^{j}Q^{j(j-1)}}{(Q;Q)_j(\eta;Q)_j}. V ( η ) = j ≥ 0 ∑ f j ( η ) , f j ( η ) = ( Q ; Q ) j ( η ; Q ) j η j Q j ( j − 1 ) .
For ξ = Q j \xi=Q^j ξ = Q j , its two rational ratios are
r ( ξ ) = f j + 1 ( η ) f j ( η ) = η ξ 2 ( 1 − Q ξ ) ( 1 − η ξ ) , s ( ξ ) = f j ( Q η ) f j ( η ) = ξ ( 1 − η ) 1 − η ξ . r(\xi)=\frac{f_{j+1}(\eta)}{f_j(\eta)}=\frac{\eta\xi^{2}}{(1-Q\xi)(1-\eta\xi)},\qquad s(\xi)=\frac{f_j(Q\eta)}{f_j(\eta)}=\frac{\xi(1-\eta)}{1-\eta\xi}. r ( ξ ) = f j ( η ) f j + 1 ( η ) = ( 1 − Q ξ ) ( 1 − η ξ ) η ξ 2 , s ( ξ ) = f j ( η ) f j ( Q η ) = 1 − η ξ ξ ( 1 − η ) .
The telescoping certificate Ψ ( ξ ) = ( 1 − η ) ( ξ − 1 ) \Psi(\xi)=(1-\eta)(\xi-1) Ψ ( ξ ) = ( 1 − η ) ( ξ − 1 ) satisfies
r ( ξ ) Ψ ( Q ξ ) − Ψ ( ξ ) = ( 1 − η ) − s ( ξ ) (40) r(\xi)\Psi(Q\xi)-\Psi(\xi)=(1-\eta)-s(\xi) \tag*{(40)} r ( ξ ) Ψ ( Q ξ ) − Ψ ( ξ ) = ( 1 − η ) − s ( ξ ) (40)
Its lower boundary is Ψ ( 1 ) = 0 \Psi(1)=0 Ψ ( 1 ) = 0 , and the upper flux at j = N + 1 j=N+1 j = N + 1 has η \eta η -valuation at least N + 1 N+1 N + 1 ; after η = z q \eta=zq η = z q , Q = q 2 Q=q^{2} Q = q 2 , it has q q q -valuation at least 2 ( N + 1 ) 2 − ( N + 1 ) 2(N+1)^{2}-(N+1) 2 ( N + 1 ) 2 − ( N + 1 ) . Summing (40) therefore gives
V ( Q η ) = ( 1 − η ) V ( η ) , V ( 0 ) = 1. V(Q\eta)=(1-\eta)V(\eta),\qquad V(0)=1. V ( Q η ) = ( 1 − η ) V ( η ) , V ( 0 ) = 1.
These equations determine V V V uniquely in Q ( Q ) [ [ η ] ] \mathbb{Q}(Q)[[\eta]] Q ( Q ) [[ η ]] : its coefficient of η j \eta^j η j is 1 / ( Q ; Q ) j 1/(Q;Q)_j 1/ ( Q ; Q ) j . The product ( η ; Q ) ∞ − 1 (\eta;Q)_{\infty}^{-1} ( η ; Q ) ∞ − 1 , interpreted coefficientwise Q Q Q -adically, satisfies the same equations and initial value, proving (38) . The proof is uniform in the statistic variable z z z .
For the quantum-plane certificate, put S = H d z − H o z S=H^{z}_{d}-H^{z}_{o} S = H d z − H o z and let D ( x , R ) D(x,R) D ( x , R ) be its normal symbol. Set
L d = ( 1 0 0 − q Υ 0 1 − q Υ 0 − q R 0 1 − q R 0 0 − z q R 0 1 − z q R ) , L = diag ( L d , 1 − z q Υ − q 2 R 2 ) , L_d=\begin{pmatrix}
1&0&0&-q\Upsilon\\
0&1&-q\Upsilon&0\\
-qR&0&1-qR&0\\
0&-zqR&0&1-zqR
\end{pmatrix},\qquad L=\operatorname{diag}(L_d,1-zq\Upsilon-q^{2}R^{2}), L d = 1 0 − q R 0 0 1 0 − z q R 0 − q Υ 1 − q R 0 − q Υ 0 0 1 − z q R , L = diag ( L d , 1 − z q Υ − q 2 R 2 ) ,
f = ( 0 , 0 , q , z q , q ) T , d = ( 0 , 0 , 0 , R , − z R ) . f=(0,0,q,zq,q)^{\mathsf T},\qquad d=(0,0,0,R,-zR). f = ( 0 , 0 , q , z q , q ) T , d = ( 0 , 0 , 0 , R , − z R ) .
The first four states are: after a U U U ending a B 1 B_1 B 1 block, after a U U U ending a B 0 B_0 B 0 block, inside a B 0 B_0 B 0 block, inside a B 1 B_1 B 1 block. The normalized solution of L m = f Lm=f L m = f gives the marginal C = d m C=dm C = d m of S S S ; its initial data at R = 0 R=0 R = 0 are
m ( 0 ) = ( z q 2 , q 2 , q , z q , q / ( 1 − z q ) ) T , m(0)=(zq^{2},q^{2},q,zq,q/(1-zq))^{\mathsf T}, m ( 0 ) = ( z q 2 , q 2 , q , z q , q / ( 1 − z q ) ) T ,
and all subsequent coefficients are fixed by L L L , whose constant-coefficient matrix at every R R R -degree is invertible over Q ( q , z ) \mathbb{Q}(q,z) Q ( q , z ) . The marked system is obtained by replacing Υ \Upsilon Υ by x Υ x\Upsilon x Υ in L L L and f f f by x f xf x f ; the output row d d d is unchanged. The general construction gives G b G_b G b of size 7 and G a G_a G a of size 12.
The constant terms 1 cancel in S S S , all other accepted words lie in R { R , U } ∗ U R\{R,U\}^{*}U R { R , U } ∗ U , and the statistic weights lie in Q [ z ] \mathbb{Q}[z] Q [ z ] , so Lemma 3.5 and Corollary 3.6 apply unchanged: by (38) , ν ( a i , j ) ≥ i + j \nu(a_{i,j})\ge i+j ν ( a i , j ) ≥ i + j and ν ( b n ) ≥ n + 1 \nu(b_n)\ge n+1 ν ( b n ) ≥ n + 1 . This is a convergent holonomic certificate for every coefficient identity in (38) .
Residues of parts at prescribed indices
For a partition λ 1 > ⋯ > λ ℓ > 0 \lambda_1 > \cdots> \lambda_\ell> 0 λ 1 > ⋯ > λ ℓ > 0 into distinct parts, let o ( λ ) o(\lambda) o ( λ ) and e ( λ ) e(\lambda) e ( λ ) count its odd parts at odd and even indices, respectively. Berkovich and Uncu proved that these two statistics have the same joint distribution as the numbers of parts congruent to 1 and 3 modulo 4, respectively [4 Theorem 1.1] . Thus
∑ λ d i s t i n c t u o ( λ ) v e ( λ ) q ∣ λ ∣ = ( − q 2 ; q 2 ) ∞ ( − u q ; q 4 ) ∞ ( − v q 3 ; q 4 ) ∞ . (41) \sum_{\lambda\ \mathrm{distinct}} u^{o(\lambda)}v^{e(\lambda)}q^{\lvert\lambda\rvert}=(-q^2;q^2)_\infty(-uq;q^4)_\infty(-vq^3;q^4)_\infty. \tag*{(41)} λ distinct ∑ u o ( λ ) v e ( λ ) q ∣ λ ∣ = ( − q 2 ; q 2 ) ∞ ( − u q ; q 4 ) ∞ ( − v q 3 ; q 4 ) ∞ . (41)
Both sides count partitions into distinct parts. The empty partition contributes 1 on each side and cancels in their difference. We give the automaton kernels, and certify convergence by a recurrence for the bounded refinement of [4 Theorem 2.1] .
Automata and holonomic kernels. The boundary word reads parts in increasing order. To recover their indices in decreasing order, the first automaton records the parity of the number of U U U ’s still to be read. Its eight interior states are ( h , k , r ) ∈ { 0 , 1 } 3 (h,k,r) \in\{0,1\}^3 ( h , k , r ) ∈ { 0 , 1 } 3 : h h h is the current part size modulo 2, k k k is this remaining parity, and r = 1 r=1 r = 1 means that an R R R has occurred since the preceding U U U . Its transitions and weights are
Table 1.
w ( h , k ) = { 1 , h = 0 , u , h = 1 , k = 1 , v , h = 1 , k = 0. w(h,k)=
\begin{cases}
1, & h=0,\\
u, & h=1,\ k=1,\\
v, & h=1,\ k=0.
\end{cases} w ( h , k ) = ⎩ ⎨ ⎧ 1 , u , v , h = 0 , h = 1 , k = 1 , h = 1 , k = 0.
The initial R R R enters either ( 1 , 0 , 1 ) (1,0,1) ( 1 , 0 , 1 ) or ( 1 , 1 , 1 ) (1,1,1) ( 1 , 1 , 1 ) ; a final U U U leaves ( h , 1 , 1 ) (h,1,1) ( h , 1 , 1 ) with weight w ( h , 1 ) w(h,1) w ( h , 1 ) . Exactly one initial choice accepts any given boundary word. The flag r r r forbids consecutive U U U ’s, hence repeated parts.
The second automaton has eight states ( h , r ) ∈ Z / 4 Z × { 0 , 1 } (h,r) \in\mathbb{Z}/4\mathbb{Z} \times\{0,1\} ( h , r ) ∈ Z /4 Z × { 0 , 1 } and transitions
Table 2.
w ~ ( h ) = { u , h = 1 , v , h = 3 , 1 , h = 0 , 2. \widetilde{w}(h)=
\begin{cases}
u, & h=1,\\
v, & h=3,\\
1, & h=0,2.
\end{cases} w ( h ) = ⎩ ⎨ ⎧ u , v , 1 , h = 1 , h = 3 , h = 0 , 2.
Its initial R R R enters ( 1 , 1 ) (1,1) ( 1 , 1 ) ; its final U U U leaves any ( h , 1 ) (h,1) ( h , 1 ) with weight w ~ ( h ) \widetilde{w}(h) w ( h ) .
Let A R , A U A_R,A_U A R , A U be the incoming transition matrices of their disjoint union: an edge j → i j \to i j → i contributes its weight to entry ( i , j ) (i,j) ( i , j ) . Let r 0 r_0 r 0 be the column recording the three initial entries just specified, and d d d the row of final weights, positive for the first automaton and negative for the second. These tables give the marked kernel
( I − q R A R − q x Υ A U ) m x = x r 0 , D ( x , R ) = q R d m x . (42) (\mathbb{I}-qR A_R-qx\Upsilon A_U)m_x=xr_0,\qquad D(x,R)=qRdm_x. \tag*{(42)} ( I − q R A R − q x Υ A U ) m x = x r 0 , D ( x , R ) = q R d m x . (42)
The matrix is congruent to I \mathbb{I} I modulo q q q , which selects its solution in Z [ u , v , x , R ] [ [ q ] ] \mathbb{Z}[u,v,x,R][[q]] Z [ u , v , x , R ] [[ q ]] , and (7) gives the canonical holonomic a , b a,b a , b . Equivalently, the following is a kernel for the preimage ϕ \phi ϕ of Theorem A(3d):
( I − R A R − q Υ A U ) m ~ = r 0 , ( q − R ) ϕ = R d m ~ , b = q Υ ϕ . (43) \begin{aligned}
(\mathbb{I}-R A_R-q\Upsilon A_U)\widetilde{m}&=r_0,\\
(q-R)\phi&=Rd\widetilde{m},\qquad b=q\Upsilon\phi. \tag*{(43)}
\end{aligned} ( I − R A R − q Υ A U ) m ( q − R ) ϕ = r 0 , = R d m , b = q Υ ϕ . (43)
Here m ~ = m 1 ( R / q ) \widetilde{m}=m_1(R/q) m = m 1 ( R / q ) is selected coefficientwise in R R R ; the constant matrix at degree n n n is I − q n + 1 A U \mathbb{I}-q^{n+1}A_U I − q n + 1 A U , with inverse I + q n + 1 A U \mathbb{I}+q^{n+1}A_U I + q n + 1 A U since A U 2 = 0 A_U^2=0 A U 2 = 0 . The second equation forces ϕ ( 0 ) = 0 \phi(0)=0 ϕ ( 0 ) = 0 . The argument below proves that ϕ ∈ Z [ u , v ] [ [ q , R ] ] \phi\in\mathbb{Z}[u,v][[q,R]] ϕ ∈ Z [ u , v ] [[ q , R ]] .
A recurrence certificate. Write E N ( i , j ) E_N(i,j) E N ( i , j ) for the coefficient of u i v j u^iv^j u i v j on the left of (41) with largest part at most 2 N 2N 2 N . Choosing which of 2 N , 2 N − 1 2N,2N-1 2 N , 2 N − 1 occur gives
E N ( i , j ) = E N − 1 ( i , j ) + q 2 N E N − 1 ( i , j ) + q 2 N − 1 E N − 1 ( j , i − 1 ) + q 4 N − 1 E N − 1 ( i , j − 1 ) . (44) \begin{aligned}
E_N(i,j)&=E_{N-1}(i,j)+q^{2N}E_{N-1}(i,j)\\
&\quad+q^{2N-1}E_{N-1}(j,i-1)+q^{4N-1}E_{N-1}(i,j-1). \tag*{(44)}
\end{aligned} E N ( i , j ) = E N − 1 ( i , j ) + q 2 N E N − 1 ( i , j ) + q 2 N − 1 E N − 1 ( j , i − 1 ) + q 4 N − 1 E N − 1 ( i , j − 1 ) . (44)
The initial data are E 0 ( i , j ) = δ i , 0 δ j , 0 E_0(i,j)=\delta_{i,0}\delta_{j,0} E 0 ( i , j ) = δ i , 0 δ j , 0 ; negative indices give zero. Its proposed solution is
T N ( i , j ) = q 2 i 2 − i + 2 j 2 + j ( q 4 ; q 4 ) N ( q 4 ; q 4 ) i ( q 4 ; q 4 ) j ( q 2 ; q 2 ) N − i − j ( i , j ≥ 0 , i + j ≤ N ) . (45) T_N(i,j)=\frac{q^{2i^2-i+2j^2+j}(q^4;q^4)_N}{(q^4;q^4)_i(q^4;q^4)_j(q^2;q^2)_{N-i-j}}\qquad(i,j\ge0,\ i+j\le N). \tag*{(45)} T N ( i , j ) = ( q 4 ; q 4 ) i ( q 4 ; q 4 ) j ( q 2 ; q 2 ) N − i − j q 2 i 2 − i + 2 j 2 + j ( q 4 ; q 4 ) N ( i , j ≥ 0 , i + j ≤ N ) . (45)
and zero elsewhere. To verify it, set X = q 2 N X=q^{2N} X = q 2 N , Y = q 2 i Y=q^{2i} Y = q 2 i , Z = q 2 j Z=q^{2j} Z = q 2 j , and divide each term on the right of (44) by T N ( i , j ) T_N(i,j) T N ( i , j ) :
Table 3.
Their sum is 1 by the polynomial identity
( Y Z − X ) ( Z + X Y ) + X Z ( 1 − Y 2 ) + X 2 Y ( 1 − Z 2 ) = Y Z 2 ( 1 − X 2 ) . (46) (YZ-X)(Z+XY)+XZ(1-Y^2)+X^2Y(1-Z^2)=YZ^2(1-X^2). \tag*{(46)} ( Y Z − X ) ( Z + X Y ) + X Z ( 1 − Y 2 ) + X 2 Y ( 1 − Z 2 ) = Y Z 2 ( 1 − X 2 ) . (46)
On the support boundary i + j = N i+j=N i + j = N , the factor Y Z − X YZ-X Y Z − X makes the first two terms zero; at i = 0 i=0 i = 0 or j = 0 j=0 j = 0 , the factors 1 − Y 2 1-Y^2 1 − Y 2 or 1 − Z 2 1-Z^2 1 − Z 2 give the required zeros. Outside the support all four terms vanish. Together with T 0 ( 0 , 0 ) = 1 T_0(0,0)=1 T 0 ( 0 , 0 ) = 1 , this proves E N = T N E_N=T_N E N = T N .
Convergence data. Put σ i , j = 2 i 2 − i + 2 j 2 + j \sigma_{i,j}=2i^2-i+2j^2+j σ i , j = 2 i 2 − i + 2 j 2 + j and
G i , j = q σ i , j ( − q 2 ; q 2 ) ∞ ( q 4 ; q 4 ) i ( q 4 ; q 4 ) j . G_{i,j}=\frac{q^{\sigma_{i,j}}(-q^2;q^2)_\infty}{(q^4;q^4)_i(q^4;q^4)_j}. G i , j = ( q 4 ; q 4 ) i ( q 4 ; q 4 ) j q σ i , j ( − q 2 ; q 2 ) ∞ .
These are the coefficients of u i v j u^iv^j u i v j in the product in (41) : the expansion of ( − z ; q 4 ) ∞ (-z;q^4)_\infty ( − z ; q 4 ) ∞ is selected by the first-order equation g ( z ) = ( 1 + z ) g ( q 4 z ) g(z)=(1+z)g(q^4z) g ( z ) = ( 1 + z ) g ( q 4 z ) and g ( 0 ) = 1 g(0)=1 g ( 0 ) = 1 . For i + j ≤ N i+j\leq N i + j ≤ N , the quotient is
T N ( i , j ) G i , j = ( q 2 ( N − i − j + 1 ) ; q 2 ) ∞ ( q 4 N + 4 ; q 4 ) ∞ ≡ 1 ( m o d q 2 ( N − i − j + 1 ) ) . \frac{T_N(i,j)}{G_{i,j}}=\frac{(q^{2(N-i-j+1)};q^2)_\infty}{(q^{4N+4};q^4)_\infty}\equiv1\pmod{q^{2(N-i-j+1)}}. G i , j T N ( i , j ) = ( q 4 N + 4 ; q 4 ) ∞ ( q 2 ( N − i − j + 1 ) ; q 2 ) ∞ ≡ 1 ( mod q 2 ( N − i − j + 1 ) ) .
The valuation check is
σ i , j − 2 ( i + j ) + 1 = ( i − 1 ) ( 2 i − 1 ) + j ( 2 j − 1 ) ≥ 0 ( i , j ∈ Z ≥ 0 ) , (47) \sigma_{i,j}-2(i+j)+1=(i-1)(2i-1)+j(2j-1)\geq0\qquad(i,j\in\mathbb{Z}_{\geq0}), \tag*{(47)} σ i , j − 2 ( i + j ) + 1 = ( i − 1 ) ( 2 i − 1 ) + j ( 2 j − 1 ) ≥ 0 ( i , j ∈ Z ≥ 0 ) , (47)
each summand being a product of two integers of the same sign. Consequently ν ( T N ( i , j ) − G i , j ) ≥ σ i , j + 2 ( N − i − j + 1 ) ≥ 2 N + 1 \nu(T_N(i,j)-G_{i,j})\geq\sigma_{i,j}+2(N-i-j+1)\geq2N+1 ν ( T N ( i , j ) − G i , j ) ≥ σ i , j + 2 ( N − i − j + 1 ) ≥ 2 N + 1 . When i + j > N i+j>N i + j > N , the same inequality gives ν ( G i , j ) = σ i , j ≥ 2 N + 1 \nu(G_{i,j})=\sigma_{i,j}\geq2N+1 ν ( G i , j ) = σ i , j ≥ 2 N + 1 . The estimate is therefore uniform in the marking variables:
ν ( ∑ i , j u i v j E N ( i , j ) − ∑ i , j u i v j G i , j ) ≥ 2 N + 1. (48) \nu\left(\sum_{i,j}u^iv^jE_N(i,j)-\sum_{i,j}u^iv^jG_{i,j}\right)\geq2N+1. \tag*{(48)} ν ( i , j ∑ u i v j E N ( i , j ) − i , j ∑ u i v j G i , j ) ≥ 2 N + 1. (48)
By the remark after (11) , b n b_n b n is the difference between the two weighted counts with largest part at most n n n . The second count, truncated at 2 N 2N 2 N , differs from its infinite product by terms of valuation at least 2 N + 1 2N+1 2 N + 1 . Hence (48) gives $\nu(b_{2N})\geq 2N+1$. Also $\nu(b_{2N+1}-b_{2N})\geq 2N+1$, since the newly added partitions contain a part $2N+1$. Thus $b_n\to 0$, and Theorem A(3)`, whose proof is unchanged for weights in Z [ u , v ] \mathbb{Z}[u,v] Z [ u , v ] , gives the sharper bound and the integral preimage:
ν ( b n ) ≥ n + 1 , b = q Υ ϕ , ϕ ∈ Z [ u , v ] [ [ q , R ] ] . \nu(b_n)\geq n+1,\qquad b=q^\Upsilon\phi,\quad\phi\in\mathbb{Z}[u,v][[q,R]]. ν ( b n ) ≥ n + 1 , b = q Υ ϕ , ϕ ∈ Z [ u , v ] [[ q , R ]] .
The same theorem supplies ν ( a i , j ) ≥ i + j \nu(a_{i,j})\geq i+j ν ( a i , j ) ≥ i + j . This completes the convergent holonomic certificate. Its data are the two transition tables, the rational ratios, the support and initial value, and (47) ; its main check is (46) .
Larger moduli. Fu and Zeng extended (41) as follows. Fix L ≥ 2 L\geq2 L ≥ 2 and require distinct parts, with at most one part from each block
{ L j + 1 , … , L j + L − 1 } , j ≥ 0. \{Lj+1,\ldots,Lj+L-1\},\qquad j\geq0. { L j + 1 , … , L j + L − 1 } , j ≥ 0.
Let o r , e r o_r,e_r o r , e r count the parts congruent to r r r (mod L L L ) at odd and even indices. With z = ∑ r = 1 L − 1 u r q r z=\sum_{r=1}^{L-1}u_rq^r z = ∑ r = 1 L − 1 u r q r and w = ∑ r = 1 L − 1 v r q r w=\sum_{r=1}^{L-1}v_rq^r w = ∑ r = 1 L − 1 v r q r , the one-variable case of [9 Theorem 2.6] reads
∑ λ q ∣ λ ∣ ∏ r = 1 L − 1 u r o r ( λ ) v r e r ( λ ) = ( − q L ; q L ) ∞ ( − z ; q 2 L ) ∞ ( − q L w ; q 2 L ) ∞ . (49) \sum_{\lambda}q^{|\lambda|}\prod_{r=1}^{L-1}u_r^{o_r(\lambda)}v_r^{e_r(\lambda)}=(-q^L;q^L)_\infty(-z;q^{2L})_\infty(-q^Lw;q^{2L})_\infty. \tag*{(49)} λ ∑ q ∣ λ ∣ r = 1 ∏ L − 1 u r o r ( λ ) v r e r ( λ ) = ( − q L ; q L ) ∞ ( − z ; q 2 L ) ∞ ( − q L w ; q 2 L ) ∞ . (49)
For L = 2 L=2 L = 2 the block condition adds nothing to distinctness, so this is exactly (41) .
The certificate extends to all L L L . For largest part at most L N LN L N , the generating polynomial, expressed in z , w z,w z , w , satisfies
F N ( z , w ) = ( 1 + q L ( 2 N − 1 ) w ) F N − 1 ( z , w ) + q L ( N − 1 ) ( q L + z ) F N − 1 ( w , z ) , F 0 = 1. (50) F_N(z,w)=(1+q^{L(2N-1)}w)F_{N-1}(z,w)+q^{L(N-1)}(q^L+z)F_{N-1}(w,z),\qquad F_0=1. \tag*{(50)} F N ( z , w ) = ( 1 + q L ( 2 N − 1 ) w ) F N − 1 ( z , w ) + q L ( N − 1 ) ( q L + z ) F N − 1 ( w , z ) , F 0 = 1. (50)
There are four possibilities in the last block: neither part, its multiple of L L L , one nonmultiple, or both. They give these four terms, with the swap z ↔ w z\leftrightarrow w z ↔ w whenever one part is added. The coefficient of z i w j z^iw^j z i w j is certified by
T N ( L ) ( i , j ) = q L i ( i − 1 ) + L j 2 ( q 2 L ; q 2 L ) N ( q 2 L ; q 2 L ) i ( q 2 L ; q 2 L ) j ( q L ; q L ) N − i − j , (51) T_N^{(L)}(i,j)=\frac{q^{Li(i-1)+Lj^2}(q^{2L};q^{2L})_N}{(q^{2L};q^{2L})_i(q^{2L};q^{2L})_j(q^L;q^L)_{N-i-j}}, \tag*{(51)} T N ( L ) ( i , j ) = ( q 2 L ; q 2 L ) i ( q 2 L ; q 2 L ) j ( q L ; q L ) N − i − j q L i ( i − 1 ) + L j 2 ( q 2 L ; q 2 L ) N , (51)
on the same support. With X = q L N X=q^{LN} X = q L N , Y = q L i Y=q^{Li} Y = q L i , Z = q L j Z=q^{Lj} Z = q L j , the four ratios are exactly those in the preceding table. Thus (46) verifies (50) for every L L L . The limit of (51) , summed against z i w j z^iw^j z i w j , is the product in (49) , by the same first-order product equations.
Since ν ( z ) , ν ( w ) ≥ 1 \nu(z),\nu(w)\geq1 ν ( z ) , ν ( w ) ≥ 1 , each summand has valuation at least σ i , j ( L ) = L i ( i − 1 ) + L j 2 + i + j \sigma_{i,j}^{(L)}=Li(i-1)+Lj^2+i+j σ i , j ( L ) = L i ( i − 1 ) + L j 2 + i + j . The convergence check becomes
σ i , j ( L ) − L ( i + j ) + L = L ( i − 1 ) 2 + L j ( j − 1 ) + i + j ≥ 1. \sigma_{i,j}^{(L)}-L(i+j)+L=L(i-1)^2+Lj(j-1)+i+j\geq1. σ i , j ( L ) − L ( i + j ) + L = L ( i − 1 ) 2 + L j ( j − 1 ) + i + j ≥ 1.
With the tail quotient
( q L ( N − i − j + 1 ) ; q L ) ∞ ( q 2 L ( N + 1 ) ; q 2 L ) ∞ , \frac{(q^{L(N-i-j+1)};q^L)_\infty}{(q^{2L(N+1)};q^{2L})_\infty}, ( q 2 L ( N + 1 ) ; q 2 L ) ∞ ( q L ( N − i − j + 1 ) ; q L ) ∞ ,
this gives the uniform error ν ( F N − F ∞ ) ≥ L N + 1 \nu(F_N-F_\infty)\geq LN+1 ν ( F N − F ∞ ) ≥ L N + 1 , and the same argument proves convergence of the canonical boundary whenever the marking weights have nonnegative valuation.
Both sides have finite automata: on the index side, the first automaton, with h h h read modulo L L L , needs one extra bit, which records whether a part not divisible by L L L has been used in the current block and is reset when an R R R crosses a multiple of L L L ; on the product side, each factor 1 + z q 2 L j 1+zq^{2Lj} 1 + z q 2 L j or 1 + w q 2 L j + L 1+wq^{2Lj+L} 1 + w q 2 L j + L selects at most one part from its block.
For example, retain only residues 0 , a 0,a 0 , a at odd indices and 0 , c 0,c 0 , c at even indices, where 1 ≤ a , c < L 1\leq a,c<L 1 ≤ a , c < L . Within the stated block condition, these partitions are equinumerous with distinct partitions whose parts lie in
0 , a , L , L + c ( m o d 2 L ) . 0,\ a,\ L,\ L+c\pmod{2L}. 0 , a , L , L + c ( mod 2 L ) .
Taking L = 4 L=4 L = 4 , a = 1 a=1 a = 1 , c = 3 c=3 c = 3 gives an index-dependent mod-4 condition equal to the residue classes 0 , 1 , 4 , 7 ( m o d 8 ) 0,1,4,7\pmod8 0 , 1 , 4 , 7 ( mod 8 ) .
Arbitrary index periods. The automata do not depend on a product formula. For L , M ≥ 1 L,M\geq1 L , M ≥ 1 , give a part of residue h h h mod L L L at index k k k mod M M M the weight z k , h z_{k,h} z k , h ; weights 0 0 0 or 1 1 1 impose arbitrary index-dependent residue restrictions. The states
( h , k , r ) ∈ Z / L Z × Z / M Z × { 0 , 1 } (h,k,r)\in\mathbb{Z}/L\mathbb{Z}\times\mathbb{Z}/M\mathbb{Z}\times\{0,1\} ( h , k , r ) ∈ Z / L Z × Z / M Z × { 0 , 1 }
have transitions
( h , k , r ) → R ( h + 1 , k , 1 ) , ( h , k , 1 ) → z k , h U ( h , k − 1 , 0 ) . (h,k,r)\xrightarrow{R}(h+1,k,1),\qquad(h,k,1)\xrightarrow{z_{k,h}U}(h,k-1,0). ( h , k , r ) R ( h + 1 , k , 1 ) , ( h , k , 1 ) z k , h U ( h , k − 1 , 0 ) .
The initial R R R enters ( 1 , k , 1 ) (1,k,1) ( 1 , k , 1 ) for every k k k , and the final U U U leaves ( h , 1 , 1 ) (h,1,1) ( h , 1 , 1 ) with weight z 1 , h z_{1,h} z 1 , h , all indices being read modulo their periods. As before, exactly one initial choice accepts each word, and the kernel (42) applies unchanged, with at most 2 L M 2LM 2 L M interior states.
Bounded counts also have a small kernel. Let v N ( k ) v_N^{(k)} v N ( k ) count distinct partitions with largest part at most N N N , assigning index k k k to the largest part and then increasing indices modulo M M M . Adding N N N gives the recurrence
v N ( k ) = v N − 1 ( k ) + q N z k , N m o d L v N − 1 ( k + 1 ) , v 0 ( k ) = 1. (52) v_N^{(k)}=v_{N-1}^{(k)}+q^N z_{k,N\ \mathrm{mod}\ L}v_{N-1}^{(k+1)},\qquad v_0^{(k)}=1. \tag*{(52)} v N ( k ) = v N − 1 ( k ) + q N z k , N mod L v N − 1 ( k + 1 ) , v 0 ( k ) = 1. (52)
In vector form put C h = diag ( z k , h ) C_h=\operatorname{diag}(z_{k,h}) C h = diag ( z k , h ) , and let S S S be the cyclic shift ( S v ) k = v k + 1 (Sv)_k=v_{k+1} ( S v ) k = v k + 1 . Expand the ordered matrix product
B ( Z ) = ( I + q L Z C 0 S ) ( I + q L − 1 Z C L − 1 S ) ⋯ ( I + q Z C 1 S ) = ∑ d = 0 L B d Z d , B(Z)=(\mathbb{I}+q^L ZC_0S)(\mathbb{I}+q^{L-1}ZC_{L-1}S)\cdots(\mathbb{I}+qZC_1S)=\sum_{d=0}^{L}B_dZ^d, B ( Z ) = ( I + q L Z C 0 S ) ( I + q L − 1 Z C L − 1 S ) ⋯ ( I + q Z C 1 S ) = d = 0 ∑ L B d Z d ,
where for L = 1 L=1 L = 1 there is just the factor I + q Z C 0 S \mathbb{I}+qZC_0S I + q Z C 0 S . Then f ( R ) = ∑ N ≥ 0 v L N R N f(R)=\sum_{N\geq0}v_LNR^N f ( R ) = ∑ N ≥ 0 v L N R N is selected by
( I − R ∑ d = 0 L B d R L d ) f = 1 , f ( 0 ) = 1 , (53) \left(\mathbb{I}-R\sum_{d=0}^{L}B_dR^{Ld}\right)f=1,\qquad f(0)=1, \tag*{(53)} ( I − R d = 0 ∑ L B d R L d ) f = 1 , f ( 0 ) = 1 , (53)
since v L ( N + 1 ) = B ( q L N ) v L N v_{L(N+1)}=B(q^{LN})v_{LN} v L ( N + 1 ) = B ( q L N ) v L N . This is an M M M -component holonomic kernel.
A matching product, and the convergence of the difference between two proposed classes, remain separate assertions. For L > 2 L>2 L > 2 the Fu–Zeng block condition cannot simply be dropped: without it, marking the parts ≡ 1 ( m o d 4 ) \equiv1\pmod4 ≡ 1 ( mod 4 ) by their index parity gives the coefficient 1 + v 1+v 1 + v at q 3 q^3 q 3 , from (3) and (2,1), whereas marking the parts ≡ 1 ( m o d 8 ) \equiv1\pmod8 ≡ 1 ( mod 8 ) by u u u and those ≡ 5 ( m o d 8 ) \equiv5\pmod8 ≡ 5 ( mod 8 ) by v v v gives 1 + u 1+u 1 + u .
An identity between five-state languages
Consider the two partition languages recognized by the automata of Figure 3 . They were found by enumerating the 14 614 080 minimal deterministic interior automata with at most five states (§ 3.8 ) and grouping them by their series to order q 50 q^{50} q 50 . Neither is defined by conditions of the form (5), for any M M M , W W W , W W W and initial vectors, even with unbounded multiplicities. Indeed, for n ≥ W n\geq W n ≥ W and k ≡ − 2 ( m o d M ) k\equiv-2\pmod M k ≡ − 2 ( mod M ) , the words R R U ( R U ) n U RRU(RU)^nU R R U ( R U ) n U and R R U R R U ( R U ) n + k U R R U RRURRU(RU)^{n+k}URRU R R U R R U ( R U ) n + k U R R U lie in L A \mathcal{L}_A L A but R R U ( R U ) n U R R U RRU(RU)^nURRU R R U ( R U ) n U R R U does not, although each window of the latter, with its height modulo M M M , is a window of one of the former; similarly, R ( R U ) n R U R(RU)^nRU R ( R U ) n R U and R R ( R U ) n + M − 1 U RR(RU)^{n+M-1}U R R ( R U ) n + M − 1 U lie in L B \mathcal{L}_B L B but R ( R U ) n U R(RU)^nU R ( R U ) n U does not. Yet both languages count partitions of each size equally. We give algebraic certificates of this equality.
Figure 3. The two automata for the identity. Each has one initial state s s s and one final state f f f ; the numbered states are interior states.
Write H A , H B H_A,H_B H A , H B for their path series. Reversing a boundary word and interchanging R , U R,U R , U conjugates its partition, so it preserves area. Let H B † H_B^{\dagger} H B † be the series of these conjugate partitions. We certify
S = H A − H B † = ( 1 − U ) a + b ( 1 − R ) . S=H_A-H_B^{\dagger}=(1-U)a+b(1-R). S = H A − H B † = ( 1 − U ) a + b ( 1 − R ) .
The normal symbol D ( x , R ) D(x,R) D ( x , R ) of S S S has the following kernel. Every term involving an unknown on the right has a factor q q q , so the equations determine unique solutions in Q [ x , R ] [ [ q ] ] \mathbb{Q}[x,R][[q]] Q [ x , R ] [[ q ]] :
m 1 = x + q x Υ m 3 , c 1 = q x Υ c 2 , m 2 = q R m 1 , c 2 = x + q R c 1 + q x Υ c 3 , m 3 = q x Υ ( m 2 + m 4 ) , c 3 = q R c 2 + q x Υ c 4 , m 4 = q R m 3 + q x Υ m 5 , c 4 = x + q R c 5 , m 5 = q R m 4 , c 5 = q x Υ c 3 . (54) \begin{aligned}
m_1 &= x+qx\Upsilon m_3, & c_1 &= qx\Upsilon c_2,\\
m_2 &= qRm_1, & c_2 &= x+qRc_1+qx\Upsilon c_3,\\
m_3 &= qx\Upsilon(m_2+m_4), & c_3 &= qRc_2+qx\Upsilon c_4,\\
m_4 &= qRm_3+qx\Upsilon m_5, & c_4 &= x+qRc_5,\\
m_5 &= qRm_4, & c_5 &= qx\Upsilon c_3. \tag*{(54)}
\end{aligned} m 1 m 2 m 3 m 4 m 5 = x + q x Υ m 3 , = q R m 1 , = q x Υ ( m 2 + m 4 ) , = q R m 3 + q x Υ m 5 , = q R m 4 , c 1 c 2 c 3 c 4 c 5 = q x Υ c 2 , = x + q R c 1 + q x Υ c 3 , = q R c 2 + q x Υ c 4 , = x + q R c 5 , = q x Υ c 3 . (54)
They are read directly from incoming edges: m m m uses the first automaton, and c c c uses the second with all edges reversed and labels interchanged, whose initial interior states are therefore 2 and 4, and whose final interior state is 1. The source terms x x x and the factors q q q account for the two outer letters, as in §3.5. Thus
D ( x , R ) = q R ( m 2 + m 3 + m 4 + m 5 − c 1 ) , a ( x , R ) = D ( x , R ) − D ( 1 , R ) 1 − x . (55) D(x,R)=qR(m_2+m_3+m_4+m_5-c_1),\qquad a(x,R)=\frac{D(x,R)-D(1,R)}{1-x}. \tag*{(55)} D ( x , R ) = q R ( m 2 + m 3 + m 4 + m 5 − c 1 ) , a ( x , R ) = 1 − x D ( x , R ) − D ( 1 , R ) . (55)
Equations (54)–(55) give a holonomic definition of a a a .
For b b b , in the operator algebra Υ R = q R Υ \Upsilon R=qR\Upsilon Υ R = q R Υ , put
P = 1 − q 2 ( 1 + q ) R Υ − q 3 R Υ R 2 + q 8 R 2 Υ 3 , Q = − 1 − ( 1 + q ) R − R 2 + q ( q 2 − 1 ) R Υ + q 5 R 2 Υ 2 . (56) P=1-q^2(1+q)R\Upsilon-q^3R\Upsilon R^2+q^8R^2\Upsilon^3,
\qquad
Q=-1-(1+q)R-R^2+q(q^2-1)R\Upsilon+q^5R^2\Upsilon^2. \tag*{(56)} P = 1 − q 2 ( 1 + q ) R Υ − q 3 R Υ R 2 + q 8 R 2 Υ 3 , Q = − 1 − ( 1 + q ) R − R 2 + q ( q 2 − 1 ) R Υ + q 5 R 2 Υ 2 . (56)
Define
P ψ = q 2 R , ψ ( 0 ) = 0 , ϕ = Q ψ − R , b = q Υ ϕ . (57) P\psi=q^2R,\qquad\psi(0)=0,\qquad\phi=Q\psi-R,\qquad b=q\Upsilon\phi. \tag*{(57)} P ψ = q 2 R , ψ ( 0 ) = 0 , ϕ = Q ψ − R , b = q Υ ϕ . (57)
Since P ≡ 1 ( m o d R ) P\equiv1\pmod R P ≡ 1 ( mod R ) , the first equation determines ψ \psi ψ uniquely; the displayed initial value is forced by it.
We now bind these data to the automata; from now on, put x = 1 x=1 x = 1 in (54). Elimination gives
q R ( m 2 + m 3 + m 4 + m 5 ) = Θ ψ , q R c 1 = ( Υ − q 5 R Υ 3 ) c 3 − q , Ω c 3 = q + ( q − q 2 ) R . (58) \begin{aligned}
qR(m_2+m_3+m_4+m_5)&=\Theta\psi,\\
qRc_1&=(\Upsilon-q^5R\Upsilon^3)c_3-q,\\
\Omega c_3&=q+(q-q^2)R. \tag*{(58)}
\end{aligned} q R ( m 2 + m 3 + m 4 + m 5 ) q R c 1 Ω c 3 = Θ ψ , = ( Υ − q 5 R Υ 3 ) c 3 − q , = q + ( q − q 2 ) R . (58)
where
Θ = R + q R ( 1 − q 2 R + q 2 R 2 ) Υ − q 5 R 2 Υ 2 , Ω = 1 − q ( 1 + q ) R Υ − q 4 R Υ 2 + q 6 R 2 Υ 3 . \Theta=R+qR(1-q^2R+q^2R^2)\Upsilon-q^5R^2\Upsilon^2,
\qquad
\Omega=1-q(1+q)R\Upsilon-q^4R\Upsilon^2+q^6R^2\Upsilon^3. Θ = R + q R ( 1 − q 2 R + q 2 R 2 ) Υ − q 5 R 2 Υ 2 , Ω = 1 − q ( 1 + q ) R Υ − q 4 R Υ 2 + q 6 R 2 Υ 3 .
For the first equality, note that the first five equations are satisfied by
q 2 R ( m 1 m 2 m 3 m 4 m 5 ) = ( 1 − q 2 ( 1 + q ) R Υ q R − q 3 ( 1 + q ) R 2 Υ q 2 R Υ − q 6 R 2 Υ 2 q 3 R 2 Υ q 4 R 3 Υ ) ψ . q^2R
\begin{pmatrix}
m_1\\
m_2\\
m_3\\
m_4\\
m_5
\end{pmatrix}
=
\begin{pmatrix}
1-q^2(1+q)R\Upsilon\\
qR-q^3(1+q)R^2\Upsilon\\
q^2R\Upsilon-q^6R^2\Upsilon^2\\
q^3R^2\Upsilon\\
q^4R^3\Upsilon
\end{pmatrix}
\psi. q 2 R m 1 m 2 m 3 m 4 m 5 = 1 − q 2 ( 1 + q ) R Υ q R − q 3 ( 1 + q ) R 2 Υ q 2 R Υ − q 6 R 2 Υ 2 q 3 R 2 Υ q 4 R 3 Υ ψ .
substituting, the first equation of (54) becomes P ψ = q 2 R P\psi=q^2R P ψ = q 2 R and the other four hold identically. Each entry on the right side is divisible by R R R because ψ ( 0 ) = 0 \psi(0)=0 ψ ( 0 ) = 0 .
Put
N = q + q ( 1 + q ) R + q 2 ( 1 + R ) Υ − q 6 R ( 1 + R ) Υ 2 , J = q + q ( 1 + q ) R + q 2 ( 1 + q 2 R ) Υ − q 3 R ( 1 + q 3 R ) Υ 2 . N=q+q(1+q)R+q^2(1+R)\Upsilon-q^6R(1+R)\Upsilon^2,
\qquad
J=q+q(1+q)R+q^2(1+q^2R)\Upsilon-q^3R(1+q^3R)\Upsilon^2. N = q + q ( 1 + q ) R + q 2 ( 1 + R ) Υ − q 6 R ( 1 + R ) Υ 2 , J = q + q ( 1 + q ) R + q 2 ( 1 + q 2 R ) Υ − q 3 R ( 1 + q 3 R ) Υ 2 .
The two operator checks
Ω N = J P , Ω ( q ( 1 + R ) ) + J ( q 2 R ) = q + ( q − q 2 ) R (59) \Omega N=JP,\qquad\Omega\bigl(q(1+R)\bigr)+J(q^2R)=q+(q-q^2)R \tag*{(59)} Ω N = J P , Ω ( q ( 1 + R ) ) + J ( q 2 R ) = q + ( q − q 2 ) R (59)
show, by uniqueness for Ω ≡ 1 ( m o d R ) \Omega\equiv1\pmod R Ω ≡ 1 ( mod R ) , that c 3 = q ( 1 + R ) + N ψ c_3=q(1+R)+N\psi c 3 = q ( 1 + R ) + N ψ . Finally, direct multiplication gives
Θ − ( Υ − q 5 R Υ 3 ) N − ( 1 − R ) q Υ Q = ( R − q 2 ( 1 + q 3 R ) R 2 ) P , \Theta-(\Upsilon-q^5R\Upsilon^3)N-(1-R)q\Upsilon Q
=(R-q^2(1+q^3R)R^2)P, Θ − ( Υ − q 5 R Υ 3 ) N − ( 1 − R ) q Υ Q = ( R − q 2 ( 1 + q 3 R ) R 2 ) P ,
q − ( Υ − q 5 R Υ 3 ) ( q ( 1 + R ) ) + ( 1 − R ) q 2 R + ( R − q 2 ( 1 + q 3 R ) R 2 ) ( q 2 R ) = 0. (60) q-(\Upsilon-q^5R\Upsilon^3)\bigl(q(1+R)\bigr)+(1-R)q^2R+(R-q^2(1+q^3R)R^2)(q^2R)=0. \tag*{(60)} q − ( Υ − q 5 R Υ 3 ) ( q ( 1 + R ) ) + ( 1 − R ) q 2 R + ( R − q 2 ( 1 + q 3 R ) R 2 ) ( q 2 R ) = 0. (60)
Applying these identities to ψ \psi ψ and using (58) proves
D ( 1 , R ) = ( 1 − R ) q Υ ( Q ψ − R ) = ( 1 − R ) b . D(1,R)=(1-R)q\Upsilon(Q\psi-R)=(1-R)b. D ( 1 , R ) = ( 1 − R ) q Υ ( Q ψ − R ) = ( 1 − R ) b .
Together with (55), this verifies the boundary identity. All checks are polynomial calculations using Υ R = q R Υ \Upsilon R=qR\Upsilon Υ R = q R Υ ; the operators N N N , J J J are verification data, not additional unknown series.
It remains to certify convergence. Every term of 1 − P 1-P 1 − P has positive R R R -degree and coefficients in Z [ q ] \mathbb{Z}[q] Z [ q ] , so P P P is an integral kernel, and induction on the R R R -degree gives
ψ ∈ Z [ q ] [ [ R ] ] , ϕ = Q ψ − R ∈ Z [ q ] [ [ R ] ] ⊂ Q [ [ q , R ] ] . \psi\in\mathbb{Z}[q][[R]],\qquad\phi=Q\psi-R\in\mathbb{Z}[q][[R]]\subset\mathbb{Q}[[q,R]]. ψ ∈ Z [ q ] [[ R ]] , ϕ = Q ψ − R ∈ Z [ q ] [[ R ]] ⊂ Q [[ q , R ]] .
Together with the checks above, this is a direct certificate: ϕ \phi ϕ is the integral preimage of Theorem A(3d), since b = q Υ ϕ b=q\Upsilon\phi b = q Υ ϕ . Hence
ν ( [ R n ] b ) = n + 1 + ν ( [ R n ] ϕ ) ≥ n + 1. \nu([R^n]b)=n+1+\nu([R^n]\phi)\ge n+1. ν ([ R n ] b ) = n + 1 + ν ([ R n ] ϕ ) ≥ n + 1.
Theorem A(2) gives ν ( a i , j ) ≥ i + j \nu(a_{i,j})\ge i+j ν ( a i , j ) ≥ i + j for the canonical a a a in (55). Both boundaries consequently converge q q q -adically, proving the partition identity.
Implementation
The accompanying Julia package [3 ] , built on Oscar [17 ] , enumerates automata, computes their series in A A A and their marginal systems (9), and searches for certificates as in §3.8 and Appendix A. Its directory paper lists the automata of the figures as transition tables, and contains one script for each of §1 and §§4.1–4.6. The scripts check exactly the polynomial, rational and operator identities of the certificates, the sizes and kernels of the searches that find them, and, in §§1, 4.3 and 4.6, that the marginal systems compiled from the automata are those displayed at x = 1 x=1 x = 1 . They also compare, to a fixed order, the series of the automata with the stated sums and products, and the coefficients b n b_n b n with the differences of bounded counts given after (11); these comparisons are consistency checks, on which no proof relies. With Julia 1.13.1 and Oscar 1.8.2 the checks run in under a minute.
Further convergence tests
The following constructions give further sufficient tests for the convergence of b b b , by finite checks on holonomic data; the examples of §4 do not need them.
A degree bound for auxiliary normalized systems
A degree bound for auxiliary normalized systems. For L m = f Lm=f L m = f with L ≡ I ( m o d q ) L\equiv I\pmod q L ≡ I ( mod q ) and C = d m C=dm C = d m as in (9), choose
κ ≥ max ( 1 , deg R L ) , κ 0 ≥ deg R f + deg R d . \kappa\ge\max(1,\deg_R L),\qquad\kappa_0\ge\deg_R f+\deg_R d. κ ≥ max ( 1 , deg R L ) , κ 0 ≥ deg R f + deg R d .
Each factor I − L I-L I − L raises q q q -valuation by at least one and R R R -degree by at most κ \kappa κ . Thus
deg R ( [ q t ] C ) ≤ κ 0 + κ t , ν ( [ R n ] C ) ≥ max ( 0 , ⌈ n − κ 0 κ ⌉ ) . \deg_R([q^t]C)\le\kappa_0+\kappa t,\qquad
\nu([R^n]C)\ge\max\left(0,\left\lceil\frac{n-\kappa_0}{\kappa}\right\rceil\right). deg R ([ q t ] C ) ≤ κ 0 + κ t , ν ([ R n ] C ) ≥ max ( 0 , ⌈ κ n − κ 0 ⌉ ) .
The same estimate applies to each component of m m m , omitting deg R d \deg_R d deg R d . This supplies convergent inputs to the searches below. For the marginal of a partition difference, Lemma 3.5 gives the sharper bound ν ( [ R n ] C ) ≥ n \nu([R^n]C)\ge n ν ([ R n ] C ) ≥ n .
Polynomial and rational operator witnesses
Polynomial and rational operator witnesses. In (9), seek operator rows Λ 1 , Λ 2 \Lambda_1,\Lambda_2 Λ 1 , Λ 2 with
d = ( 1 − R ) Λ 1 + Λ 2 L , Λ 2 f = ( 1 − R ) b 0 , b 0 ∈ K [ R ] . d=(1-R)\Lambda_1+\Lambda_2L,\qquad\Lambda_2f=(1-R)b_0,\qquad b_0\in\mathbb{K}[R]. d = ( 1 − R ) Λ 1 + Λ 2 L , Λ 2 f = ( 1 − R ) b 0 , b 0 ∈ K [ R ] .
Then C = ( 1 − R ) Λ 1 m + Λ 2 f C=(1-R)\Lambda_1m+\Lambda_2f C = ( 1 − R ) Λ 1 m + Λ 2 f , so b = Λ 1 m + b 0 b=\Lambda_1m+b_0 b = Λ 1 m + b 0 . If Λ 1 \Lambda_1 Λ 1 has polynomial entries with coefficients regular at q = 0 q=0 q = 0 , the preceding degree bound, applied to m m m , proves convergence. An ansatz of bounded T \mathbb{T} T -degree for Λ 2 ∈ K [ R ] 1 × s \Lambda_2\in\mathbb{K}[R]^{1\times s} Λ 2 ∈ K [ R ] 1 × s gives a finite linear system: reduction modulo the right ideal ( 1 − R ) O (1-R)\mathcal{O} ( 1 − R ) O is evaluation at R = 1 R=1 R = 1 in normal order, and the remaining condition is ( Λ 2 f ) ( 1 ) = 0 (\Lambda_2f)(1)=0 ( Λ 2 f ) ( 1 ) = 0 . Rational or Laurent functions of T \mathbb{T} T enlarge the search, but an inverse shift costs valuation, T − k ∑ q n R n = ∑ q − k n q n R n \mathbb{T}^{-k}\sum q_nR^n=\sum q^{-kn}q_nR^n T − k ∑ q n R n = ∑ q − k n q n R n , and poles at T = q n \mathbb{T}=q^n T = q n require separately specified coefficients at finitely many exceptional indices; such witnesses need a corresponding additional convergence proof.
Contracting coefficient recurrences
Coefficient extraction from a holonomic system gives, beyond finitely many singular indices, a finite state recurrence v n + 1 = M ( q n ) v n v_{n+1}=M(q^n)v_n v n + 1 = M ( q n ) v n and an output b n = π ( q n ) v n b_n=\pi(q^n)v_n b n = π ( q n ) v n . A diagonal gauge with entries powers of q n q^n q n , and a fixed block of ℓ \ell ℓ steps, may make every transition entry have valuation at least one and every output entry have valuation bounded below. This proves linear decay along that progression by induction. Rational functions evaluated at q n q^n q n have eventually affine valuations, by Lemma 3.4 (2); consequently the required eventual inequalities reduce to finitely many checks. This sufficient test does not capture cancellation special to one initial vector, such as the one recorded by Proposition 3.7 .
Spans of convergent holonomic series
Call a holonomic series g = ∑ g n R n ∈ K [ [ R ] ] g=\sum g_nR^n\in K[[R]] g = ∑ g n R n ∈ K [[ R ]] a certified convergent input when it comes with selecting data and a proved bound ν ( g n ) ≥ δ n − c \nu(g_n)\ge\delta n-c ν ( g n ) ≥ δ n − c , δ > 0 \delta>0 δ > 0 . Normalized state systems supply such inputs. Sums, products, sections and polynomial operator images have effective holonomic constructions and explicit bounds. Another useful constructor is
∇ λ g ( R ) = g ( R ) − g ( R λ ) 1 − R , λ ≥ 2 , for which [ R n ] ∇ λ g = ∑ j = ⌈ ( n + 1 ) / λ ⌉ n g j . \nabla_\lambda g(R)=\frac{g(R)-g(R^\lambda)}{1-R},\qquad\lambda\ge2,\qquad\text{for which}\qquad[R^n]\nabla_\lambda g=\sum_{j=\lceil(n+1)/\lambda\rceil}^{n}g_j. ∇ λ g ( R ) = 1 − R g ( R ) − g ( R λ ) , λ ≥ 2 , for which [ R n ] ∇ λ g = j = ⌈( n + 1 ) / λ ⌉ ∑ n g j .
so that its slope is at least δ / λ \delta/\lambda δ / λ . It is holonomic by closure under monomial pullback and rational multiplication. For monomial pullback, first eliminate to an equation involving only powers of Υ λ \Upsilon^\lambda Υ λ , then substitute R ↦ R λ R\mapsto R^\lambda R ↦ R λ ; this turns Υ λ \Upsilon^\lambda Υ λ into the q q q -shift of the pulled-back series.
For a finite family g 1 , … , g t g_1,\ldots,g_t g 1 , … , g t and a degree bound, seek
C = ∑ i , k , ℓ c i , k , ℓ ( 1 − R ) R ℓ Υ k g i , c i , k , ℓ ∈ K . C=\sum_{i,k,\ell}c_{i,k,\ell}(1-R)R^\ell\Upsilon^k g_i,\qquad c_{i,k,\ell}\in K. C = i , k , ℓ ∑ c i , k , ℓ ( 1 − R ) R ℓ Υ k g i , c i , k , ℓ ∈ K .
Take a common scalar annihilator Ω \Omega Ω of C C C and of the finitely many terms on the right. By Lemma 3.4 (1), a linear combination vanishes identically if and only if its coefficients at E Ω E_\Omega E Ω vanish, so the relation space is the kernel of the finite matrix of those exceptional coefficients. A relation with nonzero coefficient of C C C expresses b b b through certified convergent inputs, and so proves its convergence.
Enumerating finite expressions in these constructors, with bounds on operator degrees, gives a sound search, which terminates whenever the canonical b b b belongs to the span of finitely many enumerated certified convergent inputs. In parallel, the criterion (15) eventually detects every false identity.
Acknowledgments
I am deeply grateful to Jehanne Dousse for stimulating discussions on partition identities, and for pointing me to many examples to which the methods of this article apply.
ChatGPT 6 Astra and Claude Opus 5.5 helped write the code that searches for and checks the certificates, and Claude helped prepare the manuscript. No AI was used to set up the underlying mathematical theory here.
Julien David generously contributed the enumerator of finite automata used in §§§ 3.8 and § 4.7 (exhaustivautomata in [3 ] ).