research/lab/py/restricted-franel
0 directories and 2 files in research/lab/py/restricted-franel.
restricted-franel
- Franel's identity on a digit design: an exact formula for the Fourier and the rank
L^2discrepancy of the denominator-restricted Farey set, written as a gcd-weighted double sum over DILATED Mertens sums. - The object is the Farey page's THE METER ON A DIGIT DESIGN, denominator convention, written
F_Q^d(S_F)here because the Farey page bindsF_Q(S_F)to the STRICT set:F_Q^d(S_F) = {a/b reduced, b in S_F, b <= Q, 1 <= a <= b}, withS_Fthe whole numbers whose every base digit lies in a digit setF. - Five results: the frequency-
mexponential sum (Theorem 1), the identity in both its Fourier and its rank form (Theorem 2), the one-line inequality that carries the restricted Mertens sum back out of it (Theorem 3), the strict convention's divisor-sum identity with the refutation that kills every Mertens-type closed form for it (Theorem 4), and the dilate automaton with the converse it closes (Theorem 5). - The template is gaussian-franel, which is the same transcription one field up. The enumeration is this study's own: farey-discrepancy is a Rust crate whose
Designtype and mediant walk are not exported to Python, so the digit test and the design list are rebuilt here in eight lines and the node counts are checked against that study's published table.
THE OBJECT
- Design:
base >= 2, digit setF,S_Fthe whole numbers whose every base digit lies inF, mass exponentalpha = log|F|/log base. The two designs run here are base 3 withF = {0,1}and base 10 with the digit 9 missing, the pair the Farey page meters. - Denominator set:
F_Q^d(S_F) = {a/b : b in S_F, b <= Q, 1 <= a <= b, gcd(a,b) = 1}, of sizem_F(Q) = sum_{b in S_F, b <= Q} phi(b). - Strict set:
F_Q^s(S_F) = {a/b : a in S_F, b in S_F, b <= Q, 1 <= a <= b, gcd(a,b) = 1}. - The dilated Mertens sums.
M_F(x; d) = sum_{c <= x, dc in S_F} mu(c), the Mertens function of the dilated designd^{-1} S_F = {c : dc in S_F}. Atd = 1it is the design's own meterM_F(x) = sum_{n in S_F, n <= x} mu(n), the Mobius page's. The dilates are the whole content of the restriction:d^{-1} S_Fis notS_F, is not a digit design in general, and carries no digit test, so everyd > 1brings a new function rather than a rescaling of the old one. On the full set every dilate isZand all of them collapse toM. - Exponential sums at frequency
k:S_F(k, Q) = sum_{r in F_Q^d(S_F)} e(kr)withe(x) = exp(2 pi i x), andS_F^s(k, Q)the same over the strict set.
THEOREM 1, THE EXPONENTIAL SUM
- Ramanujan's sum over a digit design. For every
m >= 1andQ >= 1,S_F(m, Q) = sum_{d | m} d M_F(Q/d; d). Proved. - Proof. Partition by denominator:
S_F(m, Q) = sum_{b in S_F, b <= Q} c_b(m)withc_b(m) = sum_{a mod b, gcd(a,b) = 1} e(ma/b)Ramanujan's sum, which depends onbalone and not on the design. Kluyver's formula isc_b(m) = sum_{d | gcd(m, b)} d mu(b/d). Exchanging the two sums givesS_F(m, Q) = sum_{d | m} d sum_{b in S_F, b <= Q, d | b} mu(b/d), and substitutingb = dcturns the inner sum intoM_F(Q/d; d). - At frequency 1 the exponential sum IS the design's Mertens meter.
S_F(1, Q) = M_F(Q). Proved, them = 1case, onlyd = 1surviving. It needs no Ramanujan input:sum_{a mod b, gcd(a,b) = 1} e(a/b) = sum_{d | b} mu(d) sum_{c mod b/d} e(c/(b/d)) = mu(b), the inner complete sum vanishing unlessb/d = 1. - The classical case is
S_F = Z, whereM_F(Q/d; d) = M(Q/d)for everydand the formula readssum_{d | m} d M(Q/d), the divisor-shifted Mertens sums. - Verified. Every denominator
b in S_Fup toQ = 10^5at base 3{0,1}, and up toQ = 10^4at base 10 without 9 and on the control, has its literal sum ofphi(b)roots of unity equal tomu(b): worst deviation1.09e-11atb = 86293on the base 3 design and1.36e-12atb = 7247on base 10 without 9,0roundings to the wrong integer at either. TheQ = 10^5rung on base 10 without 9 costs2.6e9roots of unity, past the machine budget, so it is not walked and nothing is claimed at it. The aggregated literal sum agrees with the exact integerM_F(Q)to2.02e-09at the top rung, floating-point accumulation over5.4e7roots of unity and not a mismatch. The frequency-mformula is checked atm = 1, 2, 3, 4, 5, 6, 12on both designs and the control (verb_denominator).
THEOREM 2, FRANEL ON A DIGIT DESIGN
- The kernel.
G_F(Q) = sum_{d, e >= 1} (gcd(d,e)^2/(d e)) M_F(Q/d; d) M_F(Q/e; e), a finite sum of exact rationals, every term withd > Qore > Qvanishing. - The Fourier form.
sum_{k != 0} |S_F(k, Q)|^2 / k^2 = (pi^2/3) G_F(Q). Proved. Substitute Theorem 1, expand the square and exchange: for fixedd, ethe inner sum issum_{k != 0, lcm(d,e) | k} k^-2 = 2 zeta(2)/lcm(d,e)^2, andd e/lcm(d,e)^2 = gcd(d,e)^2/(d e). - The rank form. If the node set is closed under
r -> 1 - raway from the node1, then withrho_1 < ... < rho_mthe nodes ofF_Q^d(S_F)ascending,m = m_F(Q)anddelta_j = rho_j - j/m,G_F(Q) - 1 = 12 m sum_j delta_j^2, an identity of exact rationals. Proved. Closure is used only to fixsum_r r = (m+1)/2, so THAT mean-value condition is the surviving hypothesis and closure is one sufficient condition for it; the piecewise integration below needs nothing of the node set beyond a top node of1. Closure holds for every denominator-restricted set, sincea/breduced withb in S_Fgives(b-a)/breduced with the samebanda = bonly at the node1; it fails for every proper strict set. - Proof of the rank form. With
B1barthe sawtooth andA(v)the counting function of the nodes,sum_r B1bar(u + r) = D(1 - u) + cwithD(v) = A(v) - m vandc = sum_r r - m/2, so Parseval givessum_{k != 0} |S_F(k,Q)|^2/k^2 = 4 pi^2 (int_0^1 D^2 - c^2). IntegratingD^2piecewise between consecutive nodes and telescoping givesint_0^1 D^2 = m sum_j delta_j^2 - (m+1)/2 + sum_r r + 1/3, the two boundary cubes vanishing becauserho_m = 1. Under the reflection hypothesis the nodes other than1pair intom - 1values summing to(m-1)/2, sosum_r r = (m+1)/2andc = 1/2, and the remainder collapses to1/12. - The kernel is the gcd matrix
gcd(d,e)^2/(d e), the same Smith kernel that carries the moire correlation law on the stack page and the Gaussian identity in gaussian-franel. Digit restriction moves the entries, never the kernel. - Verified. The rank form is an exact rational identity at base 3
{0,1}Q = 81andQ = 243, base 10 without 9Q = 40, and the full-set controlQ = 40:Trueat all four. The Fourier form is checked against the literal Fourier side truncated at|k| <= 200000with the printed tail bound2 m^2/K, the gap inside the bound in every row. The control atQ = 40regenerates Edwards section 12.2,G_F(40) - 1 = 12 m sum delta^2withm = 490andsum delta^2 = 0.0104270117, which is the number gaussian-franel prints (verb_identity).
THEOREM 3, WHAT COMES BACK OUT
- The inequality.
2 M_F(Q)^2 <= (pi^2/3) G_F(Q), and by the rank form(pi^2/3) G_F(Q) = 4 pi^2 m_F(Q) sum_j delta_j^2 + pi^2/3. Proved, in one line:k = 1andk = -1contribute2|S_F(1,Q)|^2 = 2 M_F(Q)^2to a sum of nonnegative terms. The constant is not decorative: the weaker-looking2 M_F(Q)^2 <= G_F(Q)is FALSE, failing already at the full-set controlQ = 5(2 M(5)^2 = 8againstG(5) = 64/15) and at base 3{0,1}Q = 37(18againstG_F(37) = 14.230517). - So a square-root bound on the denominator-restricted Farey discrepancy forces a square-root bound on the design's Mertens meter. The converse needs the whole
ksum and is not claimed here. - The implication.
A_F(Q) = #{n in S_F : n <= Q} << Q^alpha, the block count|F|^levelatQ = base^levelwith a constant depending on the design alone, som_F(Q) = sum_{b in S_F, b <= Q} phi(b) <= Q A_F(Q) << Q^(1+alpha). Hence the denominator lane's CONJECTURED shapesum_j delta_j^2 = O(Q^(-1+eps))on the Farey page givesG_F(Q) = O(Q^(alpha+eps))through the rank form, and the inequality then gives|M_F(Q)| = O(Q^(alpha/2+eps)), which is the square-root ceiling for the design's Mertens meter on the Mobius page. Proved, no unproved input entering, the constantpi^2/3being absorbed. What yields the ceiling is the conjectured exponent-1and not the measured one: the lane's measurede_2sits at-0.959and-0.899at the top rungs andS2 Qis still climbing there, and a proof of only the measured shape would give no ceiling. Digit restriction of the denominator is invisible to the SHAPE of the meter, and this is what that invisibility costs. - The converse is open and is where the road ends. It needs the whole
ksum controlled from Mertens bounds, so it needsM_F(x; d)ford > 1, and no bound for a dilated design's Mertens function exists on the tree. - Verified, at every integer
Qrather than at a sample. Both sides step only atQ in S_F, sinceM_F(Q)and everyM_F(Q/d; d)change only whenQitself joinsS_F, so scanningS_Fcovers every integerQbelow the bound. The inequality is asserted at0violations overQ <= 2187at base 3{0,1}(128 jumps) andQ <= 400at base 10 without 9 (324) and on the control (400). The ratio2 M_F(Q)^2/((pi^2/3) G_F(Q))peaks at0.607927at the trivialQ = 1on all three; overQ >= 100its maximum is0.340071atQ = 253on base 3{0,1},0.137645atQ = 221on base 10 without 9 and0.086385atQ = 114on the control (scan_backward). - That the
k = +-1terms carry a third of the whole functional atQ = 253on a thin design is worth its own look: if the fraction does not fall withQ, the inequality is near sharp there and the converse is nearer than OBSTRUCTIONS claims. Not measured here.
THEOREM 4, THE STRICT SET
- The divisor-sum identity.
S_F^s(1, Q) = sum_{b in S_F, b <= Q} sum_{d | b} mu(d) sum_{a <= b/d, da in S_F} e(a/(b/d)). Proved, Mobius inversion of the coprimality condition followed bya -> da. It is exact and it is where the road stops: the inner sum is a digit-restricted exponential sum over an arithmetic progression, the Type II object the Mobius page has no bound for. - Verified at every denominator: literal summation over
{a in S_F : a <= b, gcd(a,b) = 1}against the divisor route agrees to6.28e-15over the 64 denominators of base 3{0,1}below 729 and to1.95e-14over the 162 denominators of base 10 without 9 below 200. - No Mertens-type sum over
S_Fequals it. Refuted, two ways, the first exact. - First, the strict sum is not real. Base 3
{0,1}atQ = 3gives1 + e(1/3) = 0.5 + (sqrt 3/2) i, and base 10 without 9 atQ = 10gives-1.809016994 + 0.587785252 i; every Mertens-type sum overS_Fis a real integer or a real rational. The two witnesses carry the whole refutation. Beside them sits an observation and not a mechanism: the strict set also fails the pairinga -> b - athat makes the denominator set's sum real, sincea in S_Fdoes not giveb - a in S_F. Failure of that pairing is not shown to force a non-real sum, and nothing here rests on it. - Second, it rides the mass.
|S_F^s(1,Q)|/cardsettles at0.335693, 0.343837, 0.345905, 0.346338overQ = 3^5, 3^7, 3^9, 3^11at base 3{0,1}and at0.015138, 0.012250, 0.011561overQ = 10^2, 10^3, 10^4at base 10 without 9, so|S_F^s(1,Q)|grows like the node countQ^(2 alpha)while every Mertens-type sum overS_Fis bounded by#{n in S_F : n <= Q} = O(Q^alpha). AtQ = 3^11the strict sum has modulus374203.231againstM_F(Q) = -10. - The three candidates the refutation names, each killed at the smallest
Qprinted:M_F(Q), the count-weightedsum_{b in S_F, b <= Q} mu(b) phi_F(b)and the normalisedsum_{b in S_F, b <= Q} mu(b) phi_F(b)/phi(b), withphi_F(b) = #{a in S_F : a <= b, gcd(a,b) = 1}. All three are real, so the first refutation kills all three at once, atQ = 3andQ = 10respectively. - This is the Fourier face of the Farey page's reading that the strict lane's
S1andS2ride the mass: a set whose frequency-1 sum is proportional to its own count has no cancellation at frequency 1 and does not equidistribute.
THEOREM 5, THE DILATE AUTOMATON AND THE CONVERSE
- The dilate is a regular language.
d^(-1) S_Fis recognised least-significant-digit-first by a deterministic automaton whosedstates are the carries of long multiplication byd: from carryrthe digitewrites the output digit(de + r) mod base, which must lie inF, and moves to the carryfloor((de + r)/base), which stays belowdbecausefloor((d(base-1) + d - 1)/base) = d - 1. Afterleveldigitsdcis theleveloutput digits with the terminal carryr_levelwritten above them, so the run accepts exactly whenr_levellies inAcc_d = {0} union (S_F intersect [1, d)), of sizeA_F(d-1) + 1. Proved for a design carrying the digit0. The dilated Mertens sums therefore run over regular sets and not over digit designs. Without0 in Fthe run tests every one of thelevelpadded output digits, a leading output digit0is not a digit ofdc, and the automaton recognises the padded set of mobius instead: at base 3 withF = {1,2},d = 1andlevel = 3it reads8where the true count is14. - The transfer matrix.
T_d(r, r') = #{e < base : (de + r) mod base in F, floor((de + r)/base) = r'}and#{c < base^level : dc in S_F} = e_0 T_d^level 1_(Acc_d). Proved for a design carrying0, path counting. Overbase = 3, 4, 5, everyF, everyd <= 6and everylevel <= 5the identity holds in all840cases carrying0and fails in399of the750without it. - Every column of
T_dsums to|F|, at every base, digit set andd. Proved. The pairs(e, r)in[0, base) x [0, d)are in bijection withv = de + rin[0, d base)by the division algorithm; the column atr'counts thevwithv - base r'inF, and[base r', base r' + base)lies inside[0, d base)for everyr' < d. So the all-ones vector is a positive left eigenvector and the spectral radius ofT_dis|F|for everyd: a dilate carries the design's own mass exponent as its Perron root. The rows sum togtimes#(F intersect (r + gZ))withg = gcd(d, base), hence equal|F|whengcd(d, base) = 1, where#{c < base^level : dc in S_F} <= |F|^levelwith constant1. - Base powers are free. If
0 in FthenM_F(x; base^j d) = M_F(x; d), since appending zero digits neither enters nor leavesS_F. Proved. At base 3{0,1}thed = 3column is thed = 1column,M_F(3^12; 3) = 56with peak61. - The matrix digit transform.
sum of e(ct)overc < base^levelwithdc in S_Fise_0 M(t) M(base t) ... M(base^(level-1) t) 1_(Acc_d)withM(t)(r, r') = sum of e(et)over the digitsecarryingrtor', andM(0) = T_d. Proved. This is the Mobius page's ladderprod_j g_F(base^j t)with the scalar digit symbol replaced by a matrix, and the replacement is the price of the route: an ordered product of non-commuting matrices does not factor, so the sup-over-shiftl^1exponent that carries a Type I estimate for a digit design has no scalar analogue here. The matrix gives the exact count, the exact mass constant and exact evaluation at anyt, and no cancellation inmu. - The mass constant is the accepting set. Conjecture under
gcd(d, base Delta_F) = 1, withDelta_Fthe gcd of the differences of the digits inF. Forgcd(d, base) = 1both sums makeT_d/|F|doubly stochastic, so the stationary law is uniform on each closed class, and where the carry chain is irreducibleA_d(base^level)/|F|^levelconverges to#Acc_d/d = (A_F(d-1) + 1)/d = O(d^(alpha-1)), which is exactly the saving a level of distribution forS_Fat the modulusdwould give. Coprimality to the base alone is not enough: at base 3 withF = {0,2}andDelta_F = 2the dilated = 2hasT_2 = [[2,0],[0,2]], carry1unreachable from carry0, counts2, 4, 8, 16, 32, 64, 128, 256atlevel = 1to8and a constant of1against#Acc_2/2 = 1/2. Verified to three decimals atlevel = 24at every coprimedprinted, on both designs, both havingDelta_F = 1; and the hypothesis is verified as a split over everybase <= 7, everyFcarrying0and every2 <= d <= 24coprime tobaseatlevel = 400,1747agreements and0failures atgcd(d, Delta_F) = 1against0agreements and148failures atgcd(d, Delta_F) > 1. - The saving is not uniform in
d. Refuted, by the base-power ladder:0 in Fmakes(base^j)^(-1) S_Fequal toS_F, so the constant atd = base^jis1at everyjwhile the ceilingbase^(j(alpha-1))tends to0, andA_d(x)/(d^(alpha-1) x^alpha)is at leastbase^(j(1-alpha)), unbounded. Off the ladder the base-smooth dilates are denser than the design in the same way: base 10 without 9 reads1.1111, 1.1358, 1.1111, 1.1413, 1.0700, 1.0343atd = 2, 4, 5, 8, 16, 32against the ceilings0.968781, 0.938537, 0.929003, 0.909237, 0.880851, 0.853352, and the accepting-set law fails there too. The last digit of an element ofS_Fis uniform onFandFis unbalanced modulo a prime dividing the base. It costs the converse its first hypothesis and not its conclusion: over the29base-smoothd <= 1000at base 10 without 9 the constant lies in[0.9273, 1.1637]and over everyd <= 200its inflation over the value at the coprime part ofdlies in[0.9375, 1.1413], bounded on the metered range and unmeasured past it. - (U), the
d-uniform bound|M_F(x; d)| = O_eps(d^((alpha-1)/2) x^(alpha/2 + eps)), is Refuted. For any design carrying0and1, base powers being free givesM_F(x; base^j) = M_F(x), so (U) atd = base^jdemands|M_F(x)| <= C_eps base^(j(alpha-1)/2) x^(alpha/2+eps)at everyj, andalpha < 1drives the right side to0at fixedx, forcingM_Fidentically zero againstM_F(1) = 1. At base 3{0,1}andx = 3^12the left side is56at everyj = 0to12while the bound falls64.0000to5.6187overd = 3^0to3^12, the ratio climbing0.875to9.967. - The converse. Proved from (U'),
|M_F(x; d)| = O_eps(d_base^((alpha-1)/2) x^(alpha/2 + eps))uniform ind >= 1andx >= 1, withd_basethe part ofdcoprime tobase: the same bound with the dilate's true mass constant, whosed = 1case is exactly the square-root ceiling Theorem 3 delivers. ThenG_F(Q) = O_eps(Q^(alpha + eps))andsum_j delta_j^2 = O_eps(Q^(-1+eps)), the denominator lane's conjecture. Writed = a d_baseande = b e_basewithaandbsupported on the primes dividingbase; the parts have disjoint prime support, sogcd(d,e) = gcd(a,b) gcd(d_base,e_base)and the kernel sum factors. Each term is at mostgcd(d,e)^2 (d e)^(-1-alpha/2-eps) (d_base e_base)^((alpha-1)/2) Q^(alpha+2eps); the coprime factor carries exponent-3/2-epsand, withd_base = g uande_base = g v,gcd(u, v) = 1, is at mostzeta(1 + 2eps) zeta(3/2 + eps)^2; the base factor isprod over p | base of sum over i, j >= 0 of p^(2 min(i,j) - (i+j)s)withs = 1 + alpha/2 + eps, equal toprod over p | base of (1 + p^(-s))/((1 - p^(-s))(1 - p^(-alpha-2eps)))and finite becausealpha > 0. Andm_F(Q) >> Q^(1+alpha)/log log Q, the>> Q^alphamembers ofS_Fin the top block belowQeach exceedingQ/basewithphi(b) >> b/log log b, so the rank form divides that down. - The exponent
(alpha-1)/2ond_baseis critical. Proved. At(alpha-1)/2 + deltawithdelta > 0the coprimegsum becomessum of g^(-1+2delta), of sizeQ^(2delta), and the conclusion weakens toG_F(Q) = O(Q^(alpha + 2delta + eps))with no threshold; atdelta = 0it is the harmonic sum and only theepscloses it. The base factor never sees the exponent, so the refutation above does not touch this. - It is one implication and not an equivalence. The reverse is unproved and fails by the natural route: the threshold gives
|S_F(k, Q)| <= k (pi^2 G_F(Q)/6)^(1/2)termwise, and Mobius inversion of Theorem 1 givesd M_F(Q/d; d) = sum over c | d of mu(d/c) S_F(c, Q), hence only|M_F(Q/d; d)| <= (sigma(d)/d)(pi^2 G_F(Q)/6)^(1/2) << log log dtimesQ^(alpha/2+eps), which grows indwhere (U') asks ford_base^(-1/2-eps)decay. - The Mobius-free surrogate. Conjecture.
M_F(Q/d; d)sums over exactly them in S_FbelowQdivisible byd, so square-root cancellation in that mass is|M_F(Q/d; d)| <= N_F(Q; d)^(1/2+eps)withN_F(Q; d) = #{m in S_F : m <= Q, d divides m}, and under that hypothesis the converse reduces toB(Q) = sum over d, e of (gcd(d,e)^2/(d e)) sqrt(N_F(Q; d) N_F(Q; e))beingO(Q^(alpha+eps)), a divisor statement with nomuin it. It stays consistent where (U) does not, reading|M_F(Q/base^j)| <= A_F(Q/base^j)^(1/2+eps)atd = base^j, which is thed = 1ceiling again.B(Q)/Q^alphareads12.5146, 17.8640, 24.7369, 31.5935, 39.0671atQ = 3^4to3^8with local exponents0.955, 0.927, 0.854, 0.824falling towardalpha = 0.630930; no exponent is claimed. - Verified. The transfer matrix
[[1,1],[1,1]]at base 3{0,1}andd = 2, with both carries accepting, counts2^level - 1against literal enumeration of{c : 2c in S_F}at everylevel <= 12, both reading4095atx = 3^12 = 531441. Overd <= 64on both designs no column ofT_dis off|F|; rows are off|F|at 21 and 38 of the 64 and every one of thosedshares a factor with the base. The falsification run meters (U') itself, by the ratio ofmax |M_F(y;d)|overy <= xtod_base^((alpha-1)/2) x^(alpha/2), which (U') asks to stay bounded ind, and it does not fire: at base 3{0,1}andx = 3^12the ratio reads0.9531, 0.7991, 0.9531, 0.6054, 0.9883, 0.3580, 0.5504, 0.8269, 1.0535, 0.4431, 0.5528, 0.3828atd = 1, 2, 3, 4, 5, 7, 8, 11, 13, 16, 22, 31, and at base 10 without 9 andx = 10^7it reads1.1276, 0.9264, 0.5523, 0.6173, 0.8583, 1.2702, 0.4594atd = 1, 2, 3, 4, 5, 7, 11. Under (U)'s exponent the same ratio reads1.1673atd = 3on base 3 against1.0535as the maximum over the coprimed. The raw readingslog max |M_F(x;d)|overlog xadd nothing:0.311823atd = 1and at most0.292046over the ten base 3 dilates, thed = 3row being thed = 1row rather than an independent reading, and0.484570atd = 1against0.489199atd = 7at base 10 with thex = 10^6crossing reversing byx = 10^7; the local exponents swing over0.24to0.845, so none of these readings is an exponent (verb_dilate).
WHAT IT PRINTS
verb_denominator: the frequency-mtable atm = 1, 2, 3, 4, 5, 6, 12for base 3{0,1}atQ = 2187, base 10 without 9 atQ = 1000and the control atQ = 300, exact integer against literal sum; then the frequency-1 ladder toQ = 10^5on base 3,Q = 10^4on base 10 and on the control, with the per-denominator worst deviation, the denominator it sits at, and the count of denominators rounding to the wrong integer.verb_identity: the Fourier form against the truncated Fourier side with its tail bound at four settings; the rank form as an exact rational identity at the same four; the Edwards control atQ = 40; the backward inequality's ratio at six settings.verb_dilate: thed = 2automaton at base 3{0,1}with its transfer matrix, accepting carries and row and column sums, then the count against brute-force enumeration at everylevel <= 12; the column-sum law overd <= 64on both designs with the count of rows off|F|and whether each suchdshares a factor with the base; the accept-mass tableK_dagainst#Acc_d/dandd^(alpha-1)atlevel = 24; the dilated meter at base 3 to3^12and base 10 to10^7with the mass, the endpoint, the peak, the exponent reading atxand atx/base, the ratio against the refuted bound and against (U), and the local exponents; thed = 1againstd = 2ladder at base 3; and the Smith reductionB(Q)atQ = 3^4to3^8with its three normalisations and its local exponent.verb_converse: the scope of the automaton, automaton against brute force overbase = 3, 4, 5, everyF,d <= 6,level <= 5, split on whetherFcarries0; the mass constant against#Acc_d/datlevel = 400over everybase <= 7, everyFcarrying0and every2 <= d <= 24coprime tobase, split ongcd(d, Delta_F); the refutation table for thed-uniform bound at base 3{0,1}andx = 3^12overd = 3^0to3^12; and the base-smooth mass constants at base 10 without 9 overd <= 1000with the inflation ratio overd <= 200.verb_strict: the divisor-sum identity's worst deviation at base 3 to 729 and base 10 to 200; the smallestQwith a nonzero imaginary part on each design; and the refutation ladder toQ = 3^11andQ = 10^4carryingcard,Re,Im,|S_F^s|,|S_F^s|/cardand all three candidates.
| set | Q | card | Re | Im | abs | abs/card | M_F(Q) |
|---|---|---|---|---|---|---|---|
base 3 {0,1} | 243 | 278 | 75.524 | 54.820 | 93.323 | 0.335693 | 4 |
base 3 {0,1} | 2187 | 4286 | 1251.588 | 777.991 | 1473.683 | 0.343837 | -7 |
base 3 {0,1} | 19683 | 67561 | 19889.522 | 12269.876 | 23369.702 | 0.345905 | -4 |
base 3 {0,1} | 177147 | 1080458 | 320313.430 | 193461.533 | 374203.231 | 0.346338 | -10 |
| base 10 without 9 | 100 | 1830 | 24.736 | 12.471 | 27.702 | 0.015138 | 1 |
| base 10 without 9 | 1000 | 147096 | 1730.230 | 503.091 | 1801.887 | 0.012250 | -1 |
| base 10 without 9 | 10000 | 11890654 | 130648.220 | 42778.695 | 137473.540 | 0.011561 | 17 |
- The
cardcolumn is the independent control:278, 4286, 67561, 1080458and1830, 147096, 11890654are the strict counts farey-discrepancy prints from a Mobius sieve that enumerates no fraction, and this study reaches them by literal enumeration.
RUN
uv run python research/lab/py/restricted-franel/restricted_franel.py- From the repository root. One core, 32.4 s summed over the five verbs:
denominator6.7 s,identity10.8 s,strict1.6 s,dilate5.5 s,converse7.8 s. Appenddenominator,identity,strict,dilateorconversefor one. - Domain: the frequency-1 identity to
Q = 10^5at base 3{0,1}andQ = 10^4at base 10 without 9 and the control; the frequency-midentity at seven frequencies; the Franel identity atQ <= 243; the strict refutation toQ = 3^11andQ = 10^4; the dilate automaton overd <= 64, the dilated meter tox = 3^12andx = 10^7, and the Smith reduction toQ = 3^8. - The wall: the literal frequency-1 check at
Q = 10^5on base 10 without 9 costs2.6e9roots of unity, past the machine budget, so that rung is not walked. - Nothing is written to disk.
WITNESSES
- Theorem 1:
ladder_checkandmertens_dilated, withliteral_denominatoras the multi-frequency route. - Theorem 2:
kernel_sumagainstfourier_sidewith the printed tail bound, and againstfarey_delta_squareas an exact rational identity. - Theorem 3: the ratio column of
verb_identity. - Theorem 4:
strict_ramanujan_literalagainststrict_ramanujan_divisor, thenstrict_literalandstrict_weightson the ladder. - Theorem 5:
dilate_matrixanddilate_acceptbuilt intoautomaton_count, againstdilate_maskas the literal enumeration, withsmith_bilinearoverdivisor_countsfor the surrogate. - Theorem 5, the hypotheses:
brute_dilateagainstpositive_countfor the scope,digit_gapfor the mass constant's split, andqfreefor the coprime part in (U) and in the falsification ratio. - The node counts:
strict_literal'scardagainst the published strict counts of farey-discrepancy. - The Farey page THE RESTRICTED FRANEL IDENTITY, the objects and the identity: the dilated sums
M_F(x; d), the frequency-mformula and its frequency-1 case, with the deviations1.09e-11atb = 86293toQ = 10^5at base 3{0,1}and1.36e-12atb = 7247toQ = 10^4at base 10 without 9 and on the control, theQ = 10^5rung on base 10 not walked,0wrong roundings, and the seven frequencies (Theorem 1). - The Farey page THE RESTRICTED FRANEL IDENTITY, the two
L^2forms:G_F(Q), the Fourier form inside its tail bound and the rank form as exact rationals at the four settings, with the Edwards controlG_F(40) = 62.310829atm = 490andsum delta^2 = 0.0104270117(Theorem 2). - The Farey page THE RESTRICTED FRANEL IDENTITY, the inequality and the implication: the false display's witnesses
Q = 5andQ = 37,0violations at every integerQ <= 2187andQ <= 400, the peak0.607927atQ = 1and the maxima0.340071atQ = 253,0.137645atQ = 221and0.086385atQ = 114(Theorem 3). - The Farey page THE RESTRICTED FRANEL IDENTITY, the strict set: the divisor identity, the non-real sums at
Q = 3andQ = 10, the three killed candidates, the ratios0.335693to0.346338and0.015138to0.011561, and the modulus374203.231againstM_F(3^11) = -10(Theorem 4). - The Farey page THE RESTRICTED FRANEL IDENTITY, the dilates: the automaton and its accepting set, the column-sum law and the Perron root
|F|, the matrix digit transform, the accept-mass readings, the base-smooth refutation, the converse under (U) with its critical exponent, the surrogateB(Q)and the falsification run (Theorem 5).
- README.md29.5 kB
- restricted_franel.py23.6 kB