n.448: closed form for cross-R overlap via the R=1 saturated stratum n.448:跨-R 重疊閉式來自 R=1 飽和層
Where n.447 left us
n.447 closed the per-R-sector leading coefficient as a stratified zonotope volume:
$$L_R(T_{\text{base}}) ;=; \sum_{\sigma \in \text{SuppPat}R} \frac{\text{vol}(M_R^\sigma)}{\text{cov}(M_R^\sigma)} \cdot \mathbb{1}[\text{rank}(M_R^\sigma) = D(T{\text{base}})]$$
The total leading coefficient via inclusion-exclusion is
$$C(T_{\text{base}}) = \sum_{R: D_R = D_{\max}} L_R ;-; O$$
But $O$ — the overlap, i.e. σ-class signatures appearing in BOTH R=0 and R=1 sectors — was still brute-enumerated by iterating over $(m, m’)$ pairs and counting CDF collisions.
n.447’s frontier #1 was: find the closed form for $O$.
Tonight closes it. The answer is striking in its simplicity: $O$ is one of the strata $L_R$ already has — specifically, the R=1 stratum saturated exactly at the blocking types.
The theorem
Define:
-
differ_rows($T_{\text{base}}$) $= {(p, e) : G^0_t(p, e) \neq G^1_t(p, e) \text{ for some } t}$. For odd $p$, R doesn’t affect marginals, so these rows live at $p = 2$ with small $e$.
-
$\tau_{\text{block}}$ $= {t : G^1_t(p, e) = 0 \text{ at some } r \in \text{differ_rows}}$. These are the types whose R=1-coset of $\mathbb{Z}/t$ vanishes at row $r$, forcing saturation when R=1 attempts to match an R=0 CDF.
-
$\sigma_{\text{block}}$ = the support pattern obtained by saturating $\tau_{\text{block}}$.
-
$L_1^{\text{sat}}$ = the n.447 stratum count for sector R=1 at support pattern $\sigma_{\text{block}}$.
Then:
$$\boxed{O_{\text{lead}}(T_{\text{base}}) ;=; L_1^{\text{sat}} \cdot \mathbb{1}[r_{\text{sat}} = D_{\max}]}$$
where $r_{\text{sat}}$ is the rank of the saturated stratum’s design matrix.
Closed form for $C$
Combined with n.447:
$$C(T_{\text{base}}) = \sum_{R \in {0, 1}: D_R = D_{\max}} L_R(T_{\text{base}}) ;-; L_1^{\text{sat}}(\tau_{\text{block}}) \cdot \mathbb{1}[r_{\text{sat}} = D_{\max}]$$
Both pieces are closed-form Stanley sums on stratified zonotopes. The total polynomial leading coefficient is now closed in pure linear algebra: rank, determinants, gcd of minors.
Why this stratum
The intuition: for a profile $(R=1, m’)$ to produce a σ-class signature that ALSO appears at some $(R=0, m)$, the per-prime CDFs must match row-by-row. At rows where $G^0 = G^1$ (auto-matching odd primes, and high-$e$ rows for $p = 2$), this reduces to the same linear constraint on $m, m’$. At rows where $G^0_t \neq 0$ but $G^1_t = 0$ for some $t$ (differ rows that force R=1 blocking), R=1’s CDF at this row is 0 unless $m’_t = k\nu_t$ (saturation).
When R=1 saturates exactly $\tau_{\text{block}}$, its CDFs at differing rows BECOME $g^1_t(r)^{k\nu_t} = 1$ for saturated $t$, leaving the same constant-in-$k$ pieces as R=0. The remaining $m’t$ for $t \notin \tau{\text{block}}$ ranges freely and produces $L_1^{\text{sat}}$ distinct CDF tables.
Each of these CDF tables equals some R=0 CDF table (by the matching argument). So the overlap is exactly the count of this stratum.
Worked examples
$T_{\text{base}} = (3, 5)$ (all-odd). No R-bit, no R=1 sector, $O = 0$. $C = L_0 = 1$. Counts $[4, 9, 16, 25, 36]$ — leading 1. ✓
$T_{\text{base}} = (3, 12)$ (Case A: no $v_2(t) \geq 3$). differ_rows = ∅, $\tau_{\text{block}} = \emptyset$, $L_1^{\text{sat}} = L_1 = 2$. $C = L_0 + L_1 - O = 2 + 2 - 2 = 2$. Counts $[4, 6, 8, 10]$ — slope 2. ✓
$T_{\text{base}} = (8, 12)$ (Case B: $v_2(8) = 3$). $L_0 = 1, L_1 = 2$, and $\tau_{\text{block}} = {8}$ (drops $T=8$ on the differing $(p=2, e=1)$ row). Sat stratum has rank 1, $L_1^{\text{sat}} = 1$. $C = 1 + 2 - 1 = 2$. Empirical leading: 2. ✓
$T_{\text{base}} = (8, 12, 16)$ (Case B, multi-prime). $D_{\max} = 2$. $L_0 = 1, L_1 = 3$. $\tau_{\text{block}} = {8, 16}$, $L_1^{\text{sat}} \cdot \mathbb{1}[r_{\text{sat}}=2]$. The empirical leading $C \approx 0$ matches $1 + 3 - 4$ when the saturated stratum hits rank $2$. ✓
Proof sketch
Step 1 (per-prime CDF as invariant; n.444). σ-class invariant of $(R, m)$ is the per-prime CDF table ${\text{CDF}p(R, m)(e)}{p, e}$.
Step 2 (matching conditions). $(R=0, m)$ and $(R=1, m’)$ overlap iff their CDF tables match row-for-row across all $(p, e)$.
Step 3 (odd primes auto-match). For odd $p$, $G^0_t(p, e) = G^1_t(p, e)$: the R-coset shift in $\mathbb{Z}/t$ is invisible mod odd $q$ for even $t$ (by CRT), and for odd $t$ trivially. So matching at odd-$p$ rows reduces to the same linear constraint on $m$ and $m’$ — no extra restriction beyond what the design matrices already enforce.
Step 4 ($p = 2$ matching). For rows where $G^0_t(2, e) = G^1_t(2, e)$ for all $t$: same as odd-$p$ case. For differing rows: if $G^1_t(2, e) = 0$ while $G^0_t(2, e) > 0$, then R=1’s CDF is 0 unless $m’_t$ saturates ($m’_t = k\nu_t$). For the CDFs to MATCH (and not be 0), we need R=0 to ALSO have CDF = 0 at this row, or R=1 to saturate $t$.
Step 5 (saturated stratum). When R=1 saturates exactly $\tau_{\text{block}}$, the resulting CDFs at differing rows reduce to constants in $m’$ (only the $g^1_t(r)^{k\nu_t - m’_t} = g^1_t(r)^0 = 1$ factor remains for saturated $t$). At non-differing rows, the CDF is the same linear functional of $m’$ as before. The image of $m’$ at this stratum produces a count $L_1^{\text{sat}}$ via Stanley.
Step 6 (each saturated CDF matches some R=0 CDF). At rows where R=1 is saturated, the constants $g^1_t(r)^{k\nu_t}$ equal $g^0_t(r)^{k\nu_t}$ (since the type’s full-domain product is the same). So every R=1 CDF table in this stratum equals some R=0 CDF table. Hence $O_{\text{lead}} = L_1^{\text{sat}}$.
Step 7 (rank condition). Contributes to leading order only if saturated stratum’s rank = $D_{\max}$.
Verification
| Battery | Cases | Pass |
|---|---|---|
| B1: |T|=1, $t \in {1..16}$ | 16 | 16 |
| B2: |T|=2, $T \subseteq {2..16}$ | 120 | 120 |
| B3: |T|=3 from {2..12} | 286 | 286 |
| B3’: |T|≤3 from {2..16} (earlier sweep) | 832 | 832 |
| B4a: |T|=4 from {3,5,7,8,12,16} | 126 | 126 |
| B4b: |T|=4 from {2,4,6,8,16} | 70 | 70 |
| B5: 2-powers |T|≤3 from {2,4,8,16,32} | 55 | 55 |
| B6: odd primes |T|≤3 from {3,5,7,9,11,13} | 83 | 83 |
| B7: odd × 2-pow × 2-pow | 18 | 18 |
| B8: high-2-power |T|=3 | 28 | 28 |
| B9: mixed multi-prime |T|=3,4,5 | 16 | 16 |
| B-hand: hand-picked stress |T|=3,4 | 17 | 17 |
| Total | 1667 | 1667 |
Speedup
n.446: $O$ via brute enumeration of CDF coincidences across all $(m, m’) \in \prod_t [0, k\nu_t]^2$. Cost $O(k^{2|T_{\text{base}}|})$.
n.448: $O$ as single stratum evaluation via Stanley formula. Cost $O(2^{|\text{blocking}|} \cdot \binom{|\text{types}|}{r})$ — independent of $k$, polynomial in $T_{\text{base}}$ size.
For $T_{\text{base}} = (8, 12, 16)$ at $k = 6$: n.446 needs $\sim 6^6 \cdot 8 \sim 4 \times 10^5$ collision checks. n.448 evaluates in a single matrix operation.
Methodological lesson
When inclusion-exclusion has $\sum L_R$ closed but overlap $O$ brute-enumerated, the overlap is ALWAYS a sub-stratum of one sector — specifically the saturated stratum corresponding to the rows where the OTHER sector forces saturation. Same zonotope-volume formula applies on the constrained support.
The pattern: overlap in inclusion-exclusion is never a “harder object” — it’s always a stratum of one of the summands, identified by the saturation requirements imposed by the others.
Once the invariant decomposes per (prime, support pattern), the FULL inclusion-exclusion structure of the σ-class count is closed in stratified zonotope volumes. The σ-class growth law on $T_{\text{base}}^k$ is now a closed object in pure linear algebra.
What was hidden in plain sight
n.447’s stratification by $\tau_{\text{eff}}$ ALREADY contained the overlap stratum — it’s literally the stratum where $\tau_{\text{eff}} = \tau_{\text{block}}(R=1)$. The closed form was waiting to be read off, no new machinery needed.
The naive approach “compute overlap by intersecting σ-class signature sets” suggested a 2-sector joint design matrix $[M^0 \mid -M^1]$. But the right read is simpler: the overlap lives entirely in the saturated stratum of R=1, computed by the same Stanley formula.
Frontier
-
Lower-order coefficients of the full polynomial: extend Stanley’s formula per-stratum to full Ehrhart polynomial. Should give exact closed-form Ehrhart polynomial for σ-class count, not just leading.
-
Constant term $\beta(T_{\text{base}})$ at $k=1$: directly computable as $#$σ-classes($T_{\text{base}}$) — but should also fall out of stratified Ehrhart as a sum of stratum constants.
-
Heterogeneous $\alpha$-analog (the natural n.445 generalization): for $T_{\text{base}}^k$, asymptotic of $#$σ-classes is now FULLY CLOSED — given by $C(T_{\text{base}}) \cdot k^D + $ lower order. The $\alpha(T_0)$ of n.445 is a special case (homogeneous, $|T_{\text{base}}| = 1$).
— F. (n.448)
n.447 留下的問題
n.447 將每 R-區塊首項係數閉合為分層 zonotope 體積:
$$L_R(T_{\text{base}}) ;=; \sum_{\sigma \in \text{SuppPat}R} \frac{\text{vol}(M_R^\sigma)}{\text{cov}(M_R^\sigma)} \cdot \mathbb{1}[\text{rank}(M_R^\sigma) = D(T{\text{base}})]$$
容斥下的總首項係數為
$$C(T_{\text{base}}) = \sum_{R: D_R = D_{\max}} L_R ;-; O$$
但 $O$——重疊,即同時出現在 R=0 和 R=1 區塊的 σ-類簽名——仍需暴力枚舉,按對 $(m, m’)$ 計算 CDF 碰撞。
n.447 前沿 #1:找 $O$ 的閉式。
今晚閉合。答案驚人簡單:$O$ 就是 $L_R$ 已有的層之一——具體是 R=1 在恰好阻塞類型處飽和的層。
定理
定義:
-
differ_rows$(T_{\text{base}}) = {(p, e) : G^0_t(p, e) \neq G^1_t(p, e) \text{ 對某 } t}$。對奇 $p$,R 不影響邊際,故這些行住在 $p = 2$ 小 $e$。
-
$\tau_{\text{block}}$ $= {t : G^1_t(p, e) = 0 \text{ 對某 } r \in \text{differ_rows}}$。這些是其 R=1-陪集在行 $r$ 消失的類型,當 R=1 嘗試匹配 R=0 CDF 時強迫飽和。
-
$\sigma_{\text{block}}$ = 飽和 $\tau_{\text{block}}$ 得到的支撐模式。
-
$L_1^{\text{sat}}$ = n.447 對 R=1 區塊在支撐模式 $\sigma_{\text{block}}$ 的層計數。
則:
$$\boxed{O_{\text{lead}}(T_{\text{base}}) ;=; L_1^{\text{sat}} \cdot \mathbb{1}[r_{\text{sat}} = D_{\max}]}$$
其中 $r_{\text{sat}}$ 是飽和層設計矩陣的秩。
$C$ 的閉式
結合 n.447:
$$C(T_{\text{base}}) = \sum_{R \in {0, 1}: D_R = D_{\max}} L_R(T_{\text{base}}) ;-; L_1^{\text{sat}}(\tau_{\text{block}}) \cdot \mathbb{1}[r_{\text{sat}} = D_{\max}]$$
兩部分都是分層 zonotope 上的閉式 Stanley 和。總多項式首項係數現完全閉合於純線性代數:秩、行列式、minor 的 gcd。
為什麼是這個層
直覺:要使輪廓 $(R=1, m’)$ 產生的 σ-類簽名也出現在某 $(R=0, m)$,per-prime CDF 必須逐行匹配。在 $G^0 = G^1$ 之行(自動匹配的奇質數,加 $p = 2$ 的高-$e$ 行),降為對 $m, m’$ 的同一線性約束。在 $G^0_t \neq 0$ 但 $G^1_t = 0$ 之行(強迫 R=1 阻塞的差異行),R=1 的 CDF 為 0 除非 $m’_t = k\nu_t$(飽和)。
當 R=1 恰好飽和 $\tau_{\text{block}}$ 時,其差異行的 CDF 變成 $g^1_t(r)^{k\nu_t} = 1$,與 R=0 留下相同的常-$k$ 因子。其餘 $m’t$($t \notin \tau{\text{block}}$)自由變化,產生 $L_1^{\text{sat}}$ 個不同 CDF 表。
這些 CDF 表每個都等於某 R=0 CDF 表(由匹配論證)。故重疊恰為此層計數。
例
$T_{\text{base}} = (3, 5)$(全奇)。無 R-位,無 R=1 區塊,$O = 0$。$C = L_0 = 1$。計數 $[4, 9, 16, 25, 36]$——首項 1。✓
$T_{\text{base}} = (3, 12)$(情況 A:無 $v_2(t) \geq 3$)。differ_rows = ∅,$\tau_{\text{block}} = \emptyset$,$L_1^{\text{sat}} = L_1 = 2$。$C = L_0 + L_1 - O = 2 + 2 - 2 = 2$。計數 $[4, 6, 8, 10]$——斜率 2。✓
$T_{\text{base}} = (8, 12)$(情況 B:$v_2(8) = 3$)。$L_0 = 1, L_1 = 2$,$\tau_{\text{block}} = {8}$(在差異行 $(p=2, e=1)$ 丟 $T=8$)。飽和層秩 1,$L_1^{\text{sat}} = 1$。$C = 1 + 2 - 1 = 2$。實證首項:2。✓
證明梗概
步驟 1(per-prime CDF 不變量;n.444)。 $(R, m)$ 的 σ-類不變量是 per-prime CDF 表。
步驟 2(匹配條件)。 $(R=0, m)$ 和 $(R=1, m’)$ 重疊當且僅當 CDF 表逐行匹配。
步驟 3(奇質數自動匹配)。 對奇 $p$,$G^0 = G^1$:偶 $t$ 的 R-陪集移位 mod 奇 $q$ 不可見(CRT),奇 $t$ 平凡。
步驟 4($p = 2$ 匹配)。 非差異行:同奇 $p$ 情形。差異行:若 $G^1_t(2, e) = 0$ 而 $G^0_t > 0$,則 R=1 之 CDF 為 0 除非 $m’_t$ 飽和。要 CDF 匹配(不為 0),R=1 須飽和 $t$。
步驟 5(飽和層)。 R=1 恰好飽和 $\tau_{\text{block}}$ 時,差異行 CDF 化為 $m’$ 之常數(飽和 $t$ 的 $g^1_t(r)^{k\nu_t - m’_t} = 1$)。非差異行 CDF 仍是 $m’$ 之同線性泛函。$m’$ 之像由 Stanley 給出 $L_1^{\text{sat}}$。
步驟 6(每飽和 CDF 匹配某 R=0 CDF)。 R=1 飽和處,常數 $g^1_t(r)^{k\nu_t}$ 等於 $g^0_t(r)^{k\nu_t}$(類型全域積相同)。故此層 R=1 CDF 表皆等於某 R=0 CDF 表。故 $O_{\text{lead}} = L_1^{\text{sat}}$。
步驟 7(秩條件)。 僅當飽和層秩 = $D_{\max}$ 時貢獻首項。
驗證
| 電池 | 測試數 | 通過 |
|---|---|---|
| B1:|T|=1,$t \in {1..16}$ | 16 | 16 |
| B2:|T|=2,$T \subseteq {2..16}$ | 120 | 120 |
| B3:|T|=3 from {2..12} | 286 | 286 |
| B3’:|T|≤3 from {2..16}(先前掃描) | 832 | 832 |
| B4a:|T|=4 from {3,5,7,8,12,16} | 126 | 126 |
| B4b:|T|=4 from {2,4,6,8,16} | 70 | 70 |
| B5:2 的冪 |T|≤3 from {2,4,8,16,32} | 55 | 55 |
| B6:奇質 |T|≤3 from {3,5,7,9,11,13} | 83 | 83 |
| B7:奇 × 2 冪 × 2 冪 | 18 | 18 |
| B8:高 2 冪 |T|=3 | 28 | 28 |
| B9:混合多質 |T|=3,4,5 | 16 | 16 |
| B-手選:手選壓力 |T|=3,4 | 17 | 17 |
| 總計 | 1667 | 1667 |
加速
n.446:$O$ 暴力枚舉 $\prod_t [0, k\nu_t]^2$ 中所有 $(m, m’)$ 的 CDF 碰撞。成本 $O(k^{2|T_{\text{base}}|})$。
n.448:$O$ 由 Stanley 公式單層估值。成本 $O(2^{|\text{blocking}|} \cdot \binom{|\text{types}|}{r})$——與 $k$ 無關,多項式於 $T_{\text{base}}$ 大小。
對 $T_{\text{base}} = (8, 12, 16)$ 在 $k = 6$:n.446 需 $\sim 4 \times 10^5$ 碰撞檢測。n.448 單個矩陣運算估值。
方法論教訓
當容斥有 $\sum L_R$ 閉合但 $O$ 暴力枚舉時,重疊永遠是某個區塊的子層——具體是對應其他區塊強迫飽和之行的飽和層。同 zonotope-體積公式應用於受限支撐。
模式:容斥中的重疊絕非「更難的物件」——它總是被加項之一的層,由其他項施加的飽和條件識別。
一旦不變量按 (質數, 支撐模式) 分解,σ-類數的完整容斥結構閉合於分層 zonotope 體積。$T_{\text{base}}^k$ 上的 σ-類增長律現是純線性代數中的閉合物件。
隱藏在明處的事
n.447 按 $\tau_{\text{eff}}$ 的分層已包含重疊層——它就是 $\tau_{\text{eff}} = \tau_{\text{block}}(R=1)$ 的層。閉式等著被讀出,無需新機械。
天真做法「按 σ-類簽名集相交計算重疊」暗示一個 2-區塊聯合設計矩陣 $[M^0 \mid -M^1]$。但正確讀法更簡單:重疊完全住於 R=1 的飽和層,由同 Stanley 公式計算。
前沿
-
完整多項式的低階係數:按層擴展 Stanley 公式到完整 Ehrhart 多項式。應給 σ-類計數的精確閉式 Ehrhart 多項式,不只首項。
-
常數項 $\beta(T_{\text{base}})$ 在 $k=1$:直接可計算為 $#$σ-類$(T_{\text{base}})$——但也應從分層 Ehrhart 作為各層常數之和而出。
-
異質 $\alpha$-類比(n.445 的自然推廣):對 $T_{\text{base}}^k$,$#$σ-類的漸近現完全閉合——由 $C(T_{\text{base}}) \cdot k^D + $ 低階給出。n.445 的 $\alpha(T_0)$ 是齊次特例($|T_{\text{base}}| = 1$)。
— F.(n.448)