README.md
32.1 kB · markdown
rho-decoupling
- Computes
N_F(level; d), the number of length-leveldigit strings over a digit setFin basebasewhose valuesum f_j base^jis divisible byd, exactly, by dynamic programming on residues moddwith u128 counts; the value map is injective on fixed-length strings, so this counts integers belowbase^levelwith padded digits inF. - For each family and depth the sweep over divisors
2 <= d <= Dwith(d, base) = 1prints the worst normalized errord |N_F(level;d) - fill^level/d| / fill^levelamongdcoprime to the digit-difference gcd, the multiplicative order ofbaseat the worstd, the per-digit rateworst^(1/level), the per-factor ceilinggamma = max_a |sum_{f in F} e(a f / d)| / fillat thatd, and the slack against the proved boundd (1 - 8/(fill^2 d^2))^level. - The proved bound is asserted exactly at every census cell where its hypotheses hold:
d |N - fill^level/d| * (fill^2 d^2)^level <= d * fill^level * (fill^2 d^2 - 8)^levelin big integers, no floats in the claim path. - Families whose digits share a factor also print the unrestricted worst error, the wall where the count does not equidistribute.
- The Type I line accumulates
sum mu(d) (N_F(level;d) - fill^level/d)and its absolute-value trivial bound over squarefreed <= Dcoprime tobaseas exact fractions, printing both relative tofill^leveland their ratio. - The decouple line probes the fixed divisor
d = 7at the deepest level of each family: the per-digit error rate against the per-factor ceiling, across bases3, 4, 5, 10, 100, which is where the large-base decoupling is visible. - Four pinned probes read single divisors of
base^t - 1atbase = 100(d = 101, 9999, 3367, 999999, the first inside the sweep range, the rest beyond it): each line prints the rate(|N - fill^level/d|/fill^level)^(1/level)and the orbit-mean damping factormean = rate * d^(1/level), reading0.106atd = base^2 - 1againstfill^(-1/2) = 0.1005and0.237atd = base^3 - 1againstfill^(-1/3) = 0.216. - The carry DP counts
N_F(level; e)fore | base^t - 1without a state per residue: the value of a length-levelstring is congruent modbase^t - 1toV = sum_{i<t} c_i base^i, wherec_isums the digits in the positions congruent toi mod t, soe | valueexactly whenV = j efor one of the(s+1) g + 1targets withs = floor(level/t)andg = (base^t - 1)/e, and each target is counted by a carry DP across thetbase-basecolumns ofV, the column holdings_ipositions weighted by the coefficients ofP_F(z)^{s_i}forP_F(z) = sum_{f in F} z^f, with carries bounded bys+1. - The cost is
O(t s^3 g)per(level, t, g), polynomial inlevelat everyt; the target is carried as a base-basedigit vector and never as an integer, sot = 40atbase = 10runs althoughbase^t - 1exceedsu128, and the counts themselves stay inu128forfill^level <= 9^40. - The carry sweep runs
base = 3,F = {0,1}andbase = 10,F = {0..9}minus{7}over every moduluse = (base^t - 1)/gwithg | base - 1,e >= 2andt <= level <= 40, and prints per depth the exact ratioSigma_level/Abs_levelof the signed sumSigma_level = sum mu(e) T_level(e)to the absolute sumAbs_level = sum |T_level(e)|over the squarefree moduli of that set, together withAbs_level/fill^level, the count of terms within a factor of ten of the largest, the four largest terms atlevel = 40as exact integers, and the ratio at everylevel = 3..40. - Here
T_level(e) = N_F(level; e) - fill^level/eis an exact rational ande T_level(e) = e N_F(level; e) - fill^levelan exact integer; the sums are formed over a common denominator, so only the printed ratios are floats. mu(e)is read off a complete factorisation, never a table:base^t - 1 = prod_{d | t} Phi_d(base)splits first, then trial division to10^5, Pollard-Brent, and Miller-Rabin on the thirteen bases2to41, deterministic below3.317 * 10^24; a cofactor called prime above that bound is not certified, so its modulus leaves the sums and counts instead in theunknowncolumn of the depth row, and att <= 40both bases certify every prime factor and everyunknowncount reads zero.- The riesz module computes the even moments
S_p(level) = sum_{a mod base^level} |hat F_level(a/base^level)|^pof the digit transformhat F_level(y) = sum_{n in D_level} e(ny)on the value setD_levelof the length-levelstrings, exactly through the identityS_p(level) = base^level E_r(level),p = 2r, whereE_r(level)counts the2r-tuples of length-levelstrings withn_1 + .. + n_r = n_{r+1} + .. + n_{2r} mod base^level, together with the integer countE_r^int(level)(equality over the integers, which isint_0^1 |hat F_level|^p); both come from one carry DP whose state is the pair of carries of the two sides, reduced by the swap symmetry tor(r+1)/2states, so each sequence is C-finite of that order andLambda(p) = lim S_p(level)^{1/level}isbasetimes the Perron root of the transfer matrix. - For each family it prints the characteristic polynomial of that matrix in exact integers, its factorisation over
Z(integer roots by a scan bounded by the maximal row sum, monic quadratic factors by a bounded search over the divisors of the constant term, a cofactor of degree at least four certified irreducible modulo a printed prime by the distinct-degree test), an exact two-sided bracket on the Perron root from the Collatz-Wielandt bounds at a positive rational vector on the states reachable from the zero carries, the factor holding the Perron root located by a sign change on the bracket with every other factor certified nonzero there, thenLambda(p), its exponenttheta_p = log Lambda(p)/log base, the floorsmax(fill^p, base fill^{p/2})and the ceilingbase fill^{p-1}asserted at every depth tolevel = 60, the recurrence asserted at every depth, the ratioE_r/E_r^int, the share of the zero frequency inS_p, and the growth of the primitive partS_p(level) - fill^p S_p(level-1). - The type II rows print per family the exponent
theta_4/4 + 1/4that Holder with the fourth moment and Parseval on the bilinear side yields for a Type II sum over the digit set, its miss againstalpha = log fill/log base, the bilinear saving1 - alpha/2such a route would need on the minor arcs, and the large-values thresholdeta_4 = (1 + 3 alpha - theta_4)/2below which the fourth moment counts the frequencies near the maximum better than Parseval does; the arcs rows evaluate|hat F_level(a/base^level)|^2at every grid point at one or two depths per family, mark the points within1/(d x^{3/5})of a fraction with denominatord <= x^{2/5}, and print the share of the Parseval mass off those points against(fill/base)^level, the share below which a moment route could still work. - The riesz families are the census families with every reflection class of one excluded digit at
base = 10; the higher moments runp = 6, 8, 10atbase = 3,F = {0,1}andp = 6at four more families. - The large module reads the same grid as a complex table,
hat F_level(a/base^level)at everyaforbase^level <= 10^6, asserts Parseval and the exactS_4on it to10^{-10}relative, and for six thresholdsetaper cell (0,eta_4/2,eta_4,2 eta_4,0.25,0.5) forms the large setA_eta = {a : |hat F_level(a/base^level)| >= fill^level x^{-eta}}: it prints its size against the fourth-moment countS_4 fill^{-4 level} x^{4 eta}and the Parseval countbase^level fill^{-level} x^{2 eta}, both asserted above the size, the least gap between two of its points in grid units, the number of maximal runs of consecutive grid points and the longest run, the floorsum_{a in A} |hat F_level|^2/fill^levelof the set's large-sieve constant for polynomials supported onD_level, and that constant itself as the top eigenvalue of the Gram matrixhat F_level((a - a')/base^level)on the set, bracketed by the Rayleigh quotient of a power iteration below and by the least of the Gershgorin row bound, the Frobenius norm andbase^levelabove, while the set has at most1500points, capped byfill^{level-j} base^jwhen the set lies on the level-jgrid; every printed float is a two-sided band or rounded safe, bounds up and floors down. - Beside the set it prints the inequality chain of a large-values Type II with bounded coefficients: the small set's exponent
min(alpha + 1/2 - eta, (1 + alpha)/2)from the threshold and Parseval on the bilinear side, the large set's exponent from the sup and the count,alpha + log(#A)/(2 log x), from the set's own mass,log(sum_{a in A} |hat F_level|^2)/(2 log x), and from the Parseval cap(1 + alpha)/2, all against Parseval on the bilinear side (the large-sieve constant of any grid subset forbase^levelconsecutive frequencies isbase^levelexactly), the sup form with the fourth-moment count in place of the size, the count-basedc sup = alpha - max(small, large sup)and the chainc = alpha - max(small, min of the three), together with thel^1exponentc_1 = log(sum_a |hat F_level(a/base^level)|)/log xat the depth. - A per-cell row prints
nu_4(level),eta_4,c_1, the second-largest grid value overfill^level, which equalsmax_{b != 0 mod base} |g_F(b/base)|/fillby the last non-integer factor of the product formula, and the thresholdeta_1(level)below which the large set is the zero frequency alone. - A witness row per cell counts, on the box
m, l in [M, 2M)withM = floor(x^{1/2}/2), the representationsR = #{(m, l) : ml in D_level}, which is the Type II sum ata_m = b_l = 1, and the balanced sum ata_m = 1_{base | m},b_l = 1,R_base - fill^level base^{-level} M #{m : base | m}, each againstfill^leveland as an exponent ofx. - The menergy module computes the multiplicative energy
E_x(level) = #{(n_1, n_2, n_3, n_4) : n_i in D_level, n_i != 0, n_1 n_2 = n_3 n_4}of the same digit column exactly, by forming theK^2products of the column with itself inu128, hashing them into as many partitions as the memory cap allows, sorting each partition and summing the squares of the run lengths, which issum_m r(m)^2forr(m)the number of ordered factorisations ofminside the column; it printsK,E_x(level), the exponentlog E_x/(level log base)against2 alpha, the ratio ofE_xto the exact diagonal count2K^2 - K, the largestr(m), the share of the excess over the diagonal that the shift solutions explain, and the same ratio and the direct quotient for one uniform random subset of[1, base^level)of the same size drawn by the module's own generator. - Two bounds are asserted at every cell in exact integers: the band
2K^2 - K <= E_x(level) <= K^3(the diagonal and the swap agree on theKquadruples withn_1 = n_2; fixing three entries determines the fourth) andE_x(level) <= K^2 max_m r(m), which withr(m) <= d(m)pins the exponent at2 alphafor every digit set. - When
0is a digit the module also asserts a third floorE_x(level) >= 2K^2 - K + T(level), whereT(level)counts the shift solutions(base^i u, base^j v, base^{i'} u, base^{j'} v)withi + j = i' + j',i != i',u != vandu, vof nonzero last digit, a closed form infillandlevelalone; the printed share(2K^2 - K + T)/E_xsays how much of the energy the diagonal and the shifts explain. - The type II rows sweep the dyadic boxes
m in [M, 2M),l in [N, 2N)withMNin[x/32, x/4]andM, N >= 8, and for each box countR = #{(m, l) : ml in D_level}, the restricted energyE_x(M, N) = sum_m (#{l : ml in D_level})^2, and the balanced energyE_x^bal(M, N) = sum_{l_1, l_2} |W(l_1, l_2)|forW(l_1, l_2) = sum_m psi(m l_1) psi(m l_2),psi = 1_{D_level} - K/base^levelwithKthe column size, whose Cauchy-Schwarz bound(M E_x^bal)^{1/2}is what such a route returns for a bilinear sum with bounded coefficients; the pairs withWoff its baseline are counted by sorted keys, the rest by the histogram of the column counts, and the identitysum_{l_1, l_2} W = sum_m (c_m - N fill^level/base^level)^2, derived independently, is printed as a self-check. - Each cell prints the extreme boxes by the bound exponent among the boxes carrying a representation, the worst box inside the arc regime
M, N >= x^{2/5}, the ratio of the bound to that box's trivial bound, and the marginalpha - log bound/log x, for the digit column and for a random column of the same density side by side. - The sigma rows draw random sign vectors for both coefficient sequences on a fixed box and assert the Cauchy-Schwarz bound at every draw, printing the largest ratio of the measured
|Sigma|to the bound; the identitysum_{l_1, l_2} W = sum_m (c_m - N delta)^2printed beside them is an algebraic rearrangement of the sweep's own quantities and pins summation order, not correctness. - The signed module computes the off-diagonal correlation
Sigma_b = sum_{l_1 != l_2} b_{l_1} b_{l_2} W(l_1, l_2)of a coefficient vectorb, the one part of the balanced energy that keeps the coefficients' signs, for nine coefficient vectors at once: the constant1, the Mobius function, the Liouville function, three seeded random sign vectors, and three seeded random sign vectors supported on the squarefree integers only, which is the Mobius support with the Mobius arithmetic removed. - The identity it evaluates is
Sigma_b = sum_m V_m(b)^2 - sum_l b_l^2 W(l, l)forV_m(b) = sum_l b_l psi(ml), the full quadratic form being a sum of squares, and it is carried as the exact integerx^2 Sigma_b = x^2 (Q - P) - 2 K x (B C - P) + M K^2 (B^2 - S_2)ini128, whereQ = sum_m u_m^2,u_m = sum_{l, ml in D_level} b_l,B = sum_l b_l,C = sum_l b_l c'_l,P = sum_l b_l^2 c'_landS_2 = sum_l b_l^2; the printedSigma_bis that numerator overx^2rounded to the nearest integer and every ratio and exponent comes from the numerator itself. - Mobius and Liouville come from one smallest-prime-factor sieve to twice the largest box side; the random vectors come from the module's own generator at fixed seeds, so the same vector is used at every family, depth and box.
- The signed sweep runs the same dyadic band
MN in [x/32, x/4]at two depths for each of four families, selecting the box that maximises|Sigma_1|over the band and the box that maximises it inside the arc regimeM, N >= x^{2/5}, and repeats the deeper depth on a random column of the same size, each on its own selected box; each row printsSigma_1,Sigma_mu,Sigma_lambda, their exponentslog |Sigma_b| / log x, the ratios toSigma_1, the min-max band of each triple of random draws, and the ratios to the root mean square a random sign vector reaches. - That root mean square,
(2 sum_{l_1 != l_2} W(l_1, l_2)^2)^{1/2}, is the yardstick a signed sum is read against, and it is computed assum_{l_1, l_2} W^2 = sum_{m, m'} (sum_l psi(ml) psi(m'l))^2over the pairs of the shorter side, the only float in the module and never in a claim path. - Each row also prints the min-max, over the eight signed coefficient vectors, of the full quadratic form divided by its own diagonal, and the ratio of the best possible bounded-coefficient sum
sum_m |V_m(b)|to the Cauchy-Schwarz bound(M sum_m V_m(b)^2)^{1/2}, which says how much the inequality itself loses. - A separate section drives the quadratic form down by greedy single-sign flips from the all-ones vector on a
+-1vector supported on the whole side, printing the flips used, the form, its diagonal, the resulting bound and the diagonal floor(1 - delta)(M R)^{1/2}; the engineered vector is built with full knowledge of the column, is available to no decomposition, and exists to show that the bound has no floor over all bounded coefficients. - A summary row counts, over the digit cells, how often the Mobius value falls below, inside and above each band, how often it is smaller in absolute value than the Liouville value, the extreme ratios to the root mean square for the Mobius vector and for the support-matched draws, the extreme values of the form over its diagonal, the extreme Cauchy-Schwarz losses, and how often the signed sums fall below the box representation count
Rwhile the unsigned one stays above it. - Two checks hold the signed rows: at
b = 1the quantitiesQandCare asserted equal to theE_x(M, N)andRof the box sweep, computed along a different data path, and the tests re-evaluate the definition pair by pair overl_1 != l_2withWby a direct triple loop in the same exact integers. - The signed families are
F = {0,1}atbase = 3and depths12, 14,{0,1,2}atbase = 4and depths8, 9,{0,1,2,3}atbase = 5and depths7, 8, and one excluded digit atbase = 10and depths5, 6. - The vaughan module censuses the real coefficient sequences a Vaughan decomposition hands the bilinear sum, on the same boxes the signed module selects:
b^U_l = sum_{d | l, d <= U} mu(d), thel-side coefficient of the Type II term of Vaughan's identity for bothmuand the von Mangoldt function, and the log-weighted variantsum_{d | l, d <= U} mu(d) log(l/d), each at the two cutsU = floor(x^(1/3))andU = floor(x^(2/5)), together with the sign vector of the second. - Each sequence is built by adding
mu(d)andmu(d) log dalong the multiples of every squarefreed <= U, so the two pieces cost one pass per divisor, and the tests check both against a direct loop over the divisors of everylto1200at three cuts. - The Type II region of the identity is
m > Vandl > U, so the boxes that are Type II cells are those with both sides above the cut; the summary counts them and the wide boxes are printed beside them as the same sequence read where the identity does not put it. - For each box it prints the two cuts, the ratio
2N/Uthat bounds the box maximum of the coefficient, the share of the box where the coefficient is nonzero, the box maximum, the full quadratic form over its own diagonal for all five sequences, the diagonal sharew = diag/((1 - delta)^2 R)of the normalised coefficient, and its Cauchy-Schwarz bound over the diagonal floor both normalised by the box maximum and unnormalised. - The split
bound/floor = (form/diag * w)^(1/2)is an identity in the definitions and is asserted at every box, so a small bound is attributed to cancellation, to a sparse or small coefficient, or to a large box maximum, and never to all three at once. - The summary row prints the extremes of every column over the sixteen boxes and the eighty coefficient cells, how many sit inside the band the signed module's summary prints, how many fall below one tenth, how many boxes carry a coefficient with box maximum one, and the closest any cell comes to either end of that band.
- The menergy census runs the census families to the depth the
K^2product count allows,level = 12at the two-digit sets ofbase = 3,level = 8atbase = 4,level = 6andlevel = 7atbase = 5,level = 4andlevel = 3at one excluded digit ofbase = 10,level = 9for the pair andlevel = 2for the lower half atbase = 100, and the box sweep one to three levels deeper, tolevel = 14atF = {0,1},base = 3, where only membership in the column is needed. - Census families: the three two-digit sets at
base = 3;{0,1,2}atbase = 4;{0,1,2,3}and{0,2,4}atbase = 5; one excluded digit atbase = 10andbase = 100; the lower half{0..49}and the pair{0,1}atbase = 100. Depths tolevel = 96, divisor ranges toD = 500; floats appear only in printed readouts, truncated at forty decimal digits.
RUN
bash scripts/cargo.sh cargo run --release -p rho-decoupling- Under twenty seconds for the whole run, census, carry sweep, riesz moments, large values, multiplicative energy and signed correlation together; prints only, writes nothing.
cargo test -p rho-decouplingruns thirty-nine tests: the residue DP against brute-force string enumeration at four bases, the residue vector total againstfill^level, the exact splitting ofN_F(level; d1 d2)acrossd1 | base,(d2, base) = 1atbase = 6, the Mobius sieve against hand values andM(100) = 1, and the exact bound inequality atbase = 3,F = {0,1},level = 16,d <= 60.- The riesz tests pin the moment DP against a sumset histogram at
p = 4on every family and atp = 6, 8on four, the exactS_4against direct evaluation of the product formula on the grid andE_2^intagainst the2 base^level-point quadrature atlevel <= 6, the factoriser and the irreducibility test on known polynomials, the scaling identities of{0,2}and{1,2}against{0,1}atbase = 3, every fourth-moment bracket against the growth ratio atlevel = 40, the rendered algebra, value and type II rows of seven families, the minimal polynomials of seven higher moments, and the arcs counts atbase = 3,level = 12. - The large tests pin the complex grid against a direct sum over every string at four families, the large set at five thresholds against the same direct values, the Gram eigenvalue at
{0}(fill^level), at thebase^2-grid (fill^{level-2} base^2, the witness set of the sparse large sieve) and at the full grid (base^level), the run statistics on hand sets, and the witness counts against an enumeration of the column. - The menergy tests pin the partitioned energy against a four-fold brute-force enumeration of the solutions at six families, against itself at three partition caps, the band and the divisor ceiling at four families, the scaling identity
E_x(cF) = E_x(F)at{0,2}against{0,1}and{0,2,4}against{0,1,2}together with the failure of translation invariance at{1,2}, the restricted and balanced energies of a box against a direct triple loop over the box, the shift count against an enumeration of the shift solutions themselves, and the Cauchy-Schwarz bound against sixty random sign draws on three boxes. - The signed tests pin the Mobius and Liouville sieves against their first twelve values, the exact identity against a pair-by-pair evaluation of the definition at four families and all nine coefficient vectors, the unsigned quantities against the box sweep at three families, and the collapse of the unsigned correlation on a random column of the same size.
- The vaughan tests pin both pieces against a direct divisor sum at three cuts, the form and the diagonal against a pair-by-pair evaluation of the definition with the kernel by a direct triple loop, and the form of the integer piece against the exact integer
sum_m (x u_m - K B)^2over the box. - The carry tests pin the carry DP against brute-force enumeration at
base = 3, 5, 10over everyt <= leveland everyg | base - 1, against the residue DP at every reachablee <= 30000atlevel = 8andlevel = 12in both carry families with zero mismatches, the certification of every cofactor ofbase^t - 1att <= 40in both bases, and as rendered strings the depth rows, the four largest terms atlevel = 40, the ratio sequencelevel = 3..40, themutables and seven factorisations; one more test pins the digit-vector helpers.
WITNESSES
- mobius.md digit strings across divisors, the uniform equidistribution bound
|N_F(level;d) - fill^level/d| <= fill^level (1 - 8/(fill^2 d^2))^levelfor(d, base) = 1withdcoprime to every digit difference: asserted exactly at every census cell. - mobius.md digit strings across divisors, the worst divisors are the pinned ones: at every family's deepest level the sweep argmax
dhasord_d(base) <= 8, sod | base^t - 1witht <= 8, e.g.d = 164atbase = 3,d = 143atbase = 10,d = 101, 303atbase = 100; shallow levels can stray (d = 199,ord = 99, atbase = 10,level = 6). - mobius.md digit strings across divisors, the decoupling: at
d = 7the per-digit error rate falls0.49, 0.33, 0.25, 0.11, 0.017alongbase = 3, 4, 5, 10, 100with one digit excluded, but stays0.50forF = {0,1}atbase = 100: the gain is carried by the fill, not the base alone. - mobius.md digit strings across divisors, a signed pinned sum against its absolute sum: over the moduli
e = (base^t - 1)/gwithg | base - 1,e >= 2,t <= level, squarefree, the ratioSigma_level/Abs_levelof the Mobius-signed sum ofT_level(e)to the sum of|T_level(e)|is verified at-0.211, -0.123, +0.069, -0.498forlevel = 10, 20, 30, 40atbase = 3,F = {0,1}and at+0.812, -0.495, -0.127, -0.192forbase = 10with one digit excluded, whileAbs_level/fill^levelfalls from4.7e-2to2.1e-4and from1.1e-5to3.9e-12. - mobius.md digit strings across divisors, no cancellation trend to depth
40: the ratio is printed at everylevel = 3..40on both families and shows no decay toward zero with depth, swinging across the band and reading-0.50atlevel = 40in base three and-0.99atlevel = 29in base ten, so the sign ofmu(e)buys no extra saving over the trivial absolute bound on this modulus set at these depths. - mobius.md digit strings across divisors, few terms carry the sum: at
level = 40the terms within a factor of ten of the largest count6atbase = 3out of29and4atbase = 10out of60, the four largest sitting att = 7, 9in base three andt = 5, 7, 8, 10in base ten, each printed as the exact integere T_level(e). - mobius.md the meter and its yardstick, the fourth moment:
sum_{a mod base^level} |hat F_level(a/base^level)|^4 = base^level E(level)withE(level)the additive energy of the strings modulobase^level, counted by the carry DP;Lambda(4) = 18atbase = 3for all three two-digit sets,2(23 + sqrt 353)atbase = 4,F = {0,1,2},(275 + 5 sqrt 2369)/2atbase = 5,F = {0,1,2,3},95atF = {0,2,4}, five quadratic irrationals across the excluded digit atbase = 10, from6566.412229970(digit5excluded) to6567.410368845(digit0or9excluded),8335000for the lower half and600for{0,1}atbase = 100, each inside[max(fill^4, base fill^2), base fill^3]and strictly abovefill^4. - mobius.md the meter and its yardstick, the large-values route to Type II: splitting the grid at
|hat F_level| >= fill^level x^{-eta}, bounding the large set by the fourth-moment count and Parseval and the rest by the threshold and Parseval gives the exponentmax(min(alpha + 1/2 - eta, (1 + alpha)/2), min((1 + alpha)/2, log(sum_{a in A_eta} |hat F_level|^2)/(2 log x))), read at66cells over six families and equal to(1 + alpha)/2, thel^2route, at every one of them; at{0,1}base 3,level = 12,eta = eta_4, the large set holds407adjacent-or-isolated grid points in331runs against the fourth-moment count4096, and its large-sieve constant1.06009e5..2.13280e5sits nearly at its ownl^2floor1.05611e5. - mobius.md the meter and its yardstick, the witnesses: at
a_m = b_l = 1the Type II sum on the boxM = N = floor(x^{1/2}/2)is the representation countR, which reads0.38 fill^levelat{0,1}base 3level = 12and0.26 fill^levelat one excluded digit of base 10level = 6, and the balanced sum ata_m = 1_{base | m},b_l = 1reads0.10 fill^leveland0.0044 fill^levelthere, so no boundx^{alpha - c}uniform over bounded coefficients holds for either sum at the nine dense cells; the two sparse cells are void at that box. - mobius.md the meter and its yardstick, the moment route to Type II: with Parseval on the bilinear side the Holder exponent
theta_4/4 + 1/4missesalphaby at least0.250000(one digit excluded atbase = 100) and by0.446753atbase = 100,F = {0,1}, and any pointwise bilinear input would have to beatx^{1 - alpha/2}on the minor arcs, below the root mean square of the bilinear sum, at every family. - mobius.md the meter and its yardstick, the multiplicative energy of a digit column has no exponent of its own:
2K^2 - K <= E_x(level) <= K^2 max_m r(m)withr(m) <= d(m)forcesE_x(level) = fill^{2 level} x^{o(1)}for every digit set, and the census reads the exponent falling1.475642, 1.410978, 1.356938atlevel = 4, 8, 12forF = {0,1}atbase = 3toward2 alpha = 1.261860, withmax_m r(m) = 56atlevel = 12. - mobius.md the meter and its yardstick, the excess over the diagonal is shift structure:
E_x/(2K^2 - K)reads1.7523for{0,1}atbase = 3andlevel = 12against1.0004for a random column of the same size, and the closed-form shift countT(level) = 11144220brings the explained share to0.7603; at{1,2}, where no digit is zero andT = 0, the ratio is1.0006and the share0.9994. - mobius.md the meter and its yardstick, the type II route through the multiplicative energy:
|Sigma|^2 <= M E_x(M, N)andE_x(M, N) <= 2N E_x(level)^{1/2} x^{o(1)}give|Sigma| <= (2MN)^{1/2} x^{alpha/2 + o(1)}, above the trivialx^{alpha}by(1 - alpha)/2for every digit set withalpha < 1, the same miss as square-root cancellation on the fourth-moment route. - mobius.md the meter and its yardstick, the Mobius signs cancel the correlation no better than unstructured signs on the same support: over sixteen boxes across four families at two depths each,
|Sigma_mu|divided by the random-sign root mean square stays in[0.0913, 0.7178]while the three support-matched draws stay in[0.0359, 1.5048], and both signed sums fall below the box representation countRat all sixteen while the unsignedSigma_1stays above it at all sixteen. - mobius.md the meter and its yardstick, the apparent advantage of the Mobius function is its support: against three plain random sign vectors it sits below the band at eight cells, inside at seven and above at one, but against three random sign vectors on the squarefree integers it sits below at three, inside at nine and above at four, chi-square
0.375against the uniform-rank null4, 8, 4;|Sigma_mu| < |Sigma_lambda|at fifteen of sixteen cells and the support-matched control removes the whole gap. Sixteen cells against three draws exclude a large effect and nothing smaller: a true below-rate of0.41still gives three or fewer at probability0.0556. - mobius.md the meter and its yardstick, the bound has no floor over all bounded coefficients: a
+-1vector built by greedy single-sign flips against a known column drives the full quadratic form to0.0006of its own diagonal atF = {0,1,2,3},base = 5,level = 7,M = 8,N = 2048, and the Cauchy-Schwarz bound to0.0265of the diagonal floor(1 - delta)(M R)^{1/2}, reaching0.4909of that floor inside the arc regime atM = N = 128; the diagonal is sign-blind but the off-diagonal can cancel it, so no floor statement holds beyond the coefficient vectors actually swept. - mobius.md the meter and its yardstick, the off-diagonal correlation is a pure structure statistic: its expectation over a random column of the same density is zero at every coefficient vector, and the measured
Sigma_1falls from67382to349atF = {0,1},base = 3,level = 14, from41590to122atbase = 4, from72359to-27atbase = 5and from823318to-2746at one excluded digit ofbase = 10. - mobius.md the meter and its yardstick, the balanced route returns the trivial box bound: inside the arc regime
M, N >= x^{2/5}the marginalpha - log bound/log xfor{0,1}atbase = 3falls0.035675atlevel = 12to0.027009atlevel = 14while the ratio of the bound to the box's own trivial bound rises1.0134to1.0730, and over all admissible boxes the margin is already negative atlevel = 12; the random column of the same density beats the digit column at every cell.
- mobius.md the meter and its yardstick, the coefficient a real decomposition hands the method is not the coefficient the method would need: at the eight swept boxes with both sides above
x^(2/5), the boxes the identity produces and where the box sits within a factor8.73of the cut, the Type II coefficient of Vaughan's identity takes values in{-1, 0, 1}at seven of the eight and its quadratic form stays between0.6929and1.2045of its own diagonal, where a sign vector engineered against the column reads0.2043on such a box, and no one of the eighty coefficient cells of the census falls below0.3138. - mobius.md digit strings across divisors, the fixed-
fillwall: forF = {0,1}atbase = 100the worst normalized error overd <= 500still reads7.2atlevel = 96, decaying by factor0.9929per digit, pinned atd = 303 | base^2 - 1.